lightoj 1007 - Mathematically Hard 欧拉函数

求a,b之间的欧拉函数的平方和。

直接筛选法求出欧拉函数，然后来个前缀和。

注意会爆ll，开成ull

#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;
}