常见数论函数模板

常见数论函数

1）$u(n),e(n),I(n)$

$u(n)\equiv 1,n\ge 1$

$e(n)=n,n\ge 1$

$I(n)=\{\begin{array}{l}1,n=1\\ 0,n>1\end{array}$

2）$d(n)$

$n$的所有正除数的个数，含1和$n$本身：$d(n)={\sum }_{d\mid n}1$

完全分解： n = p 1 k 1 p 2 k 2 ⋯ p s k s ,   { n > 1 } n=p_1^{k_1}p_2^{k_2}\cdots p_s^{k_s},\ \{n>1\} n=p1k1​​p2k2​​⋯psks​​, {n>1} ，$d(1)=1$

$d(n)={\prod }_{i=1}^s(k_i+1)$

可乘函数可使用线性筛：

$p$存放质数，用d表示约数个数，num表示最小质因子出现次数。考虑$n>1$。

$d(p)=2,num(p)=1$

设$d=\frac{n}{p}$， $p$为$n$的最小质因子：

$p\mid d:num(n)=num(d)+1,d(n)=\frac{d(d)\times (num(n)+1)}{num(d)+1}$

考虑$d,p$互质，因为为可乘函数：

$p\nmid d:num(n)=1,d(n)=d(d)\times d(p)$

//数组 d, num, p
//bool v   int tot, n
void pre(){
	d[1]=1, tot=0;
	for(int i=2; i<=n; i++){
		if(!v[i]) p[++tot]=i, d[i]=2, num[i]=1;
		for(int j=1; j<=tot && i<=n/p[j]; j++){
			v[p[j]*i]=1;
			if(i%p[j]==0){
				num[i*p[j]]=num[i]+1;
				d[i*p[j]]=d[i]/(num[i]+1)*(num[i*p[j]]+1);
				break;
			}
			num[i*p[j]]=1;
			d[i*p[j]]=d[i]*d[p[j]];
		}
	}
}

3）$\Omega (n)$

$n$的全部素因子的个数$\Omega (n)$。

$\Omega (1)=0$

$\Omega (n)=$${\sum }_{i=1}^s{a}_i$

int b[maxn];
int cnt[maxn];
void pre() {
	for(int i=2; i<maxn; i++) {
		if(b[i]!=0) continue;
		for(int j=i; j<maxn; j+=i) {
			if(b[j]==0) b[j]=i;
		}
	}
	cnt[2]=1;
	for(int i=3; i<maxn; i++) 
		cnt[i]=cnt[i/b[i]]+1;
	return;
}

4）$\omega (n)$

$n$的不同素因子个数$\omega (n)$。

$\omega (1)=0,\omega (n)=s$

bool vis[maxn];
int cnt[maxn];
void pre() {
	for(int i=2; i<maxn; i++) {
		if(vis[i]!=0) continue;
		cnt[i]=1;
		for(int j=i+i; j<maxn; j+=i) {
			cnt[j]++, vis[j]=1;
		}
	}
	return;
}

5）${\sigma }_{\lambda }(n)$

$n$的正除数的幂和函数${\sigma }_{\lambda }(n),\lambda$为实数。

${\sigma }_{\lambda }(n)={\sum }_{d\mid n}{d}^{\lambda }$

特殊的，约数和：$\lambda =1$时，$\sigma (n)={\sum }_{d\mid n}d$

完全分解： n = p 1 k 1 p 2 k 2 ⋯ p s k s ,   { n > 1 } n=p_1^{k_1}p_2^{k_2}\cdots p_s^{k_s},\ \{n>1\} n=p1k1​​p2k2​​⋯psks​​, {n>1} ，$d(1)=1$

$\sigma (n)={\prod }_{i=1}^s(1+p_i+p_i^2+\cdots +p_i^{k_i})={\prod }_{i=1}^s{\sum }_{j=0}^{k_i}{p}_i^j$

又因：${\sum }_{j=0}^{k_i}{p}_i^j=\frac{p_i^{k_i+1}-1}{p_i-1}$

$\sigma (n)={\prod }_{i=1}^s(1+p_i+p_i^2+\cdots +p_i^{k_i})={\prod }_{i=1}^s\frac{p_i^{k_i+1}-1}{p_i-1}$

以下计算均采用前者推导。

可乘函数可使用线性筛：

$p$存放质数，$f_i$表示$i$的约数和，$g_i$表示$i$的最小质因子$p$的${\sum }_{j=0}^k{p}^j$ 。考虑$n>1$。

$g(p)=p+1,f(p)=i+1$

设$d=\frac{n}{p}$，$p$为$n$的最小质因子：

$p\mid d:g(n)=g(d)\times p+1,f(n)=\frac{f(d)\times g(n)}{g(d)}$

考虑$d,p$互质，因为为可乘函数：

$p\nmid d:g(n)=p+1,f(n)=f(d)\times f(p)$

void pre(){
	g[1]=f[1]=1, tot=0;
	for(int i=2; i<=n; i++){
		if(!v[i]) p[++tot]=i, g[i]=i+1, f[i]=i+1;
		for(int j=1; j<=tot && i<=n/p[j]; j++){
			v[p[j]*i]=1;
			if(i%p[j]==0){
				g[i*p[j]]=g[i]*p[j]+1;
				f[i*p[j]]=f[i]/g[i]*g[i*p[j]];
				break;
			}
            g[i*p[j]]=p[j]+1;
			f[i*p[j]]=f[i]*f[p[j]];
		}
	}
}

注：求$\sigma (n)$是否为奇数，根据第一个公式，$p_i=2$时，$k_i$为任意数可，对其他$p_i,k_i$必须为偶数，所以$n$为平方数或$n/2$为平方数为充要条件。求$1$到$n$中，$\sigma (i)$为奇数的个数，$ans=sqrt(n)+sqrt(n/2)$即可。

6）$\phi (n)$

欧拉函数 $\phi (n)$：小于等于$n$，和$n$互质的个数。

1与任何正整数互质，$\phi (1)=1$

求解：完全分解：$n=p_1^{k_1}{p}_2^{k_2}\cdots p_s^{k_s}$

推导表达式，由容斥原理：$\phi (n)=n+(-\frac{n}{p_1}-\frac{n}{p_2}-\cdots -\frac{n}{p_s})+(\frac{n}{p_1{p}_2}+\frac{n}{p_1{p}_3}+\cdots +\frac{n}{p_{s-1}{p}_s})\dots$ φ ( n ) = ∑ S ⊆ { p 1 , p 2 … p s } ( − 1 ) ∣ S ∣ n ∏ p i ∈ S p i = n ∏ i = 1 s ( 1 − 1 p i ) \varphi(n)=\sum_{S\subseteq\{p_1,p_2\dots p_s\}}(-1)^{\mid S\mid}\frac{n}{\prod_{p_i\in S} p_i}=n\prod_{i=1}^s(1-\frac{1}{p_i}) φ(n)=∑S⊆{p1​,p2​…ps​}​(−1)∣S∣∏pi​∈S​pi​n​=n∏i=1s​(1−pi​1​)

$\phi (n)=n\times {\prod }_{i=1}^s\frac{p_i-1}{p_i}$

int phi(int n) {
	int m=(int)sqrt(n+0.5);
	int ans=n;
	for (int i=2; i<=m; i++)
		if (n%i==0) {
			ans=ans/i*(i-1);
			while(n%i==0) n/=i;
		}
	if (n>1) ans=ans/n*(n-1);
	return ans;
}

性质：

（1）为可乘函数，$gcd(a,b)=1$：$\phi (a\times b)=\phi (a)\times \phi (b)$。

（2）特殊的，$n$为奇数：$\phi (2n)=\phi (n)$，$\because 2\times \frac{2-1}{2}=1$

（3）$n={\sum }_{d\mid n}\phi (d)$ ，含$1$和$n$，证明：

法一： 要证明$e(n)=n(n\ge 1)$为$\phi (n)$的莫比乌斯变换，

只需证明$\phi (n)$为$e(n)$的反莫比乌斯变换即可。

$e(n)$的反莫比乌斯变换：$\phi (n)={\sum }_{d\mid n}\mu (d)\frac{n}{d}$

$\phi (n)$从定义式正经变换：$\phi (n)={\sum }_{1\le d\le n,(d,n)=1}1={\sum }_{1\le d\le n}{\sum }_{l\mid (d,n)}\mu (l)={\sum }_{l\mid n}\mu (l){\sum }_{1\le d\le n,l\mid d}1={\sum }_{l\mid n}\mu (l)\frac{n}{l}$

可见上下两式相等，故$\phi (n)$为$e(n)$的反莫比乌斯变换。

法二：$gcd(a,n)=d$，按$d$对$a$进行分类，令$n=dk$则$gcd(\frac{a}{d},k)=1$。固定的$d$对应$a$的个数为

$\phi (k)=\phi (\frac{n}{d})$，所以$n={\sum }_{d\mid n}\phi (\frac{n}{d})={\sum }_{d\mid n}\phi (d)$ 。

（4）特殊的$n=p^k$：$\phi (n)=p^k-p^{k-1}$

$=(p-1)\times p^{k-1}=n\times \frac{p-1}{p}$，由定义可得。

可用于求${\sum }_{i=1}^n{\sum }_{j=i+1}^{n}gcd(i,j)$ ：

显然$ans[n]=ans[n-1]+{\sum }_{i=1}^{n-1}gcd(i,n)$，仅需证明此处${\sum }_{i=1}^{n-1}gcd(i,n)={\sum }_{d\mid n}\phi (\frac{n}{d})\times d$ ：

按$d$对$a$进行分类，令$n=dk$则$gcd(\frac{a}{d},k)=1$。固定的$d$对应$a$的个数为

$\phi (k)=\phi (\frac{n}{d})$，和为$\phi (\frac{n}{d})\times d$ 。$f(n)={\sum }_{d\mid n}\phi (\frac{n}{d})\times d=(\phi \ast e)(n)$，解题关键是$f(n)$的表达式也可用线性筛，所以不会超时。因为此题$d<n$而不是$d\le n$，故内层循环从$i\ast 2$开始，而不是从$i$开始。$(\phi \ast e)(n)$筛法：

for(int i=1; i<=n; i++){
	for(int j=i+i; j<=n; j+=i){
		f[j]+=(ll)i*(ll)phi[j/i];
	}
}

筛法：

void phi_table(int n){
	memset(phi, 0, sizeof(phi));
	phi[1]=1;
	for(int i=2; i<=n; i++)
		if(!phi[i])
			for(int j=i; j<=n; j+=i){
				if(!phi[j]) phi[j]=j;
				phi[j]=phi[j]/i*(i-1);
			}
}

7）$\mu (n)$

莫比乌斯函数$\mu (n)$

$\mu (n)=\{\begin{array}{ll}1,& n=1\\ (-1{)}^s,& n=p_1{p}_2\dots p_s,p_1<p_2<\cdots <p_s\\ 0,& 其他\end{array}$

线性筛：

//mu p v
//tot
void pre() {
	mu[1]=1;
	for(int i=2; i<=maxn-13; i++) {//注意这里的上界
		if(!v[i]) mu[i]=-1, p[++tot]=i;
		for(int j=1; j<=tot && i*p[j]<=maxn-13; j++) {//改上界
			v[i*p[j]]=1;
			if(i%p[j]==0) {
				mu[i*p[j]]=0;
				break;
			}
			mu[i*p[j]]=-mu[i];
		}
	}
}

8）$\Lambda (n)$

Mangoldt函数$\Lambda (n)$：

$\Lambda (n)=\{\begin{array}{ll}\log p,& n=p^k,k\ge 1\\ 0& 其他\end{array}$

9）$\lambda (n)$

Liouville函数$\lambda (n)$：

$\lambda (n)=(-1{)}^{\Omega (n)}$

可乘函数（积性函数）

$f(mn)=f(m)\times f(n),\gcd (m,n)=1$

欧拉函数$\phi (n)$，莫比乌斯函数$\mu (n)$，正因子和$\sigma (n)$，正因子数$d(n)$等。

均可采用线性筛法。

可乘函数的卷积仍是可乘函数，可用线性筛：

求$f=g\ast h$预先计算$g,h$后，如下：

for(int i=1; i<=n; i++){
	for(int j=i+i; j<=n; j+=i){
		f[j]+=g[i]*h[j/i];
	}
}

完全（绝对）可乘函数

$f(mn)=f(m)f(n)$

$u(n),e(n),I(n)$