欧拉函数之HDU2824 The Euler function

#include <stdio.h>
#include <ctype.h>
#include <string.h>
#include <stdlib.h>
#include <limits.h>
#include <math.h>
#include <algorithm>
#include <stack>
#include <queue>
#include <vector>
#include <map>
#include <set>
#include <string>
#include <sstream>
using namespace std;

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

typedef long long LL;
const double pi=4.0*atan(1.0);
const int MAXN=3000005;
int euler_phi(int n)
{
int ans=n;
int k=n;
for(int i=2;i*i<=k;i++)
{
if(n%i==0)
{
ans=ans/i*(i-1);
while(n%i==0)
n/=i;
}
}
if(n>1)
ans=ans/n*(n-1);
return ans;
}

LL phi[MAXN];
void phi_table(int n)
{
for(int i=2;i<=n;i++)
phi[i]=0;
phi[1]=1;
for(int i=2;i<=n;i++)
if(!phi[i])
for(int j=i;j<=n;j+=i)
{
if(!phi[j])
phi[j]=j;
phi[j]=phi[j]/i*(i-1);
}

}
int main()
{
int T;
int n,m;
phi_table(3000000);
phi[0]=0;
for(int i=1;i<=3000000;i++)
phi[i]=phi[i]+phi[i-1];
while(scanf("%d%d",&n,&m)!=EOF)
{
printf("%I64d\n",phi[m]-phi[n-1]);//也可以直接求n到m的欧拉函数之和
}
return 0;
}