扩展欧几里得+求解逆元

在为扩展前，什么是欧几里得算法？？

欧几里得又称辗转相除法，用于计算两个整数a,b的最大公约数（最大公约数：指两个或多个整数共有约数中最大的一个）其计算原理依赖于下面的定理：

gcd函数就是用来求(a,b)的最大公约数的。

gcd函数的基本性质：

gcd(a,b)=gcd(b,a)=gcd(-a,b)=gcd(|a|,|b|)

公式表述

gcd(a,b)=gcd(b,a mod b)

证明：a可以表示成a = kb + r，则r = a mod b

假设d是a,b的一个公约数，则有

d|a, d|b，而r = a - kb，因此d|r

因此d是(b,a mod b)的公约数

假设d 是(b,a mod b)的公约数，则

d | b , d |r ，但是a = kb +r

因此d也是(a,b)的公约数

因此(a,b)和(b,a mod b)的公约数是一样的，其最大公约数也必然相等，得证

代码实现：

int gcd(int a,int b)
{
    int c;
    while(b!=0)
    {
        c=a%b;
        a=b;
        b=c;
    }
    return a;
}

扩展欧几里得：对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，必然存在整数对 x，y ，使得 gcd（a，b）=ax+by。

验证：我们知道gcd(a,b)=gcd(b,a%b)；所以我们可以得到：bx1+a%by1=ax+by;

化简后可以到bx1+(a-a/b*b)y1=ay1+b(x1-a/by1)=ax+by;

这里我们令x=y1;所以可以得到解：y=x1-a/by1;

对于关于x，y的方程ax+by=g来说，让x增加b/g，让y减少a/g，等式两边还相等

所以我们可以计算出所有的解，但是b到最后可能会为0，所以到b=0的时候停止求解

那么对于方程a * x + b * y = c 来说，它的解怎么求呢？ 
答：如果c % gcd(a,b) != 0，即c不是gcd的整数倍，则无解。 
如果c % gcd(a,b) == 0 且 c / gcd(a,b) = t，那么求出方程 a * x + b * y = gcd(a,b)的所有解x,y，将x,y乘上t，对应的x’,y’即是方程a * x + b * y = t * gcd(a,b)的解 

直接说吧,假设d=gcd(a,b). 那么x=x0+b/d*t; y=y0-a/d*t;其中t为任意常整数。

推导过程：由题意可得a*x0+ab/d+b*y0-ab/d=d;

得到:a(x0+b/d)+b(y0-a/d)=d;

因为这是不等式所以得到：x=x0+b/d*t,y=y0-a/d*t,t为常数

大概就是这样

代码：

int gcd(int a,int b,int &x,int &y)
{
    if(a==0&&b==0)
        return -1;//无最大公约数
    if(b==0)
    {
        x=1;
        y=0;
        return a;//返回最大公约数
    }
    int d=gcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}

可能有的人对于为什么x，y要加取地址符不是很清楚，这里给出一个代码对比应该就能理解了

样例1：没有取地址符的

样例2：有取地址符

那么扩展欧几里得有什么用呢？？？可以用它求如下的问题：

形如ax+by=c这样的数或者让你求逆元，逆元就是ax同余（1modb）这里面x的最小正整数即为逆元

a同余（1%b）就等价于ax+by=1解这个方程即可

注意这里的gcd=1，也就是说当两数互质的时候存在解，然后还需要注意这里b可能为非正整数所以需要取正。

代码：

ll cal(ll a,ll b)///求逆元
{
    ll x,y;
    ll gcd=exgcd(a,b,x,y);
    if(1%gcd)
        return -1;///无解的情况
    x*=1/gcd;
    b=abs(b);
    ll ans=x%b;
    if(ans<=0)
        ans+=b;
    return ans;
}