lightoj 1007 - Mathematically Hard 欧拉函数

#include<bits/stdc++.h>
using namespace std;
#define N 5010000
#define ll long long
#define ull unsigned long long
int phi[N+20];
ll ans[N+20];
void euler()
{
for(int i=1;i<=N;i++)
phi[i]=i;
phi[1]=1;
for(int i=2;i<=N;i++)
{
if(phi[i]==i)
{
phi[i]=i-1;
for(int j=i+i;j<=N;j+=i)
{
phi[j]=phi[j]/i*(i-1);
}
}
}
ans[0]=0;
for(int i=1;i<=N;i++)
{
ans[i]=ans[i-1]+(ll)phi[i]*phi[i];
}
}
int main()
{
euler();
int t;
scanf("%d",&t);
for(int cas=1;cas<=t;cas++)
{
int a,b;
scanf("%d %d",&a,&b);
printf("Case %d: %llu\n",cas,ans[b]-ans[a-1]);
}
return 0;
}