ACM 数论（一些定理及应用）

常用定理

1. 欧拉定理（费马-欧拉定理）

若$n,a$为正整数且$n,a$互素，则有如下公式： $a^{\phi (n)}\equiv 1$ (mod n)
含义： $a^{\phi (n)}$与1在模n下 同余，其中$\phi (n)$为欧拉函数，代表小于或等于$n$的正整数中与$n$互质的数的数目。
举个栗子：$a=3,n=5$，小于5的正整数中与5互质的数为1,2,3,4，所以$\phi (5)=4$，那么$3^4\equiv 1$ mode 5，成立。
应用：本定理可大程度地简化幂的模运算。${}^{[1]}$比如计算$7^{222}$的个位数时，可将此命题视为求$7^{222}$被10除的余数。因为7和10互质，而$\phi (10)=4$，所以$7^{222}=7^{4\cdot 55+2}=(7^4{)}^{55}\cdot 7^2\equiv 1^{55}\cdot 7^2\equiv 49\equiv 9(\mathrm{m}\mathrm{o}\mathrm{d}10)$。

2. 费马小定理

假如a是一个整数，p是一个质数，那么$a^p-a$是p的倍数，可以表示为$a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$。如果a不是p的倍数，这个定理也可以写成$a^{p-1}\equiv 1$ (mod p)。
应用：计算$2^{100}$除以13的余数:
$2^{100}\equiv 2^{12\times 8+4}(\mathrm{m}\mathrm{o}\mathrm{d}13)\equiv (2^{12}{)}^8\cdot 2^4(\mathrm{m}\mathrm{o}\mathrm{d}13)\equiv 1^8\cdot 16(\mathrm{m}\mathrm{o}\mathrm{d}13)\equiv 16(\mathrm{m}\mathrm{o}\mathrm{d}13)\equiv 3(\mathrm{m}\mathrm{o}\mathrm{d}13)$

3.威尔逊定理

当且仅当p为质数时：$(p-1)!\equiv -1$ mod p
含义：给出了判定一个自然数是否为质数的充分必要条件。但是由于阶乘是呈爆炸增长的，其结论对于实际操作意义不大。

4.卢卡斯定理

A、B是非负整数，$p$是质数。A、B写成p进制：$A=a[n]a[n-1]...a[0]，B=b[n]b[n-1]...b[0]$。则组合数$C(A,B)$与 $C(a[n],b[n])\ast C(a[n-1],b[n-1])\ast ...\ast C(a[0],b[0])$ mod p同余，即
$Lucas(n,m,p)=c(n\%p,m\%p)\ast Lucas(n/p,m/p,p)$
应用：（一道模板题）用于处理组合数取模，$C(n,m)$ mod p，其中p为素数。

附上模板题AC代码：

#include<bits/stdc++.h>
#define N 100010
using namespace std;
typedef long long ll;
ll a[N];  //代表N的阶乘mod p后的结果 
int p;
ll pow(ll y, int z, int p) {
	y %= p;
	ll ans = 1;
	for(int i=z; i; i>>=1, y=y*y%p)  //快速幂 
		if(i&1)  //奇数 
			ans = ans*y%p;
	return ans;
}

ll C(ll n, ll m) {
	if(m>n)	return 0;
	return ((a[n]*pow(a[m], p-2, p))%p * pow(a[n-m], p-2, p)%p);  //逆元，所以是乘法 
}

ll Lucas(ll n, ll m) {
	if(!m)	return 1;
	return C(n%p, m%p) * Lucas(n/p, m/p)%p;
}

int main() {
	int T, n, m;
	cin>>T; 
	while(T--) {
		cin>>n>>m>>p; 
		a[0] = 1;
		for(int i=1; i<=p; i++)
			a[i] = (a[i-1]*i)%p;
		cout<<Lucas(n+m, n)<<endl;
	}
	return 0; 
}

---

参考资料：
[1] https://zh.wikipedia.org/wiki/%E6%AC%A7%E6%8B%89%E5%AE%9A%E7%90%86_(%E6%95%B0%E8%AE%BA)