AcWing.873.欧拉函数

给定 n n n 个正整数 ai,请你求出每个数的欧拉函数。

欧拉函数的定义1∼ N N N 中与 N N N 互质的数的个数被称为欧拉函数,记为 ϕ(N)。
若在算数基本定理中, N N N=p1a1p2a2…pmam, 则: ϕ ( N ) ϕ(N) ϕ(N) = N N N× p 1 − 1 p 1 \frac{p1 −1}{p1} p1p11× p 2 − 1 p 2 \frac{p2−1}{p2} p2p21×…× p m − 1 p m \frac{pm−1}{pm} pmpm1

输入格式
第一行包含整数 n n n

接下来 n n n 行,每行包含一个正整数 ai

输出格式
输出共 n n n 行,每行输出一个正整数 ai 的欧拉函数。

数据范围
1 ≤ n ≤ 100 , 1 ≤ a i ≤ 2 × 109 1≤n≤100,1≤ai≤2×109 1n100,1ai2×109

输入样例:

3
3
6
8

输出样例:

2
2
4

对于 n n n,分解质因数后为: N N N=p1a1p2a2…pmam,那么欧拉函数的结果就是 ϕ ( N ) ϕ(N) ϕ(N) = N N N× p 1 − 1 p 1 \frac{p1 −1}{p1} p1p11× p 2 − 1 p 2 \frac{p2−1}{p2} p2p21×…× p m − 1 p m \frac{pm−1}{pm} pmpm1,即运用此公式进行求欧拉函数的过程

代码:

#include<iostream>
using namespace std;

int main() {
   
	int n; cin >> n;
	
	while (n--) {
   
		int a; cin >> a;

		int res = a;	//答案初始为a
		for (int i = 2; i <= a / i; i++) {
   	//分解质因数
			if (a % i == 0) {
   	//如果a整除i,那么是质因子
				res = res / i * (i - 1);	//( res = res * (1 - 1/i) )化简避免小数
				while (a % i == 0) a /= i;	//消倍数
			}
		}
		if (a > 1) res = res / a * (a - 1);	//如果有余下的质因子

		cout << res << endl;
	}
	return 0;
}

筛法求欧拉函数代码:

#include<iostream>
using namespace std;
const int N = 1e6 + 10;

int primes[N], cnt;	//存质数
bool st[N];			//判断哪个数被筛掉
int phi[N];			//存欧拉函数值

long long get_eulers(int n) {
   
	phi[1] = 1;	//1的欧拉函数值为1
	for (int i = 2; i <= n; i++) {
   	
		if (!st[i]) {
   	//如果没有被筛过
			primes[cnt++] = i;	//说明i为质数
			phi[i] = i - 1;		//如果i是质数,那么欧拉函数值为i-1,因为1~i-1都与他互质
		}
		for (int j = 0; primes[j] <= n / i; j++) {
   
			st[primes[j] * i] = 1;
			if (i % primes[j] == 0) {
   	//primes[j]是i的质因子
				//phi(pj * i) = pj * phi(i)
				phi[primes[j] * i] = primes[j] * phi[i];
				break;
			}
			//如果不是,那么为( phi(pj * i) = phi(i) * ( pj - 1 ) )
			phi[primes[j] * i] = phi[i] * (primes[j] - 1);
		}
	}

	long long res = 0;
	for (int i = 1; i <= n; i++)res += phi[i];
	return res;

}

int main() {
   
	int n; cin >> n;
	
	cout << get_eulers(n) << endl;

	return 0;
}

相关推荐

  1. AcWing.873.函数

    2024-01-30 04:14:01       29 阅读
  2. 浅谈函数

    2024-01-30 04:14:01       27 阅读
  3. 【算法基础 & 数学】函数

    2024-01-30 04:14:01       28 阅读

最近更新

  1. TCP协议是安全的吗?

    2024-01-30 04:14:01       18 阅读
  2. 阿里云服务器执行yum,一直下载docker-ce-stable失败

    2024-01-30 04:14:01       19 阅读
  3. 【Python教程】压缩PDF文件大小

    2024-01-30 04:14:01       18 阅读
  4. 通过文章id递归查询所有评论(xml)

    2024-01-30 04:14:01       20 阅读

热门阅读

  1. VUE就是最强!

    2024-01-30 04:14:01       34 阅读
  2. 十个鼠标事件

    2024-01-30 04:14:01       40 阅读
  3. 1.基于C#的Dbf读写(文件结构概述)

    2024-01-30 04:14:01       32 阅读
  4. cpp-stub 打桩失败

    2024-01-30 04:14:01       40 阅读
  5. 题解:CF1922C(Closest Cities)

    2024-01-30 04:14:01       35 阅读
  6. 面试 CSS 框架八股文十问十答第一期

    2024-01-30 04:14:01       42 阅读
  7. C++算法学习心得七.贪心算法(2)

    2024-01-30 04:14:01       21 阅读
  8. Quartz在spring boot项目中重启后不能继续执行问题

    2024-01-30 04:14:01       32 阅读
  9. OpenSSH 9.6/9.6p1 (2023-12-18)的发布说明(中文译文)

    2024-01-30 04:14:01       41 阅读
  10. Apache孵化器领路人与导师的职责

    2024-01-30 04:14:01       46 阅读
  11. jQuery 和 Zepto 的区别? 各自的使用场景?

    2024-01-30 04:14:01       36 阅读