acwing算法基础之数学知识--求数a的欧拉函数值phi(a)

目录

  • [1 基础知识](#1 基础知识)
  • [2 模板](#2 模板)
  • [3 工程化](#3 工程化)

1 基础知识

数a的欧拉函数 ϕ ( a ) \phi(a) ϕ(a):表示1~n中与n互质的数的个数。其中两个数互质,是指这两个数的最大公约数为1。

根据定义,我们可以写出如下方法,

cpp 复制代码
int gcd(int a, int b) {
	return b ? gcd(b, a % b) : a;
}

int phi(int a) {
	int res = 0;
	for (int i = 1; i <= a; ++i) {
		if (gcd(i, a) == 1) {
			res += 1;
		}
	}
	return res;
}

但存在更快的求解方法,见如下关键步骤:

  1. 对数a进行分解质因子操作。
    a = p 1 α 1 ⋅ p 2 α 2 ⋯ p k α k a=p_1^{\alpha_1} \cdot p_2^{\alpha_2}\cdots p_k^{\alpha_k} a=p1α1⋅p2α2⋯pkαk
cpp 复制代码
unordered_map<int,int> get_prime_divisors(int a) {
	unordered_map<int,int> mp;
	for (int i = 2; i <= a / i; ++i) {
		if (a % i == 0) {
			int s = 0;
			while (a % i == 0) {
				a /= i;
				s++;
			}
			mp[i] = s;
		}
	}
	if (a > 1) mp[a] = 1;
	return mp;
}
  1. 计算数a的欧拉函数,
    ϕ ( a ) = a ⋅ ( 1 − 1 p 1 ) ⋅ ( 1 − 1 p 2 ) ⋯ ( 1 − 1 p k ) \phi(a)=a\cdot (1-\frac{1}{p_1}) \cdot (1-\frac{1}{p_2}) \cdots (1-\frac{1}{p_k}) ϕ(a)=a⋅(1−p11)⋅(1−p21)⋯(1−pk1)
cpp 复制代码
int phi(int a, unordered_map<int,int> mp) {
	int res = a;
	for (auto [x, y] : mp) {
		res = res / x * (x - 1);
	}
	return res;
} 

可以将以上两步合并,请看如下代码,

cpp 复制代码
int phi(int a) {
	int res = a;
	for (int i = 2; i <= a / i; ++i) {
		if (a % i == 0) {
			res = res / i * (i - 1);
			while (a % i == 0) {
				a /= i;
			}
		}
	}
	if (a > 1) {
		res = res / a * (a - 1);
	}
	return res;
}

2 模板

cpp 复制代码
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;
}

3 工程化

题目1:输入n个数,请分别求出它们的欧拉函数值。

cpp 复制代码
#include <iostream>

using namespace std;

int main() {
    int n;
    cin >> n;
    
    while (n--) {
        int x;
        cin >> 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);
        cout << res << endl;
    }
    
    return 0;
}
相关推荐
凌肖战2 小时前
力扣网编程55题:跳跃游戏之逆向思维
算法·leetcode
88号技师3 小时前
2025年6月一区-田忌赛马优化算法Tianji’s horse racing optimization-附Matlab免费代码
开发语言·算法·matlab·优化算法
ゞ 正在缓冲99%…3 小时前
leetcode918.环形子数组的最大和
数据结构·算法·leetcode·动态规划
Kaltistss4 小时前
98.验证二叉搜索树
算法·leetcode·职场和发展
知己如祭4 小时前
图论基础(DFS、BFS、拓扑排序)
算法
mit6.8244 小时前
[Cyclone] 哈希算法 | SIMD优化哈希计算 | 大数运算 (Int类)
算法·哈希算法
c++bug4 小时前
动态规划VS记忆化搜索(2)
算法·动态规划
哪 吒4 小时前
2025B卷 - 华为OD机试七日集训第5期 - 按算法分类,由易到难,循序渐进,玩转OD(Python/JS/C/C++)
python·算法·华为od·华为od机试·2025b卷
军训猫猫头5 小时前
1.如何对多个控件进行高效的绑定 C#例子 WPF例子
开发语言·算法·c#·.net
success5 小时前
【爆刷力扣-数组】二分查找 及 衍生题型
算法