数论-常见定理的探究

关键词：威尔逊定理 欧拉定理 费马小定理 中国剩余定理

威尔逊定理

当且仅当p为素数时，(p - 1)! ≡ -1 (mod p)

• 上式也可写作 (p - 1)! ≡ p-1 (mod p) 或 (p-1)! + 1 ≡ 0 (mod p)
• 上式可作素数的判别条件，但复杂度较大，实用性较差。

欧拉定理

若n,a为正整数，且n,a互质，即gcd(a,n) = 1，则
a^φ(n) ≡ 1 (mod n)

• φ(n) 称欧拉函数，欧拉函数是求小于n的数中与n互质的数的数目

• 如果n是质数，那么1到n-1所有数都是与n互质的，所以φ(n) = n-1

• 如果n是合数，例如φ(8)=4，因为1,3,5,7均和8互质

• 欧拉函数求法：

  • 基本原理：1到n-1中找到与n不互质的数，然后把他们减掉
  • 具体操作*：找到n的质因子，把它们和它们的倍数都删掉
  • 优化技巧：利用容斥原理避免重复处理
  • 例子：φ(12)
    • 找质因数：12 $\to$ {2 3}
    • 找倍数：[2] ≡ {2 4 6 8 10} 、[3] ≡ {3 6 9}
    • 做统计：φ(12) = 12 - (12 / 2) - (12 / 3 ) + 1 = 12 * (1 - 1/2 - 1/3 + 1/6)
    • 容斥优化：φ(12) = 12 * (1 - 1/2) * (1 - 1/3) = 12 * (1 - 1/2 - 1/3 + 1/6)

• 欧拉定理可用于求逆元，显然 inv(a, n) = a ^ (φ(n) - 1)

实现代码

• 求单个欧拉函数（O(n^0.5）

//欧拉函数
int phi(int x){
    int ans = x;
    for(int i = 2; i*i <= x; i++){
        if(x % i == 0){
            ans = ans / i * (i-1);
            while(x % i == 0) x /= i;
        }
    }
    if(x > 1) ans = ans / x * (x-1);
    return ans;
}

• 求n个欧拉函数（O(n)）

  注：类似于素数筛

#include<cstdio>
const int N = 100000 + 5;
int phi[N];
void Euler(){
    phi[1] = 1;
    for(int i = 2; i < N; i ++){
        if(!phi[i]){
            for(int j = i; j < N; j += i){
                if(!phi[j]) phi[j] = j;
                phi[j] = phi[j] / i * (i-1);
            }
        }
    }
}
int main(){
    Euler();
}

费马小定理

如果p是一个质数，而整数a不是p的倍数，则有a ^ (p - 1) ≡ 1(mod p)。

• 该定理可用于求逆元，显然 inv(a, p) = a ^ (p-2)

中国剩余定理

目标问题

求解线性同余方程组的最小整数解

\begin{cases}
    x ≡ a_1 (mod m_1) \\
    x ≡ a_2 (mod m_2) \\
    ...\\
    x ≡ a_n (mod m_n) \\
\end{cases}    

其中，$m_1$、$m_2$、… 、$m_n$为两两互质的整数

定理内容

令$M={\prod }_{i=1}^n{m}_i$,$M_i=M/m_i$,$t_i=inv(M_i,m_i)$,则方程组通解为
$$x=\sum _{i=1}^n{a}_i{t}_i{M}_i+kM,k\in Z$$
在模M意义下，方程唯一解为
$$x=\sum _{i=1}^n{a}_i{t}_i{M}_i$$

简单证明

由逆元定义知
$$\begin{array}{rl}{x}_i& =a_i\ast t_i{M}_i\\ & =a_i\ast inv(M_i)\ast M_i\\ & =a_i\ast (m_1...m_{i-1}{m}_{i+1}...m_n)\ast inv((m_1...m_{i-1}{m}_{i+1}...m_n))\\ & \equiv a_i\ast 1(modmi)\\ & \equiv a_i(mod{m}_i)\end{array}$$
故
$$\begin{array}{rl}{a}_i& |x_i\\ m_k& |x_i,k\ne i\end{array}$$
即
$$\begin{array}{rl}{x}_i& \equiv a_i(mod{m}_i)\\ x_i& \equiv 0(mod{M}_k),k\ne i\end{array}$$
令
$$x=\sum _{i=1}^n{x}_i$$
则
$$\begin{array}{rl}x\%m_i& =(x_1+...+x_i+...+x_n)\%m_i\\ & =(x_1\%m_i+...+x_i\%m_i+x_n\%m_i)\%m_i\\ & =(0+...+a_i+...+0)\\ & =a_i\end{array}$$
即
$$x\equiv a_i(mod{m}_i)$$

实现代码

//n个方程：x=a[i](mod m[i]) (0<=i<n)
LL china(int n, LL *a, LL *m){
    LL M = 1, ret = 0;
    for(int i = 0; i < n; i ++) M *= m[i];
    for(int i = 0; i < n; i ++){
        LL w = M / m[i];
        ret = (ret + w * inv(w, m[i]) * a[i]) % M;
    }
    return (ret + M) % M;
}

参考博客

• https://www.cnblogs.com/linyujun/p/5194142.html