(扩展欧几里得:存在 x,y 使得 gcd(a,b)=ax+by )
int exgcd(int a, int b, int &x, int &y)
//int函数中设置x,y=0,可求得a,b公约数以及一组x,y
//(求出来的x,y可能为负,对其做以下操作可使其非负:(x%9973+9973)%9973
{
if (b == 0)
{
x = 1;
y = 0;
return a;
}
int gcd = exgcd(b, a%b, x, y);//调用自身
int x2 = x, y2 = y;
x = y2;
y = x2 - (a / b)*y2;
return gcd;//返回公约数