Java实现算法导论中求解模线性方程解（基于最大公约数欧几里得扩展算法）

package cn.ansj;

/*
假设方程ax=b(mod n)有解，且x0是方程的任意一个解，
则该方程对模n恰有d个不同的解(d=gcd(a,n))，分别为：xi=x0+i*(n/d) mod n 。
*/
public class ModularLinearEquation {
static int x,y;
public static int extend_gcd(int a, int b){
if (b == 0){
x = 1;
y = 0;
return a;
}
int re = extend_gcd(b, a % b);
int tmp = x;
x = y;
y = tmp - a / b * y;
return re;
}

public static void main(String args[]) {
int a=777;
int b=66;
int n=6;
int d=extend_gcd(a,n);
System.out.println("最大公约数d="+d+" And x="+x+" And y="+y);
if(b%d!=0) System.out.println("no solution");
else{
int x0=x*(b/d)%n;
System.out.println("x0="+x0);
for(int i=1;i<d;i++)
System.out.println("solution:"+(x0+i*(n/d))%n);
}
}
}

最大公约数d=3 And x=1 And y=-129
x0=4
solution:0
solution:2