(AcWing)筛法求欧拉函数

给定一个正整数 n，求 1∼n 中每个数的欧拉函数之和。

输入格式

共一行，包含一个整数 n。

输出格式

共一行，包含一个整数，表示 1∼n 中每个数的欧拉函数之和。

数据范围

1≤n≤10^6

输入样例：

6

输出样例：

12

#include<iostream>
using namespace std;

typedef long long LL; // 定义长整型LL

const int N = 1000010; // 定义常量N，表示数组的大小上限

int prime[N], cnt = 0; // 存储所有素数的数组prime和计数器cnt
int pri[N]; // 存储n以内的欧拉函数值的数组pri
bool str[N]; // 记录sy[x]中的x是否被筛掉的布尔数组str

LL find(int n)
{
    pri[1] = 1; // 欧拉函数定义：phi(1) = 1
    for(int i=2;i<=n;i++)
    {
        if(!str[i]) // 如果i没有被筛掉，则i是素数
        {
            prime[cnt++] = i; // 将i存入素数数组prime，并增加计数器cnt
            pri[i] = i-1; // 欧拉函数定义：如果i是素数，则phi(i) = i - 1
            str[i] = true;
        }
        for(int j=0;prime[j]<=n/i;j++) // 遍历素数数组prime
        {
            str[prime[j]*i] = true; // 筛掉合数
            if(i%prime[j]==0) // 如果i能整除prime[j]
            {
                pri[prime[j]*i] = pri[i]*prime[j]; // 根据欧拉函数性质计算phi(prime[j]*i)
                break;
            }
            pri[prime[j]*i] = pri[i]*(prime[j]-1); // 根据欧拉函数性质计算phi(prime[j]*i)
        }
    }
    LL res = 0; // 定义结果变量res，并初始化为0
    for(int i=1;i<=n;i++) res+=pri[i]; // 将数组pri中的所有欧拉函数值累加到res中

    return res; // 返回结果res
}

int main()
{
    int n;
    cin>>n; // 输入一个整数n
    cout<<find(n)<<endl; // 调用find函数计算n以内所有数的欧拉函数值之和，并输出结果
}