2013寒假练习 1052 青蛙的约会

地址：http://acm.bit.edu.cn/mod/programming/view.php?a=538

扩展欧几里得算法解ax+by=c，经典题。注意最后要输出x的最小的正整数解，故求出通解应该先%(b/gcd),再把负的加成正的。直接加减会超时。

#include<iostream>#define ll __int64using namespace std;ll x0,y0;ll exgcd(ll a,ll b){if(b==0){x0=1,y0=0;return a;}ll gcd=exgcd(b,a%b);ll temp=x0;x0=y0,y0=temp-a/b*y0;return gcd;}int main(){ll x,y,m,n,l,a,b,c;while(~scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&l)){a=m-n,b=(-1)*l,c=y-x;ll gcd=exgcd(a,b);if(c%gcd!=0) printf("Impossible\n");else{x0*=c/gcd;x0%=b/gcd;  //不然TLEif(b/gcd>0)while(x0<0)  x0+=b/gcd;printf("%I64d\n",x0);}}return 0;}