Acwing_874（筛法求欧拉函数）

笨笨sg

已于 2024-03-13 20:45:36 修改

阅读量382

点赞数 5

分类专栏： Acwing算法基础课文章标签：线性筛欧拉函数

于 2024-03-13 20:45:17 首次发布

本文链接：https://blog.csdn.net/a131529/article/details/136690995

版权

Acwing算法基础课专栏收录该内容

107 篇文章 0 订阅

订阅专栏

本文介绍了如何使用筛法求欧拉函数的C++代码示例，通过OláPrimes算法计算1到n范围内的欧拉函数总和，适用于AcWing题库中的相关问题。

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

原题链接：

874. 筛法求欧拉函数 - AcWing题库

题解：

欧拉函数共有以下三种情况：

参考讲解：

AcWing 874. $\Huge\color{gold}{筛法求欧拉函数}$ - AcWing

代码：

#include<bits/stdc++.h>
using namespace std;
int n;
using ULL = unsigned long long;
ULL res;
const int N = 1e6 + 10;
int p[N], s[N], phi[N], cnt;
//p：下标-质数（从小到大）
//s：数-是不是质数
//phi：质数-欧拉函数值
//cnt：已经找到的质数个数

void ola_primes(int x) {
	phi[1] = 1;
	for (int i = 2;i <= x;i++) {
		if (!s[i]) {//当前数为素数的话
			p[cnt++] = i;
			phi[i] = i - 1;
		}
		for (int j = 0;p[j] <= x / i;j++) {
			s[i * p[j]] = 1;
			if (i % p[j] == 0) {//两者不互质时，p[j]可以直接拎出来
				phi[i * p[j]] = phi[i] * p[j];
				break;
			}
			phi[i * p[j]] = phi[i] * (p[j] - 1);//两者互质时，直接拆分
		}
	}
}

int main() {
	cin >> n;
	ola_primes(n);
	for (int i = 1;i <= n;i++) res += phi[i];
	cout << res;
}