Java解二元一次不定方程[ax+by=c].docVIP

Java解二元一次不定方程[ax+by c] 什么是GCD？GCD是最大公约数的简称。在开头，我们先下几个定义：①a|b表示a能整除b（a是b的约数）②a mod b表示a-[a/b]b（[a/b]在java中相当于a/b）③ a,b 表示a和b的最大公约数④a和b的线性组合表示ax+by（x,y为整数）。我们有：若d|a且d|b，则d|ax+by（这很重要！） 线性组合与GCD现在我们证明一个重要的定理：gcd a,b 是a和b的最小的正线性组合。证明：设gcd a,b 为d，a和b的最小的正线性组合为s∵d|a且d|b，∴d|s。而a mod s a-[a/s]s a-[a/s] ax+by a 1-[a/s]x -b[a/s]y亦为a和b的线性组合 ∵a mod s s，a mod s不能是a和b的最小的正线性组合 因为已经假设s 是 a和b的最小的正线性组合 ∴a mod s 0，即s|a同理由s|b∴s为a,b的公约数∴s d∵d|s∴d s。证毕。 由这条定理易推知：若d|a且d|b，则d|gcd a,b Euclid（欧几里德）算法现在的问题是如何快速的求gcd a,b 。穷举明显不是一个好方法（O n ），所以需要一个更好的方法。首先我们先提出一个定理：gcd a,b gcd b,a-bx （x为正整数）。 证明：设gcd a,b d,gcd b,a-bx e,则∵d|a,d|b∴d|a-bx∴d|gcd b,a-bx ，即d|e∵e|b,e|a-bx∴e|bx+ a-bx ，即e|a∴e|gcd a,b ，即e|d∴d e。证毕。 这个定理非常有用，因为它能快速地降低数据规模。当x 1时，gcd a,b gcd b,a-b 。这就是辗转相减法。当x达到最大时，即x [a/b]时，gcd a,b gcd b,a mod b 。这个就是Euclid算法。它是不是Euclid提出的我不知道，但听说是在Euclid时代形成的，所以就叫Euclid算法了。程序非常的简单： public int Gcd int a, int b if b 0 return a; return Gcd b, a % b ; 当然你也可以写成迭代形式： 　public 　int Gcd int a, int b while b ! 0 int r b; b a % b; a r; return a; Euclid算法比辗转相减法好，不仅好在速度快，而且用起来也方便。两种算法都有一个隐含的限制：a b。用辗转相减法时，必须先判断大小，而Euclid算法不然。若a b，则一次递归就会转为gcd b,a ，接着就能正常运行了。 扩展Euclid前面我们说过，gcd a,b 可以表示为a和b的最小的正线性组合。现在我们就要求这个最小的正线性组合ax+by中的x和y。这个可以利用我们的Euclid算法。从最简单的情况开始。当b 0时，我们取x 1,y 0。当b≠0时呢？假设gcd a,b d，则gcd b,a mod b d。若我们已经求出了gcd b,a mod b 的线性组合表示bx+ a mod b y，则gcd a,b d bx+ a mod b y bx+ a-[a/b]b y ay+b x-[a/b]y 那么取x y，y x-[a/b]y就可以将gcd a,b 表示为a和b的最小的正线性组合，这样就可以在Euclid的递归过程中求出x和y。程序：? import java.util.Scanner; public class GED static long x0, y0; public static void main String[] args Scanner scan new Scanner System.in ; long a scan.nextInt ; long b scan.nextInt ; long d gcd a , b ; //这个地方太麻烦了，将就吧 System.out.println gcd +a+,+b+ +a+*+x0+++b+*+y0 ; //将gcd a,b 表示为a,b的线性组合。 public static long gcd long a, long b long t, d; if b 0 x0 1; y0 0; return a; d gcd b, a % b ; t x0; x0 y0; y0 t - a / b * y0; // System.out.println a +a+ b +b+ x0 +x0