[Acwing] 878. 线性同余方程

题意

求出一个$x_i$满足$a_i\ast x_i\equiv b_i(mo{m}_i)$

思路

($m|x$表示$m$是$x$的因子,或者$x$能被$m$整除)

首先什么是同余 ?

设整数$m\ne 0$

若$m|(a-b)$,称$m$位模数,$a$同余$b$模$m$,$b$是$a$对模$m$的剩余记作$a\equiv b(modm)$

则我们可以知道 $a-b$是$m$的倍数,因此上面的式子可以化为
$$\begin{array}{rl}& a\ast x+y\ast m=b\end{array}$$

但是显然这里并不能使用$扩展欧几里得$

因为扩展欧几里得需要满足
$$\begin{array}{rl}& a\ast x+y\ast m=gcd(a,m)\end{array}$$

那么上面那个定理什么时候有解呢 ?

    a x + b y = c 设 d = g c d ( a , b )         d ∣ a , d ∣ b     d ∣ ( a x ) , d ∣ ( b y ) c = m ∗ d = m ∗ g c d ( a , b ) \begin{aligned} &\ \ \ ax+by=c\\ &设d=gcd(a,b)\\ &\ \ \ \ \ \ \ d|a,d|b\\ &\ \ \ d|(ax),d|(by)\\ &c=m*d=m*gcd(a,b) \end{aligned} ​   ax+by=c设d=gcd(a,b)       d∣a,d∣b   d∣(ax),d∣(by)c=m∗d=m∗gcd(a,b)​
则显然$c$是$gcd(a,b)$的时候有解,这就是贝祖定理

因此我们又回去考虑
$$\begin{array}{rl}& a\ast x+y\ast m=b\end{array}$$

显然当$gcd(a,m)|b$的时候有解

因此我们不妨先用扩展欧几里得求解出
$$\begin{array}{rl}& a\ast x_0+y_0\ast m=gcd(a,m)\end{array}$$

然后我们将式子左右两边同时乘上$\frac{b}{gcd(a,m)}$
$$\begin{array}{rl}& (a\ast x_0+y_0\ast m)\ast \frac{b}{gcd(a,m)}=gcd(a,m)\ast \frac{b}{gcd(a,m)}\\ & a\ast x_0\ast \frac{b}{gcd(a,m)}+y_0\ast \frac{b}{gcd(a,m)}=b\end{array}$$

因此最后对于上面由扩展欧几里得求出的$x_0$ 可以用来推导出所求式的$x$

$x=x_0\ast gcd(a,m)/b\%m$

Code :

#include <iostream>
#include <algorithm>

using namespace std;

typedef long long LL;


int exgcd(int a, int b, int &x, int &y)
{
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}


int main()
{
    int n;
    scanf("%d", &n);
    while (n -- )
    {
        int a, b, m;
        scanf("%d%d%d", &a, &b, &m);

        int x, y;
        int d = exgcd(a, m, x, y);
        if (b % d) puts("impossible");
        else printf("%d\n", (LL)b / d * x % m);
    }

    return 0;
}