扩展欧几里德

对扩展欧几里德的理解，首先我们需要知道欧几里得，

扩展欧几里德的出现来自于a*x+b*y=gcd(a,b) ，这个是一定会有解的，而扩展欧几里德就是求出一组x，y。

而我们要做的一般是a*x+b*y=c；

令r=gcd(a,b)；为了求出一组x,y；a*x+b*y=r;   同时乘以c/r;    c/r*a*x+c/r*b*y=r*c/r;得到a*(c/r*x)+b*(c/r*y)=c;

得到了x,y；那也就可以得到c/r*x;

怎么求x,y；扩张欧几里德算法的返回条件是b==0；当b==0时返回a就是最大公约数,而这时a*x+b*y==gcd(a,b)=r;

假设a*x+b*y=gcd(a,b)     b*x1+(a%b)*y1=gcd(a,b)   a*x+b*y=b*x1+(a-a/b*b)y1;

a*x+b*y=b*x1+a*y1-a/b*b*y1

a*x+b*y=a*y1+b*(x1-a/b*y1)     x=y1,  y=x1-a/b*y1;当到递归函数的返回条件时，就知道了x和y的值，就能一层一层的返回

当满足返回条件时b==0；a*x+b*y=gcd(a,b);b=0  gcd(a,b)=a; x=1，y=0;

注意:

         当我们使用exgcd()的时候，我们需要注意，求的是两个正整数，当a*x+b*y=c时，把符号转移到c，x，y上面去

ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
      {
          x=1;
          y=0;
          return a;
      }
    ll r=exgcd(b,a%b,x,y);
    ll tem=y;
    y=x-a/b*y;
    x=tem;
    return r;
}

菜鸟的理解

补充，平时我们用扩展欧几里德我们经常需要求最小正整数解

这是 当x变化时，y也要变化来抵消x的变化，x的加减b/gcd(a,b) y的就是 a/gcd(a,b)

令b/gcd(a,b)=t;   x的最小正整数解  (x%t+t)%t