题目:

给出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算法相关推荐

  1. 数列互质(莫队算法)

    数列互质 题目描述 给出一个长度为 n 的数列 { a[1] , a[2] , a[3] , - , a[n] },以及 m 组询问 ( l[i] , r[i] , k[i]). 求数列下标区间在 [ ...

  2. python判断两个数是否互质_《算法》第一章——判断两个整数是否互质

    判断两个整数互质的方法 概念:公约数只有1的两个数叫做互质数.根据互质数的概念可以对一组数是否互质进行判断.如:9和11的公约数只有1,则它们是互质数. 求商判断法:用大数除以小数,如果除得的余数与其 ...

  3. 容斥原理 —— 求1~n有多少个数与k互质(二进制算法详细解释模板)

    这里有一道经典的例题,可以看一下:点击打开链接 这里的n可能要大于k的,所以不能用欧拉函数去做. 我们首先把k分解质因数,储存到p数组中,num表示质因子的数量. void pr(int k) //求 ...

  4. Discrete Logging ZOJ - 1898 (模板题大小步算法)

    就是求Ax三B(mod C)当C为素数时 #include<cstdio> #include<cstring> #include<cmath> #include&l ...

  5. c语言 判断互质,[经典算法] 找出某数以内与其互质的数

    起因:前几天做了一个这样的题目,感觉自己变成了一个智障,来写写算法: 相关定义: ①质数定义:没有除 '1' 以外的因子的数就是质数.如:1.3.5.7.11等等. ②互质定义:两个数之间除 '1'  ...

  6. 【noi 2.5_7834】分成互质组(dfs)

    有2种dfs的方法: 1.存下每个组的各个数和其质因数,每次对于新的一个数,与各组比对是否互质,再添加或不添加入该组. 2.不存质因数了,直接用gcd,更加快.P.S.然而我不知道为什么RE,若有好心 ...

  7. #6164. 「美团 CodeM 初赛 Round A」数列互质-莫队

    #6164. 「美团 CodeM 初赛 Round A」数列互质 思路 : 对这个题来言,莫队可以 n*根号n 离线处理出各个数出现个的次数 ,同时可以得到每个次数出现的次数 , 但是还要处理有多少 ...

  8. 中国剩余定理matlab非互质,中国剩余定理模板(互质版和非互质版)

    互质版: #include #include #include using namespace std; typedef __int64 int64; int64 a[15],b[15]; int64 ...

  9. bzoj 4921: [Lydsy六月月赛]互质序列

    4921: [Lydsy六月月赛]互质序列 Time Limit: 1 Sec  Memory Limit: 256 MB Submit: 188  Solved: 110 [Submit][Stat ...

最新文章

  1. 哪个更值得入手?荣耀8与小米6对比评测
  2. java面试题2 牛客:定义类中成员变量时不可能用到的修饰是
  3. vs2005 生成安装包的时候提示“不可恢复的生成错误”
  4. 在SQL Server 2005中实现异步触发器架构
  5. 2.4.PHP7.1 狐教程-【PHP常量】
  6. Transformer论文阅读(一):CoTr: Efficiently Bridging CNN and Transformer for 3D Medical Image Segmentation
  7. mac 设置网页字体
  8. 思迈特软件完成C轮融资,让BI真正实现“普惠化”
  9. VMware虚拟机不能识别U盘
  10. uni-app云打包失败
  11. 腾讯云---语音识别--一句话识别
  12. CATT eCATT
  13. Rust vs Pandas 数据处理
  14. 100+文档格式预览,私有化部署,10分钟搭建专属IM
  15. 【2022年的STM32】 04-GPIO特性、使用及与NXP GPIO比较
  16. Wordpress后台寻找新的主题,缩略图无法显示(问题解决方案)
  17. 示波器的触发设置详谈
  18. 计算机主机如何启动,教你电脑关机后自动开机怎么办
  19. 全国计算机一级上机考试试题,全国计算机等级考试一级上机WORD题(六)
  20. python弹球游戏实验报告_弹珠小游戏课程设计实验报告

热门文章

  1. 【CyberSecurityLearning 30】Linux操作系统的用户和组、文件及目录权限
  2. Oracle 根据字符串的长度排序
  3. int指令---汇编学习笔记
  4. @RequestParam和@RequestBody的区别 (结合 Get/Post )
  5. JVM_03 运行时数据区 [ 方法区]
  6. 操作系统(二十七)管程
  7. FFmpeg再学习 -- FFmpeg解码知识
  8. UNIX再学习 -- 信号
  9. Hi3516A开发--视频输入和输出接口
  10. 吴恩达机器学习笔记——第一章