扩展欧几里得刷题合集

P1516 青蛙的约会

分析：

• $x+km\equiv y+kn(modl)$

  转换一下：$k(m-n)+lz=y-x$

• 然后直接套 $exgcd$ 求解即可

• 注意：$m-n<0$ 时，等号两边要同时取反

• 这题有助于理解取最小正整数解时的操作，方程 $ax+by=gcd$ 的解不唯一

  • 设 $x_0,y_0$ 是解，则 $x_0+t\ast b/gcd和y_0-t\ast a/gcd$ 也是解，$t\in N$ （要理解）
  • 在此题还有一个倍数 $k(b/d)$ , 同上：$k{x}_0+t\ast b/gcd$ 为方程解

#include <bits/stdc++.h>
#define int long long 
using namespace std;

void exgcd(int a,int b,int &x,int &y,int &d)
{
    if(!b) { x=1; y=0; d=a; return; } // 边界
    exgcd(b,a%b,x,y,d); // 当前层
    int t=y; 
    y=x-a/b*y; x=t; // 解出上一层的
}

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int sx,sy,m,n,l;
    cin>>sx>>sy>>m>>n>>l;
    int a=m-n, b=sy-sx;
    if(a<0) a=-a, b=-b;
    int x,y,d;
    exgcd(a,l,x,y,d);
    if(b%d!=0) { cout<<"Impossible"<<endl; } // 无法遇见
    else {
        int ans=(x*(b/d)%(l/d)+(l/d)-1)%(l/d)+1; // 取最小正整数解
        cout<<ans<<endl;
    }

    return 0;
}

Modular Inverse ZOJ (pintia.cn)

板子题

#include <bits/stdc++.h>
#define int long long 
using namespace std;

void exgcd(int a,int b,int &x,int &y,int &d)
{
    if(!b) { x=1; y=0; d=a; return; } // 边界
    exgcd(b,a%b,x,y,d); // 当前层
    int t=y; 
    y=x-a/b*y; x=t; // 解出上一层的
}

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int a,b,T;
    cin>>T;
    while(T--)
    {
        cin>>a>>b;
        int x,y,d;
        exgcd(a,b,x,y,d);
        if(d!=1) { cout<<"Not Exist"<<endl; continue; }
        x=(x%b+b-1)%b+1; // 最小正整数解，因为gcd=1，故b/gcd=b
        cout<<x<<endl;
    }
    return 0;
}

One Person Game ZOJ (pintia.cn)

分析：

• 求解 $ax+by=|B-A|$ ,

• 设 $x_0,y_0$ 是解，则通解为 $x_0+k\ast b/gcd和y_0-k\ast a/gcd$ ，$k\in N$

• 根据题意，要求的是 $|x|+|y|$ , 且 $x$ 和 $y$ 同号时，可以共有一部分（可以 $a+b=c$ ），此时 $ans=max(|x|,|y|)$

  要使 $ans$ 最小，即公共部分最大，最优解是 $x=y$ 的情况

• 若 $x=y$ ，即 $x_0+k\ast b/gcd=y_0-k\ast a/gcd$

  故，$k=(y0-x0)/(a/gcd+b/gcd)$ , 再偏移一下，遍历求出答案

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N=1e5+5;

void exgcd(ll a,ll b,ll &x,ll &y,ll &d)
{
    if(!b) { x=1; y=0; d=a; return; }
    exgcd(b,a%b,x,y,d);
    int t=y;
    y=x-a/b*y; x=t; 
}
void solve(ll a,ll b,ll c)
{
    ll x,y,d;
    exgcd(a,b,x,y,d);
    if(c%d!=0) 
    {
        cout<<"-1"<<endl; 
        return ;
    }
    x*=c/d; y*=c/d;
    a/=d; b/=d;
    ll mid=(y-x)/(a+b), ans=1e18;
    for(ll i=mid-1;i<=mid+1;i++)
    {
        ll xn=x+i*b, yn=y-i*a;
        if(xn*yn>0) // 同号
        {
            ans=min(ans,max(abs(xn),abs(yn)));
        }
        else ans=min(ans,abs(xn)+abs(yn)); // 异号
    }
    cout<<ans<<endl;
}

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int T;
    cin>>T;
    while(T--)
    {
        ll a,b,an,bn;
        cin>>a>>b>>an>>bn;
        ll c=abs(a-b);
        solve(an,bn,c);
    }

    return 0;
}

hdu 2669 Romantic

板子题

#include <bits/stdc++.h>
#define int long long 
using namespace std;

void exgcd(int a,int b,int &x,int &y,int &d)
{
    if(!b) { x=1; y=0; d=a; return; } // 边界
    exgcd(b,a%b,x,y,d); // 当前层
    int t=y; 
    y=x-a/b*y; x=t; // 解出上一层的
}

signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int a,b;
    while(cin>>a>>b)
    {
        int x,y,d;
        exgcd(a,b,x,y,d);
        if(d!=1) { cout<<"sorry"<<endl; continue; }
        x=(x%b+b)%b; // 最小正整数解，因为gcd=1，故b/gcd=b
        y=(1-x*a)/b;
        cout<<x<<' '<<y<<endl;
    }
    return 0;
}