[积性函数、杜教筛、莫比乌斯函数入门]学习总结

什么是积性函数？

1.数论函数

在数论上，算术函数（或称数论函数）指定义域为正整数、陪域为复数的函数，每个算术函数都可视为复数的序列。
最重要的算术函数是积性及加性函数。算术函数的最重要操作为狄利克雷卷积，对于算术函数集，以它为乘法，一般函数加法为加法，可以得到一个阿贝尔环。

2.积性函数

积性函数指对于所有互质的整数a和b有性质f(ab)=f(a)f(b)的数论函数。

• 定理1.1：积性函数的和函数也是积性函数，如果f是积性函数，那么f的和函数$F(n)={\sum }_{d|n}f(d)$也是积性函数。
• 常见的积性函数：
  $I(x)恒等函数，I(n)=1$
  $id(n):单位函数，id(n)=n$
  $I_k(n):幂函数，I_k(n)=n^k$
  $ϵ(n):原函数，ϵ(n)=\{\begin{array}{l}1(n=1)\\ 0(n>1)\end{array}$
  $\sigma (n):因子和函数，\sigma （n）={\sum }_{d|n}d$
  $d(n)：约数个数，d(n)={\sum }_{d|n}1$
  $\mu (n):莫比乌斯函数$
  $gcd(n,k)：最大公因子，当k固定的情况$

狄利克雷卷积

设 $f,g$是算数函数,记$f,g$的狄利克雷卷积是$f\ast g$, 定义为 $$(f\ast g)(n)=\sum _{d|n}f(d)g(\frac{n}{d})$$

狄利克雷卷积满足:交换律,结合律,分配律

莫比乌斯函数

1.莫比乌斯函数($\mu (x)$)的定义:
    $$\mu (x)=\{\begin{array}{l}0(\exists {\alpha }_i⩾2)\\ (-1{)}^k(\forall {\alpha }_i=1)\end{array}$$
2.莫比乌斯函数($\mu (x)$)的性质:
    $$若有：ϵ(x)=\sum _{d|n}\mu (d)，那么:ϵ(x)=\{\begin{array}{l}1(n=1)\\ 0(n>1)\end{array}$$

莫比乌斯反演

提到莫比乌斯函数，那莫比乌斯反演必不可少。
   公式一:
$$若F(n)=\sum _{d|n}f(d),则f(n)=\sum _{d|n}\mu (d)\ast F(\frac{n}{d})$$
   公式二:
$$若F(n)=\sum _{n|d}f(d),则f(n)=\sum _{n|d}\mu (\frac{d}{n})\ast F(d)$$
关于莫比乌斯反演的证明，会在不久后更新。

杜教筛公式

杜教筛是用于解决数论函数f(x)的前缀和的问题.
如:$s(n)={\sum }_{i=1}^{n}f(i)$

构造$$h=g\ast f$$
得
$$\sum _{i=1}^{n}h(i)=\sum _{i=1}^n\sum _{d|i}g(d)f(\frac{i}{d})$$
$$=\sum _{d=1}^{n}g(d)\sum _{d|i}^{n}f(\frac{i}{d})$$
$$=\sum _{d=1}^{n}g(d)\sum _{i=1}^{[\frac{n}{d}]}f(i)$$
$$=\sum _{d=1}^{n}g(d)s([\frac{n}{d}])$$
得到
$$\sum _{i=1}^{n}h(i)=\sum _{d=1}^{n}g(d)s([\frac{n}{d}]）$$
提出第一项
$$\sum _{i=1}^{n}h(i)=g(1)s(n)+\sum _{d=2}^{n}g(d)s([\frac{n}{d}]）$$
得到杜教筛公式：
$$g(1)s(n)=\sum _{i=1}^{n}h(i)-\sum _{d=2}^{n}g(d)s([\frac{n}{d}]）$$

理解：翻遍CSDN,鲜有人详细证明步骤二,但第二步的变换却是杜教筛的精华,这边记录下自己的理解:1式枚举了每一个i,对于每一个i,又枚举了每一个d,我们目光关注在g(d)这个东西上,发现:d是可以取遍所有n的,对于每一个g(d),都有$f(\frac{i}{d})$其中,$i\times d<=n$,并且全部能取到,所以第二步"杜教筛变换"是对的.


栗子：求欧拉函数的前n项。
可以直接套用杜教筛公式：
$$由于id(n)=\sum _{d|n}\varphi (x)$$
$$构造h(x)=id(x),g(x)=\varphi (x),f(x)=I(x)$$
代入杜教筛公式得：
$$s(n)=\sum _{i=1}^{n}i-\sum _{i=2}^{n}s([\frac{n}{i}])$$
代码

#include <unordered_map>
#include <iostream>

using namespace std;

typedef long long LL;

const int N = 5e6 + 10;
const int mod = 1e9 + 7;
const int inv2 = 500000004;

unordered_map<LL , int> mp;
int primes[N] , cnt , phi[N] , sum[N];
bool st[N];

void init(int n)
{
	phi[1] = 1;
	for(int i = 2;i <= n;i ++)
	{
		if(!st[i])
		{
			primes[cnt ++ ] = i;
			phi[i] = i - 1;	
		}
		
		for(int j = 0;primes[j] <= n / i;j ++)
		{
			st[primes[j] * i] = true;
			if(i % primes[j] == 0)
			{
				phi[i * primes[j]] = phi[i] * primes[j];
				break;
			}
			phi[i * primes[j]] = phi[i] * (primes[j] - 1);
		}
	}
	
	for(int i = 1;i < n;i ++) sum[i] = (sum[i - 1] + phi[i]) % mod;
}

LL get(LL k ,LL x)
{
	return k / (k / x);
}

LL sol(LL n)
{
	if(n < N - 1) return sum[n];
	if(mp[n]) return mp[n];
	
	LL ans = ((n % mod) * ((n + 1) % mod)) % mod * inv2 % mod;
//	LL ans = n % mod * (n + 1) % mod * inv2 % mod;

	for(LL l = 2 , r;l <= n;l = r + 1)
	{
		r = get(n , l);
//		ans = ans - (r - l + 1) * sol(n / i);
		ans = (ans - (r - l + 1) % mod * sol(n / l) % mod + mod) % mod;
	}
	
	return mp[n] = ans;
}

int main()
{
	init(N - 1);
	LL n;cin >> n;cout << sol(n) << endl;
	return 0;
}

技巧

• $\mu \ast I=ϵ$
• $\varphi \ast I=id$
• $\mu \ast id=\varphi$

参考文献：

1.杜教筛 以及积性函数的前世今生 --算法竞赛专题解析（4)
2.百度百科
------数论函数
------积性函数