4804: 欧拉心算莫比乌斯反演线性筛

最新推荐文章于 2020-02-22 00:48:50 发布

200815147

最新推荐文章于 2020-02-22 00:48:50 发布

阅读量203

点赞数

分类专栏：莫比乌斯反演

本文链接：https://blog.csdn.net/baidu_36797646/article/details/81214193

版权

莫比乌斯反演专栏收录该内容

5 篇文章 0 订阅

订阅专栏

题解：

最后推出来是这个式子：（过程很常规，我就懒得写了）

\sum T = 1 n ⌊ n T ⌋ 2 \sum d | T ϕ (d) μ (T d)

$\sum_{T=1}^n{\lfloor {n\over T} \rfloor}^2\sum_{d|T}\phi(d)\mu({T\over d})$ 后面的东西为两个积性函数的狄利克雷卷积，是可以线性筛的，筛的方法就是对于每个数，假如它的质因数分解为

p1k1p2k2……pnkn p 1 k 1 p 2 k 2 … … p n k n ${p_1}^{k_1}{p_2}^{k_2}……{p_n}^{k_n}$ (

p p $p$ 是递增的)，记录它的

{p_{1}}^{k_{1}}

${p_1}^{k_1}$ ，对于

f[pk] f [ p k ] $f[p^k]$ 、

f[p] f [ p ] $f[p]$ 、

f[1] f [ 1 ] $f[1]$ 这些特殊的特判，然后其他的通过这个来分成两个互质的数求。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=10000010;
const int inf=10000000;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int phi[Maxn],prime[Maxn],len=0,f[Maxn],g[Maxn];LL sum[Maxn];
bool mark[Maxn];
void pre()
{
    phi[1]=f[1]=1;
    for(int i=2;i<=inf;i++)
    {
        if(!mark[i])phi[i]=i-1,prime[++len]=i,f[i]=i-2,g[i]=i;
        for(int j=1;j<=len&&prime[j]*i<=inf;j++)
        {
            mark[prime[j]*i]=true;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i]=phi[i]*prime[j];
                g[prime[j]*i]=g[i]*prime[j];
                if(i==g[i])f[prime[j]*i]=phi[prime[j]*i]-phi[i];
                else f[prime[j]*i]=f[i/g[i]]*f[g[prime[j]*i]];
                break;
            }
            phi[prime[j]*i]=phi[i]*(prime[j]-1);
            f[prime[j]*i]=f[prime[j]]*f[i];
            g[prime[j]*i]=prime[j];
        }
    }
    sum[0]=0;
    for(int i=1;i<=inf;i++)sum[i]=(LL)sum[i-1]+f[i];
}
int main()
{
    pre();
    int T=read();
    while(T--)
    {
        int n=read(),pos;LL ans=0;
        for(int i=1;i<=n;i=pos+1)
        {
            pos=n/(n/i);
            ans+=(LL)(n/i)*(n/i)*(sum[pos]-sum[i-1]);
        }
        printf("%lld\n",ans);
    }
}