1 /**
2 大意： 求[a,b] 之间 phi(a) + phi(a+1)...+ phi(b)；
3 思路： 快速求欧拉函数
4 **/
5
6 #include <iostream>
7 #include <cstring>
8 using namespace std;
9 #define Max 3000000
10
11 long long phi[Max+5];
12 int prime[Max/10];
13 bool flag[Max+5];
14
15 void init()
16 {
17 int i,j,num=0;
18 memset(flag,1,sizeof(flag));
19 phi[1]=1;
20 for(i=2;i<=Max;i++)//欧拉筛选
21 {
22 if(flag[i])
23 {
24 prime[num++]=i;
25 phi[i]=i-1;
26 }
27 for(j=0;j<num && prime[j]*i<=Max;j++)
28 {
29 flag[i*prime[j]]=false;
30 if(i%prime[j]==0)
31 {
32 phi[i*prime[j]]=phi[i]*prime[j];
33 break;
34 }
35 else phi[i*prime[j]]=phi[i]*(prime[j]-1);
36 }
37 }
38 for(i=2;i<=Max;i++)
39 phi[i] = phi[i-1]+phi[i];
40 }
41
42 int main(){
43 init();
44 long long a,b;
45 while(cin>>a>>b){
46 long long res = phi[b]-phi[a-1];
47 cout<<res<<endl;
48 }
49 }
hdu 2814 快速求欧拉函数
时间： 2024-12-22 07:41:47