线性筛——莫比乌斯函数

如果不会线性筛素数的话,建议先看这篇博客了解一下线性筛素数。
莫比乌斯函数函数(积性函数都可以线性筛)主要是在线性筛素数的基础上得到的

我们知道:
n = ∏ i = 1 n p i t i n=\prod_{i=1}^n p_i^{t_i} n=i=1npiti

μ ( n ) = { ( − 1 ) k k = ∑ i = 1 n t i a n d max ⁡ ( t 1 , t 2 , ⋯   , t n ) ≤ 1 0 o t h e r s μ(n)=\begin{cases} (-1)^k &k=\sum_{i=1}^n t_i \quad and\quad \max(t_1,t_2,\cdots,t_n)\leq1\\ 0 & others \end{cases} μ(n)={(1)k0k=i=1ntiandmax(t1,t2,,tn)1others

所以:
1、当 n n n 是质数时, μ ( n ) = − 1 μ(n)=-1 μ(n)=1
对于 2和3 设 d = n p d=\frac{n}{p} d=pn 其中 p p p n n n 的最小质因子
2、当 p p p d d d 的某个质因子时, 则 μ ( n ) = 0 μ(n)=0 μ(n)=0
3、当 p p p d d d 互质时, μ ( n ) = − μ ( d ) μ(n)=-μ(d) μ(n)=μ(d)

good luck and have fun!!!
附上代码:

int mu[MAXN];
void Mobius(int n)
{
	memset(vis,0,sizeof(vis));
	mu[1]=1;
	prime[0]=0;
	for(int i=2;i<=n;i++)
	{
		if(!vis[i])
		{
			prime[++prime[0]]=i;
			mu[i]=-1;
		}
		for(int j=1;j<=prime[0]&&i<=n/prime[j];j++)
		{
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)
			{
				mu[i*prime[j]]=0;
				break;
			}
			mu[i*prime[j]]=-mu[i];
		}
	}
}

你可能感兴趣的:(数论)