连分数因子分解法——C语言实现

2024-09-07 08:04

本文主要是介绍连分数因子分解法——C语言实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

参考网址:连分数分解法寻找整数的因子(Python)-CSDN博客

大数运算:C语言实现 大数运算 加减乘除模运算 超详细_64编程 加减乘除取模 复杂运算-CSDN博客

连分数因子分解法‌是一种用于大整数因子分解的算法,它是计算数论中的一个重要方法。连分数因子分解法通过寻找x2≡y2 (mod p)x2≡y2 (mod p)的形式来分解N。具体来说,这种方法涉及到计算N的简单连分数展开,并通过组合得到的同余式来寻找或组合得到因子。

连分数因子分解法的步骤大致如下:

  1. 计算NN的简单连分数展开‌,同时计算PkPk​和QkQk​序列,其中P0=a0,P1=a0a1+1,Pk=akPk−1+Pk−2P0​=a0​,P1​=a0​a1​+1,Pk​=ak​Pk−1​+Pk−2​,Q0=1,Q1=a1,Qk=akQk−1Q0​=1,Q1​=a1​,Qk​=ak​Qk−1​。
  2. 通过连分数展开寻找或组合得到同余式‌:x2≡y2 (mod p)x2≡y2 (mod p),其中pp是素数,xx和yy是整数。
  3. 利用找到的同余式‌,通过计算和组合这些同余式,可以进一步分解给定的整数NN。

这种方法在处理大整数时特别有用,因为它提供了一种有效的手段来分解那些无法通过简单试除法直接分解的大数。连分数因子分解法与二次筛法、数域筛法等方法一起,构成了处理大数因子分解的一系列工具‌1。

此外,连分数因子分解法在实际应用中也有其局限性,尤其是在处理非常大的数时,其计算复杂度和所需资源可能会成为限制因素。因此,尽管这种方法在理论和实践上都有其价值,但在实际应用中还需要考虑其他更高效的算法或技术的结合使用‌

我使用C语言实现的连分数因子分解法部分代码,由于是测试代码,目前还没有优化,存在一些问题:

/*
-----------------------连分数因子分解法----------------------实现时间:2024年9月5日作者:赵良军
*/
//求完全平方数
int is_square(unsigned char* seq_Q){unsigned char* p = my_sqrt(seq_Q);p = mul(p, p);//printf( "s = %s\n",p );if(maxer( p,seq_Q ) == 0){return 1;	}	  return 0;
}
//求最大公约数
unsigned char* my_gcd(unsigned char* a,unsigned char* b){unsigned char* c = NULL;while (a){   	c = mol(b,a);   b = a;a = c;if( maxer( a, "0" ) == 0 || maxer( c, "0" ) == 0 ){return b;}     	} return NULL;   
}
void cfrac(char* n){unsigned char seq_P[3][BUFLEN] = {{'\0'},{'\0'},{'\0'}};unsigned char seq_Q[3][BUFLEN] = {{'\0'},{'\0'},{'\0'}};unsigned char seq_a[3][BUFLEN] = {{'\0'},{'\0'},{'\0'}};unsigned char seq_p[3][BUFLEN] = {{'\0'},{'\0'},{'\0'}};strncpy(seq_P[0],"0",1);//strncpy(seq_P[1],"0",1);//strncpy(seq_P[2],"0",1);strncpy(seq_Q[0],"1",1);//strncpy(seq_Q[1],"1",1);//strncpy(seq_Q[2],"1",1);unsigned char* p = NULL;unsigned char buf[BUFLEN] = {'\0'};printf("N = %s\n",n);unsigned char* a0 = my_sqrt(n);strcpy(seq_a[0],a0);   strncpy(seq_p[0],"0",1);strcpy(seq_p[1],a0 ); unsigned char i[BUFLEN] = {'\0'};strncpy(i,"1",1);unsigned char t[BUFLEN] = {'\0'};while(1){//P_{k},Q_{k}序列printf("Start------------------------->\n");//P序列的模数p = mul(seq_a[0], seq_Q[0]);strcpy(buf,p);p = sub(buf,seq_P[0]);strcpy( seq_P[0],p );printf("P_M = %s\n",p );p = mul(seq_P[0], seq_P[0]);strcpy(buf,p);//printf("M0 = %s\n",p );p = sub(n,buf);strcpy(buf,p);//printf("M1 = %s\n",p );p = divi( buf,seq_Q[0]);strcpy(seq_Q[0],p);//printf("M2 = %s\n",p );//p = mol(buf,seq_Q[0]);//strcpy(buf,p);//printf("M3 = %s\n",p );p = add(seq_P[0],my_sqrt(n));strcpy(buf,p);//printf("M3 = %s\n",p );p = divi( buf,seq_Q[0]);strcpy(seq_a[0],p);//printf("M4 = %s\n",seq_a[0] );//p_{k}序列if( maxer(i,"1") == 0 ){p = mul(a0, seq_a[0]);strcpy(buf,p);p = add(buf,"1");strcpy(seq_p[2],p);}else{strcpy(seq_p[0],seq_p[1]);strcpy(seq_p[1],seq_p[2]);p = mul(seq_a[0], seq_p[1]);strcpy(buf,p);p = add( buf,seq_p[0] );strcpy(seq_p[2],p);}if( maxer( mol( i, "2" ), "0" ) == 0 && is_square(seq_Q[0]) ){unsigned char* s = my_sqrt(seq_Q[0]);unsigned char* p_p = NULL;printf("S = %s\n",s );p = add(seq_p[1],s);p_p = my_gcd(p,n);     	  printf("P+ = %s\n",p_p);p = sub(seq_p[1],s);p_p = my_gcd(p,n);     	  printf("P- = %s\n",p_p);if(p_p != NULL && maxer( p_p, "1" ) != 0){	      	  	  break;}} if( maxer( mol( i, "1000" ), "0" ) == 0 ){printf("Run to %s\n",i);}p = add(i,"1");strcpy(i,p);}
}int main(){//unsigned char* n = "70191551";//2**67-1unsigned char* n = "147573952589676412927";cfrac(n);}

这篇关于连分数因子分解法——C语言实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

pandas中位数填充空值的实现示例

《pandas中位数填充空值的实现示例》中位数填充是一种简单而有效的方法,用于填充数据集中缺失的值,本文就来介绍一下pandas中位数填充空值的实现,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是中位数填充?为什么选择中位数填充?示例数据结果分析完整代码总结在数据分析和机器学习过程中,处理缺失数

Golang HashMap实现原理解析

《GolangHashMap实现原理解析》HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持高效的插入、查找和删除操作,:本文主要介绍GolangH... 目录HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持

Pandas使用AdaBoost进行分类的实现

《Pandas使用AdaBoost进行分类的实现》Pandas和AdaBoost分类算法,可以高效地进行数据预处理和分类任务,本文主要介绍了Pandas使用AdaBoost进行分类的实现,具有一定的参... 目录什么是 AdaBoost?使用 AdaBoost 的步骤安装必要的库步骤一:数据准备步骤二:模型

使用Pandas进行均值填充的实现

《使用Pandas进行均值填充的实现》缺失数据(NaN值)是一个常见的问题,我们可以通过多种方法来处理缺失数据,其中一种常用的方法是均值填充,本文主要介绍了使用Pandas进行均值填充的实现,感兴趣的... 目录什么是均值填充?为什么选择均值填充?均值填充的步骤实际代码示例总结在数据分析和处理过程中,缺失数

Java对象转换的实现方式汇总

《Java对象转换的实现方式汇总》:本文主要介绍Java对象转换的多种实现方式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录Java对象转换的多种实现方式1. 手动映射(Manual Mapping)2. Builder模式3. 工具类辅助映

C语言中位操作的实际应用举例

《C语言中位操作的实际应用举例》:本文主要介绍C语言中位操作的实际应用,总结了位操作的使用场景,并指出了需要注意的问题,如可读性、平台依赖性和溢出风险,文中通过代码介绍的非常详细,需要的朋友可以参... 目录1. 嵌入式系统与硬件寄存器操作2. 网络协议解析3. 图像处理与颜色编码4. 高效处理布尔标志集合

Go语言开发实现查询IP信息的MCP服务器

《Go语言开发实现查询IP信息的MCP服务器》随着MCP的快速普及和广泛应用,MCP服务器也层出不穷,本文将详细介绍如何在Go语言中使用go-mcp库来开发一个查询IP信息的MCP... 目录前言mcp-ip-geo 服务器目录结构说明查询 IP 信息功能实现工具实现工具管理查询单个 IP 信息工具的实现服

SpringBoot基于配置实现短信服务策略的动态切换

《SpringBoot基于配置实现短信服务策略的动态切换》这篇文章主要为大家详细介绍了SpringBoot在接入多个短信服务商(如阿里云、腾讯云、华为云)后,如何根据配置或环境切换使用不同的服务商,需... 目录目标功能示例配置(application.yml)配置类绑定短信发送策略接口示例:阿里云 & 腾

python实现svg图片转换为png和gif

《python实现svg图片转换为png和gif》这篇文章主要为大家详细介绍了python如何实现将svg图片格式转换为png和gif,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录python实现svg图片转换为png和gifpython实现图片格式之间的相互转换延展:基于Py

Python利用ElementTree实现快速解析XML文件

《Python利用ElementTree实现快速解析XML文件》ElementTree是Python标准库的一部分,而且是Python标准库中用于解析和操作XML数据的模块,下面小编就来和大家详细讲讲... 目录一、XML文件解析到底有多重要二、ElementTree快速入门1. 加载XML的两种方式2.