博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
日渐成熟的Pollard_Rho板子
阅读量:5250 次
发布时间:2019-06-14

本文共 1419 字,大约阅读时间需要 4 分钟。

namespace PR{    map
book; ll mul(ll a,ll b,ll MOD) {return (a*b-(ll)((long double)a/MOD*b)*MOD+MOD)%MOD;} ll pp[]={2,3,5,7,11,13,17}; ll ksm(ll a,ll b,ll MOD) { ll ans=1; for(;b;b>>=1,a=mul(a,a,MOD)) if(b&1LL)ans=mul(ans,a,MOD); return ans; } bool ip(ll x) { if(x==1)return 0; for(int i=0;i<7;i++) if(x==pp[i])return 1; else if(x%pp[i]==0)return 0; ll tmp=x-1; int cnt=0; while((tmp&1)^1)++cnt,tmp>>=1; for(int i=0;i<7;i++){ ll val=ksm(pp[i],tmp,x),nxt; for(int op=1;op<=cnt;op++){ nxt=mul(val,val,x); if(nxt==1&&val!=1&&val!=x-1)return 0; val=nxt; } if(val!=1)return 0; } for(int i=0;i<3;i++){ ll val=ksm(rand()%(x-1)+1,tmp,x),nxt; for(int op=1;op<=cnt;op++){ nxt=mul(val,val,x); if(nxt==1&&val!=1&&val!=x-1)return 0; val=nxt; } if(val!=1)return 0; } return 1; } ll Abs(ll x) {return x<0?-x:x;} ll Found(ll x,ll c) { ll a=rand()%(x-1)+1,b=a,res; for(int i=2,j=2;;i++){ a=(mul(a,a,x)+c)%x; if(a==b)return -1; if((res=__gcd(x,Abs(a-b)))>1&&res

转载于:https://www.cnblogs.com/CQNKinclude/articles/10853390.html

你可能感兴趣的文章
JSON跨域解决方案收集
查看>>
图的深度优先遍历
查看>>
C# 之 提高WebService性能大数据量网络传输处理
查看>>
[bzoj1004] [HNOI2008] Cards
查看>>
原生HttpClient详细使用示例
查看>>
几道面试题
查看>>
Factory Design Pattern
查看>>
python中贪婪与非贪婪
查看>>
guava API整理
查看>>
无锁编程笔记
查看>>
jquery mobile
查看>>
如何在vue单页应用中使用百度地图
查看>>
Springboot使用步骤
查看>>
Spring属性注入
查看>>
Springboot-配置文件
查看>>
Springboot-日志框架
查看>>
P1192-台阶问题
查看>>
一、使用pip安装Python包
查看>>
spring与quartz整合
查看>>
Kattis之旅——Eight Queens
查看>>