题目描述
给定 n n n 个正整数 a i a_i ai,请你求出每个数的欧拉函数。
输入格式
第一行包含整数 n n n 。
接下来 n n n 行,每行包含一个正整数 a i a_i ai。
输出格式
输出共 n n n行,每行输出一个正整数 a i a_i ai 的欧拉函数。
数据范围
1
≤
n
≤
100
1≤n≤100
1≤n≤100,
1
≤
a
i
≤
2
×
1
0
9
1≤a_i≤2×10^9
1≤ai≤2×109
样例
输入样例:
3
3
6
8
输出样例:
2
2
4
定义
φ ( n ) \varphi(n) φ(n): 1 1 1~ n n n中与 n n n互质的数的个数
求法
若一个数
N
N
N可以分解质因子
N
=
p
1
α
1
p
2
α
2
.
.
.
p
m
α
m
N=p_1^{\alpha_1}p_2^{\alpha_2}...p_m^{\alpha_m}
N=p1α1p2α2...pmαm
则
φ
(
N
)
=
N
×
p
1
−
1
p
1
×
p
2
−
1
p
2
×
.
.
.
×
p
m
−
1
p
m
\varphi(N) = N×\frac{p_1-1}{p_1}×\frac{p_2-1}{p_2}×...×\frac{p_m-1}{p_m}
φ(N)=N×p1p1−1×p2p2−1×...×pmpm−1
推导过程(利用容斥原理)
利用分解后的质因子
- 从1~n中去掉 p 1 , p 2 , . . . , p k p_1,p_2,...,p_k p1,p2,...,pk的所有的倍数
r e s = N − N p 1 − N p 2 − . . . res=N-\frac{N}{p_1}-\frac{N}{p_2}-... res=N−p1N−p2N−...
- 加上所有 p i × p j p_i×p_j pi×pj的倍数(公共的倍数在前面被减了两次)
r e s + = N p 1 p 2 + N p 1 p 3 . . . res += \frac{N}{p_1p_2} + \frac{N}{p_1p_3}... res+=p1p2N+p1p3N...
- 减去所有 p i × p j × p k p_i×p_j×p_k pi×pj×pk的倍数
r e s − = N p 1 p 2 p 3 − . . . res -= \frac{N}{p_1p_2p_3}-... res−=p1p2p3N−...
显然,将上面的公式展开后就是容斥原理的式子
时间复杂度
每个数的判断复杂度为
O
(
n
)
O(\sqrt n)
O(n)
一共100个数据,最坏情况下达到
1
0
6
10^6
106
代码
#include<iostream>
using namespace std;
int phi(int x){
int res = x;
for(int i = 2; i <= x / i; i++){
if(x % i == 0){
res = res / i * (i - 1);
while(x % i == 0) x /= i;
}
}
if(x > 1) res = res / x * (x - 1);
return res;
}
int main(){
int n;
cin >> n;
while(n--){
int x;
cin >> x;
cout << phi(x) << endl;
}
return 0;
}