POJ1061青蛙的约会

Description

两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止。可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置。不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的。但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的。为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面。 我们把这两只青蛙分别叫做青蛙A和青蛙B,并且规定纬度线上东经0度处为原点,由东往西为正方向,单位长度1米,这样我们就得到了一条首尾相接的数轴。设青蛙A的出发点坐标是x,青蛙B的出发点坐标是y。青蛙A一次能跳m米,青蛙B一次能跳n米,两只青蛙跳一次所花费的时间相同。纬度线总长L米。现在要你求出它们跳了几次以后才会碰面。

Input

输入只包括一行5个整数x,y,m,n,L,其中x≠y < 2000000000,0 < m、n < 2000000000,0 < L < 2100000000。

Output

输出碰面所需要的跳跃次数,如果永远不可能碰面则输出一行"Impossible"

这道题需要用到数论的同余思想进行求解,对于ax与b在mod L情况下恒等时,当且仅当gcd(a,L)能被b整除时才有解,所以不能整除时,直接输出Impossible

能整除时利用同余条件下的定理,和欧几里得扩展定理,来得到符合条件的x,而对于这一组解来说,所有与此x在mod L条件下恒等的都是符合条件的值,我们在里面找一个最小的正整数作为结果输出即可

题目中可轻易得到公式(m-n)*k=q-p(mod L) 这里的k即为上面所讲的所要求的x

代码如下:

 1 #include <iostream>
 2 #include <iostream>
 3 #include <iomanip>
 4 #include<string>
 5 #include<cstring>
 6 #define LL long long
 7 using namespace std;
 8
 9 LL a,b,x,y,d;
10
11 LL gcd(LL a,LL b)
12 {
13     if(b==0) return a;
14     else return gcd(b,a%b);
15 }
16 void ex_gcd(LL a,LL b,LL &x,LL &y,LL &d)
17 {
18     LL t;
19     if(b==0){
20         d=a,x=1,y=0;
21     }
22     else{
23         ex_gcd(b,a%b,x,y,d);
24         t=x,x=y,y=t-a/b*y;
25     }
26 }
27 int main()
28 {
29     LL p,q,m,n,L;
30     cin>>p>>q>>m>>n>>L;
31     a=m-n,b=q-p;
32
33     if(a<0) a=-a,b=-b;
34     LL Gcd=gcd(a,L);
35
36     if(b%Gcd!=0) cout<<"Impossible"<<endl;
37     else{
38         LL temp;
39         a=a/Gcd,temp=L/Gcd,b=b/Gcd;
40         ex_gcd(a,temp,x,y,d);
41         b=b*x;
42
43         b%=L;
44         if(b<0) cout<<b+L<<endl;
45         else cout<<b<<endl;
46     }
47
48
49     return 0;
50 }

POJ1061青蛙的约会,布布扣,bubuko.com

时间: 2024-10-13 02:03:14

POJ1061青蛙的约会的相关文章

POJ1061——青蛙的约会(扩展欧几里德)

青蛙的约会 Description两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面.

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

题目链接: 啊哈哈,点我点我 这道题是扩展欧几里得问题...哎,数学太弱了,看了半天才看懂.... 如果要相遇的话,则(n-m)*T+p*c=x-y成立,那么进行代换得到a*x+b*y=c,那么就转换成小白上面讲的了,所以用扩展欧几里得算法求得一组解,那么最后得到解的通式为x=x0+k*b/gcd(a,b),那么直接另右式子等于0及可..还有就是没有解的情况就是c%gcd(a,b)不等于0,那么就没有整数解...那么这个问题就得到了解决.... 题目: 青蛙的约会 Time Limit: 100

解题报告 之 POJ1061 青蛙的约会

解题报告 之 POJ1061 青蛙的约会 Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙

POJ-1061 青蛙的约会-数论扩展欧几里德算法入门及推导

Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面. 我们把这

POJ1061青蛙的约会[扩展欧几里得]

青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 108911   Accepted: 21866 Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总

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

题目链接:http://poj.org/problem?id=1061 青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 131879   Accepted: 29100 Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的

POJ1061 青蛙的约会 __一维世界的爱情

由于今天上午在做数论知识的笔记,发现那时候赵老师讲的线性丢番图(求ax+by=c的特解)部分完全搞不懂,后来网上查了一下才发现这个公式就是求同余方程,所用方法就是扩展欧几里得算法.正好红皮书上有这么一个模板,直接敲了下来然后稍作修改.后来发现POJ上的1061题就是这样一个类型,用了三个小时,荒废无数精力才AC.毕竟无聊,写个题解记录一下. 题目描述: 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是

[poj1061]青蛙的约会&lt;扩展欧几里得&gt;

题目链接:http://poj.org/problem?id=1061 其实欧几里得我一直都知道,只是扩展欧几里得有点蒙,所以写了一道扩展欧几里得裸题. 欧几里得算法就是辗转相除法,求两个数的最大公约数,算法是,a,b的最大公约数是gcd(b,a%b)然后不断递归下去,直到b=0 转换成c++语言就是 1 int ex_gcd(int a,int b) 2 { 3 if(b==0)return a; 4 return ex_gcd(b,a%b); 5 } 扩展欧几里得就是假设c=gcd(a,b)

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

本文出自:http://blog.csdn.net/svitter 题意:青蛙绕圈跳, 初始位置X,Y,速度M,N,方向相反,L为模.最后能否相遇?相遇时间是什么? 本题目为扩展欧几里德,扩展欧几里德介绍: 关于扩展欧几里德方程 ax + by = c (1) 可以用来求是否有解.即是否存在c满足这个方程. exgcd(a, b, x, y)是用来求ax + by = gcd(a, b)中x的值和y的值的.如果仅仅只是判断(1)是否有解,直接看gcd(a, b)能否整除c即可. 然后开始分析本题