LightOJ1007 Mathematically Hard（欧拉函数打表+前缀和（附：欧拉函数模板（2套）））

题目链接：https://cn.vjudge.net/contest/319236#problem/A
Sample Input

3
6 6
8 8
2 20

Sample Output

Case 1: 4
Case 2: 16
Case 3: 1237

给定一个区间，求区间内的每个数的欧拉函数的平方进行累加。

知识点：

求单个数的欧拉函数

欧拉函数的定义：phi(n)=n*(1-p1)*(1-p2)*....*(1-pn)//p1,p2,,,pn为n的质因数
void eular(int x)
{
    int res=x;//res的值即为x的欧拉函数
    for(int i=2; i*i<=x; i++) //下面要找x%i==0的i，在sqrt(x)之后就没有意义了
    {
        if(x%i==0)//找到一个质因数
        {
            res-=res/i;//相当于n*(1-p1)的展开
            while(x%i==0)
                x/=i;//保证这个因子完全除去
        }
    }
    if(x>1)//还剩一个一次幂的因子
        res-=res/x;
    return res;
}

欧拉函数打表（求某个区间各个数的欧拉函数）

void solve()
{
    memset(phi,0,sizeof(phi));
    phi[1]=1;//1的欧拉函数为本身
    for(int i=2; i<N; i++) //区间
    {
        if(phi[i]==0)//找到每一个素数
        {
            for(int j=i; j<N; j=+i)//
            {
                if(phi[j]==0)
                    phi[j]=j;//相当于公式中的第一个数字n
                phi[j]=phi[j]/i*(i-1);//如果j只能分解为i，此时的book[j]即为j的欧拉函数，但往往不是
            }                           //例如i=2是，book[6]=3(实则6的欧拉函数为2)，
                                        //当执行i=3时，再次执行，book[6]=3/3*(3-1)=2
            //book[j]先/i再*(i-1)的目的，避免book[j]*(i-1)超数据范围
        }
    }
}

此题就是欧拉函数打表+前缀和：
代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define LL unsigned long long
#define N 5000001
LL book[N];
void init()
{
    memset(book,0,sizeof(book));
    book[1]=1;
    for(int i=2; i<N; i++)
    {
        if(book[i]==0)
        {
            for(int j=i; j<N; j+=i)
            {
                if(book[j]==0)
                    book[j]=j;
                book[j]=book[j]/i*(i-1);
            }
        }
    }
    for(int i=2; i<N; i++)
        book[i]=book[i-1]+book[i]*book[i];
}
int main()
{
    init();
    int T,f=1;
    scanf("%d",&T);
    while(T--)
    {
        int a,b;
        scanf("%d %d",&a,&b);
        printf("Case %d: ",f++);
        printf("%llu\n",book[b]-book[a-1]);
    }
    return 0;
}