莫比乌斯反演公式
若对于数论函数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;
}
}