洛谷P2568 GCD【欧拉函数】

时空限制 1000ms / 256MB

题目描述

给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对.

说明

1<=N<=10^7

---

题目分析

第一眼想到的是莫比乌斯反演，然而愉快地T飞了
仔细一想单词询问莫比乌斯反演的预处理确实挺费时的

考虑转化题意为对每个不大于n的素数p
求$1\le x,y\le \lfloor \frac{n}{p}\rfloor$内有多少个$(x,y)$满足$gcd(x,y)==1$
即$ans={\sum }_{p\in prime}{\sum }_{i=1}^{\lfloor \frac{n}{p}\rfloor }{\sum }_{j=1}^{\lfloor \frac{n}{p}\rfloor }[gcd(i,j)=1]={\sum }_{p\in prime}(sum[\lfloor \frac{n}{p}\rfloor ]\ast 2-1)$

#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
typedef long long lt;
 
int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=1e7+10;
int vis[maxn],prim[maxn],cnt;
lt phi[maxn],sum[maxn],ans;

void Phi(int n)
{
    phi[1]=1;
    for(int i=2;i<=n;++i)
    {
        if(!vis[i]) prim[++cnt]=i,phi[i]=i-1;
        for(int j=1;j<=cnt;++j)
        {
            if(i*prim[j]>n) break;
            vis[i*prim[j]]=1;
            if(i%prim[j]==0){ phi[i*prim[j]]=phi[i]*prim[j]; break;}
            else phi[i*prim[j]]=phi[i]*phi[prim[j]];
        }
    }
    for(int i=1;i<=n;++i)
    sum[i]=sum[i-1]+phi[i];
}

int main()
{
    int n=read();
    if(n==1){cout<<0;return 0;}
    
    Phi(n);
    for(int i=1;i<=cnt;++i)
    ans+=sum[n/prim[i]]*2-1;
    
    printf("%lld",ans);
    return 0;
}