【数论】——同余

【数论】——同余

定义

• 若存在整数 a.b，除以 m 的余数相同，则称 a，b mod m 同余，记为：
  $$a\equiv b(modm)$$

同余系 & 剩余系

1. 同余系：对区间 [1,m-1], 集合 {a+km} 的所有数 mod m，余数都是 a，称集合为一个 mod m 的同余类，记做：$\overline{a}$。
2. 剩余系：mod m 的同余类有 m 个：
   $$\overline{0},\overline{1}...\overline{m-1}$$
   构成 m 的完全剩余系；
   [1,m] 中与 m 互质的数代表的同余类有：$\varphi (m)$ 个，构成 m 的简化剩余系。

• 简化剩余系关于 mod m 乘法封闭：
  任取其中 a，b 与 m 互质，则 a*b 也与 m 互质（反证法证明）。

扩展欧几里得

• 对于任意整数 a，b，存在一对整数 x，y，有：
  $$a\ast x+b\ast y=\gcd (a,b)$$

证明

1. 若 b = 0，则当：
   $$x=1,y=0$$
   公式成立。
2. 若 b>0, 根据欧几里得算法：
   $$\gcd (a,b)=\gcd (b,amodb)$$
   此时不妨假设存在：$x_0,y_0$ 使得
   $$b\ast x_0+(amodb)\ast y_0=\gcd (b,amodb)$$
   对 $(amodb)$ 展开有：
   $$amodb=a-b\ast [\frac{a}{b}]$$
   代入有：
   $$b\ast x_0+(a-b\ast [\frac{a}{b}])\ast y_0=\gcd (b,amodb)$$
   分理出 x，y：
   $$a\ast y_0+b\ast (x_0-[\frac{a}{b}]\ast y_0)=\gcd (b,amodb)=\gcd (a,b)$$
   令，$x^{\prime }=y_0,y^{\prime }=x_0-[\frac{a}{b}]\ast y_0$, 代入上式：
   $$a\ast x^{\prime }+b\ast y^{\prime }=\gcd (a,b)$$
   通过以上方式，可以运用数学归纳法，得证扩展欧几里得。

代码实现

• 在欧几里得的基础上，传入 x，y，运用：
  $$x^{\prime }=y,y^{\prime }=x-[\frac{a}{b}]\ast y$$
  即可递推得的到系数 x，y。
• 代码

int extend_gcd(int a, int b, int& x, int& y)// 函数返回 a，最大公约数
{if (!b) {
        x = 1, y = 0;
        return a;
    }
    int d = gcd(b, a % b, x, y);
    int t = x;// 临时记录 x
   // 利用：x'= y,y' = x-[a/b]*y
    x = y;
    y = t - y * (a / b);
    return d;
}

线性同余方程

定义

• 对于整数 a，b，m，满足：
  $$a\ast x\equiv b(modm)$$
  求整数 x（有解或无解）。由于未知数最高次为 1，因此又被称为一次同余方程。

求解

• 由定义式易知：$a\ast x-b$ 是 m 的倍数，不妨设商为 $-y$
  对定义式去模运算：
  $$a\ast x+m\ast y=b$$
  则原方程的解等价于求满足上式的 x 的值。
• 方程有解等价于：
  $$b\mid gcd(a,m)$$
• 因此可以借助扩展欧几里得进行求解

代码

int solve(int a, int m, int& x, int& y)
{if (!m) {
        x = 1, y = 0;
        return a;
    }
    int d = solve(m, a % m, x, y);
    int t = x;
    y = x;
    x = t - (a / m) * y;
    return d;
}
bool check(int m,int d)
{return m % d == 0;}
void answer(int a,int m)
{
    int x, y;
    int d = solve(a, m, x, y);
    int answer;
    if (check(m, d))
        answer = (x % m + m) % m;
    else
        puts("impossible");
}

线性同余方程组

定义

• 对于 n 对模数与被模数 $a_i,m_i$, 有：
  $$x\equiv a_1(mod{m}_1)$$
  $$x\equiv a_2(mod{m}_2)$$
  $$...$$
  $$x\equiv a_n(mod{m}_n)$$

求解

1. 中国剩余定理
   设 $m,M_i$ 为 $$m=\prod _{i=1}^n{m}_i$$
   $$M_i=\frac{m}{m_i}$$
   若 $m_i$ 之间两两互质，设 $t_i$ 为：（$M_i^{-1}$）
   $$M_i\ast t_i\equiv 1(mod{m}_i)$$
   解 $x$ 为：
   $$x=\sum _{i=1}^n{a}_i\cdot M_i\cdot t_i$$

• 证明——（来自 李煜东 《算法竞赛进阶指南》）
  

2. m 不满足互余下的通解

• 取出前两个式子：
  $$x\equiv a_1(mod{m}_1)$$
  $$x\equiv a_2(mod{m}_2)$$
  去模处理：
  $$x=k_1\ast a_1+m_1$$
  $$x=k_2\ast a_2+m_2$$
  联立并化简：
  $$k_1\ast a_1-k_2\ast a_2=m_2-m_1$$
  根据扩展欧几里得, 求解 $k_1,k_2$
  不定方程 $x=k_1\ast a_1+m_1$ 的解为：
  $$x=(k_1+k\frac{a_1}{gcd(a_1,a_2)})\ast a_1+m_1$$
  化简：
  $$x=a_1\ast k_1+m_1+k\ast lcm(a_1,a_2)$$
  令 $m=a_1\ast k_1+m_1,a=lcm(a_1,a_2)$
  $$x=k\ast a+m$$
  因此方程组中所有式子两两合并，n-1 次后最终可以化为一个形如 $x=k\ast a^{\prime }+m^{\prime }$ 的等式。
  最终求解：
  $$xmod{a}^{\prime }\equiv m^{\prime }$$
  中 x 的解，等价于：$mmoda$的最小正整数。
• 代码

typedef long long LL;
LL exgcd(LL a, LL b, LL& x, LL& y)
{
    if (!b) {
        x = 1, y = 0;
        return a;
    }
    LL d = exgcd(b, a % b, x, y);
    LL t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
}
LL mod(LL a, LL b)
{
    return ((a % b) + b) % b;
}
int main()
{
    int n;
    scanf("%d", &n);
    LL a, m;
    bool flag = true;
    scanf("%lld%lld", &a, &m);//先读入第一个式子
    n--;
    //读入剩余的并合并
    while (n--) {
        LL a0, m0;
        scanf("%lld%lld", &a0, &m0);
        LL k1, k2;
        LL d = exgcd(a, -a0, k1, k2);
        if ((m0 - m) % d) {//如果同余方程无解
            flag = false;
            break;
        }
        k1 = mod(k1 * (m0 - m) / d, abs(a0 / d));//防止溢出
        m = a * k1 + m;//合并后的m = a*k1+m
        a = abs(a / d * a0);//合并后的a = lcm(a,a0）
    }
    if (flag)
        printf("%lld", mod(m, a));
    else
        puts("-1");
    return 0;
}