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;
}

 

posted @ 2017-04-11 20:45  IIIIIIIIIU  阅读(247)  评论(0编辑  收藏  举报