poj 1061 青蛙的约会扩展gcd

最新推荐文章于 2020-03-02 18:06:33 发布

Jeromiewn

最新推荐文章于 2020-03-02 18:06:33 发布

阅读量334

点赞数

分类专栏： gcd 经典题模板题文章标签： acm hdu 模板算法 gcd

本文链接：https://blog.csdn.net/Jeromiewn/article/details/54695489

版权

模板题同时被 3 个专栏收录

42 篇文章 0 订阅

订阅专栏

经典题

24 篇文章 0 订阅

订阅专栏

gcd

6 篇文章 0 订阅

订阅专栏

青蛙相遇可列方程

(x+mt)%L=(y+nt)%L，即(m-n)*t+k*L=y-x

可转化为ax+by=c

先求出ax+by=gcd，再将解得的x*（c/gcd）

求出最小整数解即可:

用x=(x*b+b)%b....(这里的b为b/gcd

#include <iostream>
#define LL long long
using namespace std;
LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    LL r=exgcd(b,a%b,x,y);
    LL t=x;
    x=y;
    y=t-a/b*y;
    return r;
}


int main()
{
    LL x,y,m,n,L;
    cin>>x>>y>>m>>n>>L;
    LL a,b,c,gcd;
    a=n-m;
    b=L;
    c=x-y;
    if(a<0)
    {
        a=-a;
        c=-c;
    }
    gcd=exgcd(a,b,x,y);
    if(m==n||c%gcd)
    {
        cout<<"Impossible"<<endl;
    }
    else
    {
        c/=gcd;
        b/=gcd;
        long long ans=x*c;//x是ax+by=gcd的解，ax+by=c的解只需x*(c/gcd)
//        long long ans=t%b;
//        if(ans<=0)
//            ans+=b;
           ans=(ans%b+b)%b;//使解>0,且最小
        cout<<ans<<endl;


    }
    return 0;
}