洛谷P1516 青蛙的约会 题解

题目内容

题目描述

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

我们把这两只青蛙分别叫做青蛙 A 和青蛙 B，并且规定纬度线上东经 $0$ 度处为原点，由东往西为正方向，单位长度 $1$ 米，这样我们就得到了一条首尾相接的数轴。设青蛙 A 的出发点坐标是 $x$，青蛙 B 的出发点坐标是 $y$。青蛙 A 一次能跳 $m$ 米，青蛙 B 一次能跳 $n$ 米，两只青蛙跳一次所花费的时间相同。纬度线总长 $L$ 米。现在要你求出它们跳了几次以后才会碰面。

输入格式

输入只包括一行五个整数 $x,y,m,n,L$。

输出格式

输出碰面所需要的次数，如果永远不可能碰面则输出一行一个字符串 Impossible。

样例 #1

样例输入 #1

1 2 3 4 5

样例输出 #1

4

提示

对于 $100\%$ 的数据，$1\le x\ne y\le 2\times 10^9$，$1\le m,n\le 2\times 10^9$，$1\le L\le 2.1\times 10^9$。

解题思路

本题实际上需要的是让我们求解一个同余方程：
$$x+km\equiv y+kn(modL)$$
我们想不出直接求解同余方程的解法，但是我们可以把同余方程转化为线性方程来求解：
$$x+km-y-kn=aL(a\in \mathbb{Z})$$
使用乘法分配律的逆定理得到：
$$x-y+k(m-n)=aL$$
移项之后得到：
$$x-y=aL+k(n-m)$$
这里为了更直观地观察，我们令 $S\leftarrow x-y,T\leftarrow n-m$ 。再交换一下等式的左右两边，于是就有：
$$aL+kT=S$$
其中 $L,T,S$ 为定值，满足求二元一次不定方程的算法要求。

题目要求的是求出满足条件的最小正整数 $k$ 的值，若不存在则输出 Impossible 。对有解或无解的判断，只需要使用裴蜀定理即可，即 $ax+by=z$ 有解当且仅当 $\gcd (a,b)|z$ 。

如果存在解，只需要使用解线性不定方程的模板求出最小解即可。如果不太了解的同学可以看一看这篇文章。

另外，如果 $T$ 为负数需要特殊处理。具体请看代码。

代码

代码很丑，请见谅。

#include<cstdio>
#define f_inline inline __attribute__((always_inline))
using namespace std;
using UL=unsigned long;//用于存储非负整数
using ULL=unsigned long long;//用于存储大型非负整数
using LL=long long;//用于存储大型带符号整数
f_inline char gc()//快速读入一个字符
{
#define MAXN 100000
    static char buf[MAXN],*l,*r;
    return (l==r)&&(r=(l=buf)+fread(buf,1,MAXN,stdin),l==r)?EOF:*l++;
#undef MAXN
}
f_inline UL read()//快速读入一个非负整数
{
    register char ch(gc());
    register UL x(0);
    while(ch<'0'||ch>'9')
        ch=gc();
    while(ch>='0'&&ch<='9')
        x=(x<<1)+(x<<3)+(ch^48),
        ch=gc();
    return x;
}
void print(register const ULL x)//快速输出一个ULL型整数
{
    if(x>9)
        print(x/10);
    putchar(x%10|48);
}
UL exgcd(register const UL a,register const UL b,register LL &x,register LL &y)//扩展欧几里得算法求解线性不定方程
{
    if(!b)
    {
        x=1,
        y=0;
        return a;
    }
    register const UL d(exgcd(b,a%b,y,x));
    y-=x*static_cast<LL>(a/b);
    return d;
}
int main()
{
    //freopen("../stdin","r",stdin);
    register LL k,a;
    register const UL x(read()),y(read()),m(read()),n(read()),L(read());
    register int S(x-y),T(n-m);//文章中描述的S和T
    if(T<0)//如果T小于零，那么等式两边同时乘以-1，得到(-a)L+k*(-T)=—S，此时求解出的a是原本a的相反数，但是不对答案产生影响，所以不用特殊处理
        S=-S,
        T=-T;
    register const UL d(exgcd(T,L,k,a));
    if(S%d)//如果S不是d的倍数，根据裴蜀定理，该方程无解
    {
        fputs("Impossible",stdout);
    }
    else//否则有解，输出最小解
    {
        register const UL t(L/d);//将aL+kT=gcd(L,T)的解转换为aL+kT=S的解
        k*=S/static_cast<int>(d),
        print((k%static_cast<LL>(t)+t)%t);//输出答案
    }
    return 0;
}