筛法求欧拉函数（模板题）

chycjq

于 2024-04-16 11:08:23 发布

阅读量276

点赞数 12

文章标签：算法 c++ 数据结构

本文链接：https://blog.csdn.net/chq66666/article/details/137817470

版权

本文介绍了一种计算1到n范围内每个数的欧拉函数之和的方法，利用质数分解的思想，通过编程实现，主要关注了欧拉函数的性质和计算效率。

摘要由CSDN通过智能技术生成

给定一个正整数 n，求 1∼n 中每个数的欧拉函数之和。

输入格式

共一行，包含一个整数 n。

输出格式

共一行，包含一个整数，表示 1∼n 中每个数的欧拉函数之和。

数据范围

1≤n≤106

输入样例：

输出样例：

思路：

质数i的欧拉函数即为phi[i] = i - 1。因为1 ~ i−1均与i互质，所以共i−1个。
phi[primes[j] * i]分为两种情况：
① i % primes[j] == 0时：primes[j]是i的最小质因子，也是primes[j] * i的最小质因子，因此1 - 1 / primes[j]这一项在phi[i]中计算过了，只需将基数N修正为primes[j]倍，最终结果为phi[i] *primes[j]。
② i % primes[j] != 0：primes[j]不是i的质因子，只是primes[j] * i的最小质因子，因此不仅需要将基数N修正为primes[j]倍，还需要补上1 - 1 / primes[j]这一项，因此最终结果phi[i] * (primes[j] - 1)。

代码：

#include<iostream>

using namespace std;
const int N = 1e6 + 10;
int cnt = 0, prime[N], n, x;
int phi[N];
bool st[N];
long long get_prime(int x) {
	phi[1] = 1;
	for (int i = 2; i <= x; i++) {
		if (!st[i]) {
			prime[cnt++] = i;
			phi[i] = i - 1;
		}
		for (int j = 0; prime[j] <= x / i; j++) {
			st[prime[j] * i] = true;
			
			if (i % prime[j] == 0) {
				phi[prime[j] * i] = phi[i] * prime[j];
				break;
			}
			phi[prime[j] * i] = phi[i] * (prime[j] - 1);
		}
	}
	long long res = 0;
	for (int i = 1; i <= x; i++)res += phi[i];
	return res;
}


int main() {
	cin >> n;
	cout << get_prime(n) << endl;
	return 0;
}