poj 1061 青蛙的约会 扩展gcd

24 篇文章 0 订阅
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;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值