【欧拉函数+线性筛】筛法求欧拉函数

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

输入格式

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

输出格式

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

数据范围

1≤n<=1e6

输入样例：

6

输出样例：

12

 

思路：

首先要知道欧拉函数的定义：1 ~ N 中与 N 互质的数的个数被称为欧拉函数，记为ϕ(N)

然后就是一个计算欧拉函数的公式：若将N分解质因数，得到N=(p1^a1)*(p2^a2)*...*(pk^ak),则ϕ(N)=N*(1-1/p1)*(1-1/p2)*...*(1-1/pk)

然后就可以在处理线性筛的过程中顺便求出欧拉函数之和啦~

关于线性筛的原理昨天刚写了一篇可供参考：https://blog.csdn.net/qq_44661312/article/details/103301213，这里的代码就不写有关线性筛原理的注释了，只写有关求欧拉函数的注释。

 

 

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 1e6+10;
int prime[maxn],cnt=0,ph[maxn];
bool st[maxn];

ll getprime(int n)
{
    //根据欧拉函数的定义可得1的欧拉函数为1
    ph[1]=1;
    for(int i = 2; i <= n; i ++)
    {
        if(!st[i])
        {
            prime[++cnt]=i;
            //若i为质数，则所有小于i的正整数都与i互质
            ph[i]=i-1;
        }
        for(int j = 1; i*prime[j]<=n;j++)
        {
            st[prime[j]*i]=true;
            if(i%prime[j]==0)
            {
                //prime[j]为i的质因子，所以i的欧拉函数包括了prime[j]的欧拉函数
                //所以公式中后面所乘的大串因子是相同的，只有前面的整数需要多乘prime[j]
                ph[i*prime[j]]=ph[i]*prime[j];
                break;
            }
            //prime[j]不是i的质因子，因此除了前面的整数需要多乘prime[j]
            //后面的大串因子还需乘上代表prime[j]的因子
            //即ph[i*prime[j]]=ph[i]*prime[j]*(1-1/prime[j])=ph[i]*(prime[j]-1)
            ph[i*prime[j]]=ph[i]*(prime[j]-1);
        }
    }
    //最后统计所有欧拉函数之和
    ll ans=0;
    for(int i = 1; i <= n; i++)ans+=ph[i];
    return ans;
}
int main()
{
    int n;
    cin>>n;
    printf("%lld\n",getprime(n));
}