BZOJ 2818 Gcd

Posted Wisdom+.+

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ 2818 Gcd相关的知识,希望对你有一定的参考价值。

BZOJ 2818 Gcd

思路:枚举每个质数pi,求出∑phi[n/pi],对phi函数前缀和处理一下。

代码:

 

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pb push_back
#define mem(a,b) memset(a,b,sizeof(a))

const int N=1e7+5;
bool not_prime[N]={false};
int phi[N];
int prime[N];
ll sum[N];
int k,n;
void Euler(int n)
{
    phi[1]=1;
    k=0;
    for(int i=2;i<=n;i++)
    {
        if(!not_prime[i])
        {
            phi[i]=i-1;
            prime[k++]=i;
        }
        for(int j=0;i*prime[j]<=n;j++)
        {
            not_prime[i*prime[j]]=true;
            if(i%prime[j]==0)//根据欧拉函数的性质。 
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    int n;
    cin>>n;
    Euler(n);
    sum[1]=0;
    ll ans=0;
    for(int i=2;i<=n;i++)sum[i]=sum[i-1]+phi[i];
    for(int i=0;i<k;i++)ans+=1+sum[n/prime[i]]*2;
    cout<<ans<<endl;
    return 0;
}

 

线性筛速度大约是普通筛法的5倍

 

以上是关于BZOJ 2818 Gcd的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 2818 Gcd(莫比乌斯反演)

bzoj 2818: Gcd

BZOJ2818 Gcd

BZOJ-2818Gcd 线性筛

BZOJ 2818: Gcd

bzoj 2818 Gcd