青蛙的约会【同余】【扩欧】

>Link

luogu P1516

---

>Description

求关于$t$的同余方程 $x+tn\equiv y+tm(modL)$

---

>解题思路

把常数和未知数分开：$tn-tm\equiv y-x(modL)$
然后同余转换成等式：
$tn-tm+sL=y-x$
$t(n-m)+sL=y-x$
$t,s$都是未知数，现在就是如 $ax+by=c$ 的式子

为使得 $a,b$ 互质（就可以用扩欧求了），$a,b$ 同除以 $gcd(a,b)$，那等式右边的 $c$ 也需要除以 $(a,b)$。如果 $c$ 不能整除 $(a,b)$，就说明无解了。
设 $g=gcd(a,b)$，现在求 $\frac{a}{g}x+\frac{b}{g}y=\frac{c}{g}$。
$\frac{a}{g},\frac{b}{g}$ 已互质，用扩展欧几里得可以求出 $\frac{a}{g}x+\frac{b}{g}y=1$ 的解，两边再同乘 $\frac{c}{g}$ 就可以得到题目要求的 $t$ 了。

要注意负数，要把它变成非负数。

---

>代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL long long
using namespace std;

LL gcd (LL a, LL b)
{
	if (!b) return a;
	return gcd (b, a % b);
}
void ex_gcd (LL a, LL b, LL &x, LL &y)
{
	if (!b) {x = 1; y = 0; return;}
	LL X = x, Y = y;
	ex_gcd (b, a % b, X, Y);
	x = Y;
	y = X - (a / b) * Y;
}

int main()
{
	LL p, q, n, m, a, b, c, x, y, g;
	scanf ("%lld%lld%lld%lld%lld", &p, &q, &n, &m, &b);
	a = ((n - m) % b + b) % b;
	c = ((q - p) % b + b) % b;
	g = gcd (a, b);
	if (c % g != 0)
	{
		printf ("Impossible");
		return 0;
	}
	a /= g, b /= g, c /= g;
	ex_gcd (a, b, x, y);
	x = (x % b + b) % b;
	x = x * c % b;
	printf ("%lld", x);
	return 0;
}