hdu 2814 快速求欧拉函数

 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 }