积性函数（欧拉，莫比乌斯函数），狄利克雷卷积，莫比乌斯反演，杜教筛

今天终于能将这几者的关系依自己的理解好好写一写了。

积性函数

积性函数：对于任意互质的整数a，b有f(a×b) = f(a)×f(b)的数论函数。
完全积性函数：对于任意的整数a，b有f(a×b) = f(a) × f(b)。
常见的积性函数：
$$\begin{gathered} \mu :莫比乌斯函数 \\ \phi :欧拉函数 \\ ϵ:单位元函数 \\ id:单位函数 \\ I:恒等函数 \\ \sigma :约数和函数 \\ d:约数个数函数（除数函数） \end{gathered}$$

• 较常用的积性函数：
• 

狄利克雷卷积

• 先定义数论函数的的加法和数乘：

加法：$(F+G)(i)=F(i)+G(i)$;数乘：$(xF)(i)=xF(i)$

其中F+G表示新函数H(i) = F(i) + G(i)

两个数论函数$F=G$的充要条件是$\forall n\in {\mathbb{N}}^{+}，F(n)=G(n)$。

• 定义数论函数的狄利克雷卷积*如下：

设$G(n),F(n)$是两个数论函数,则其狄利克雷卷积为:$$T(n)=F(n)\ast G(n)=\sum _{d|n}F(n)G(\frac{n}{d})$$

等价形式：$T(n)={\sum }_{i\times j}F(i)G(j)$，这个形式常用于证明某些性质。

• 定义狄利克雷卷积的逆元：

对于$F(1)\ne 0$,都存在$F$*$G=ϵ$，称$G$是$F$的逆元，记作:$F^{-1}$

$F(1)\ne 0$的原因是:$(F\ast F^{-1})(1)={\sum }_{ij=1}F(i)F(j{)}^{-1}=F(1)F(1{)}^{-1}=1$，因此$F(1)\ne 0$。

• 可得${\mu }^{-1}=I$

证: $\mu \ast I={\sum }_{d|n}\mu (\frac{n}{d})={\sum }_{d|n}\mu (d)=[(n==1)]=ϵ$

• 狄利克雷卷积的性质：

设$F,G$是两个积性函数,有:

• 单位函数的运算:$ϵ\ast F=F$
• $F\ast G$是积性函数
• $F^{-1}$也是积性函数
• 交换律:$F\ast G=G\ast F$
• 结合律:$(F\ast G)\ast H=F\ast (G\ast H)$
• 分配率:$(F+G)\ast H=F\ast H+G\ast H$

莫比乌斯反演

可由狄利克雷卷积得到:

• $G(n)={\sum }_{d|n}F(d)\Rightarrow F(n)={\sum }_{d|n}G(d)\mu (\frac{n}{d})$

因为有$G=F\ast I(\Leftarrow G(n)={\sum }_{d|n}F(d))$
又因为$F=F\ast ϵ,ϵ=I\ast \mu$,
可得:$F=F\ast ϵ=F\ast I\ast \mu =G\ast u={\sum }_{d|n}G(d)u(\frac{n}{d})$

• $G(d)={\sum }_{d|n}F(n)\Rightarrow F(d)={\sum }_{d|n}G(n)\mu (\frac{n}{d})$

也可通过上面得证,具体见:数学/数论专题-学习笔记：莫比乌斯反演

欧拉函数与莫比乌斯函数的关系

• $$\begin{gathered} \phi =id\ast u \\ ={\sum }_{d|n}d\mu (\frac{n}{d}) \\ ={\sum }_{d|n}\frac{n}{d}\mu (d) \\ (两边除n) \\ \Rightarrow \frac{\phi (n)}{n}={\sum }_{d|n}\frac{\mu (d)}{d} \\ (个人感觉还是\phi (n)={\sum }_{d|n}d\mu (\frac{n}{d})好记) \end{gathered}$$
• $$\begin{gathered} \phi =id\ast \mu \\ id=\phi \ast {\mu }^{-1}=\phi \ast I \\ \Rightarrow n={\sum }_{d|n}\phi (d) \\ \Rightarrow \phi (n)=n-{\sum }_{d|n,d<n}\phi (d) \end{gathered}$$

总结:$\phi (n)={\sum }_{d|n}d\mu (\frac{n}{d})=n-{\sum }_{d|n,d<n}\phi (d)$
and$n={\sum }_{d|n}\phi (d)$

杜教筛

用途：可以在低于线性的时间复杂度内求解积性函数的前缀和。
时间复杂度：$O(n^{\frac{2}{3}})$

证明：P4213题解有杜教筛时间复杂度证明

• 公式(具体推导见数论学习笔记2之杜教筛初探（含例题练习）)：
  $$\begin{gathered} 令h=f\ast g，G为要求的函数 \\ 适用条件：\sum _{i=1}^{n}h(i),\sum _{i=1}^{n}f(i)好求 \\ 则\sum _{i=1}^{n}g(i)就可套该公式求出 \\ 令:H(n)=\sum _{i=1}^{n}h(i)=\sum _{i=1}^n(f\ast g)(i) \\ 令:S(n)=\sum _{i=1}^{n}g(i) \\ 可得:S(n)=\frac{H(n)-{\sum }_{d=2}^{n}f(d)g(\frac{n}{d}) \\ }{f(1)} \end{gathered}$$

• 应用

• 求${\sum }_{i=1}^n\phi (i)$

因为:$id=\phi \ast I$
则令$g=\phi ,f=I,h=id$
带入上式得:
$$\begin{gathered} S(n)=\frac{H(n)-{\sum }_{d=2}^{n}f(d)g(\frac{n}{d}) \\ }{f(1)} \\ =\frac{\frac{n(n+1)}{2}-S(\frac{n}{d})}{1} \\ =\frac{n(n+1)}{2}-S(\frac{n}{d}) \end{gathered}$$

• 求${\sum }_{i=1}^n\mu (i)$

因为:$ϵ=\mu \ast I$
则令$g=\mu ,f=I,h=ϵ$
带入上式得:
$$\begin{gathered} S(n)=\frac{H(n)-{\sum }_{d=2}^{n}f(d)g(\frac{n}{d}) \\ }{f(1)} \\ =\frac{1-S(\frac{n}{d})}{1} \\ =1-S(\frac{n}{d}) \end{gathered}$$

Code板子

线性筛法求积性函数

1. 积性函数
//d(n):n正因子个数 -- d(nm) = d(n) * d(m)；n，m需要互质;
const int N = 1e7 + 10;
int primes[N], tot, d[N], cnt[N];
bool st[N];

int calc_d(int x, int k)
{
	int res = 1;
	while(k)
	{
		if(k & 1) res = res * x;
		x *= x;
		k >>= 1;
	}
	
	int ans = 0;
	for(int i = 1; i <= res / i; i ++)
		if(res % i == 0) 
		{
			ans ++;
			if(res / i != i) ans ++;
		}
	return ans;
}

void init()
{
	d[1] = 1;
	for(int i = 2; i < N; i ++)
	{
		if(!st[i]) primes[tot ++] = i, d[i] = 2, cnt[i] = 1;
		for(int j = 0; primes[j] * i < N; j ++)
		{
			st[primes[j] * i] = true;
			if(i % primes[j] == 0)//pj之前已经是i的质因子了
			{
				cnt[i * primes[j]] = cnt[i] + 1;
				//先将pj对d[i]的贡献去掉，此时i(去掉pj贡献)与pj互质，可以使用积性函数的性质求d[i*pj]
				d[i * primes[j]] = d[i] / calc_d(primes[j], cnt[i]) * calc_d(primes[j], cnt[i] + 1);
				break;
			}
			cnt[i * primes[j]] = 1;
			d[i * primes[j]] = d[i] * 2;
		}
	}
}

2.完全积性函数
/*
f(n) = n ^ N --> f(nm) = (nm)^ N = n^N * m^N = f(n)(m)
故f(n)为积性函数（n,m不需要互质）

int(1e6) = 4M char(1e6) = 1M 
*/

const int N = 2e7, mod = 1e9 + 7;
LL primes[N], tot, f[N], n;
bool st[N];

LL calc_f(LL x)
{
	LL res = 1, k = n;
	while(k)
	{
		if(k & 1) res = res * x % mod;
		x = (LL)x * x % mod;
		k >>= 1;
	}
	return res % mod;
}

void init()
{
	f[1] = 1;
	for(int i = 2; i <= n; i ++)
	{
		if(!st[i]) primes[tot ++] = i, f[i] = calc_f(i);
		for(int j = 0; (LL)primes[j] * i <= n; j ++)
		{
			st[i * primes[j]] = true;
			//f(nm) = f(n)f(m)不需要互质，就不用分开计算f(i * pj)
			f[i * primes[j]] = (LL)f[i] * f[primes[j]] % mod;
			
			if(i % primes[j] == 0) break;
		}
	}
}

杜教筛板子

const int N = 1e7 + 10;
LL primes[N], tot, mu[N], phi[N];//积性函数都可以用线性筛算
bool st[N];
unordered_map<int, LL> ans_mu, ans_phi;

void init()
{
	mu[1] = phi[1] = 1;
	for(int i = 2; i < N; i ++) 
	{
		if(!st[i]) primes[tot ++] = i, mu[i] = -1, phi[i] = i - 1;
		for(int j = 0; primes[j] * i < N; j ++)
		{
			st[i * primes[j]] = true;
			if(i % primes[j] == 0) 
			{
				phi[i * primes[j]] = phi[i] * primes[j];
				mu[i * primes[j]] = 0;
				break;
			}
			mu[i * primes[j]] = -mu[i];
			phi[i * primes[j]] = phi[i] * (primes[j] - 1);
		}
	}
	for(int i = 2; i < N; i ++) mu[i] += mu[i - 1], phi[i] += phi[i - 1];
}

LL calc_mu(LL x)//求莫比乌斯函数前缀和
{
	if(x < N) return mu[x];//预处理
	if(ans_mu[x]) return ans_mu[x];//记忆化
	
	LL res = 1;
	for(LL l = 2, r; l <= x; l = r + 1)
	{
		r = x / (x / l);
		res -= (r - l + 1) * calc_mu(x / l);
	}
	ans_mu[x] = res;
	return res;
}

LL calc_phi(LL x)//求欧拉函数前缀和
{
	if(x < N) return phi[x];//预处理
	if(ans_phi[x]) return ans_phi[x];//记忆化
	
	LL res = 1ll * x * (x + 1) / 2;
	for(LL l = 2, r; l <= x; l = r + 1)
	{
		r = x / (x / l);
		res -= (r - l + 1) * calc_phi(x / l);
	}
	ans_phi[x] = res;
	return res;
}

参考

数论学习笔记2之杜教筛初探（含例题练习）
浅谈一类积性函数的前缀和
数学/数论专题-学习笔记：狄利克雷卷积
数学/数论专题-学习笔记：莫比乌斯反演
数学/数论专题：莫比乌斯函数与欧拉函数