欧拉函数的学习

欧拉函数的定义
1 ~ N 中与 N 互质的数的个数被称为欧拉函数，记为ϕ(N)。
若在算数基本定理中，N = p1 ^ a1 * p2 ^ a2 * … * pk ^ ak
那么ϕ(N) = N * （1 - 1 / p1） * (1 - 1 / p2) * … * (1 - 1 / pk).

证明过程：
所用到的原理：容斥原理
1.从1 - N 中去掉 p1, p2, …, pk的所有的倍数
2. 加上所有 pi * pj 的倍数
…

所以有：
ϕ(N) = N - N / P1 - N / P2 - N / P3 - … - N / pk + N / p1 * p2 + N / p1 * p3 + … N / pk * p(k - 1) - N / p1 * p2 * p3 …

所有得到：ϕ(N) = N * (1 - 1 / P1) * (1 - 1/ P2) * … * (1 - 1 / Pk).

模板题：873. 欧拉函数
题目链接：https://www.acwing.com/problem/content/875/
代码:

#include <iostream>
#include <cstdio>

using namespace std;

int main(void) {
    int n, x;
    scanf("%d", &n);
    while(n --) {
        scanf("%d", &x);
        int res = x;
        
        for(int i = 2; i * i <= x; i ++) {
            if(x % i == 0) {
                res = res / i * (i - 1);
                while(x % i == 0) x/= i;
            }
        }
        if(x > 1) res = res / x * (x - 1);
    
        printf("%d\n", res);
        
    }
    
    return 0;
}

线性筛欧拉函数
题目：AcWing 874. 筛法求欧拉函数
题目链接:https://www.acwing.com/problem/content/876/
做法:线筛
代码：

#include <iostream>
#include <cstdio>

using namespace std;

const int maxn = 1e6 + 5;

bool vis[maxn];
int prime[maxn], phi[maxn], cnt;

inline void get_phi(int n) {
    long long res = 0;
    phi[1] = 1;
    for(int i = 2; i <= n; i ++) {
        if(!vis[i]) phi[i] = i - 1, prime[cnt ++] = i;
        for(int j = 0; prime[j] <= n / i; j ++) {
            vis[i * prime[j]] = true;
            if(i % prime[j] == 0) {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            } else phi[i * prime[j]] = phi[i] * (prime[j] - 1);
        }
    }
    
    for(int i = 1; i <= n; i ++)
        res += phi[i];
    
    printf("%lld\n", res);
}

int main(void) {
    int n; 
    scanf("%d", &n);
    get_phi(n);
    return 0;
}

证明：
当 i 是质数的时候，phi[i] = i - 1;
当 i 是和数的时候：
1.当 i % prime[j] == 0 时：phi[i * prime[j]] = phi[i] * prime[j]
2.当 i % prime[j] != 0 时, phi[i * prime[j]] = phi[i] * primep[j] * (1 - 1 / prime[j])
= phi[i] * (prime[j] - 1)