首页 > 代码库 > AC日记——青蛙的约会 poj 1061

AC日记——青蛙的约会 poj 1061

青蛙的约会

 POJ - 1061 
 
思路:
  扩展欧几里得;
  设青蛙们要跳k步,我们可以得出式子
    m*k+a≡n*k+b(mod l)
   式子变形得到
      m*k+a-n*k-b=t*l
      (m-n)*k-t*l=b-a
    然后,exgcd函数求出k
    然后输出刚刚大于0的k即可
 
来,上代码:
#include <cmath>#include <cstdio>#include <cstring>#include <iostream>using namespace std;long long n,m,a1,b1,l,xx,yy;long long exgcd(long long a,long long b,long long &x,long long &y){    if(b==0)    {        x=1,y=0;        return a;    }    long long r=exgcd(b,a%b,x,y),tmp=x;    x=y,y=tmp-a/b*y;    return r;}int main(){    scanf("%lld%lld%lld%lld%lld",&a1,&b1,&m,&n,&l);    long long a=m-n,b=l,c=b1-a1;    if(a<0) a=-a,c=-c;    long long r=exgcd(a,b,xx,yy);    if(c/r*r==c)    {        xx=xx*c/r,l=l/r;        if(xx>0) xx=xx%l;        else xx=xx%l+l;        cout<<xx;    }    else cout<<"Impossible";    return 0;}

 

AC日记——青蛙的约会 poj 1061