【刷题】数学知识——欧拉函数

欧拉函数$\varphi (n)$：表示1~n中与n互质的数的个数
例如$\varphi (6)=2$，因为1 2 3 4 5 6中只有1和5跟6互质。
注意1和任何数互质，1的欧拉函数是1。

求$\varphi (n)$方法如下：
当$n=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_k^{a_k}$时，其中$p$均为质数
$\varphi (n)=n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ...\ast (1-\frac{1}{p_k})$
例如$\varphi (6)=6\ast (1-\frac{1}{2})\ast (1-\frac{1}{3})=6\ast \frac{1}{2}\ast \frac{2}{3}=2$

---

下面证明为什么 $\varphi (n)=n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ...\ast (1-\frac{1}{p_k})$：
求1 ~ n中与n互质的数的个数，相当于求从1 ~ n中去掉所有因数有$p_1,p_2,...,p_k$的数后剩下多少个数。从一个集合里去掉若干元素，可以想到容斥原理。
1、先从1 ~ n中去掉$p_1,p_2,...,p_k$的所有倍数，即$n-\frac{n}{p_1}-\frac{n}{p_2}-...-\frac{n}{p_k}$。此处所有$\frac{n}{p_k}$都是向下取整。
在这个步骤中，当一个数同时是多个$p$的倍数时，它会被去除很多次。

2、加上所有$p_i\ast p_j$的倍数，即$n-\frac{n}{p_1}-\frac{n}{p_2}-...-\frac{n}{p_k}+\frac{n}{p_1{p}_2}+\frac{n}{p_1{p}_3}+...+\frac{n}{p_{k-1}{p}_k}$
这一步加回去了同时是2个$p$的倍数的数，那么还要接着减去同时是3个$p$的倍数的数。
不妨设这个数是$p_1\ast p_2\ast p_3$，在步骤1，这个数被减了3次，在步骤2，这个数又被$p_1{p}_2,p_1{p}_3,p_2{p}_3$加了三次，不减也不加，但这个数最终应该被减1次。

3、加上所有$p_i\ast p_j\ast p_k$的倍数$n-\frac{n}{p_1}-\frac{n}{p_2}-...-\frac{n}{p_k}+\frac{n}{p_1{p}_2}+\frac{n}{p_1{p}_3}+...+\frac{n}{p_{k-1}{p}_k}-\frac{n}{p_1{p}_2{p}_3}-\frac{n}{p_1{p}_2{p}_4}-...-\frac{n}{p_{k-2}{p}_{k-1}{p}_k}$。同理还要接着加上同时是4个$p$的倍数的数。

依次类推，最终
$n-\frac{n}{p_1}-\frac{n}{p_2}-...-\frac{n}{p_k}+\frac{n}{p_1{p}_2}+\frac{n}{p_1{p}_3}+...+\frac{n}{p_{k-1}{p}_k}-\frac{n}{p_1{p}_2{p}_3}-\frac{n}{p_1{p}_2{p}_4}-...-\frac{n}{p_{k-2}{p}_{k-1}{p}_k}+...$
$=n(1-\frac{1}{p_1}-\frac{1}{p_2}-...-\frac{1}{p_k}+\frac{1}{p_1{p}_2}+\frac{1}{p_1{p}_3}+...+\frac{1}{p_{k-1}{p}_k}-\frac{1}{p_1{p}_2{p}_3}-\frac{1}{p_1{p}_2{p}_4}-...-\frac{1}{p_{k-2}{p}_{k-1}{p}_k}+...\frac{1}{p_1{p}_2...p_k})$
这个式子就是$n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ...\ast (1-\frac{1}{p_k})$展开后的式子，通过对比各项之前的系数即可知道，如$\frac{1}{p_1{p}_2}$系数是$+1$；$n(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_k})$中$-\frac{1}{p_1}$和$-\frac{1}{p_2}$相乘，其他项全取1，系数也是$+1$。

---

该题用公式求即可，分解质因数时间复杂度$\sqrt{a}=10^4$，再乘上n=100是$10^6$，符合要求。

#include <iostream>
using namespace std;

int n, a;

int main() {
    scanf("%d", &n);
    while(n -- ) {
        scanf("%d", &a);
        double ans = a;
        for (int i = 2; i * i <= a; i ++ ) {
            int prime = 0;
            while(a % i == 0) {
                a /= i;
                prime = i;
            }
            if (prime) ans *= 1 - 1.0 / i; // 只有是质数时才进行运算 
        }
        if (a > 1) ans *= 1 - 1.0 / a; // a本身是质数 
        printf("%.0lf\n", ans);
    }
    return 0;
}

实际上不需要小数运算，只要将$ans\ast (1-1/p_k)$改成$ans/p_k\ast (p_k-1)$即可避免小数。

    while(n -- ) {
        scanf("%d", &a);
        int ans = a;
        for (int i = 2; i * i <= a; i ++ ) {
            if (a % i == 0) {
            	ans = ans / i * (i - 1);
            	while (a % i == 0) a /= i;
			} 
        }
        if (a > 1) ans = ans / a * (a - 1); // a本身是质数 
        printf("%d\n", ans);
    }