欧拉函数实现与分析

问题

https://www.acwing.com/problem/content/875/

分析

对于一个正整数 $n$ 欧拉函数指小于或等于 $n$ 的正整数中与 $n$ 互质的数的数目，记作 $\phi (n)$ ，其中有 $\phi (1)=1$

欧拉函数本身是一个积性函数，$\phi (m\ast n)=\phi (m)\ast \phi (n)$

已知对于每一个正整数 $n$ 而言都有 $n=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \cdots \ast p_k^{{\alpha }_k}$

则可以得到 $\phi (n)=\phi (p_1^{{\alpha }_1})\ast \phi (p_2^{{\alpha }_2})\ast \cdots \ast \phi (p_k^{{\alpha }_k})$

对于 $\phi (p_1^{{\alpha }_1})$ 与 $p_1^{{\alpha }_1}$ 不互质的数有 $p_1$、$2\ast p_1$、$3\ast p_1$、$\cdots$ 、$p_1^{{\alpha }_1-1}\ast p_1$ 共 $p_1^{{\alpha }_{1-1}}$ 个

可以得到 $\phi (p_1^{{\alpha }_1})=p_1^{{\alpha }_1}-p_1^{{\alpha }_1-1}=p_1^{{\alpha }_1}(1-\frac{1}{p_1})$

最终可以得到
$$\phi (n)=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \cdots \ast p_k^{{\alpha }_k}\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast \cdots \ast (1-\frac{1}{p_k})$$
最终化简可得
$$\phi (n)=n\ast \prod _{i=1}^k(1-\frac{1}{p_i})$$

代码

#include<iostream>

using namespace std ;

int main()
{
    int n ;
    cin >> n ;
    while(n --)
    {
        int a ;
        cin >> a ;
        int res = a ;
        for(int i = 2 ; i <= a / i ; i ++)
        {
            if(a % i == 0)
            {
                // 此处需要先除 i 再乘 i - 1 要不然可能无法整除
                res = res / i * (i - 1) ;
                while(a % i == 0)
                    a /= i ;
            }
            if(a > 1)
                res = res / a * (a - 1) ;
        }
        
        cout << res << endl ;
    }
    
    return 0 ;
}