洛谷 P1516 青蛙的约会 题解

题目传送门

题目大意： 给一个长度为 $L$ 的圆环，有两只青蛙，给出两者的起点和速度，问何时相遇。

题解

显然是个追击问题，先列出柿子：
$$x+tn\equiv y+tm(\mathrm{m}\mathrm{o}\mathrm{d}L)$$

其中 $t$ 即为我们要求的答案，其他的变量意义如题所述。

那么我们大力变换：
$$\begin{array}{rl}x+tn& \equiv y+tm(\mathrm{m}\mathrm{o}\mathrm{d}L)\\ x-y+t(n-m)& \equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}L)\end{array}$$

不妨设 $A=x-y,B=m-n$，那么有：
$$A-tB\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}L)$$

这个$(\mathrm{m}\mathrm{o}\mathrm{d}L)$ 有点棘手，不妨设多一个 $kL$ 带进去，有
$$\begin{array}{rl}A-tB& =kL\\ kL+tB& =A\end{array}$$

根据裴蜀定理，如果 $A$ 不是 $(L,B)$ 的倍数，那么该方程无解，此时输出 ${}^{\prime \prime }Impossibl{e}^{\prime \prime }$。

如果有解，还要注意，$B$ 不能是负数，所以如果 $B$ 是负数的话就要将 $A,B$ 同时取反，这是不影响答案的。（回到 $A-tB\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}L)$ 这个柿子就能明白）

那么我们可以用 $exgcd$ 求出下面这个方程的一组合法解：
$$kL+tB=\gcd (L,B)$$

对于求出来的一组合法的解 $k_1,t_1$，并不一定是最后的答案，我们要的是最小的 $t_1$，所以我们还得搞一搞。

设 $g=\gcd (L,B)$，我们发现，让 $t_1$ 加 $\frac{L}{g}$ 再让 $k_1$ 减 $\frac{B}{g}$，这样又能得到一组另外的合法的解，那么我们不妨让 $t_1$ 对 $\frac{L}{g}$ 取模，那么就能得到最小解了。

但是我们求出的是 $kL+tB=\gcd (L,B)$ 这个方程的最小解，要让整个方程乘以 $\frac{A}{\gcd (L,B)}$，才能得到原方程的最小解。（注意这个解也需要对 $\frac{L}{g}$ 取模）

于是代码如下：

#include <cstdio>
#include <cstring>
#define ll long long

ll x,y,n,m,L;
ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
void exgcd(ll a,ll b,ll &xx,ll &yy)
{
	if(b==0){xx=1;yy=0;return;}
	exgcd(b,a%b,yy,xx);
	yy-=a/b*xx;
}

int main()
{
	scanf("%lld %lld %lld %lld %lld",&x,&y,&n,&m,&L);
	ll A=x-y,B=m-n,k,t,g;
	if(B<0)A=-A,B=-B;
	g=gcd(L,B);//注意，要在B判负之后再求gcd
	if(A%g!=0)return printf("Impossible"),0;
	exgcd(L,B,k,t);
	t*=A/g;
	t=(t%(L/g)+L/g)%(L/g);//注意负数的情况
	printf("%lld",t);
}