数论 - 筛法求欧拉函数

数论 - 筛法求欧拉函数

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

输入格式
共一行，包含一个整数n。

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

数据范围
1≤n≤10⁶

输入样例：
6
输出样例：
12

---

分析：

$先预处理欧拉函数再求和。$

$欧拉函数：$
$$\varphi (n)=n(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_n})，其中p_1,p_2,...,p_n是n的质因子。$$

$表示的是1$~$n-1中与n互质的数的个数。$

$考虑线性筛的过程：$

$①、若i是质数，则\varphi (i)=i-1。$

$②、若i是合数，$

$$\begin{gathered} 若i\%p_j=0，则i的质因子与i\times p_j的质因子完全相同，由\varphi (i)=i(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_i})， \\ 得\varphi (p_j\times i)=p_j\times i(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_i})=p_j\times \varphi (i)。 \end{gathered}$$

$$\begin{gathered} 若i\%p_j\ne 0，则i的质因子比i\times p_j少p_j，由\varphi (i)=i(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_i})， \\ 得\varphi (p_j\times i)=p_j\times i\times (1-\frac{1}{p_j})(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_i})=p_j(1-\frac{1}{p_j})\times \varphi (i)。 \end{gathered}$$

代码：

#include<iostream>

using namespace std;

#define ll long long

const int N=1e6+10;

int n,phi[N],prime[N],cnt;
bool st[N];

ll get_euler(int x)
{
    phi[1]=1;
    for(int i=2;i<=x;i++)
    {
        if(!st[i])
        {
            prime[cnt++]=i;
            phi[i]=i-1;
        }
        for(int j=0;prime[j]<=x/i;j++)
        {
            st[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=prime[j]*phi[i];
                break;
            }
            else phi[i*prime[j]]=(prime[j]-1)*phi[i];
        }
    }
    
    ll res=0;
    for(int i=1;i<=x;i++) res+=phi[i];
    
    return res;
}

int main()
{
    cin>>n;
    cout<<get_euler(n)<<endl;
    
    return 0;
}