【数论】浅谈线性筛

前置知识

线性筛质数
积性函数定义：若一个定义域在正整数上的函数 $f(x)$ 满足 $\forall x,y\in N_{+}\text{且}\gcd (x,y)=1,f(xy)=f(x)f(y)$
凡是积性函数，都可以用线性筛。

剖析线性筛质数
回顾整个算法，我们发现每个正整数都只会被它的最小质因子筛掉，故整个算法的时间复杂度为 $O(n)$。
这里给出参考代码：

void ola(int x){
	pr[1]=1;
	for(int i=1;i<=x;i++){
		if(!pr[i])
			prv[++cnt]=i;
		for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
			pr[i*prv[j]]=1;
			if(i%prv[j]==0)
				break;
		}
	}
}

于是我们就可以以此为基础筛一些积性函数了。

欧拉函数

定义
欧拉函数 $\phi (n)$ 表示小于等于 $n$，且与 $n$ 互质的正整数的个数。特别的，$\phi (1)=1$

求单个欧拉函数
设 $n={\prod }_{i=1}^q{p}_i^{x_i}$(质因数分解，如无特殊说明，下文此形式均为质因数分解)，集合 $A_i$ 表示不大于 $n$ 且是 $p_i$ 的倍数的自然数的集合。

所有 $|A_i|=\frac{n}{p_i}$
因为 $\gcd (p_i,p_j)=1(i\ne j)$
所以 $|A_i\cap A_j|=\frac{n}{p_i{p}_j}$
再结合容斥：

参考代码：

int getphi(int x){
	int res=x;
	for(int i=2;i*i<=x;i++){
		if(x%i==0)
			res=res*(i-1)/i;
		while(x%i==0)
			x/=i;
	}
	if(x>1)
		res=res*(x-1)/x;
	return res;
}

欧拉函数是积性函数
证明：
$$\phi (n)\phi (m)=n\prod _{i=1}^x(1-\frac{1}{p_i})m\prod _{i=x+1}^{x+y}(1-\frac{1}{p_i})=\phi (nm)$$
由于 $\gcd (n,m)=1$，所以 $n,m$的每个质因子必不相同，所以得证。

性质

1. 若 $n$ 是质数，$\phi (n)=n-1$
2. 若 $p$ 为质数，且 $n=p^k$，则 $\phi (n)=(p-1)p^{k-1}$
3. 若 $q$ 为质数，且 $\exists k\in N_{+}$ 使得 $n=qk$，则 $\phi (nq)=q\phi (n)$
4. 若 $q$ 为质数，且 $\forall k\in N_{+}$，$n\ne qk$，则 $\phi (nq)=(q-1)\phi (n)$
5. $n={\sum }_{d|n}\phi (d)$
6. $\phi (ab)=\phi (a)\phi (b)\frac{d}{\phi (d)}$，其中 $d=\gcd (a,b)$

证明

1. 不证
2. $\phi (n)=n(1-\frac{1}{p})=p^k\frac{p-1}{p}=(p-1)p^{k-1}$
3. $\phi (nq)=nq{\prod }_{i=1}^x(1-\frac{1}{p_i})=qn{\prod }_{i=1}^x(1-\frac{1}{p_i})=q\phi (n)$
4. $n={\prod }_{i=1}^y{p}_i^{x_i},\phi (nq)=nq(1-\frac{1}{q}){\prod }_{i=1}^y(1-\frac{1}{p_i})=q\frac{q-1}{q}[n{\prod }_{i=1}^y(1-\frac{1}{p_i})]=(q-1)\phi (n)$
5. 如图
   
   理性理解：$n={\sum }_{d|n}{\sum }_{i=1}^n[\gcd (i,n)=d]={\sum }_{d|n}{\sum }_{i=1}^{\frac{n}{d}}[\gcd (i,\frac{n}{d})=1]={\sum }_{d|n}\phi (\frac{n}{d})={\sum }_{d|n}\phi (d)$

$$\begin{gathered} d=\prod _{i=1}^{q_1}{p_{1i}}^{x_{1i}},a=\prod _{i=1}^{q_2}{p_{2i}}^{x_{2i}},b=\prod _{i=1}^{q_3}{p_{3i}}^{x_{3i}} \\ \phi (ab)=\frac{ab{\prod }_{i=1}^{q_2}(1-\frac{1}{p_{2i}}){\prod }_{i=1}^{q_3}(1-\frac{1}{p_{3i}})}{{\prod }_{i=1}^{q_1}(1-\frac{1}{p_{1i}})}=\frac{\phi (a)\phi (b)}{\frac{1}{d}\phi (d)}=\phi (a)\phi (b)\frac{d}{\phi (d)} \end{gathered}$$

线性筛欧拉函数

回顾性质 $1,3,4$，也就是我们筛欧拉函数的基础了，即我们在筛质数的时候顺便把欧拉函数也给筛出来。
参考代码：

void ola(int x){
	pr[1]=phi[1]=1;
	for(int i=1;i<=x;i++){
		if(!pr[i])
			prv[++cnt]=i,phi[i]=i-1;
		for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
			int u=i*prv[j];
			pr[u]=1;
			if(i%prv[j]==0){
				phi[u]=phi[i]*prv[j];
				break;
			}
			else{
				phi[u]=phi[i]*(prv[j]-1);
			}
		}
	}
}

$d(x)$

设 $d(x)$ 表示 $x$ 的约数的个数。

1. $d(x)$ 是积性函数。
   证明：
   $$\begin{gathered} A=\prod _{i=1}^n{p}_i^{x_i},B=\prod _{i=1}^m{q}_i^{y_i},n=AB,gcd(A,B)=1 \\ d(A)=\prod _{i=1}^n\sum _{j=0}^{x_i}{p}_i^j,d(B)=\prod _{i=1}^m\sum _{j=0}^{y_i}{q}_i^j \\ \because \gcd (A,B)=1,\therefore n=\prod _{i=1}^n{p}_i^{x_i}\prod _{i=1}^m{q}_i^{y_i} \\ \therefore d(n)=\prod _{i=1}^n\sum _{j=0}^{x_i}{p}_i^j\prod _{i=1}^m\sum _{j=0}^{y_i}{q}_i^j \\ \therefore d(n)=d(A)d(B) \end{gathered}$$

2. 线性筛 $d(x)$：
   

• 参考代码：

void ola(int x){
	pr[1]=d[1]=1;
	for(int i=1;i<=x;i++){
		if(!pr[i])
			prv[++cnt]=i,d[i]=2,mx[i]=i,mi[i]=1;
		for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
			int u=i*prv[j];
			pr[u]=1;
			if(i%prv[j]==0){
				mi[u]=mi[i]+1;
				mx[u]=mx[i]*prv[j];
				d[u]=(mx[u]==u?mi[u]+1:d[mx[u]]*d[u/mx[u]]);
				break;
			}
			else{
				mi[u]=1;
				mx[u]=prv[j];
				d[u]=d[mx[u]]*d[u/mx[u]];
			}
		}
	}
}