poj 1060---青蛙的约会(加深对扩展GCD的理解,然而是一道水题)

连接

题意:
问最少要跳跃几次两只青蛙才能见面.
思路:
根据题意,两个青蛙跳到同一个点上才算是遇到了,所以有 (x+mt) - (y+nt) = p * k; (t是跳的次数,k是a青蛙跳的圈数跟b青蛙的圈数之差。整个就是路程差等于纬度线周长的整数倍),
转化一下: (n-m) * t + k * p = x – y;
令 a = n-m, b = k, c = gcd(a, b), d = x-y;
有 a * t + b * p = d; (1)
要求的是t的最小整数解。
用扩展的欧几里德求出其中一组解t0 ,p0, 并令c = gcd(a, b);
有 a * t0 + b * p0 = c; (2)
因为c = gcd(a, b), 所以 a * t / c是整数,b * t / c 也是整数,所以 d / c 也需要是整数,否则无解。
(2)式两边都乘(d / c) 得 a * t0 (d / c) + b * p0 * (d / c) = d;
所以t0 * (d / c)是最小的解,但有可能是负数。
因为a * ( t0 (d / c) + bn) + b * (p0 * (d / c) – a
n) = d; (n是自然数)
所以解为 (t0 * (d / c) % b + b) % b;

所以根据思路,有以下代码

#include<cstdio>
#include<cstring>
#include<iostream>
#define ll long long
#define db double
using namespace std;
ll x,y,c;

ll exgcd(int a,int b)
{
    if(!b){
        x=1;
        y=0;
        c=a;
    }
    else{
        exgcd(b,a%b);
        ll tmp=x;
        x=y;
        y=tmp-a/b*y;
    }
}

int main()
{
    ll a,b,d,x1,y1,m,n,l,flag=0;
    cin >> x1 >>y1>>m>>n>>l;
    if(m==n)
        flag=1;
    else{
        a=n-m;
        b=l;
        d=x1-y1;
        exgcd(a,b);
        if(d%c)//根据定义,如果除不尽就无解.
            flag=1;
    }
    if(flag)
        cout << "Impossible" << endl;
    else{
        b/=c;
        d/=c;
        cout << ((x*d)%b+b)%b <<endl;
    }
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 【1】7,9,-1,5,( ) A、4;B、2;C、-1;D、-3 分析:选D,7+9=16;9+(-1)=8;(...
    Alex_bingo阅读 19,081评论 1 19
  • 有些人会说再见,有些人却只能说拜拜, 多么无奈 我所做的一切只是为了让你成为更好的自己 17的我比18岁的我快乐...
    Lynn琳恩阅读 238评论 0 0
  • #DAY 56#高级框架模型之黄金思维圈-03 永澄:高级框架模型之黄金思维圈-03 1、多维度使用该模型的方法 ...
    威威单阅读 338评论 0 0
  • 看着滚圆的肚子越来越大了,有点想念在珠海跑步的日子,怎么也没想到跑步为数不多可以坚持下来的事,要不是术后有恢复期,...
    夜初微凉阅读 264评论 0 0
  • 最近身边有几个朋友都有前往印度的打算,他们问我意见或者有什么地方值得推荐。我想了想,以我两个月在印度的晃荡经历,学...
    跟着鹏叔游世界阅读 635评论 0 1