hdu 2814 快速求欧拉函数

 /**
 大意： 求[a,b] 之间 phi(a) + phi(a+1)...+ phi(b)；
 思路： 快速求欧拉函数
 **/

 #include <iostream>
 #include <cstring>
 using namespace std;
 #define Max 3000000

 long long phi[Max+];
 int prime[Max/];
 bool flag[Max+];

 void init()
 {
     int i,j,num=;
     memset(flag,,sizeof(flag));
     phi[]=;
     for(i=;i<=Max;i++)//欧拉筛选
     {
         if(flag[i])
         {
             prime[num++]=i;
             phi[i]=i-;
         }
         for(j=;j<num && prime[j]*i<=Max;j++)
         {
             flag[i*prime[j]]=false;
             if(i%prime[j]==)
             {
                 phi[i*prime[j]]=phi[i]*prime[j];
                 break;
             }
             else phi[i*prime[j]]=phi[i]*(prime[j]-);
         }
     }
     for(i=;i<=Max;i++)
         phi[i] = phi[i-]+phi[i];
 }

 int main(){
     init();
     long long a,b;
     while(cin>>a>>b){
         long long res = phi[b]-phi[a-];
         cout<<res<<endl;
     }
 }