POJ 1601 青蛙的约会

题意是两只青蛙,一开始分别在x,y点,他们分别每单位时间可以移动m,n个点位时间,路径是一个单位长度为l的圆环

通过题意很容易得到一个方程:

(n-m)*k+t*l=x-y;

其中k和t是未知数,(n-m)以及l我们可以通过输入确定

看到这个方程我们联想到ax+by=d;这个二元一次方程。d=gcd(a,b);

我们不妨另a=n-m b=l

则原式就为a*k+b*t=(x-y);。。。。1

已知1有整数解的条件为(d|(x-y)))

所以如果(x-y)%d!=0  那么必然无整数解,则输出impossible

另c= (x-y)

我们不妨给a/=d b/=d c/=d

利用欧几里得扩展 excgcd(a,b,x,y)

计算出x; 有欧几里得扩展的意义可以知道 x *=c;

然后计算最小值即可

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
long long d;//代表gcd(a,l)
long long gcd(long long  a,long long  b)
{
    if(b==0)
    {
        return a;
    }
    else{
        gcd(b,a%b);
    }
}

void  extgcd(long long a,long long b,long long &x,long long &y)
{

    if(b==0)
    {
        x=1;

        y=0;

    }
    else{
        extgcd(b,a%b,y,x);
        y -= (a/b)*x;
    }


}


long long x,y,m,n,l;
int main()
{
    while(cin>>x>>y>>m>>n>>l)
    {
        long long a = n-m;
        long long b = x-y;

        long long k,y;
        d = gcd(a,l);

        if(b%d!=0)
        {
            printf("Impossible\n");
            continue;
        }
        else{
            a/=d;
            l/=d;
            b/=d;
            extgcd(a,l,k,y);
            k *= b;
            long long t =k%l;//?????
            while(t<0)
            {
                t+=l;
            }
            cout<<t<<endl;
        }
    }
    return 0;

}



“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:大白 设计师:CSDN官方博客 返回首页
评论

打赏作者

UpCoderXH

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值