扩展欧几里德算法

最近做了两道扩展欧几里德算法的题。
写一下自己对扩展欧几里德算法的理解。
扩展欧几里德算法用于求解ax+by = c
求ax+by = c， 令d =gcd(a,b);
那么(a / d ) * x + (b / d )* y = c / d
如果c不是d的倍数，则方程无解
如有解，则方程得到其中一个解 x0 = x*c/d
方程的最小整数解 公式 (x0%(b/d)+(b/d))%(b/d)

先上模板：

void gcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b) {d=a;x=1;y=0;}
    else {gcd(b,a%b,d,y,x);y-=x*(a/b);}
}

扩展欧几里德算法主要用于求解
1.乘法逆元 ax = 1(mod n) a关于模n的逆
2.模线性方程组 ax=b(mod n)

例如2.求解模线性方程组，可以转化为ax-b=ny–>ax-ny=b

R - 青蛙的约会 poj 1061

题意就不用写了。
此题是一到扩展欧几里德的经典题。
只要推出(x+mt)-(y+nt)=SL 这个方程
转化为(n-m)t-SL=x-y; 即可。
需要注意的是当n-m是负数的时候，要特殊处理一下，因为求出的gcd也是负数。如求解最小整数解的时候x=(x%b+b)%b; b一定要取绝对值

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long LL;

LL d,x,y;

void gcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b) {d=a;x=1;y=0;}
    else {gcd(b,a%b,d,y,x);y -= x*(a/b);}
}

int main()
{
    LL p,q,m,n,l;
    while(cin>>p>>q>>m>>n>>l)
    {
          LL a = m-n;
          LL b = l;
          LL c = q-p;
          gcd(a,b,d,x,y);
          if(c%d) {printf("Impossible\n");continue;}

           //cout<<d<<" "<<c<<" "<<x<<endl;

           x=x*c/d;
           b=b/d;
           if(b<0) b=-b;

           x=(x%b+b)%b;//最小整数解

           cout<<x<<endl;
    }
    return 0;
}

S - C Looooops

此题主要在于推出方程 cx-2^k = B-A;
并且求出最小整数解
注意的是求解2^k如果不用快速幂，就必须进行转换，LL tmp=(LL)1<

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long LL;

LL d,x,y;
void gcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b) {d=a;x=1;y=0;}
    else {gcd(b,a%b,d,y,x);y-=x*(a/b);}
}


int main()
{
    LL A,B,C,K;
    while(cin>>A>>B>>C>>K)
    {
        if(A==0&&B==0&&C==0&&K==0)
            break;

        LL tmp=(LL)1<<K;
        LL q = B-A;
        gcd(C,tmp,d,x,y);
        if(q%d) {cout<<"FOREVER"<<endl;continue;}
        x=x*q/d;
        tmp=tmp/d;
        cout<<(x%tmp+tmp)%tmp<<endl;

    }
    return 0;
}

HDU 2669 romantic

求解ax+by=1，此题其实是扩展欧几里德的水题
当gcd(a,b)!=1时，方程无解。

代码:

#include <iostream>
#include <cstdio>
using namespace std;

#define ll long long

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;
}
int main()
{
    ll a,b;
    while(cin>>a>>b)
    {
        ll x,y;
        int gcd =exgcd(a,b,x,y);
        if(gcd!=1) {cout<<"sorry"<<endl;continue;}
        while(x<0)
        {
            x+=b;
            y-=a;
        }
        cout<<x<<" "<<y<<endl;
    }
    return 0;
}