LightOj1007 - Mathematically Hard（欧拉函数）

#include <stdio.h>
#include <string.h>
#include <algorithm>
typedef unsigned long long LL;
#define N 5000001
using namespace std;

int f
= {0};
LL ans
= {0};

void Init()
{
for(int i=2; i<N; i++)
{
if(f[i]) continue;
for(int j=i; j<N; j+=i)
{
if(!f[j]) f[j] = j;
f[j] = f[j]/i*(i-1);
}
}
ans[0] = 0;
for(int i=1; i<N; i++)
ans[i] = ans[i-1] + (LL)f[i]*f[i];
}

int main()
{
Init();

int T, a, b, t = 1;
scanf("%d", &T);
while(T--)
{
scanf("%d %d", &a, &b);
printf("Case %d: %llu\n", t++, ans[b]-ans[a-1]);
}
return 0;
}