【ACWing】873. 欧拉函数

题目地址：

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

给定$n$个正整数$a_i$，求$\varphi (a_i)$。$\varphi (x)$是欧拉函数，表示小于等于$x$的且与之互素的数的个数。

输入格式：
第一行包含整数$n$。接下来$n$行，每行包含一个正整数$a_i$。

输出格式：
输出共$n$行，每行输出一个正整数$a_i$的欧拉函数。

数据范围：
$1\le n\le 100$
$1\le a_i\le 2\times 10^9$

若正整数$N$有素因子分解$N=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}...p_k^{{\alpha }_k}$，那么有公式$$\varphi (N)=N(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_k})$$这可以通过容斥原理来推导。我们看哪些数与$N$不互素，显然含因子$p_1,...,p_k$的数与$N$不互素，我们只需要求小于等于$N$且含因子$p_1,...,p_k$的数的个数即可。由容斥原理，个数是$$N(\sum _i\frac{1}{p_i}-\sum _{i_1<i_2}\frac{1}{p_{i_1}{p}_{i_2}}+\sum _{i_1<i_2<i_3}\frac{1}{p_{i_1}{p}_{i_2}{p}_{i_3}}-...+(-1{)}^{k-1}\frac{1}{p_1{p}_2...p_k})$$所以$$\begin{gathered} \varphi (N)=N(1-\sum _i\frac{1}{p_i}+\sum _{i_1<i_2}\frac{1}{p_{i_1}{p}_{i_2}}-\sum _{i_1<i_2<i_3}\frac{1}{p_{i_1}{p}_{i_2}{p}_{i_3}}-...+(-1{)}^k\frac{1}{p_1{p}_2...p_k}) \\ =N(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_k}) \end{gathered}$$而代码方面，只需要套用分解质因数的代码即可。代码如下：

#include <iostream>
using namespace std;

int euler(int x) {
  int res = x;
  for (int i = 2; i <= x / i; i++) {
    if (x % i == 0) {
      // 这里i是x的素因子，所以可以先除再乘
      res = res / i * (i - 1);
      while (x % i == 0) x /= i;
    }
  }

  if (x >= 2) res = res / x * (x - 1);
  return res;
}

int main() {
  int n;
  cin >> n;
  while (n--) {
    int a;
    cin >> a;
    cout << euler(a) << endl;
  }
}

每次询问时间复杂度$O(\sqrt{x})$，空间$O(1)$。