定理:对于方程\(ax+by=c\),等价于\(a*x=c(mod b)\),有整数解的充分必要条件是c是gcd(a,b)的整数倍。
——信息学奥赛之数学一本通
避免侵权。哈哈。
两只青蛙跳到一格才行,所以说
\(x+mt=y+nt(mod l) \)
\((x-y)+(m-n)t=0(mod l)\)
\((m-n)t+ls=(y-x) s属于整数集\)
令a=n-m,b=l,c=gcd(a,b),d=x-y
则有\( at+bs=d\)
扩展欧几里得求解。
设c=gcd(a,b),若d/c 不是整数则无解。
最小解=(c*(d/c)%b+b)%b
这里是扩展欧几里得的代码。
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
long long exgcd(long long a,long long b,long long &x,long long &y){ if(b==0){ y=0;x=1; return a; } long long ret=exgcd(b,a%b,x,y); long long tmp=x;x=y;y=tmp-a/b*y; return ret; }
附上解题代码
#include<iostream> #include<cstdio> using namespace std; long long p,q; long long exgcd(long long a,long long b,long long &x,long long &y){ if(b==0){ y=0;x=1; return a; } long long ret=exgcd(b,a%b,x,y); long long tmp=x;x=y;y=tmp-a/b*y; return ret; } int main(){ long long n,m,x,y,l; cin>>x>>y>>m>>n>>l; long long a=x-y,b=n-m; if(b<0){ b=-b;a=-a; } long long c; if((a)%(c=exgcd(b,l,p,q))){ printf("Impossible"); return 0; } printf("%lld",((a)/c*p%(l/c)+(l/c))%(l/c)); return 0; }