基于最大公约数欧几里得扩展算法求解算法导论中模线性方程解。具体要结合算法导论中的有关数论算法章节理解,具体代码如下:
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
System.out.println("solution:"+(x0+i*(n/d))%n);
}
}
}
执行结果:
最大公约数d=3 And x=1 And y=-129
x0=4
solution:0
solution:2