数论——同余

数论有一部分的总括

---

设$a,b\in \mathbb{N},a̸=0$。
记$a$除$b$即$b÷a$的余数为$b\mathrm{m}\mathrm{o}\mathrm{d}a$。
若存在$m\in {\mathbb{Z}}^{+}$，使得$a\mathrm{m}\mathrm{o}\mathrm{d}m=b\mathrm{m}\mathrm{o}\mathrm{d}m$，记作$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$

1.整除

若存在整数$p̸=0$使得$aq=b$，那么称$a$是$b$的因数（约数），$b$是$a$的倍数，$a$整除$b$，$b$能被$a$整除，记作$a\mid b$，同时$b\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}a)$。

性质

$a\mid b,b\mid c\Rightarrow a\mid c$
$a\mid b,c\mid d\Rightarrow ac\mid bd$
$a\mid b,a\mid c\Rightarrow a\mid pb+qc$

整除也可以定义同余，即$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)\iff m\mid (a-b)$
所以，若$\mathbf{a}\equiv \mathbf{b}(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{m}),\mathbf{n}\mid \mathbf{m}\iff \mathbf{a}\equiv \mathbf{b}(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{n})$

2.性质

自反性：$a\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}m)$
对称性：$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)\Rightarrow b\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}m)$
传递性：$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m),b\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)\Rightarrow a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)$

$\mathbf{a}\equiv \mathbf{b}(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{m})\iff \mathbf{k}\mathbf{a}\equiv \mathbf{k}\mathbf{b}(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{k}\mathbf{m})$
证明：$m\mid (a-b)\iff km\mid k(a-b)\iff km\mid (ka-kb)\iff ka\equiv kb(\mathrm{m}\mathrm{o}\mathrm{d}km)$

$\mathbf{a}\mathbf{c}\equiv \mathbf{b}\mathbf{c}(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{m})\iff \mathbf{a}\equiv \mathbf{b}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{\mathbf{m}}{(\mathbf{c},\mathbf{m})})$
证明：（充分性）有$m\mid c(a-b)$。
若$(m,c)=1$，则$m\mid (a-b)$；
否则，有$\frac{m}{(c,m)}\mid \frac{c}{(c,m)}(a-b)$，那么$\frac{m}{(c,m)}\mid (a-b)$，因为$(\frac{m}{(c,m)},\frac{c}{(c,m)})=1$

3.辗转相除法

就是欧几里得算法，用于在$O(\log \max (a,b))$时间内求出$(a,b)$。

这是基于一个定理的——$(a,b)=(r,b)$，其中$a\equiv r(\mathrm{m}\mathrm{o}\mathrm{d}b)$
证明：
设$a=mc,b=nc,(m,n)=1,a=kb+r$，则$a=mc=kb+r=knc+r\Rightarrow r=(m-kn)c$
所以说现在需要证$((m-kn)c,n)=1$
使用反证法，设$m-kn=xd,n=yd,d>1$，则有$m=xd+kn=xd+kyd=(x+ky)d$，这样$d\mid (m,n)$，与$(m,n)=1$矛盾，所以原命题成立。

模板代码：__gcd(a,b)

inline long long gcd(long long a, long long b) { return b ? gcd(b, a % b) : a; }
inline long long lcm(long long a, long long b) { return a / gcd(a,b) * b; }

由于$a\mathrm{m}\mathrm{o}\mathrm{d}b<\frac{a}{2}$，每次递归后$a,b$中间至少有一个数缩小一半，所以该算法求$gcd$的时间复杂度为$O(\log \max (a,b))$

4.欧拉函数

4.1.完全剩余系

定义模m同余的所有整数构成一个剩余类。将整数划分为m个剩余类，同一个类中所有的数模m同余，来自不同剩余类的数模m不同余。
例如模3意义下可分为以下3个剩余类 { … , − 3 , 0 , 3 , 6 , … } , { … , − 2 , 1 , 4 , 7 , … } , { … , − 1 , 2 , 5 , 8 , … } \{…,-3,0,3,6,…\}, \{…,-2,1,4,7,…\}, \{…,-1,2,5,8,…\} {…,−3,0,3,6,…},{…,−2,1,4,7,…},{…,−1,2,5,8,…}
记号[a]表示a所在的剩余类，例如 [ 0 ] = { … , − 3 , 0 , 3 , 6 , … } [0]= \{…,-3,0,3,6,…\} [0]={…,−3,0,3,6,…}，
对于集合S，S中所有数模m的余数构成的集合称为S在模m意义下的剩余系。如果S为全体整数，此剩余系称为完全剩余系。
一般地，模m意义下的完全剩余系为 { 0 , 1 , 2 , … , m − 1 } \{0,1,2,…,m-1\} {0,1,2,…,m−1}
剩余系中与m互质的数的集合称为简化剩余系，或缩系。

完全剩余系的性质：

1. 对于$m$个整数，其构成模$m$的完系等价于其关于模$m$两两不同余
2. 若$a_i(1\le i\le m)$构成模$m$的完系$k,t\in \mathbb{Z},(k,m)=1$，则$k{a}_i+t$也构成模$m$的完系。

4.2.简化剩余系

简化剩余系又称既约剩余系，缩系。
定义：剩余系中与m互质的数的集合称为简化剩余系。
性质：

1. 设$m$为自然数，$k,t$为任意整数，$(k,m)=1$，则当x属于m的简化剩余系时，$kx+tm$亦属于模$m$的简化剩余系。
2. 设$m$为自然数，$k,t$为任意整数，$(k,m)=1$，若ai组成了模m的简化剩余系，那么$k{a}_i+tm$也组成模$m$的简化剩余系。

4.3.正题

欧拉是一个大坏人……

一个数n的简化剩余系的大小记为$\phi (n)$，称为欧拉函数，一个等价的定义是$n$以内和$n$互质的数的个数。
有$a>1$，则$a={p_1}^{w_{{}_1}}\ast {p_2}^{w_{{}_2}}\ast {p_3}^{w_{{}_3}}\ast \cdots \ast {p_m}^{w_{{}_m}}$，则$$\phi (n)=n\prod _{i=1}^m(1-p_i)=\prod _{i=1}^m{p_i}^{w_{{}_i}-1}(p_i-1)$$
（老师说：）从表达式看出，欧拉函数是积性函数。

性质（鬼火冒）：

1. 对于质数$p$，$\phi (p)=p-1$
2. 当$n>2$时，$\phi (n)$是偶数
3. 小于$n$且与$n$互质的数的总和为$\phi (n)\ast \frac{n}{2}$
4. $$n=\sum _{d\mid n}\phi (n)$$

4.4.计算

一般来说，我们可以根据上面的两个公式求出$\phi (n)$。
公式一计算($O(\sqrt{n})$)：$$\phi (n)=n\prod _{i=1}^m(1-p_i)$$

inline int phi(int n) {
    int res = n;
    for(register int i = 2; i <= sqrt(n); i++) {
        if(n % i == 0) {
            res = res / i * (i - 1);
            while(n % i == 0) n /= i;
        }
    }
    if(n > 1) res = res / n * (n - 1);
    return res;
}

该方法适合于单个计算的情况。
公式二计算：($O(n)$)
注：此方法是基于欧拉筛的，得到$\phi (1)-\phi (n)$

const int MAXN = 100001;
bool vis[MAXN];
int phi[MAXN], pr[MAXN], tot;
inline void Euler(int n) {
    phi[1] = 1;
    for(register int i = 2; i <= n; i++) {
        if(!vis[i]) {
            pr[++tot] = i;
            phi[i] = i - 1;
        }
        for(register int j = 1; j <= tot && i * pr[j] <= n; j++) {
            vis[i * pr[j]] = 1;
            if(i % pr[j]) phi[i * pr[j]] = phi[i] * (pr[j] - 1);
            else {
                phi[i * pr[j]] = phi[i] * pr[j];
                break;
            }
        }
    }
}

5.一些定理

5.1.欧拉定理

若$(\mathbf{a},\mathbf{n})=1,\mathbf{a},\mathbf{n}\in {\mathbb{Z}}^{+}$，则${\mathbf{a}}^{\phi (\mathbf{n})}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{n})$

证明：设$n$的缩系为$b_1,b_2,\cdots ,b_{\phi (n)}$
都乘以$a$，那么$a{b}_1,a{b}_2,\cdots ,a{b}_{\phi (n)}$亦构成n的缩系（缩系性质二）
$b_1{b}_2\cdots b_{\phi (n)}\equiv a{b}_{1}a{b}_2\cdots a{b}_{\phi (n)}\equiv a^{\phi (n)}{b}_1{b}_2\cdots b_{\phi (n)}(\mathrm{m}\mathrm{o}\mathrm{d}n)$
由于$n$的缩系内所有数之积与$n$是互质的，故约去$b_1{b}_2\cdots b_{\phi (n)}$
得$a^{\phi (n)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$

欧拉定理也有特殊情况——那就是费马小定理。
有$p\in \mathbb{P},p\nmid a$，则${\mathbf{a}}^{\mathbf{p}-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{p})$

欧拉定理也可以用于降幂，可以得到${\mathbf{a}}^{\mathbf{b}}\equiv {\mathbf{a}}^{\mathbf{b}\%\phi (\mathbf{m})+\phi (\mathbf{m})}(\mathrm{m}\mathrm{o}\mathrm{d}\mathbf{m}),\mathbf{b}\ge \phi (\mathbf{m})$

5.2.二次探测定理

$p$是奇素数，在$p$的缩系内解方程$x^2\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$
$x^2-1\equiv (x+1)(x-1)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$
显然有$1$和$p-1$两个解，这两个解是唯二解。

5.3.威尔逊定理

$$p\in \mathbb{P}\iff (p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$
证明：证明

5.4.Miller-Rabin

实在写不下去了 给链接
和代码

typedef long long ll;
inline ll mul(ll a, ll b, ll mod) {
    ll res = 0;
    for (; b ; a = (a + a) % mod, b >>= 1)
        if (b & 1)
            res = (res + a) % mod;
    return res;
}
inline ll ksm(ll a, ll n, ll mod) {
    ll res = 1;
    for (; n ; a = mul(a, a, mod), n >>= 1)
        if (n & 1)
            res = mul(res, a, mod);
    return res;
}
inline bool Miller_Rabin(ll n) {
    if(n == 2 ||  n == 3 ||  n == 5 ||  n == 7) return 1;
    if(n % 2 == 0|| n % 3 == 0 ||  n % 5 == 0 ||  n % 7 == 0) return 0;
    ll m = n - 1, k = 0;
    while(!(m & 1)) {
        k++; m >>= 1;
        for(register int i = 1; i <= 10; i++) {
            ll a = rand() % (n - 1) + 1;
            ll x = ksm(a, m, n); ll y;
            for(register int j = 1; j <= k; j++) {
                y = mul(x, x, n);
                if (y == 1 && x != 1 && x != n - 1)
                    return 0;
                x = y;
            }
            if(y != 1) return 0;
        }
    }
    return 1;
}

另外，进行多轮素性测试，误判率将指数级收敛。

上一篇
下一篇
回到原来的地方