AcWing 874.筛法求欧拉函数(模板题)

Issue：

Thinking：

1. 在线性筛质数下同时把欧拉函数给筛出来；
2. 如果一个数 $i$ 是质数，那么欧拉函数 $\varphi (i)=i-1$；
3. if (i % primes[j] == 0)：$\varphi (i)=i(1-\frac{1}{P_1})\cdots (1-\frac{1}{P_j})\cdots (1-\frac{1}{P_k})$，//$P_j$是i的一个质因子，直接用欧拉函数的公式，同时$P_j$也是$P_j\times i$的质因子；
   $\varphi (P_{j}i)=P_{j}i(1-\frac{1}{P_1})\cdots (1-\frac{1}{P_j})\cdots (1-\frac{1}{P_k})=P_j\varphi (i)$；
4. if (i % primes[j] != 0：$\varphi (i)=i(1-\frac{1}{P_1})\cdots (1-\frac{1}{P_k})$，//$P_j$不是i的质因子，但是$P_j\times i$的质因子；
   $\varphi (P_{j}i)=P_{j}i(1-\frac{1}{P_1})\cdots (1-\frac{1}{P_k})(1-\frac{1}{P_j})=P_j\varphi (i)(1-\frac{1}{P_j})=(P_j-1)\varphi (i)$；

Code：

#include <iostream>

using namespace std;

const int N = 1e6 + 10;

int primes[N], cnt;
int phi[N];
bool st[N];

long long euler(int n)
{
    phi[1] = 1;  //1~1中和1互质的数只有1
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            primes[cnt ++ ] = i;
            phi[i] = i - 1;  //质数i的欧拉函数为i-1
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
            {
                phi[primes[j] * i] = primes[j] * phi[i];
                break;
            }
            phi[primes[j] * i] = (primes[j] - 1) * phi[i];
        }
    }
    
    long long res = 0;
    for (int i = 1; i <= n; i ++ ) 
        res += phi[i];

    return res;
}

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

Tips：

1. 欧拉筛模板：

int primes[N], cnt;
int phi[N];
bool st[N];

void euler(int n)
{
    phi[1] = 1;
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            primes[cnt ++ ] = i;
            phi[i] = i - 1;
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
            {
                phi[primes[j] * i] = primes[j] * phi[i];
                break;
            }
            phi[primes[j] * i] = (primes[j] - 1) * phi[i];
        }
    }
}