P1516 青蛙的约会 （扩展欧几里得，理解如何取解）

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;
}