POJ 1061 青蛙的约会(扩展欧几里德)

原创 2016年05月31日 01:10:16

题目链接:
POJ 1061 青蛙的约会
题意:
规定纬度线上东经0度处为原点,由东往西为正方向,单位长度1米,这样就得到了一条首尾相接的数轴。
设青蛙A的出发点坐标是x,青蛙B的出发点坐标是y。青蛙A一次能跳m米,青蛙B一次能跳n米,两只青蛙跳一次所花费的时间相同且都是从东往西跳。纬度线总长L米。现在要你求出它们跳了几次以后才会碰面。
分析:
假设两只青蛙跳了k次相遇。则:(x + k * m) % L = (y + k * n) % L.
也就是:(x + k * m) - (y + k * n) = k’ * L,
移项即为:(m - n) * k - k’ * L = y - x。
相当于求解此不定方程的最小正整数解k。利用扩展欧几里德先计算(m - n) * k - k’ * L = gcd(m, n)的最小解k.
然后判断y - x是否为gcd(m - n, L)的倍数,如果不是倍数则无解。否则先将k扩大(y - x) / gcd(m - n, L)倍。
因为(m - n) * k - k’ * L = y - x的通解(kk, kk’)满足:
kk = (k + t * (L / d)), kk’ =(k’ - t * (m - n) / d),t为整数,d = gcd(m - n, L);
那么求k的最小正整数解kk:令r = L / d, kk = (k % r + r ) % r;
还要注意需要将系数变正!即m - n > 0!

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <climits>
#include <cmath>
#include <ctime>
#include <cassert>
#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;

ll ex_gcd(ll a, ll b, ll& xx, ll& yy)
{
    if(b == 0){
        xx = 1, yy = 0;
        return a;
    }
    ll d = ex_gcd(b, a % b, yy, xx);
    yy -= a / b * xx;
    return d;
}

int main()
{
    ll x, y, n, m, L, xx, yy;
    while(~scanf("%lld%lld%lld%lld%lld", &x, &y, &m, &n, &L)){

        if(m < n){
            swap(n, m);
            swap(x, y);
        }

        ll d = ex_gcd(m - n, L, xx, yy);
        if((y - x) % d != 0) printf("Impossible\n");
        else {
            xx *= ((y - x) / d);
            //通解满足:xx' = (xx + c * (L / d)), yy' =(yy - c * (m - n) / d),c为整数; 
            ll r = L / d;
            xx = (xx % r + r) % r; //最小正整数解
            printf("%lld\n", xx);
        }
    }
    return 0;
}
版权声明:缥缈玉京人,想语然、京兆眉妩。

POJ 1061 青蛙的约会

先说一下大概题意:有两只青蛙,一只在坐标x,另一直在坐标y,青蛙x一次跳跃可以前进m单位距离,青蛙y一次跳跃可以前进n单位的距离,两青蛙都在同一纬度,该纬度长度为L。两只青蛙同方向同时跳啊跳,问你最少...
  • ilovexiaohao
  • ilovexiaohao
  • 2013年08月06日 09:37
  • 1087

pku 1061青蛙的约会

扩展欧几里得 详细介绍 : http://hi.baidu.com/%B1%BF%D0%A1%BA%A2_shw/blog/item/0676025d56a87d4afbf2c093.html 扩...
  • u013849646
  • u013849646
  • 2014年03月25日 01:40
  • 719

POJ 青蛙的约会 (扩展欧几里得)

青蛙的约会 Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 20000/10000K (Java/Other) Total S...
  • qaz135135135
  • qaz135135135
  • 2016年08月02日 19:49
  • 275

北大POJ题目---青蛙的约会

今天在北大POJ上做到一道题——青蛙的约会,连续研究了数个小时,提交了共有50多次,每次都是以wrong anwser或是Time Limit Exceeded被打了回来。后来没吃晚饭,当看到最后一次...
  • xiefubao
  • xiefubao
  • 2014年01月16日 16:38
  • 1838

poj - 1061 - 青蛙的约会(扩展欧几里德)

题意:两只青蛙在同一条纬度线上,它们约定各自朝西跳,直到碰面为止。规定纬度线上东经0度处为原点,由东往西为正方向,单位长度1米。设青蛙A的出发点坐标是x,青蛙B的出发点坐标是y。青蛙A一次能跳m米,青...
  • SCNU_Jiechao
  • SCNU_Jiechao
  • 2013年02月11日 22:28
  • 729

POJ 1061 青蛙约会(扩展欧几里德)

#include #include //VJ上用printf scanf要加上这个库 using namespace std; long long gcd(long long a,long lo...
  • mirror58229
  • mirror58229
  • 2016年07月24日 21:28
  • 253

【POJ 1061 青蛙的约会】 + 扩展欧几里德

青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 108769 Accepted: 21821...
  • WYK1823376647
  • WYK1823376647
  • 2016年10月15日 16:52
  • 196

POJ 1061 青蛙的约会 欧几里德扩展

原题: http://poj.org/problem?id=1061对于这个题,我们可以看成一个追击问题,假设其中一个人是不动的,那么另一个相对它的速度就是n-m,而最初的路径差是x-y,我们要求的...
  • qq_27508477
  • qq_27508477
  • 2015年08月14日 21:01
  • 266

poj 1061 青蛙的约会 -扩展欧几里德算法

设过s步后两青蛙相遇,则必满足以下等式:     (x+m*s)-(y+n*s)=k*l(k=0,1,2....)   稍微变一下形得:     (n-m)*s+k*l=x-y     令n-m=a,...
  • viphong
  • viphong
  • 2015年10月26日 17:56
  • 259

Poj 1061青蛙的约会 扩展欧几里德

青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 98003 Accepted: 18546...
  • FTQOOO
  • FTQOOO
  • 2015年08月18日 21:09
  • 279
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 1061 青蛙的约会(扩展欧几里德)
举报原因:
原因补充:

(最多只允许输入30个字)