POJ2417 Discrete Logging | A,C互质的bsgs算法
题目:
给出A,B,C
求最小的x使得Ax=B (mod C)
题解:
bsgs算法的模板题
bsgs 全称:Baby-step giant-step
把这种问题的规模降低到了sqrt(n)级别
首先B的种类数不超过C种,结合鸽巢原理,所以Ax具有的周期性显然不超过C
所以一般的枚举算法可以O(C)解决这个问题
但是可以考虑把长度为C的区间分为k块,每块长度为b
显然x满足x=bi-p的形式(1<=i<=k,0<=p<b),所以Ax=B (mod C)移项之后得到Abi=Ap*B (mod C)
那么这个时候可以预处理出来Ap的所有值(可以用hash表维护)
//注意!hash表在插入之前要先找有没有这个值,如果有的话直接把改了就好
处理出Ab的值,枚举i,就可以得到答案
一般来说令k=b=sqrt(C)时间复杂度最优
1 #include<cstdio> 2 #include<algorithm> 3 #include<cstring> 4 #include<map> 5 #include<cmath> 6 #define MOD 1000007 7 #define EDGE 500000 8 typedef long long ll; 9 using namespace std; 10 ll c,a,b,ok,m,tmp,t,ans,head[MOD],ecnt; 11 struct adj 12 { 13 ll nxt,w,sum; 14 }e[EDGE]; 15 void add(ll x,ll sum) 16 { 17 ll org=x; 18 e[++ecnt].w=x; 19 x%=MOD; 20 for (int i=head[x];i;i=e[i].nxt) 21 if (e[i].w==org) 22 { 23 e[i].sum=sum; 24 return ; 25 } 26 e[ecnt].sum=sum; 27 e[ecnt].nxt=head[x]; 28 head[x]=ecnt; 29 } 30 ll qow(ll x,ll y,ll P) 31 { 32 if (y==0) return 1; 33 if (y&1) return x*qow(x*x%P,y>>1,P)%P; 34 return qow(x*x%P,y>>1,P)%P; 35 } 36 ll find(ll x) 37 { 38 ll org=x; 39 x%=MOD; 40 for (int i=head[x];i;i=e[i].nxt) 41 { 42 if (e[i].w==org) 43 return e[i].sum; 44 } 45 return -1; 46 } 47 int main() 48 { 49 while (scanf("%lld%lld%lld",&c,&a,&b)!=EOF) 50 { 51 memset(head,0,sizeof(head)); 52 ecnt=0; 53 ok=0; 54 if (a%c==0) 55 { 56 puts("no solution"); 57 continue; 58 } 59 m=ceil(sqrt(c*1.0)); 60 tmp=b%c,add(tmp,0); 61 if (b==1) 62 { 63 printf("0\n"); 64 continue; 65 } 66 for (int i=1;i<m;i++) 67 { 68 tmp=tmp*a%c; 69 add(tmp,i); 70 } 71 ll base=qow(a,m,c),tmp=1; 72 for (int i=1;i<=m;i++) 73 { 74 tmp=tmp*base%c; 75 ans=find(tmp); 76 if (ans!=-1) 77 { 78 printf("%lld\n",i*m-ans); 79 ok=1; 80 break; 81 } 82 } 83 if (!ok) 84 puts("no solution"); 85 } 86 return 0; 87 }
转载于:https://www.cnblogs.com/mrsheep/p/7920967.html
POJ2417 Discrete Logging | A,C互质的bsgs算法相关推荐
- 数列互质(莫队算法)
数列互质 题目描述 给出一个长度为 n 的数列 { a[1] , a[2] , a[3] , - , a[n] },以及 m 组询问 ( l[i] , r[i] , k[i]). 求数列下标区间在 [ ...
- python判断两个数是否互质_《算法》第一章——判断两个整数是否互质
判断两个整数互质的方法 概念:公约数只有1的两个数叫做互质数.根据互质数的概念可以对一组数是否互质进行判断.如:9和11的公约数只有1,则它们是互质数. 求商判断法:用大数除以小数,如果除得的余数与其 ...
- 容斥原理 —— 求1~n有多少个数与k互质(二进制算法详细解释模板)
这里有一道经典的例题,可以看一下:点击打开链接 这里的n可能要大于k的,所以不能用欧拉函数去做. 我们首先把k分解质因数,储存到p数组中,num表示质因子的数量. void pr(int k) //求 ...
- Discrete Logging ZOJ - 1898 (模板题大小步算法)
就是求Ax三B(mod C)当C为素数时 #include<cstdio> #include<cstring> #include<cmath> #include&l ...
- c语言 判断互质,[经典算法] 找出某数以内与其互质的数
起因:前几天做了一个这样的题目,感觉自己变成了一个智障,来写写算法: 相关定义: ①质数定义:没有除 '1' 以外的因子的数就是质数.如:1.3.5.7.11等等. ②互质定义:两个数之间除 '1' ...
- 【noi 2.5_7834】分成互质组(dfs)
有2种dfs的方法: 1.存下每个组的各个数和其质因数,每次对于新的一个数,与各组比对是否互质,再添加或不添加入该组. 2.不存质因数了,直接用gcd,更加快.P.S.然而我不知道为什么RE,若有好心 ...
- #6164. 「美团 CodeM 初赛 Round A」数列互质-莫队
#6164. 「美团 CodeM 初赛 Round A」数列互质 思路 : 对这个题来言,莫队可以 n*根号n 离线处理出各个数出现个的次数 ,同时可以得到每个次数出现的次数 , 但是还要处理有多少 ...
- 中国剩余定理matlab非互质,中国剩余定理模板(互质版和非互质版)
互质版: #include #include #include using namespace std; typedef __int64 int64; int64 a[15],b[15]; int64 ...
- bzoj 4921: [Lydsy六月月赛]互质序列
4921: [Lydsy六月月赛]互质序列 Time Limit: 1 Sec Memory Limit: 256 MB Submit: 188 Solved: 110 [Submit][Stat ...
最新文章
- 哪个更值得入手?荣耀8与小米6对比评测
- java面试题2 牛客:定义类中成员变量时不可能用到的修饰是
- vs2005 生成安装包的时候提示“不可恢复的生成错误”
- 在SQL Server 2005中实现异步触发器架构
- 2.4.PHP7.1 狐教程-【PHP常量】
- Transformer论文阅读(一):CoTr: Efficiently Bridging CNN and Transformer for 3D Medical Image Segmentation
- mac 设置网页字体
- 思迈特软件完成C轮融资,让BI真正实现“普惠化”
- VMware虚拟机不能识别U盘
- uni-app云打包失败
- 腾讯云---语音识别--一句话识别
- CATT eCATT
- Rust vs Pandas 数据处理
- 100+文档格式预览,私有化部署,10分钟搭建专属IM
- 【2022年的STM32】 04-GPIO特性、使用及与NXP GPIO比较
- Wordpress后台寻找新的主题,缩略图无法显示(问题解决方案)
- 示波器的触发设置详谈
- 计算机主机如何启动,教你电脑关机后自动开机怎么办
- 全国计算机一级上机考试试题,全国计算机等级考试一级上机WORD题(六)
- python弹球游戏实验报告_弹珠小游戏课程设计实验报告