使用Pollard_rho算法分解质因数

2024-04-02 19:20

本文主要是介绍使用Pollard_rho算法分解质因数,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

分解质因数的朴素算法

最简单的算法即为从 [2, sqrt(N)] 进行遍历。

vector<int> breakdown(int N) {vector<int> result;for (int i = 2; i * i <= N; i++) {if (N % i == 0) {  // 如果 i 能够整除 N,说明 i 为 N 的一个质因子。while (N % i == 0) N /= i;result.push_back(i);}}if (N != 1) {  // 说明再经过操作之后 N 留下了一个素数result.push_back(N);}return result;
}

这个算法当 n 是质数时拥有最差时间复杂度 O(n) ,不过可以先用米勒-拉宾特判一下 n 是不是质数,这样的话最差时间复杂度就是当 n 是质数的平方时的 O(sqrt (n)) 了

Pollard_rho将一个数分解为质因数相乘的形式

如:200 = 22255
Pollard_rho算法是找到一个数的最大质因数,我参考此算法进行修改,实现了将一个数分解为质因数相乘形式的程序。

#include <iostream>
#include <stdlib.h>
#include <vector>int t;
long long max_factor, n;
std::vector<long long> factor; // 存储质因数long long gcd(long long a, long long b) {if (b == 0) return a;return gcd(b, a % b);
}long long quick_pow(long long x, long long p, long long mod) {  // 快速幂long long ans = 1;while (p) {if (p & 1) ans = (__int128)ans * x % mod;x = (__int128)x * x % mod;p >>= 1;}return ans;
}bool Miller_Rabin(long long p) {  // 判断素数if (p < 2) return 0;if (p == 2) return 1;if (p == 3) return 1;long long d = p - 1, r = 0;while (!(d & 1)) ++r, d >>= 1;  // 将d处理为奇数std::vector<long long> ud = {2, 325, 9375, 28178, 450775, 9780504, 1795265022};for (long long a:ud) {long long x = quick_pow(a, d, p);if (x == 1 || x == p - 1) continue;for (int i = 0; i < r - 1; ++i) {x = (__int128)x * x % p;if (x == p - 1) break;}if (x != p - 1) return 0;}return 1;
}long long Pollard_Rho(long long x) {long long s = 0, t = 0;long long c = (long long)rand() % (x - 1) + 1;int step = 0, goal = 1;long long val = 1;for (goal = 1;; goal *= 2, s = t, val = 1) {  // 倍增优化for (step = 1; step <= goal; ++step) {t = ((__int128)t * t + c) % x;val = (__int128)val * abs(t - s) % x;if ((step % 127) == 0) {long long d = gcd(val, x);if (d > 1) return d;}}long long d = gcd(val, x);if (d > 1) return d;}
}void fac(long long x) {if (x <= max_factor || x < 2) return;if (Miller_Rabin(x)) {              // 如果x为质数max_factor = std::max(max_factor, x);  // 更新答案return;}long long p = x;while (p >= x) p = Pollard_Rho(x);  // 使用该算法while ((x % p) == 0) x /= p;fac(x), fac(p);  // 继续向下分解x和p
}void decompose(long long n) { // 将n分解为质因数相乘的形式srand((unsigned)time(NULL));max_factor = 0;fac(n);if(max_factor == n) { // 最大的质因数是自己factor.push_back(max_factor);}else {factor.push_back(max_factor);n /= max_factor;decompose(n);}
}int main() {std::cout << "----start decompose n: [exit please input number 0 !]----" << std::endl;while(1) {std::cout << "input number: ";std::cin >> n;if(n == 0) {std::cout << "stop and exit program!" << std::endl;break;}decompose(n);std::cout << n << " = ";for(std::vector<long long>::iterator it = factor.begin(); it != factor.end()-1; ++it) {std::cout << *it << "*";}std::cout << *(factor.end()-1) << std::endl;std::vector<long long>().swap(factor);}return 0;
}

程序功能展示
在这里插入图片描述

参考文章

[1] https://oi-wiki.org/math/number-theory/pollard-rho/
[2] https://zhuanlan.zhihu.com/p/267884783
[3] Miller Rabin素数判定

这篇关于使用Pollard_rho算法分解质因数的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/870837

相关文章

postgresql使用UUID函数的方法

《postgresql使用UUID函数的方法》本文给大家介绍postgresql使用UUID函数的方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录PostgreSQL有两种生成uuid的方法。可以先通过sql查看是否已安装扩展函数,和可以安装的扩展函数

如何使用Lombok进行spring 注入

《如何使用Lombok进行spring注入》本文介绍如何用Lombok简化Spring注入,推荐优先使用setter注入,通过注解自动生成getter/setter及构造器,减少冗余代码,提升开发效... Lombok为了开发环境简化代码,好处不用多说。spring 注入方式为2种,构造器注入和setter

MySQL中比较运算符的具体使用

《MySQL中比较运算符的具体使用》本文介绍了SQL中常用的符号类型和非符号类型运算符,符号类型运算符包括等于(=)、安全等于(=)、不等于(/!=)、大小比较(,=,,=)等,感兴趣的可以了解一下... 目录符号类型运算符1. 等于运算符=2. 安全等于运算符<=>3. 不等于运算符<>或!=4. 小于运

使用zip4j实现Java中的ZIP文件加密压缩的操作方法

《使用zip4j实现Java中的ZIP文件加密压缩的操作方法》本文介绍如何通过Maven集成zip4j1.3.2库创建带密码保护的ZIP文件,涵盖依赖配置、代码示例及加密原理,确保数据安全性,感兴趣的... 目录1. zip4j库介绍和版本1.1 zip4j库概述1.2 zip4j的版本演变1.3 zip4

Python 字典 (Dictionary)使用详解

《Python字典(Dictionary)使用详解》字典是python中最重要,最常用的数据结构之一,它提供了高效的键值对存储和查找能力,:本文主要介绍Python字典(Dictionary)... 目录字典1.基本特性2.创建字典3.访问元素4.修改字典5.删除元素6.字典遍历7.字典的高级特性默认字典

使用Python构建一个高效的日志处理系统

《使用Python构建一个高效的日志处理系统》这篇文章主要为大家详细讲解了如何使用Python开发一个专业的日志分析工具,能够自动化处理、分析和可视化各类日志文件,大幅提升运维效率,需要的可以了解下... 目录环境准备工具功能概述完整代码实现代码深度解析1. 类设计与初始化2. 日志解析核心逻辑3. 文件处

一文详解如何使用Java获取PDF页面信息

《一文详解如何使用Java获取PDF页面信息》了解PDF页面属性是我们在处理文档、内容提取、打印设置或页面重组等任务时不可或缺的一环,下面我们就来看看如何使用Java语言获取这些信息吧... 目录引言一、安装和引入PDF处理库引入依赖二、获取 PDF 页数三、获取页面尺寸(宽高)四、获取页面旋转角度五、判断

C++中assign函数的使用

《C++中assign函数的使用》在C++标准模板库中,std::list等容器都提供了assign成员函数,它比操作符更灵活,支持多种初始化方式,下面就来介绍一下assign的用法,具有一定的参考价... 目录​1.assign的基本功能​​语法​2. 具体用法示例​​​(1) 填充n个相同值​​(2)

Spring StateMachine实现状态机使用示例详解

《SpringStateMachine实现状态机使用示例详解》本文介绍SpringStateMachine实现状态机的步骤,包括依赖导入、枚举定义、状态转移规则配置、上下文管理及服务调用示例,重点解... 目录什么是状态机使用示例什么是状态机状态机是计算机科学中的​​核心建模工具​​,用于描述对象在其生命

使用Python删除Excel中的行列和单元格示例详解

《使用Python删除Excel中的行列和单元格示例详解》在处理Excel数据时,删除不需要的行、列或单元格是一项常见且必要的操作,本文将使用Python脚本实现对Excel表格的高效自动化处理,感兴... 目录开发环境准备使用 python 删除 Excphpel 表格中的行删除特定行删除空白行删除含指定