/*
裴蜀定理:若ax+by = z,则 gcd(a,b)| z
如何求解 (以下讨论a>b)
显然当 b=0,gcd(a,b)=a。此时 x=1,y=0;
当a>b>0 时
设 ax1+ by1= gcd(a,b);
bx2+ (a mod b)y2= gcd(b,a mod b);
根据欧几里德原理有 gcd(a,b) = gcd(b,a mod b);
则:ax1+ by1= bx2+ (a mod b)y2;
即:ax1+ by1= bx2+ (a - [a / b] * b)y2 = ay2+ bx2- [a / b] * by2;(a mod b = a - [a / b]*b;[a / b]代表a整除b)
也就是ax1+ by1 = ay2 + b(x2- [a / b] *y2);
根据恒等定理得:x1=y2;y1=x2- [a / b] *y2;
这样我们就得到了求解 x1,y1 的方法:x1,y1 的值基于 x2,y2
由引理我们知道:ax+by = z,z为gcd(a,b)若干倍,所以我们先求解ax+by = gcd(a,b),再将求出的解乘以 z/gcd(a,b)就好了。
*/
public class 欧几里得算法 {
static long x;
static long y;
public static void main(String[] args) {
printxy(2,3,1);
System.out.println(x); //定义的x和y
System.out.println(y);
}
//辗转相除法求最大公约数,就是欧几里得算法
public int gcd(int a,int b){
if(b==0) return a;
else return gcd(b,a%b);
}
//扩展欧几里得算法的功能就更强大了,它可以用来求二元一次方程的通解.
public static long exGcd(int a,int b)//网上有很多&x,&y 的方法,在Java中没有,只能定义全局变量
{
if(b==0)
{
x = 1;
y = 0;
return a;
}
long r = exGcd(b,a%b);
long t;
t = (int) x;
x = y; //这一步其实很好理解,就是a换成b,b换成a%b
y = t-a/b*y; //这一步是a%b = a- (a/b)*b
return r;
}
//求求二元一次方程的通解
public static void printxy(int a,int b,int m){
long l = exGcd(a, b);
long beishu = m/l; //获取倍数
if(m%l != 0) System.out.println("无解");
x = x * l; //通解x
y = y * l; //通解y
System.out.println(l);
}
}
欧几里得算法求解x,y。和扩展欧几里得
最新推荐文章于 2024-08-03 21:02:21 发布