拓展欧几里得算法
//拓展欧几里得算法
/*
找出整数对(x,y),使得ax+by=gcd(a,b),(x,y)可以为负数或者0
*/
void gcd(int a,int b,int &d,int &x,int &y)
{
if(!b){d=a;x=1;y=0;}
else
{
gcd(b,a%b,d,y,x);y-=x*(a/b);
}
}
若a,b,c为任意整数. 若方程ax+by=c的一组整数解为(x0,y0), 则它的任意整数解可以写成(x+kb’,y-ka’),其中a’=a/gcd(a,b), b’=b/gcd(a,b)
有了这个结论:
移项得到ax+by=-c, 然后求出一组解即可得到所有解设a,b,c为任意整数, 方程ax+by=gcd(a,b)的一组解是(x0,y0), 则当c是gcd(a,b)的倍数时ax+by=c的一组解是(x0c/g,y0c/g); 当c不是gcd(a,b)的倍数时无整数解