hdu5605(莫比乌斯反演）

莫比乌斯反演公式
若对于数论函数f(n)f(n)和F(n)F(n)，有以下关系式：
F(n)=∑d|n f(d)
F(n)=∑d|n f(d)

则其莫比乌斯反演公式定义为：
f(n)=∑d|nμ(d)F(n/d)
f(n)=∑d|nμ(d)F(n/d)

它还有第二种形式：
F(n)=∑n∣d f(d)⇒f(n)=∑n∣dμ(d/n)F(d)

其中μ（x）函数是莫比乌斯函数
x有平方数因子时，μ=0，x=1时,μ=1,1,(−1)^k若n无平方数因数，且n=p1p2p3…pn。

对莫比乌斯反函数的进一步推导（看了别人题解得来），由于不会不会写数学符号，这里以图片方式

那么问题就变成了，求一个多项式的第i项乘以莫比乌斯函数的前n/i项和，然后累加，这里就需要用到莫比乌斯的一个重要结论结论是：
前n项和M(n)=∑n i=1 u(i)
结论是M(n)=∑i=1 nu(i)
M(n)=1−∑n i=2 M(n/i)M(n)=1−∑i=2 n M(n/i)
用递归即可
这里先用线性筛筛选出前1e7的M。

这里还有一个问题，求F（x）前n项和%(1e9+1)=x(x+1)(2x+1)/6-3x(x+1)/2+2x=(x^3-3xx+2x)/3%(1e9+7)= (xxx-3x+2x)/3(1e9+7+1)%1e9+7=((xxx-3x+2x)*((1e9+7+1)/3))%(1e9+7)

代码实现

#include<bits/stdc++.h>
const int N=5e+5;
typedef long long ll; 
const ll MOD=1e9+7;
using namespace std;
int mu[N],ji[N]; 
vector<int>p;
map<int,int>mp;
ll stil(ll y)//求F（n）d前y项和
{

//	return (((((1ll*y*(y+1))%MOD*(2*y+1))/6%MOD)-(1ll*3*y*(y+1)/2)%MOD+1ll*2*y)+MOD)%MOD;
	return (((y*y%MOD*y%MOD-3*y%MOD*y%MOD+2*y%MOD)%MOD)+MOD)%MOD*333333336%MOD;
}
void init()
{
	mu[1]=1;
	memset(ji,0,sizeof(ji));
	for(int i=2;i<N;i++)//预处理出前1e7项μ
	{
		if(!ji[i])
		{
			p.push_back(i);
			mu[i]=-1;
		}
		int m=p.size();
		for(int j=0;j<m&&i*p[j]<N;j++)
		{
			int g=i*p[j];
			ji[g]=1;					
			if(i%p[j])
				mu[g]=-mu[i];
			else 
			{
				mu[g]=0;
			break;}
		}
	}
	for(int i=2;i<N;i++)//这里的mu变成莫比乌斯前缀和了
	{
	mu[i]+=mu[i-1];
}
 } 
 int MU(int x)
 {
 	
 	if(x<N)
 	{
	 //cout<<mu[x]<<endl;
 	return mu[x];
 	}
	 if(mp[x])
 	return mp[x];
 	int ans=1,j;
 	for(int i=2;i<=x;i=j+1)
 	{
 		 j=x/(x/i);
 		ans=((ans-1ll*MU(x/i)*(j-i+1))%MOD+MOD)%MOD; //莫比乌斯前x项和递归公式
 		
 }	
 	mp[x]=ans;
 	return ans;
  }
  int main()
  {
  	int t,n;
  	cin>>t;
  	init();
  	while(t--)
  	{
  		scanf("%d",&n);
  		int j;
  		int ans=0;
  		for(int i=1;i<=n;i=j+1)
		 {
			j=n/(n/i);
			ans=(ans+((stil(j)-stil(i-1))+MOD)%MOD*MU(n/i)%MOD+MOD)%MOD; 
		//	cout<<ans<<endl;	
		 }
		 cout<<ans<<endl;
	  }
  }