同余问题共7part,我的博客链接:
一元线性同余方程
形式: a x ≡ b ( m o d m ) ax\equiv b(\mod m) ax≡b(modm) ,求解 x x x ,等价式 a x + m y = b , y ∈ N ax+my=b,y\in N ax+my=b,y∈N;引入逆元后可转化为 x ≡ b a − 1 ( m o d m ) x\equiv ba^{-1}(\mod m) x≡ba−1(modm)
有解: gcd ( a , m ) = d \gcd(a,m)=d gcd(a,m)=d , d ∣ b d\mid b d∣b 则有 d d d 个模 m m m 不同余的解,否则无解
解法:转化为 a x + m y = d ax+my=d ax+my=d ,则乘 b / d b/d b/d 即可得到原式的解。
// 求解ax+my=b的x,即ax==b(mod m),答案存在ans[0...d-1]中,返回值d为答案个数
int mod_equation(int a, int b, int m)
{
exgcd(a, m, d, x, y);
if (b % d) // 无解
return -1;
x = x * (b / d) % m;
for (int i = 0; i < d; i++)
ans[i] = (x + i * m / d) % m;
return d;
}