51nod 1594 Gcd and Phi 莫比乌斯反演

最新推荐文章于 2020-06-28 15:03:00 发布

_beginend

最新推荐文章于 2020-06-28 15:03:00 发布

阅读量291

点赞数

分类专栏：莫比乌斯反演

本文链接：https://blog.csdn.net/qq_33229466/article/details/79009750

版权

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

24 篇文章 0 订阅

订阅专栏

题意

$F(n)=\sum_{i=1}^n\sum_{j=1}^n\varphi_{(\varphi(i),\varphi(j))}$
其中 $\varphi$ 表示欧拉函数。欧拉函数 $\varphi_n$ 是不超过n的数中与n互质的数的数目。
$\varphi_{(\varphi(i),\varphi(j))}$ 表示i,j欧拉函数值的最大公约数的欧拉函数值.
给出n，求 $F(n)$ 的值。
T<=5,n<=2*10^6

分析

枚举最大公约数然后反演一波，出来的式子就可以直接预处理了。每次 $O(nln(n))$ 求答案即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int N=2000005;

int t[N],phi[N],mu[N],tot,prime[N];
bool not_prime[N];
LL f[N],g[N],ans[N];

void get_prime(int n)
{
    phi[1]=mu[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!not_prime[i]) prime[++tot]=i,phi[i]=i-1,mu[i]=-1;
        for (int j=1;j<=tot&&i*prime[j]<=n;j++)
        {
            not_prime[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
            mu[i*prime[j]]=-mu[i];
        }
    }
}

LL solve(int n)
{
    LL ans=0;
    for (int i=1;i<=n;i++) t[i]=g[i]=f[i]=0;
    for (int i=1;i<=n;i++) t[phi[i]]++;
    for (int i=1;i<=n;i++)
        for (int j=i;j<=n;j+=i)
            g[i]+=t[j];
    for (int i=1;i<=n;i++) g[i]*=g[i];
    for (int i=1;i<=n;i++)
        for (int j=i;j<=n;j+=i)
            f[i]+=(LL)g[j]*mu[j/i];
    for (int i=1;i<=n;i++) ans+=(LL)phi[i]*f[i];
    return ans;
}

int main()
{
    get_prime(2000000);
    int T;scanf("%d",&T);
    while (T--)
    {
        int x;scanf("%d",&x);
        printf("%lld\n",solve(x));
    }
    return 0;
}