真的不是我懒啊啊啊,只是人家写的太好了。
讲的扩欧的应用,题不算难(虽然我还是不会)
一点小坑就是负数的替换,花姐写的很清楚了,为的是让kw+lz=s的各符号成立,还有就是最后%l/ans不是很好理解。
实际上是找的一个素数(即b/gcd(,)与x的关系),列出了同余方程
思路出不来时就÷gcd(a,b),同余方程转二元方程,二元方程转同余方程,拿两个右侧值相同的式子联立,然后替换替换再替换。
#include <iostream> #include <cstdio> #define ll long long using namespace std; ll x,y,m,n,l; ll w,x1,y1,s,ans; ll exgcd(ll a,ll b,ll &x1,ll &y1){ if(b==0){ x1=1,y1=0; return a; } ans=exgcd(b,a%b,x1,y1); ll tmp=x1;x1=y1;y1=tmp-a/b*y1; return ans; } int main(){ scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l); s=x-y,w=n-m; if(w<0) w=-w,s=-s; exgcd(w,l,x1,y1); if(s%ans!=0) printf("Impossible"); else printf("%lld",((x1*s/ans)%(l/ans)+l/ans)%(l/ans)); return 0; }