线性筛求莫比乌斯函数前缀和

大家知道有一类问题，让你把1~N之间每一个数的莫比乌斯函数都输出来，或者把它们全部加起来再输出。这种问题应该属于求前缀和一类的啦 

用正规的方式表达：

就论求1~N间每个数的莫比乌斯函数吧。我们一般使用Eratosthenes筛来解决。（简称埃氏）大家知道埃氏求素数的过程，其实就是标记N以内所有质数小于N的倍数，这样每个数有哪些质因子都可以被标记出来，同理若未被标记就说明是质数。利用每个质因子都可以去标记这一性质，我们就可以判断这个数的质因子数量的奇偶了。

int n, v[100000], miu[100000]；
for(int i=1;i<=n;i++)miu[i]=1,v[i]=0;
	for(int i=2;i<=n;i++){
		if(v[i])continue;
		miu[i]=-1;
		for(int j=2*i;j<=n;j+=i){
			v[j]=1;
			if((j/i)%i==0)miu[j]=0;
			else miu[j]*=-1;
		}
	}

那么这个问题，能不能用线性筛来做呢？

当然也是可以的。

网上有很多线性筛求莫比乌斯函数的模板，但是由于蒟蒻比较懒，没去看dalao们的代码。于是就自己了想一个出来（保留所有版权）

大家知道线性筛有一个很重要的特征，就是每个数只会被其最小的质因子筛一次，比埃氏更省时。那么线性筛求莫比乌斯函数也是在线性筛求素数的基础上改编的。

线性筛求1~N间素数模板如下：

int v[MAXN_N],prime[MAXN_N];
void primes(int n){
	memset(v,0,sizeof(v));
	int m=0;
	for(int i=2;i<=n;i++){
		if(v[i]==0){
			v[i]=i;
			prime[++m]=i;
		}
		for(int j=1;j<=m;j++){
			if(prime[j]>v[i]||prime[j]>n/i)break;
			v[i*prime[j]]=prime[j];
		}
	} 
	for(int i=1;i<=m;i++)cout<<prime[i]<<endl;
} 

这个大家都十分熟悉。把所有的素数装到一个数组里，然后用每个数i去依次乘以这些素数prime，得到的i*prime就将其最小的质因子标记为所乘的prime。但前提是如果prime大于i的最小质因子就不能再继续了，因为此时的i*prime一定会有比prime更小的质因子。

那这个线性筛既然没有用每个质因子都去筛一遍，我们怎么去判断每个数质因子的数量的奇偶呢？

别急。我们很快发现，如果a*b=c，那么c分解质因子后，一定会等于a的所有质因子乘上b的所有质因子。eg: 4*6=24，4=2*2，6=2*3，24=2*2*2*3

那么既然质因子是可以累积的，质因子的数量奇偶是否可以累积呢？

奇数+奇数=偶数

偶数+偶数=偶数

奇数+偶数=奇数

与负负得正一样，相当于把异或关系反过来。

所以现在就一目了然了。若有一个因子的μ等于0，则其所有倍数的μ都等于0。

若μ(a)=-1，μ(b)=1，则μ(a*b)=-1*1=-1。

然鹅还有一点需要注意：a和b需要互质。

比如6和21，μ(6)!=0，μ(21)!=0，但6*21中包含两个相同的质因子3，因此μ(6*21)=0。

那么由于我们线性筛本身就是一个质数prime乘以一个正整数i，且prime不可能大于i，所以只需要判断i是否能整除prime即可。

若能，他们就不互质，则μ(i*prime)直接等于零。

若不能，他们互质，则μ(i*prime)=μ(i)*μ(prime) //由于函数性质，其实就相当于取反

int v[MAXN_N],prime[MAXN_N],mobius[MAXN_N]；
void mobius(int n){
	int m=0;
	memset(v,0,sizeof(v));
	memset(mobius,-1,sizeof(mobius));
	for(int i=2;i<=n;i++){
		if(v[i]==0){
			v[i]=i;
			prime[++m]=i;
		}
		for(int j=1;j<=m;j++){
			if(prime[j]>v[i]||prime[j]>n/i)break;
			v[i*prime[j]]=prime[j];
			if(i%prime[j]==0)mobius[i*prime[j]]=0;
			else mobius[i*prime[j]]=mobius[i]*mobius[prime[j]];
		}
	} 
	mobius[1]=1;
	for(int i=1;i<=n;i++)cout<<mobius[i]<<endl;
} 

最后献上窝丑陋的代码。感兴趣的同学可以去了解一下积性函数和完全积性函数，你会发现更多有意思的东西！