相关链接:线性筛法素数的方法见我的博客文章埃氏筛法和线性筛法求素数.
先看题目:筛法求欧拉函数
分析:
如果用分解质因数的方法求解欧拉函数之和,需要大约O(n*sqrt(n))的时间复杂度,为了更为高效的求解1~n中每个数的欧拉函数,我们可以利用欧拉函数的性质,这样O(n)的时间复杂度便能解决问题,看算法的魅力!
具体性质如下图(还是手写数学式子方便,凑合看手稿吧):
基于以上三个性质,我们可以在欧拉筛素数的过程中完成对每一个数的欧拉函数的计算。具体如代码:
//在线性欧拉筛素数的基础上计算欧拉函数,巧妙,效率高。
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1000000+7;
int primes[N],phi[N];
bool st[N];
LL getEulers(int n){
phi[1] = 1;
int cnt = 0;
for(int i = 2;i<= n; i++){
if(!st[i]){
primes[++cnt] = i;
phi[i] = i-1; //性质1
}
for(int j = 1;j<= cnt && primes[j]*i <=n; j++){
st[primes[j] * i] = true;
if(i % primes[j] == 0){
phi[primes[j] * i] = primes[j] * phi[i]; //性质2
break;
}
else{
phi[primes[j] * i] = (primes[j] - 1) * phi[i]; //性质3
}
}
}
LL ans = 0;
for(int i = 1; i<= n; i++){
ans = ans + phi[i];
}
return ans;
}
int main(){
int n;
scanf("%d",&n);
LL ans = getEulers(n);
printf("%lld",ans);
return 0;
}
好了,全文结束,如有不当之处还请看官多多批评指正。