扩展GCD

扩展gcd：已知ax+by=gcd(a,b)，求x,y。

①由gcd(a,b)==gcd(b,a%b)得ax+by==bx+(a%b)y，继续化下去，化成ax1+by1的形式即为递归转移式，因为递归过程中a,b不能改变，改的是x,y。

bx+(a%b)y == bx+ (a-(a/b)*b) y == bx+ay-(a/b)*b*y == ay + b(x-a/b*y)

即exgcd(ax+by)== exgcd( ay + b(x-a/b*y)  )

然后就可以递归求解了，共用相同的x,y，可以全局，下面代码是用引用实现。

②递归终止：当b==0时，式子为ax==gcd(a,b)，那么显然a=gcd(a,b),x=1,y=0;

 //扩展gcd求解x,y同时还把gcd(a,b)求出来了，exgcd的返回值就是gcd(a,b)

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

扩展gcd求出一个解后，那ax+by=gcd(a,b)的其他解为x=x0+k*b; y=y0-k*a;
那么ax+by=c就不难了，将ax+by=gcd(a,b)两边除以gcd可得，那么代码中就正常求exgcd，解x=c/gcd*x0,y=c/gcd*y0。(x0,y0为exgcd的解)

这里注意：x=c/gcd*x0 与 x=x0/gcd*c 不等价！有解情况下gcd一定是c的因子,而x0不一定是。

如8*3/2=12而8*(3/2)!=12。  当然x也可以=c*x0/gcd，不过要小心溢出。

 

ax+by=c无整数解的情况：c不能被gcd(a,b)整除（也即c不是gcd(a,b)的倍数）。

理由：因为ax与by均是gcd(a,b)的倍数,所以其和ax+by也是gcd(a,b)的倍数。

 

 

乘法逆元

满足a*x≡1 (mod p)的x值

通俗理解：(1/a)mod p 和 (a的逆元)mod p一样

作用：求(a/b)%p，当a过大，导致a/b无法直接求得时可以先求b对p的逆元x，再乘a去mod p，即

(a*x)%p和(a/b)%p结果是一样的

逆元求法

ll inverse(ll num,ll mod)  {  
    ll x,y;  
    exgcd(num,mod,x,y);  
    while(x<0) x+=mod,y-=num;  
    return x;  
} 

HDU 2699 Romantic

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll a, b, x, y;
ll exgcd(ll a, ll b) {
    if(!b) {
        x = 1, y = 0;
        return a;
    }
    ll d = exgcd(b, a % b);
    ll t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
}
int main() {
    while(cin >> a >> b) {
        if(exgcd(a, b) != 1)
            puts("sorry");
        else {
            while(x < 0)
                x += b, y -= a;
            cout << x << " " << y << endl;
        }
    }
    return 0;
}

HDU 1576 A/B

思路：n=A%9973=a-a/9973*9973

设x=A/B,则A=Bx，Bx-A/9973*9973=n，则Bx+9973y=n

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int t;
ll n, b, x, y;
ll exgcd(ll a, ll b) {
    if(!b) {
        x = 1, y = 0;
        return a;
    }
    ll d = exgcd(b, a % b);
    ll t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
}
int main() {
    cin >> t;
    while(t--) {
        cin >> n >> b;
        exgcd(b, 9973);
        x = n * x;
        x = (x % 9973 + 9973) % 9973;防止负数
        cout << x << endl;
    }
    return 0;
}