POJ 1061 青蛙的约会 | 同余方程和exGcd

 1 #include<cstdio>
 2 typedef long long ll;
 3 using namespace std;
 4 ll s,t,p,q,L;
 5 ll gcd (ll x,ll y)
 6 {
 7     return y==0?x:gcd(y,x%y);
 8 }
 9 ll exGcd(ll a,ll b,ll &x,ll &y)
10 {
11     if (b==0) return x=1,y=0,a;
12     ll r=exGcd(b,a%b,y,x);
13     y-=(a/b)*x;
14     return  r;
15 }
16 int main()
17 {
18     scanf("%lld%lld%lld%lld%lld",&s,&t,&p,&q,&L);
19     ll a=(p-q+L)%L,b=L,c=(t-s+L)%L,x,y,G=gcd(a,b);
20     if (c%G!=0)
21     {
22     puts("Impossible");
23     return 0;
24     }
25     a/=G,b/=G,c/=G;
26     exGcd(a,b,x,y);
27     x=(x%b+b)%b;
28     x=x*c%b;
29     printf("%lld\n",x);
30     return 0;
31 }