拓展GCD

问题引入：求直线ax+by+c=0上有多少个整点(x,y)满足x∈[x1, x2], y∈[y1, y2]。
解决这个问题之前，首先学习扩展欧几里德算法——找出一对整数(x,y)，使得ax+by=gcd(a,b)。
给出exgcd的代码

/*exgcd求解x,y同时还把gcd(a,b)求出来了，
d存的就是gcd(a,b)*/
void exgcd(int a, int b, int &d,int &x, int &y)  
{
	if (!b) { d = a; x = 1; y = 0; }
	else { exgcd(b, a%b, d, y, x); y -= x * (a / b); }
}

下面我们来证明

已知ax+by=gcd(a,b)，求x,y。
由gcd(a,b)=gcd(b,a%b)得ax+by=bx+(a%b)y
bx+(a%b)y=bx+(a-(a/b)*b)y=bx+ay-(a/b)*b*y=ay+b(x-a/b*y)
所以我们得到ax+by=ay+b(x-a/b*y)
即exgcd(ax+by)=exgcd(ay + b(x-a/b*y))
所以可以递归求解了，递归出口为b==0
显然，此时d=a，x=1，y=0

由一个解求出其他解
上面求出了ax+by=gcd(a,b)的一组解(x1,y1)，任取另外一组解(x2,y2)，则ax1 +by1 =ax2 +by2 （它们都等于gcd(a,b)），变形得a(x1-x2)=b(y2-y1)。假设gcd(a,b)=g，方程左右两边同时除以g（如果g=0，意味着a或b等于0，可以特殊判断），得a’(x1-x2)=b’ (y2-y1)，其中a’=a/g，b’=b/g。注意，此时a’和b’互素，因此x1 -x2 一定是b’的整数倍。设它为kb’，计算得y2 -y1 =ka’。注意，上面的推导过程并没有用到“ax+by的右边是什么”，因此得出如下结论。

结论
「 设a, b, c为任意整数。若方程ax+by=c的一组整数解为(x0 ,y0)，则它的任意整数解都可以写成(x0+kb’, y0-ka’)，其中a’=a/gcd(a,b) b’=b/gcd(a,b)，k取任意整数。 」

题目分析
是否还记得一开始的题目，移项得ax+by=-c，然后求出一组解即可。
例：6x+15y=9。根据欧几里德算法，已经得到了6×(-2)+15×1=3，两边同时乘以3得6×(-6)+15×3=9，即x=-6，y=3时6x+15y=9。

注意
设a, b, c为任意整数，g=gcd(a,b)，方程ax+by=g的一组解是(x0,y0)，则当c是g的倍数时ax+by=c的一组解是(x0 * c/g, y0 * c/g)；当c不是g的倍数时无整数解。