HDU 2824 The Euler function

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int N=3000001;

int phi[N];

void init(){
    for(int i=1;i<N;i++)
        phi[i]=i;

    for(int i=2;i<N;i++){
        if(i==phi[i])     //此时i是素数
            for(int j=i;j<N;j+=i)
                phi[j]=(phi[j]/i)*(i-1);  //每一个j都有素因子i
    }
}

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