质数1：质数的线性筛，欧拉函数及其线性筛

前言：本人写这篇博客的主要目的是复习，所以这些筛法直接跳到 $O(n)$的了。

质数的线性筛

对于质数的筛法，传统的是暴力枚举（废话 ），现在我们讲一讲稍微精明一点的方法。我们可以考虑用质数的乘积进行筛除，由唯一分解定理，每个合数必定是有不止一个质数因子构成的。利用已找到的质数进行筛除，找出并标记合数。当搜到的当前质数没有被标记时，证明它没有更小的素因子，将这个数加入质数中。没错，这就是埃氏筛法。

int pcnt = 0;
memset( vis , 0 , sizeof(vis) );
memset( prime , 0 , sizeof(prime) );
for ( int i = 2 ; i <= n ; ++i ) {
	if ( !vis[i] ) {
		prime[++pcnt] = i;
	}
	for ( int j = 1 ; j <= pcnt ; ++j ) {
		if ( i * prime[j] > n ) {
			break;
		}
		vis[i * prime[j]] = 1;
	}
}

不难证明这个筛法的正确性，对于任意一个合数$x$，提出它的一个质因子为$p$，则$x$可表示成$p\ast t$，$t<x$，在循环的过程中一定有一种情况为$t>p$，在满足上一个条件时一定在循环到$t$的时候通过质数$p$筛出$x$。正确性得证。
不过，稍加思考就能发现，这种筛法会出现一个数被筛除多次的情况，那么，我们怎么避免这种情况呢？最好的办法就是让每个合数只会被自己最小的质因子筛除。具体实现也不算难，在自己代码一个合适的地方加上这一句即可：

if ( i % prime[j] == 0 ) {
	break;
}

具体原理是这样的，根据唯一分解定理，对于一个合数$n$，可表示成$$n=p_1^{k_1}\ast p_1^{k_1}\ast \cdot \cdot \cdot \ast p_m^{k_m}$$其中$$p_m>\cdot \cdot \cdot >p_2>p_1$$，当我们用$p_1$筛出这个数时，要防止$p_2$也会筛出这个数，就要想办法让$\frac{n}{p_2}$未到数$p_2$时就停止筛除，也即让这个数最多筛到$p_1$，最简单的检测方法就是查看某个数是否有质因子$p_1$，当发现它为质因子$p_1$的倍数时，停止筛除。
归纳上面的思路，就是每次循环得到一个数，将这个数与已筛到的质数用乘积的方法筛去合数，当当前数的因子含有循环到的质数时，做完最后一次筛除就退出质数的循环，处理下一个数。

int pcnt = 0;
memset( vis , 0 , sizeof(vis) );
memset( prime , 0 , sizeof(prime) );
for ( int i = 2 ; i <= n ; ++i ) {
	if ( !vis[i] ) {
		prime[++pcnt] = i;
	}
	for ( int j = 1 ; j <= pcnt ; ++j ) {
		if ( i * prime[j] > n ) {
			break;
		}
		vis[i * prime[j]] = 1;
		if ( i % prime[j] == 0 ) {//不要放错位置了，必须在它上一行代码后面
			break;
		}
	}
}

欧拉函数及其线性筛

欧拉函数：
表示为$\varphi (p)$，意义是1到$p$与$p$互质（最大公因数等于1）的数的个数（$\varphi (1)=1$）。

性质1：当$p$为质数时，$\varphi (p)=p-1$。
性质2：当$p$为质数时，且$i\%p=0$，那么$\varphi (i\ast p)=p\ast \varphi (i)$。
证明：首先，考虑两个数$n$和$i$，如果两个数不互质，设$GCD(n,i)=d$，则$n=d\ast x$，$i=d\ast y$，其中$x$，$y$互质，所以$n+i=d\ast (x+y)$，即$GCD(n+i,i)=GCD(n,i)$；同时也不难看出，若$i$与$n$互质，则$i+n$与$i$也互质。
所以，区间$[1,i]$与$i$不互质的数有$i-\varphi (i)$个，区间$(i+1,i+i]$与$i$不互质的数也有$i-\varphi (i)$个。以此类推，区间$[1,i\ast p]$与$i$不互质的数有$p\ast (i-\varphi (i))$个。又因为$i\%p=0$，$i\ast p$相对$i$不会有新的因数，所以，区间$[1,i\ast p]$与$i\ast p$不互质的数有$p\ast (i-\varphi (i))$个，即$\varphi (i\ast p)=p\ast \varphi (i)$。
性质3：当$p$为质数时，且$i\%p$不等于0，那么$\varphi (i\ast p)=(p-1)\ast \varphi (i)$。
证明：因为根据欧拉函数的积性性质，当$gcd(a,b)=1$时，$\varphi (a\ast b)=\varphi (a)\ast \varphi (b)$，对于质数$p$，$\varphi (p)=p-1$，所以性质3等价于欧拉函数的积性性质。
接下来，我们证明欧拉函数的积性性质，有兴趣的继续往下看：

在证明欧拉函数的积性性质前，我先给大家安利一个概念：完全剩余系。
如果$a$，$b$关于模$m$同余，那么$a$与$b$同属一类，否则不同属于一类，这样可得$m$个类，即
M i = { i + k ∗ m ∣ k ∈ Z } , i = 0 , 1 , 2 , ⋯ &ThinSpace; , m − 1 M_i=\{i + k*m | k \in Z \} , i = 0,1,2,\cdots ,m-1 Mi​={i+k∗m∣k∈Z},i=0,1,2,⋯,m−1
它们称为模$m$的剩余类。
从每个剩余类中各取一个数作为代表，这样得到的$m$个数称为模$m$的一个完全剩余系，简称完系。
完全剩余系还有一个性质，就是如果$GCD(n,m)=1$时，那么在$a_1,a_2,\cdots ,a_m$是模$m$的一个完全剩余系时，那么$n\ast a_1+k,n\ast a_2+k,\cdots ,n\ast a_m+k$，模$m$互不同余，也是一个完全剩余系。
证明：反证法，若$n\ast a_i+k\equiv n\ast a_j+k(\mathrm{m}\mathrm{o}\mathrm{d}m)$，又因为$GCD(n,m)=1$，不难得出$a_i\equiv a_j(\mathrm{m}\mathrm{o}\mathrm{d}m)$，与定义不成立，矛盾，故这个性质成立。
下面，回到对欧拉函数积性的证明：
设$GCD(a,b)=1$且$a,b\in Z$
$\varphi (a\ast b)$就是下面的$a$行$b$列的数表中与$a\ast b$互质的数的个数：
$$\begin{array}{cccccc}1& 2& \cdots & i& \cdots & b,\\ b+1& b+2& \cdots & b+i& \cdots & 2\ast b,\\ 2\ast b+1& 2\ast b+2& \cdots & 2\ast b+i& \cdots & 3\ast b,\\ ⋮& ⋮& \ddots & ⋮& \ddots & ⋮\\ (a-1)\ast b+1& (a-1)\ast b+2& \cdots & (a-1)\ast b+i& \cdots & a\ast b\end{array}$$
在第一行中有$\varphi (b)$个数与$b$互质，其余的数与$b$不互质，根据上面证明性质2的思路不难得出第一行不与$b$互质的数所在的列中的每个数都不与$b$互质，肯定也不与$a\ast b$互质，去掉这些列，不难得出剩下的$\varphi (b)$列，同样根据上面证明性质2的思路，剩下的数都与$b$互质。

$0,1,\cdots ,a-1$构成模$a$的一个完全剩余系，由于$a$与$b$互质，根据上面已证的完系的性质，得出$0\ast b+i,1\ast b+i,\cdots ,(a-1)\ast b+i$也是模$a$的一个完系。

与$b$一样，其中有$\varphi (a)$个数与$a$互质，按照$b$的思路，不难发现每列都只有$\varphi (a)$个数与$a$互质，而一共有$\varphi (b)$个列，所以总共有$\varphi (a)\ast \varphi (b)$个数既与$a$又与$b$互质，即与$a\ast b$互质，所以欧拉函数的积性得证。

再插一句：结合性质1与性质2，对于一个质数$p$，有$\varphi (p^a)=(p-1)\ast p^{a-1}$
再结合唯一分解定理以及欧拉函数的积性，不难得出对于任意一个合数若按照下面分解成质数幂乘积表达式$$n=p_1^{a_1}\ast p_2^{a_2}\ast \cdots \ast p_m^{a_m}$$，则可得出$$\varphi (n)=n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast \cdots \ast (1-\frac{1}{p_m})$$
这可用于快速求出一个数的欧拉函数值。

鉴于性质1，2，3均与质数有关，不难发现这可以套在质数线性筛里面，保证比1大的数都能被筛到，达到线性筛欧拉函数的效果。

int pcnt = 0;
memset( vis , 0 , sizeof(vis) );
memset( prime , 0 , sizeof(prime) );
phi[1] = 0;
for ( int i = 2 ; i <= n ; ++i ) {
	if ( !vis[i] ) {
		prime[++pcnt] = i;
		phi[i] = i - 1;//性质1 
	}
	for ( int j = 1 ; j <= pcnt ; ++j ) {
		if ( i * prime[j] > n ) {
			break;
		}
		vis[i * prime[j]] = 1;
		if ( i % prime[j] == 0 ) {
			phi[i * prime[j]] = phi[i] * prime[j];//性质2 
			break;
		} else {
			phi[i * prime[j]] = phi[i] * ( prime[j] - 1 );//性质3 
		}
	}
}

本文参考资料：《初等数论的知识与问题》（单壿 著）
《信息学竞赛之数学一本通》