EXGCD - 青蛙的约会 - POJ 1061
两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。
它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止。
可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置。
不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的。
但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的。
为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面。
我们把这两只青蛙分别叫做青蛙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 。 输入只包括一行5个整数x,y,m,n,L,其中x≠y < 2000000000,0 < m、n < 2000000000,0 < L < 2100000000。 输入只包括一行5个整数x,y,m,n,L,其中x=y<2000000000,0<m、n<2000000000,0<L<2100000000。
Output
输 出 碰 面 所 需 要 的 跳 跃 次 数 , 如 果 永 远 不 可 能 碰 面 则 输 出 一 行 ′ ′ I m p o s s i b l e ′ ′ 输出碰面所需要的跳跃次数,如果永远不可能碰面则输出一行''Impossible'' 输出碰面所需要的跳跃次数,如果永远不可能碰面则输出一行′′Impossible′′
Sample Input
1 2 3 4 5
Sample Output
4
分析:
假 设 A 起 点 在 a , B 起 点 在 b 。 A 跳 一 次 m 米 , B 跳 一 次 n 米 , 当 A 、 B 相 遇 时 , 有 方 程 : 假设A起点在a,B起点在b。A跳一次m米,B跳一次n米,当A、B相遇时,有方程: 假设A起点在a,B起点在b。A跳一次m米,B跳一次n米,当A、B相遇时,有方程:
( m − n ) x = b − a + y L , 其 中 x 、 y ∈ Z , x 表 示 跳 过 的 次 数 , (m-n)x=b-a+yL,其中x、y∈Z,x表示跳过的次数, (m−n)x=b−a+yL,其中x、y∈Z,x表示跳过的次数,
求 A 、 B 何 时 相 遇 , 即 求 x 的 最 小 正 整 数 解 。 求A、B何时相遇,即求x的最小正整数解。 求A、B何时相遇,即求x的最小正整数解。
方 程 形 如 : a x + b y = C , 利 用 扩 展 欧 几 里 得 算 法 , 解 方 程 a x + b y = g c d ( a , b ) 方程形如:ax+by=C,利用扩展欧几里得算法,解方程ax+by=gcd(a,b) 方程形如:ax+by=C,利用扩展欧几里得算法,解方程ax+by=gcd(a,b)
本 题 中 , a = ( m − n ) , b = L , C = b − a 本题中,a=(m-n),b=L,C=b-a 本题中,a=(m−n),b=L,C=b−a
若 g c d ( a , b ) ∣ C , 则 有 解 , 否 则 无 解 。 若gcd(a,b)\ |\ C,则有解,否则无解。 若gcd(a,b) ∣ C,则有解,否则无解。
解 得 ( x 0 , y 0 ) 满 足 a x 0 + b y 0 = g c d ( a , b ) , 等 式 两 边 扩 大 C g c d ( a , b ) 倍 , 得 到 满 足 a x 0 ′ + b y 0 ′ = C 的 解 , 解得(x_0,y_0)满足ax_0+by_0=gcd(a,b),等式两边扩大\frac{C}{gcd(a,b)}倍,得到满足ax_0'+by_0'=C的解, 解得(x0,y0)满足ax0+by0=gcd(a,b),等式两边扩大gcd(a,b)C倍,得到满足ax0′+by0′=C的解,
x 0 ′ = C g c d ( a , b ) x 0 , 通 解 形 式 : x = x 0 ′ + k b g c d ( a , b ) x_0'=\frac{C}{gcd(a,b)}x_0,通解形式:x=x_0'+k\frac{b}{gcd(a,b)} x0′=gcd(a,b)Cx0,通解形式:x=x0′+kgcd(a,b)b
要 求 x 的 最 小 正 整 数 解 , 即 将 x 0 ′ 对 b g c d ( a , b ) 取 模 , 求 最 小 正 余 数 。 要求x的最小正整数解,即将x_0'对\frac{b}{gcd(a,b)}取模,求最小正余数。 要求x的最小正整数解,即将x0′对gcd(a,b)b取模,求最小正余数。
注意:
l o n g l o n g 类 型 求 绝 对 值 函 数 为 : long\ long类型求绝对值函数为: long long类型求绝对值函数为:
#include<algorithm>
labs(long long x);
代码:
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
ll n,m,a,b,L; // (m-n)x=b-a+yL => (m-n)x-Ly=b-a
ll exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1,y=0;
return a;
}
ll d=exgcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
int main()
{
cin>>a>>b>>m>>n>>L;
ll x,y;
ll d=exgcd(m-n,L,x,y); // x=x0+k*(L/d)
if((b-a)%d) puts("Impossible");
else
{
ll t=labs(L/d);
x*=(b-a)/d;
x=((x%t)+t)%t;
cout<<x<<endl;
}
return 0;
}