信息学竞赛中可能用到的初等数论

竞赛中会用到的初等数论

0. 线性筛

$[1,n]$的范围内找质数，复杂度线性$O(n)$。

基于一个简单的想法：每个合数只在他的最小质因数删去。

我们要删去的数$x$,
$$x=p_0\times f_{ntl}$$
其中$p_0$为$x$的最小质因数，$f_{ntl}$为$x$的次大因数。

void GetPrime(int n)
{
	memset(isPrime, 1, sizeof(isPrime));
	isPrime[1] = 0;
	
	for(int i = 2; i <= n; i++) //穷举次大因数
	{
		if(isPrime[i])
			Prime[++cnt] = i; 
			
		for(int j = 1; j <= cnt && i*Prime[j] <= n; j++) //穷举最小质因数
		{
			isPrime[i*Prime[j]] = 0;
            
			if(i % Prime[j] == 0) break; //当i里有Prime[j]这个因子的时候，Prime[j+1]一定不是i*Prime[j+1]的最小质因数，所以终止，下个不用做了，因为不满足最初的"简单的想法"
		}
	}
}

后话：如果把$isPrime[i]$看成一个函数的话，表示i是不是一个质数，那么这里其实是在用递推求值。边界条件其实是$i$为质数的时候，函数的值为1。当然这个函数不满足积性（下一节会讲到）。

1. 欧拉函数

欧拉函数 $\varphi (x)$，表示的是范围$[1,x]$中与x互质的整数的个数。举几个例子，$\varphi (1)=0,\varphi (2)=1,\varphi (3)=2,\cdots$

容易观察到，对于任意一个质数$p$, $\varphi (p)=p-1$

接下来考虑$\varphi (p^k)$，由于$p^k$只有唯一的质因数p，所以一个与其不互质的数一定包含因子p。所以$[1,p^k]$中的p的倍数都与$p^k$不互质，其他都与其互质。这个范围内p的倍数总共有$p^{k-1}$个。因此，
$$\varphi (p^k)=p^k-p^{k-1}=(p-1)p^{k-1}$$

下面证明一下$\varphi (x)$具有一下性质，
对于$(x,y)=1$，有$\varphi (x\cdot y)=\varphi (x)\cdot \varphi (y)$

这里用Dirichlet卷积，对交换代数有兴趣的话，同样也可以验证环同态，这个放在Appendix里。

考虑集合 A d = { i ∣ ( x , i ) = d } ,   d ∣ x   即 d 为 x 的 一 个 因 数 。 A_d=\{i|(x,i)=d\},\ d|x\ 即d为x的一个因数。 Ad​={i∣(x,i)=d}, d∣x 即d为x的一个因数。
$$|A_d|=\varphi (x/d)$$
显然有上式（觉得不显然的可以考虑$i=k\cdot d$的系数k，满足条件时，有几种取值），

$$x=\sum _{d|x}\varphi (x/d)$$

写成Dirichlet卷积形式

$$1\ast \varphi (x)=\sum _{d|x}\varphi (x/d)=x$$

1的卷积逆是Mobius的$\mu$, 所以$\varphi (x)=\mu \ast x$，两个积性函数的卷积还是积性函数，所以欧拉函数也是积性函数。

有了积性和唯一分解定理可以写出欧拉函数的通项

x的唯一质因数分解：
$$x=p_1^{k_1}{p}_2^{k_2}\cdots p_n^{k_n}$$
$$then\varphi (x)=(p_1-1{)}^{k_1-1}(p_2-1{)}^{k_2-1}\cdots (p_n-1{)}^{k_n-1}$$

然后可以愉快的用这三个式子来做递推

$$\begin{array}{rlr}\varphi (x)& =x-1& x为质数\\ \varphi (i\times prime(j))=& \varphi (i)\times \varphi (prime[j])& (i,prime[j])=1\\ \varphi (i\times prime(j))=& \varphi (i)\times prime[j]& (i,prime[j])=prime[j]\end{array}$$

void GetPrime(int n)
{
	memset(isPrime, 1, sizeof(isPrime));
	isPrime[1] = 0;
	
	for(int i = 2; i <= n; i++) //穷举次大因数
	{
		if(isPrime[i]){
			Prime[++cnt] = i; 
			phi[i]=i-1;//i为质数
        }	
		for(int j = 1; j <= cnt && i*Prime[j] <= n; j++) //穷举最小质因数
		{
			isPrime[i*Prime[j]] = 0;
            phi[i*Prime[j]]=phi[i]*phi[Prime[j]];//i Prime[j]互质
			if(i % Prime[j] == 0) {
				phi[i*Prime[j]]=phi[i]*Prime[j];//i Prime[j]不互质
				break; 
			}
		}
	}
}

求单个数的欧拉函数可以用通项和因式分解求。
同样利用$|A_d|=\varphi (x/d)$
可以使用分块法解决一些求和问题，比如求${\sum }_{i=1}^n(i,n)={\sum }_{d|n}d\times \varphi (n/d)$
习题：P2158

2. 费马小定理和欧拉定理

i.欧拉定理

对于任意正整数a和x，$$a^{\varphi (x)}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}x$$

ii.费马小定理

对于任意正整数a和质数p，$$a^{p-1}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$$
证明放在附录里，用的是群论的Lagrange Theorem
费马小定理可以用来判断一个数质数，用快速幂可以达到$O(\lg n)$的复杂度，也可以用来求一个数关于p的数论倒数
$$a^{(p-2)}a\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$$
$$a^{(p-2)}\equiv a^{-1}\mathrm{m}\mathrm{o}\mathrm{d}p$$

3. 裴蜀定理

两个整数$a,b$, 最大公约数为$(a,b)=d$，$\exists$整数$x,y$, 满足$$ax+by=d;$$
$\forall$整数$x,y$,
$$ax+by\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}d$$
证明很简单，这里就不展开讲了 有兴趣的同学自己去证明一下。
下面一节会去解决这个求这个a, b的问题。

4. 辗转相除和扩展欧几里得以及求逆元

辗转相除法可以看成一个递归, 递归式是$$(a,b)=(b,a\%b)$$，边界条件为$(a,0)=a$.
正确性很容易证明，复杂度在1844年被一个法国人计算出来，结果是不会大于$5\lg min(a,b)$，标志着算法复杂度理论的开端。

这里提供一种构造Fibonacci的夹逼证明，
Proof
Statement 1: 如果一组a, b在n步辗转相除计算出结果，那么有$$\begin{gathered} a>=f_{n+2} \\ b>=f_{n+1} \end{gathered}$$
$f_n$为Fibonacci数列的第n项。
这里可以用数学归纳法证明：
当第n步满足考虑(n+1)时，
$$a_{n+1}=\lfloor \frac{a_{n+1}}{b_{n+1}}\rfloor a_n+b_n>=\lfloor \frac{a_{n+1}}{b_{n+1}}\rfloor f_{n+2}+f_{n+1}>=f_{n+2}+f_{n+1}=f_{n+3}$$
下面只需要考虑Fibonacci的第n项是多少就可以了，如果学过矩阵乘法或者求矩阵特征向量特征根的话那就很容易了，如果没有的话可以去学一学，这里直接给出通项。
$$f_n=\frac{(\frac{1+\sqrt{5}}{2}{)}^n-(\frac{1-\sqrt{5}}{2}{)}^n}{\sqrt{5}}$$
显然 $f_n$和$n$是指数关系，即$f_n=a^n,a>1$
所以由$b_n\approx f_{n+1}$可以得到，$$n\approx {\log }_{a}b$$

以上为辗转相除的复习。

下面考虑这样一个问题，即解决裴蜀定理，已知a，b
$$ax+by=gcd(a,b)$$
求一组整数解x, y。

考虑辗转相除法：
边界条件：$a{x}^{(n)}+0\cdot y^{(n)}=gcd(a,0)=a$
显然有一组解：$x^{(n)}=1,y^{(n)}=0$

第一层：$a{x}^{(k)}+b{y}^{(k)}=gcd(a,b)$
第二层：$b{x}^{(k+1)}+r{y}^{(k+1)}=gcd(b,r),a\equiv r\mathrm{m}\mathrm{o}\mathrm{d}b$
有$gcd(a,b)=gcd(b,r)$, 和上面两式结合有：$$a{x}^{(k)}+b{y}^{(k)}=b{x}^{(k+1)}+r{y}^{(k+1)}$$
$$r=a-\lfloor \frac{a}{b}\rfloor b$$
$$a{x}^{(k)}+b{y}^{(k)}=b{x}^{(k+1)}+(a-\lfloor \frac{a}{b}\rfloor b)\cdot y^{(k+1)}$$
上式左右两边相等可以得到：
$$\begin{array}{rl}{x}^{(k)}& =y^{(k+1)}\\ y^{(k)}& =x^{(k+1)}-\lfloor \frac{a}{b}\rfloor y^{(k+1)}\end{array}$$
于是就得到了一个递推式，结合边界条件就可以得到解。

void Ex_gcd(int a, int b, int &x, int &y)
{
    if(b == 0)
    {
        x = 1;
        y = 0;
        return;
    }
    int x1, y1;
    Ex_gcd(b, a%b, x1, y1);
    x = y1;
    y = x1-(a/b)*y1;
}

然后求到$x^{(1)},y^{(1)}$，之后可以通过加上若干个$ap+bq=0,a>0$，那么x的最小正整数解就是$x\equiv x^{(a)}\mathrm{m}\mathrm{o}\mathrm{d}p$。

5. 中国剩余定理

韩信点兵：3个人一排剩2个，5个人一排剩3个，7个人一排剩2个，问有几个人？
答案233个或23个$\cdots$！
问题的形式化：
$$\begin{array}{rl}x& \equiv a_1\mathrm{m}\mathrm{o}\mathrm{d}{m}_1\\ x& \equiv a_2\mathrm{m}\mathrm{o}\mathrm{d}{m}_2\\ & \dots \\ x& \equiv a_k\mathrm{m}\mathrm{o}\mathrm{d}{m}_k\end{array}$$
任意两个$m_i,m_j$, 需要满足$(m_i,m_j)=1$。

构造一个$M=m_1{m}_2\dots m_k$
再定义一个
$$\begin{array}{rl}{M}_1& =\frac{M}{m_1}\\ M_2& =\frac{M}{m_2}\\ & \cdots \\ M_k& =\frac{M}{m_k}\end{array}$$
对于$M_i$, 可以求一个$t_i$，满足$M_i{t}_i\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}{m}_i$。
这时候就构造出来了！$$x=\sum _{i=1}^k{a}_i{t}_i{M}_i$$
$x$显然是上述同余方程组的解。

其中唯一有问题的就是其中求$t_i$，就是求一个整数a关于p的倒数，现有$(a,p)=1$,
$$(a\cdot x)\%p=1$$可以转化成
$$a\cdot x+p\cdot y=1=(a,p)$$
于是就可以用扩展欧几里得算法来求$x$和$y$了，当然我们只需要求$x$。

6. 威尔逊定理

$$p是质数\iff (p-1)!\equiv p-1\mathrm{m}\mathrm{o}\mathrm{d}p$$
必要性显然
充分性给一个用费马小定理的证明，$p\ge 3$且p是质数：
$$\begin{array}{rl}f(x)& =(x-1)\cdot (x-2)\dots (x-(p-1))\\ & =x^{p-1}+(1+2+\cdots +p-1)x^{p-2}+\cdots +(p-1)!\end{array}$$
注意这里满足f(x)=0的根为$1\sim (p-1)$

定义$$g(x)=x^{p-1}-1$$
根据费马小定理，满足$$g(x)\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p$$的解也是$1\sim (p-1)$
所以对于任意的$x$，满足$$h(x)=f(x)-g(x)\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p$$

但是考虑到$h(x)$的最高此项是$x^{(}p-2)$这一项，所以方程$$h(x)\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p$$最多有$(p-2)$个解，因此$h(x)$只能恒等于0。

所以他的常数项 $(p-1)!+1$必须满足$$(p-1)!+1\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}p$$，这个式子就是威尔逊定理。
题目：hdu2973

Appendix

施工中$\dots$