BZOJ 2818 Gcd
Posted Wisdom+.+
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了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的主要内容,如果未能解决你的问题,请参考以下文章