数论学习笔记2

一.素数相关算法和定理

1. 素数筛法

   -埃氏筛法

   原理：素数的整数倍是合数

   时间复杂度为：$O(nlogn)$ 。

   const int N = 1e6 + 5;
   int p[N], cnt;
   bool isp[N];
   void make_prime(int n = N - 5) {
       memset(isp, true, sizeof(isp));
       isp[0] = isp[1] = false;
       for (int i = 2; i <= n; i ++) {
           if (! isp[i]) continue;
           p[++ cnt] = i;
           for (int k = i*i; k <= n; k += i) 
               isp[k] = true;
       }
   }
   

   -欧拉筛法（也叫线性筛法）

   原理：整数的素数倍是合数

   时间复杂度为：$O(n)$ 。

   const int N = 1e6 + 5;
   int p[N], cnt;
   bool isp[N];
   void Euler_Sieve(int n = N - 5) {
       for (int i = 2; i <= n; i ++) {
           if (! isp[i]) p[++ cnt] = i;
           // 关键处 1
           for (int j = 1; j <= cnt && (ll)i*p[j] <= n; j ++) {
               isp[i*p[j]] = true;
               if (i % p[j] == 0) break; // 关键处 2
           }
       }
   }
   

   下面证明欧拉筛为什么不重不漏：

   基于一个定理：整数唯一分解定理。

   不管 $i$ 是否是素数，都会执行 “关键处 1”。

   1. 如果 $i$ 是素数的话，那么一个大的素数 $i$ 的乘以不大于 $i$ 的素数，这样筛除的数跟之前的是不会重复的。筛出的数都是 $n=p_1\ast p_2$ 的形式，$p_1$ 和 $p_2$ 不相等；

   2. 如果 $i$ 是合数，此时 $i$ 可以表示为 $i=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \dots \ast p_n^{{\alpha }_n}$ ，$p_i$ 都是素数 $(1\le i\le n),p_i\le p_j(i\le j),{\alpha }_i>0$ ，$p_1$ 是最小的素数。

      按照程序的进行，当筛到下一个数为 $nex{t}_1=i\ast p_1=p_1^{{\alpha }_1+1}\ast p_2^{{\alpha }_2}\ast \dots \ast p_n^{{\alpha }_n}$ 时停止此轮筛除。

      假设程序未停止，将会继续筛除 $nex{t}_2=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2+1}\ast \dots \ast p_n^{{\alpha }_n},nex{t}_n=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \dots \ast p_n^{{\alpha }_n+1}$ 。

      显而易见的， $nex{t}_i(2\le i\le n)$ 都将被至少筛除两次。

      为什么呢？我们可以发现 $p_1^{{\alpha }_1-1}\ast p_2^{a_2+1}\ast \dots \ast p_n^{{\alpha }_n},\dots ,p_1^{{\alpha }_1-1}\ast p_2^{a_2}\ast \dots \ast p_n^{{\alpha }_n+1}$ （这些数都是大于 $i$ 的），它们乘上 $p_1$ 后都可把 $nex{t}_i$ 筛除。

   -区间筛法

   我们什么时候会使用到区间筛法呢，当数据范围为这样的时候：$[L,R]$ ，区间长度为 $R-L+1\le 10^6，L,R\le 10^{12}$ 。

   我们直接套用线筛或者埃氏筛法都会爆炸，这时候就要换个思路：

   原理：合数可以分解为两个不为 1 且不等于本身的因子相乘，即 $n=a\ast b(n为合数)$ ，且有 $min(a,b)\le \sqrt{n}$ 。

   证明：不妨设 $a\le b$ 。
   $$\begin{gathered} \because a\le b,a\ast b=n \\ \therefore a\ast a\le a\ast b=n \\ \therefore a\le \sqrt{n} \end{gathered}$$
   所以，我们可以在 $[1,\sqrt{n}]$ 的范围内找到 $n$ 的一个因子。

   我们先筛出 $\sqrt{R}$ 内的所有素数，然后再把区间 $[L,R]$ 内的所有合数都去掉，这里去掉所有合数主要利用了埃氏筛法的思想：素数的整数倍是合数。

   if (L == 1) L = 2; // 特判，1 不是质数也不是合数
   memset(isp, false, sizeof(isp));
   for (int i = 1; i <= cnt; i ++) {
       ll P = p[i], s = (L + P - 1) / P * P; // s 为起点
       if (s < 2*P) s = 2 * P;
       for (ll j = s; j <= R; j += P) 
           isp[j-L+1] = true; // 减去 L - 1 相当于做了个简单哈希
   }
   

2. Miller-Rabin素数探测

   基于两个定理：费马小定理 和 二次探测定理。

   -费马小定理：

   如果 $p$ 是一个素数，且 $0<a<p$ ，则有 $a^{p-1}\equiv 1(modp)$。

   -二次探测定理：

   如果 $p$ 是素数， $x$ 是小于 $p$ 的正整数，且 $x^2\equiv 1(modp)$ ，则要么 $x=1$ ，要么 $x=p-1$。

   证明：
   $$\begin{gathered} \because x^2\equiv 1(modp) \\ \therefore x^2-1\equiv 0(modp) \\ \therefore p|(x+1)(x-1) \\ \because p是质数 \\ \therefore x=k_{1}p+1orx=k_{2}p-1 \\ \because 0<x<p \\ \therefore x=p-1orx=1 \end{gathered}$$
   其逆否命题为：$x$ 是小于 $p$ 的正整数，当 $x^{2}modp=1$ （触发二次探测条件）时，如果 $x\ne 1$ 且 $x\ne p-1$ ，则说明 $p$ 一定不是素数。

   选择底数：

   // 当然也是能继续加一些素数作为底数，这里给出常用的 5 个
   const int MT = 5;
   int p[MT] = {2， 3， 7， 61， 24251};
   

   一些解决乘法溢出的办法：

   // 慢速乘 (慎重使用，可能会TLE)
   inline ll fmul(ll a, ll b, ll mod) {
       if (b < 0) a = -a, b = -b;
       ll res = 0;
       while (b) {
           if (b & 1) res = (res + a) % mod;
           a = (a + a) % mod;
           b >>= 1;
       }
       return res;
   }
   // 基于一种新数据类型 __int128
   inline ll fmul(ll a, ll b, ll mod) {
       return (__int128) a * b % mod;
   }
   

   快速幂：

   inline ll fpow(ll a, ll b, ll mod) {
       ll res = 1;
       for (; b; b >>= 1, a = fmul(a, a, mod))
           if (b & 1) res = fmul(res, a, mod);
   	return res;
   }
   

   Miller_Rabin素数测试：

   （1）对于偶数和 0，1，2 可以直接判断；

   （2）设要测试的数为 $x$ ，我们取一个较小的质数 $a$ ，设 $s,t$ ，满足 $2^s\cdot t=x-1$ （其中 $t$ 是奇数）。

   （3）我们先算出 $a^t$ ，然后不断地平方并且进行二次探测（利用其逆否命题）（进行$s$ 次）。

   （4）最后我们根据费马小定理，如果最后 $a^{x-1}\not\equiv 1(modx)$ ，则说明 $x$ 是合数。

   （5）多次取不同的 $a$ 进行 $Miller-Rabin$ 素数测试，这样可以使得正确性更高。

   bool Miller_Rabin(ll n) {
   	if (n == 2) return true;
       if (n < 2 || (n&1) == 0) return false;
       ll s = 0, t = n - 1;
       while (0 == (t&1)) s ++, t >>= 1;
       for (int i = 0; i < MT && p[i] < n; i ++) {
           ll b = fpow(p[i], t, n);
           for (int j = 0; j < s; j ++) {
               ll k = fmul(b, b, n);
               if (k == 1 && b != 1 && b != n - 1) \\二次探测定理
                   return false;
               b = k;
           }
           if (b != 1) return false;  \\费马小定理
       }
       return true;
   }
   

3. Pollard-Rho 算法（待更）

   Pollard-Pho算法是John Pollard 发明的一种能快速找到大整数的一个非1、非自身的因子的算法。

二. 中国剩余定理

1. 定理说明：
   $$(S):\{\begin{array}{l}\begin{array}{rlr}x& \equiv a_1& (mod{m}_1)\\ x& \equiv a_2& (mod{m}_2)\\ & ⋮& \\ x& \equiv a_n& (mod{m}_n)\end{array}\end{array}$$

   假设整数 $m_1,m_2,\dots ,m_n$ 两两互质，则对任意的整数：$a_1,a_2,\dots ,a_n$ ，方程组 $(S)$ 有解，并且解可以通过如下方式构造得到：

   设 $M=m_1\times m_2\times \dots \times m_n=\prod _{i=1}^n{m}_i$ 是整数 $m_1,m_2,\dots ,m_n$ 的乘积，并设 $M_i=M/m_i$ ， ∀ i ∈ { 1 , 2 , … , n } \forall i\in\{1,2,…,n\} ∀i∈{1,2,…,n} 是除了 $m_i$ 以外的 $n-1$ 个整数的乘积。

   设 $t_i=M_i^{-1}$ 为 $M_i$ 模 $m_i$ 的数论倒数（$t_i$ 为 $M_i$ 模 $m_i$ 意义下的逆元） $M_i{t}_i\equiv 1(mod{m}_i)$ ， ∀ i ∈ { 1 , 2 , … , n } \forall i\in\{1,2,…,n\} ∀i∈{1,2,…,n} 。

   方程组（$S$）的通解形式为 $x=a_1{t}_1{M}_1+a_2{t}_2{M}_2+\dots +a_n{t}_n{M}_n+kM=kM+\sum _{i=1}^n{a}_i{t}_i{M}_i$。

   在模 $M$ 意义下，方程组（$S$）只有一个解：$x=(\sum _{i=1}^n{a}_i{t}_i{M}_i)modM$ 。

2. 定理证明：从假设可知，对 ∀ i ∈ { 1 , 2 , … , n } \forall i\in\{1,2,…,n\} ∀i∈{1,2,…,n} ，由于 ∀ j ∈ { 1 , 2 , … , n } ,   j ≠ i ,   G C D ( m i ,   m j ) = 1 \forall j\in\{1,2,…,n\},\ j\not=i,\ GCD(m_i,\ m_j)=1 ∀j∈{1,2,…,n}, j=i, GCD(mi​, mj​)=1，所以 $GCD(m_i,M_i)=1$。这说明存在整数 $t_i$ 使得 $t_i{M}_i\equiv 1(mod{m}_i)$ 。这样的 $t_i$ 叫做 $M_i$ 模 $m_i$ 的逆元。

   对 ∀ i ∈ { 1 , 2 , … , n } \forall i\in\{1,2,…,n\} ∀i∈{1,2,…,n}，考察乘积 $a_i{t}_i{M}_i$ 可知：
   a i t i M i ≡ a i ⋅ 1 ≡ a i   ( m o d   m i ) ∀ j ∈ { 1 , 2 , … , n } ,   j ≠ i   , a i t i M i ≡ 0   ( m o d   m j ) ∴ x = a 1 t 1 M 1 + a 2 t 2 M 2 + … + a n t n M n 满足 : ∀ i ∈ { 1 , 2 , … , n } ,   x = a i t i M i + ∑ j ≠ i a j t j M j ≡ a i + ∑ j ≠ i 0 ≡ a i   ( m o d   m i ) a_it_iM_i\equiv a_i\cdot 1\equiv a_i\ (mod\ m_i)\\ \forall j\in \{1,2,…,n\},\ j\not=i\ ,a_it_iM_i\equiv 0\ (mod\ m_j)\\ \therefore x=a_1t_1M_1+a_2t_2M_2+…+a_nt_nM_n\\ 满足: \forall i\in\{1,2,…,n\},\ x=a_it_iM_i+\displaystyle\sum_{j\not=i}a_jt_jM_j\equiv a_i+\displaystyle\sum_{j\not=i}0\equiv a_i\ (mod\ m_i) ai​ti​Mi​≡ai​⋅1≡ai​ (mod mi​)∀j∈{1,2,…,n}, j=i ,ai​ti​Mi​≡0 (mod mj​)∴x=a1​t1​M1​+a2​t2​M2​+…+an​tn​Mn​满足:∀i∈{1,2,…,n}, x=ai​ti​Mi​+j=i∑​aj​tj​Mj​≡ai​+j=i∑​0≡ai​ (mod mi​)
   这说明 $x$ 就是方程组 （$S$）的一个解。

   另外，假设 $x_1$ 和 $x_2$ 都是方程组 （$S$）的解，那么： ∀ i ∈ { 1 , 2 , … , n } ,   x 1 − x 2 ≡ 0   ( m o d   m i ) \forall i\in\{1,2,…,n\},\ x_1-x_2\equiv 0\ (mod\ m_i) ∀i∈{1,2,…,n}, x1​−x2​≡0 (mod mi​)。而 $m_1,m_2,\dots ,m_n$ 两两互质，这说明 $M=\prod _{i=1}^n{m}_i$ 整除 $x_1-x_2$ ，且 $LCM(m_1,m_2,\dots ,m_n)=M$。所以方程组的通解为： { k M + ∑ i = 1 n a i t i M i ;   k ∈ Z } \{kM+\displaystyle\sum_{i=1}^na_it_iM_i;\ k\in\Z\} {kM+i=1∑n​ai​ti​Mi​; k∈Z}。另外，在模 $M$ 意义下，方程组（$S$）只有一个解：$x=(\sum _{i=1}^n{a}_i{t}_i{M}_i)modM$ 。

   从线性代数角度：考虑 $n$ 个方程组（$1\le i\le n$），
   $$\{\begin{array}{l}\begin{array}{rlrl}x& \equiv 0& (mod& m_1)\\ & ⋮& & \\ x& \equiv 0& (mod& m_{i-1})\\ x& \equiv 1& (mod& m_i)\\ x& \equiv 0& (mod& m_{i+1})\\ & ⋮& & \\ x& \equiv 0& (mod& m_n)\end{array}\end{array}$$
   由于 $m_i(1\le i\le n)$ 两两互质，这个方程组作变量替换，令 $x=(M/m_i)\ast y$ ，方程组等价于求解同余方程：$(M/m_i)y\equiv 1(mod{m}_i)$ ，求解出 $y$ 后，若要得到特解 $x_i$ ，只要令：$x_i=(M/m_i)\ast y_i$ ，则方程组的解为：$x_0=a_1{x}_1+a_2{x}_2+\dots +a_n{x}_n(modM)$ ，在模 $M$ 的意义下唯一。

3. 代码如下：

   ll CRT(ll *a, ll *m, ll n) {
       ll M = 1, ans = 0;
       for (int i = 1; i <= n; i ++) 
           M *= m[i];
       for (int i = 1; i <= n; i ++) {
           ll Mi = M / m[i], t, y;
           exgcd(Mi, m[i], t, y); 			// Mi * t = 1 (mod m[i])
           t = (t % m[i] + m[i]) % m[i];
           // ans += a[i] * t * Mi
           ans = (ans + fmul(fmul(a[i], t, M), Mi, M)) % M;
       }
       return ans;
   }
   

三. 扩展中国剩余定理

1. 定理说明：此时的 $m_i$ （模数）不需要两两互质。

2. 算法实现及其证明：

   我们先考虑两个式子：
   $$\begin{array}{rlr}x& \equiv a_1& (mod{m}_1)\\ x& \equiv a_2& (mod{m}_2)\end{array}$$
   我们设 $k_1$ 和 $k_2$ 是满足 $x=a_1+k_1{m}_1$ 和 $x=a_2+k_2{m}_2$ 的两个整数。

   显然有 $a_1+k_1{m}_1=a_2+k_2{m}_2$ ，即 $m_1{k}_1-m_2{k}_2=a_1-a_2$ 。

   根据线性同余方程，当 $GCD(m_1,m_2)|a_2-a_1$ 时有解。我们先将 $m_1{K}_1-m_2{K}_2=GCD(m_1,m_2)$ 的解 $K_1,K_2$ 求出，接着将式子左右两边同时乘以 $\frac{a_2-a_1}{GCD(m_1,m_2)}$ 得到 $k_1,k_2$ 。从而我们可以得到 $x_0$ （一个特解）。

   接着求出通解：我们设两个解为 $x_1$ 和 $x_2$ （$x_1<x_2$），因为两个解都满足给定的两个同余方程，因此存在两个整数 ${\delta }_1$ 和 ${\delta }_2$ ，$x_2=x_1+{\delta }_1{m}_1=x_1+{\delta }_2{m}_2$ 。

   我们设 $x_2-x_1=\Delta x$ ，显然 ${\delta }_1{m}_1={\delta }_2{m}_2$ （1），设 $d=GCD(m_1,m_2),n_1=\frac{m_1}{d},n_2=\frac{m_2}{d}$，则（1）式等价于 ${\delta }_1\times n_{1}d={\delta }_2\times n_{2}d$ ，即 ${\delta }_1{n}_1={\delta }_2{n}_2$ ，所以有 $n_1|{\delta }_2\times n_2$ ，由于 $GCD(n_1,n_2)=1$ 。因此 $n_1|{\delta }_2$ 。所以 $n_1{n}_2|{\delta }_2\times n_2$ ，进而 $n_1\times n_{2}d|\Delta x$。因此 $LCM(m_1,m_2)|\Delta x$ 。

   由此得出结论：任意两个解之间的差都是模数的最小公倍数。

   通解为： { x 0 + k ⋅ L C M ( m 1 ,   m 2 ) ;   k ∈ Z } \{x_0+k\cdot LCM(m_1,\ m_2);\ k\in\Z\} {x0​+k⋅LCM(m1​, m2​); k∈Z} 。

   从另一个角度我们将两个方程合并得到了一个新的方程：
   $$x\equiv x_0(modLCM(m_1,m_2))$$

3. 代码如下：

   pair<bool, ll> exCRT(ll* r, ll* m, ll n) {
       ll x, y, ans = r[0], M = m[0];
       for (int i = 1; i < n; i ++) {
           ll a = M, b = m[i], c = (r[i] - ans%b + b) % b;
           ll d = exgcd(a, b, x, y), q = b/d;
           if (d == 0 || c % d) return make_pair(0, 0);
           x = fmul(x, c/d, q);
           ans += x * M, M *= q;
           ans = (ans%M + M) % M;
       }
       return make_pair(1, (ans%M + M) % M);
   }
   

四. 欧拉定理

1. 同余类（Residue class）:

   通过一个整数模 $n$ 的余数，我们可以把所有整数分成 $n$ 类，记 r ‾ n = { m ∈ Z   ∣   m n + r } \overline{r}_n=\{m\in \Z\ |\ mn+r\} rn​={m∈Z ∣ mn+r} 为模 $n$ 余 $r$ 的同余类（也叫剩余类）。

   举个列子： 4 ‾ 10 = { … ,   − 16 ,   − 6 ,   4 ,   14 ,   24 ,   … } \overline{4}_{10}=\{…,\ -16,\ -6,\ 4,\ 14,\ 24,\ …\} 410​={…, −16, −6, 4, 14, 24, …} 是模 10 余 4 的同余类。

2. 完全剩余系（Complete residue system）：

   从 ${\overline{0}}_n,{\overline{1}}_n,{\overline{2}}_n,\dots ,{\overline{(n-1)}}_n$ 中各挑出一个数就组成了一个模 $n$ 的完全剩余系（完系）$R_n$。

   R n = { r 0 ,   r 1 ,   … ,   r n − 1 } R_n=\{r_0,\ r_1,\ …,\ r_{n-1}\} Rn​={r0​, r1​, …, rn−1​}，其中 $r_0\in {\overline{0}}_n,r_1\in {\overline{1}}_n,r_2\in {\overline{2}}_n,\dots ,r_{n-1}\in {\overline{(n-1)}}_n$ 。

   换句话说，$n$ 个模 $n$ 互不同余的整数组成一个模 $n$ 的完全剩余系。我们称 R n = { 0 ,   1 ,   … ,   n − 1 } R_n=\{0,\ 1,\ …,\ n-1\} Rn​={0, 1, …, n−1} 为模 $n$ 的一个最小非负完全剩余系（最小非负完系）。

3. 缩剩余系（Reduced residue system）：

   取一个模 $n$ 的完全剩余系 $R_n$ ，取出里面所有和 $n$ 互质的数，这些数组成一个模 $n$ 的缩剩余系（缩系），记为 ${\mathbf{\varphi }}_{\mathbf{n}}$ 。

   ϕ n = { c 1 ,   c 2 ,   … ,   c φ ( n ) } \boldsymbol{\phi_n}=\{c_1,\ c_2,\ …,\ c_{\varphi(n)}\} ϕn​={c1​, c2​, …, cφ(n)​} ，其中 $\phi (n)$ 是欧拉函数，代表 小于 $n$ 的正整数中的和 $n$ 互质的数 的个数。

   由于 $GCD(c_i,n)=GCD(c_i+kn,n)=1$ ，每一个模 $n$ 的缩剩余系有相同数量的元素（缩剩余系中的每一个数所属的同余类是确定的，所以总有确定的 $\phi (n)$ 个同余类）

   如果缩剩余系 ϕ n = { c 1 ,   c 2 ,   … ,   c φ ( n ) } \boldsymbol\phi_n=\{c_1,\ c_2,\ …,\ c_{\varphi(n)}\} ϕn​={c1​, c2​, …, cφ(n)​} 满足 $1\le c_1,c_2,\dots ,c_{\phi (n)}\le n-1$ ，那么称其为模 $n$ 的最小正缩剩余系（最小正缩系）。

4. 欧拉函数（Euler‘s totient function）：

   对于正整数 $n$ ，$\phi (n)$ 代表 **小于 $n$ 的正整数中的和 $n$ 互质的数 **的个数，这个函数被称为欧拉函数，欧拉还告诉我们
   $$\frac{\phi (n)}{n}=\prod _{p|n}(1-\frac{1}{p})$$
   其中 $p$ 取到 $n$ 的所有质因子。

   所以我们可以很方便地计算一个整数的欧拉函数的值，比如：
   $$\phi (1926)=\phi (2\times 3^2\times 107)=1926(1-\frac{1}{2})(1-\frac{1}{3})(1-\frac{1}{107})=636$$
   欧拉函数的证明：

   首先我们考虑 $\phi (p^e)$ ，其中 $p$ 是质数，$e$ 是非负整数，如果要使 $GCD(p^e,k)\ne 1$ ，只能让 $k$ 等于 $p$ 的倍数，在 $1\le k\le p^e$ 范围内，$p$ 的倍数有 $p,2p,3p,\dots ,p^{e-1}p$ ，总共 $p^{e-1}$ 个数，所以 $\phi (p^e)=p^e-p^{e-1}=p^e(1-\frac{1}{p})$ 。

   然后我们证明对于 $GCD(m,n)=1$ 有 $\phi (mn)=\phi (m)\phi (n)$ ，我们首先构造两个集合，第一个集合是模 $mn$ 的最小正缩系 ${\mathbf{\varphi }}_{mn}$ ，第二个集合定义为：
   $$S=\{(m,n)|m\in {\mathbf{\varphi }}_m,n\in {\mathbf{\varphi }}_n\}$$
   其中 ${\mathbf{\varphi }}_{\mathbf{m}},{\mathbf{\varphi }}_{\mathbf{n}}$ 分别是模 $m,n$ 的最小正缩系，显然 $|{\mathbf{\varphi }}_{mn}|=\phi (mn)$ 和 $|S|=\phi (m)\phi (n)$ 。

   如果我们证明存在一个双射 $f$ ：${\mathbf{\varphi }}_{mn}\to S$ ，就证明了 $\phi (mn)=\phi (m)\phi (n)$ 。

   我们构造函数如下：
   $$f(a)=(amodm,amodn)$$
   首先我们用反证法证明 $f$ 是单射（即证明存在两个数 $a\ne b,f(a)=f(b)$ ，从而导出矛盾），假设 $a,b\in {\mathbf{\varphi }}_{\mathbf{mn}}$ 满足 $a\ne b且f(a)=f(b)$ ，那么就有
   $$\begin{array}{rlrl}a& \equiv b& (mod& m)\\ a& \equiv b& (mod& n)\end{array}$$
   显然因为 $GCD(m,n)=1$ 我们能得出 $a\equiv b(modmn)$ ，这与我们的假设矛盾（因为 ${\mathbf{\varphi }}_{mn}$ 是模 $mn$ 的缩系，$a,b$ 是 ${\mathbf{\varphi }}_{mn}$ 的两个不同的元素，所以他们模 $mn$ 不同余）。接下来，中国剩余定理 告诉我们

   如果整数 $r_1,r_2$ 和正整数 $GCD(n_1,n_2)=1$ ，同余方程组
   $$\begin{array}{rlr}x& \equiv r_1& (mod{n}_1)\\ x& \equiv r_2& (mod{n}_2)\end{array}$$
   在 $0\le x<n_1{n}_2$ 范围内有且只有一个解。

   通过中国剩余定理我们能够证明 $f$ 是满射，所以 $f$ 是双射。

   所以对于 $GCD(m,n)=1$ ，就有 $\phi (mn)=\phi (m)\phi (n)$ ，假设 $n$ 的唯一分解为 ：
   $$n=2^{e_1}\times 3^{e_2}\times 5^{e_3}\times \dots =\prod _{k=1}^{\infty }{p}_k^{e_k}$$
   其中 $e_1,e_2,\dots \in \mathbb{N}$ ，那么：
   $$\begin{array}{rl}\phi (n)& =\phi (\prod _{k=1}^{\infty }{p}_k^{e_k})=\prod _{k=1}^{\infty }\phi (p_k^{e_k})=\prod _{k=1}^{\infty }{p}_k^{e_k}(1-\frac{1}{p_k})\\ & =(\prod _{k=1}^{\infty }{p}_k^{e_k})\prod _{k=1}^{\infty }(1-\frac{1}{p_k})=n\times \prod _{k=1}^{\infty }(1-\frac{1}{p_k})\end{array}$$
   欧拉函数还有一些性质（证明待补）：

   • 如果正整数 $n>2$ ，那么 $\phi (n)$ 是偶数。
   • 如果 $n|N$ ，那么 $\phi (n)|\phi (N)$ 。
   • 对于正整数 $m,n$ 有 $\phi (mn)=\phi (m)\phi (n)\frac{GCD(m,n)}{\phi (GCD(m,n))}$ 。
   • 特别的，如果 $m,n$ 互质，那么 $\phi (mn)=\phi (m)\phi (n)$ 。
   • 对于正整数 $n$ 有 $\sum _{d|n}\phi (d)=n$ 。
   • 对于正整数 $n$ 有 $\sum _{1\le k\le n,GCD(k,n)=1}\frac{k}{n}=\frac{\phi (n)}{2}$ 。

   // 线性预处理欧拉函数 phi[i]
   int phi[N];
   int p[N], cnt;
   bool isp[N];
   void GetPhi(int n = N - 5) {
       for (int i = 2; i <= n; i ++) {
           if (! isp[i]) phi[i] = i - 1, p[++ cnt] = i;
           for (int j = 1; j <= cnt && i * p[j] <= n; j ++) {
               isp[i*p[j]] = true;
               if (i % p[j]) phi[i*p[j]] = phi[i] * (p[j] - 1);
               else { phi[i*p[j]] = phi[i] * p[j]; break; }
           }
       }
   }
   

接下来进入正题：欧拉定理

1. 定理描述：

   如果正整数 $n$ 和整数 $a$ 互质，那么就有 $a^{\phi (n)}\equiv 1(modn)$ ，其中 $\phi (n)$ 是欧拉函数。

   证明：考虑模 $n$ 的最小正缩系 ϕ n = { c 1 ,   c 2 ,   … ,   c φ ( n ) } \boldsymbol\phi_n=\{c_1,\ c_2,\ …,\ c_{\varphi(n)}\} ϕn​={c1​, c2​, …, cφ(n)​} ， 已知 $GCD(a,n)=1$ ，我们在 ${\mathbf{\varphi }}_n$ 的每一个元素面前都乘以一个 $a$ 。
   a ϕ n = { a c 1 ,   a c 2 ,   … ,   a c φ ( n ) } a\boldsymbol{\phi_n} = \{ac_1,\ ac_2,\ …,\ ac_{\varphi(n)}\} aϕn​={ac1​, ac2​, …, acφ(n)​}
   利用反证法可以证明 $a{\mathbf{\varphi }}_{\mathbf{n}}$ 也是一个模 $n$ 的缩系（其元素的同余类的顺序有可能会改变，但是这并没有任何影响），假设 $a{c}_i\equiv a{c}_j(modn)$ ，其中 $i\ne j$ ，因为 $a,n$ 互质可以将$a$ 消去，那么就有 $c_i\equiv c_j(modn)$ ，这是不可能的，矛盾。
   $$\begin{array}{rlr}\prod _{i=1}^{\phi (n)}{c}_i& \equiv \prod _{i=1}^{\phi (n)}a{c}_i& (modn)\\ \prod _{i=1}^{\phi (n)}{c}_i& \equiv a^{\phi (n)}\prod _{i=1}^{\phi (n)}{c}_i& (modn)\end{array}$$
   显然 $GCD(n,\prod _{i=1}^{\phi (n)}{c}_i)=1$ ，所以可以两边消去它，得到：$a^{\phi (n)}\equiv 1(modn)$ 。

   由此，当 $n=p$ 是一个质数时，我们可以得到费马小定理（Fermat’s little theorem) 。

   如果 $p$ 是一个质数，且 $GCD(a,p)=1$，那么 $a^{p-1}\equiv 1(modp)$ ，即 $p|n^p-n$ 对于任意整数 $n$ 满足 $GCD(n,p)=1$ 都成立。

   当然，费马小定理也可以用归纳法证明：假设 $p|n^p-n$ ，那么
   $$(n+1{)}^p-(n+1)\equiv \sum _{r=1}^{p-1}(\genfrac{}{}{0ex}{}{p}{r})n^r+n^p-n$$
   当 $1\le r\le p-1$ 时，二项式系数 $\left(\genfrac{}{}{0ex}{}{p}{r}\right)=\frac{p!}{(p-r)!r!}$ 的分子中有 $p$ ，分母中每一个因子都不能整除 $p$ （因为 $p$ 是质数），所以 $p$ 能够整除 $\left(\genfrac{}{}{0ex}{}{p}{r}\right)$，进而得到 $p|(n+1{)}^p-(n+1)$ 。当 $n=0$ 时显然成立，所以定理成立。特别的 $n=0$ 时显然成立。