BSGS ! x
一.引入:
若存在一个式子a^b ≡ c (mod p) (p ≡ 1000000007,且0<a,b,c<p)
已知a,b,求c.
这不就是快速幂嘛!
已知a,c,求b.
这就是我们需要研究的问题!用到了BSGS!
题目链接:poj 2417 bsgs
二.概念
BSGS:
又名大步小步算法.具体的我也不清楚啦~
那么发明来做什么事情呢?
如上所述:
就是用来求解a^x ≡ b (mod p)这样的式子
PS:已知a,b,p
求最小x
三.做法
首先,我们将x用i*m-j来表示,其中我们的m=seil(sqrt(p)),(seil为向上取整)
然后我们用i*m-j代替掉x,所以原式就变成了这个样子
a^(i * m-j) ≡ b (mod p)
又因为a^(i * m-j) => a^(i * m)/a^j
所以原式又变了模样:
a^(i * m)/a^j ≡ b (mod p)
又因为在"≡"(同余号)的两侧同时乘以一个相同的数依旧是成立的,如果不明白可以手动模拟一下,
给出个栗子~
3 ≡ 10 (mod 7),当两边同乘以5时式子变成: 15 ≡ 50 (mod 7)
因为15%7==1&&50%7==1,所以原式依旧成立
所以将式子两边同时乘以a^j
那么式子就又发生了改变:
a^(i * m) ≡ b * a^j (mod p)
= (a^m)^i ≡ b * a^j (mod p)
所以问题就简单多了!
因为a^m是常数,b是常数,p也是常数,所以只有 i 跟 j 是未知的,暴力枚举!
但是有一点是我们不能够忘记的:
m的取值范围:1 —— p(p是当b==0时),所以
首先将 j 从1 到 p-1 进行枚举一遍,求出b * a^j 的值丢到hash(大佬是用map做的)里面咯(讲真hash不会用...怪我咯?)
思路是这样的:(其实就是hash啦~)
(b*a^j)%p(如果p太小换成另外一个比较大的质数)作为下标,里面存着当前计算出来的数,以及当前的 j (可以用vector~)
然后,枚举 j ,我们大可以设a^m==c
将 i 从 1 到p-1枚举一遍,如果求出的数在hash里面找到了,则说明当前的数就是 i (能够成立的) 最小值,跳出循环,用当前的i,以及记录下来的 j 计算对应的 i 求出 x
则x就是最小的解
C++代码实现:
(map方法:)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<cmath> using namespace std;
typedef long long LL; LL a,b,c;
map<LL,LL> mp; LL qsm(LL m)
{
LL n = a; if(m == ) return ; LL t = qsm(m/); t = 1LL*t*t%c;
if(m&) t = 1LL*t*n%c; return t;
} int main()
{
//a^im=b*a^j(mod c)
while (scanf("%lld%lld%lld",&c,&a,&b)!=EOF)
{
mp.clear(); //清空
if (a%c==) //判断a,c 是否互质,因为c 是质数,所以直接判断是否整除即可
{
printf("no solution\n");
continue;
}
LL m=ceil(sqrt(c));
LL ans;
for (LL j=; j<=m; j++)
{
if (j==)
{
//当j=0时,a^j=1, b*a^j=b
ans=b%c;
mp[ans]=j;
continue;
}
ans=(ans*a)%c;
//括号里的ans指a^(j-1)*b,(a^(j-1)*b)*a=(a^j)*b
mp[ans]=j;//在((a^j)*b)%c的位置记录下j
}
LL t=qsm(m);//t=a^m
ans=;
bool p=false;
for (LL i=; i<=m; i++)
{
ans=(ans*t)%c;//括号里的ans指的是((a^m)^(i-1))*(a^m)=(a^m)^i=a^(im)
if (mp[ans])
{
LL t=i*m-mp[ans];//t=x,因为我们设的x=i*m-j
printf("%lld\n",t);
p=true;
break;
}
}
if (!p)
printf("no solution\n");
}
}
End.
BSGS ! x的更多相关文章
- 【POJ 3243】Clever Y 拓展BSGS
调了一周,我真制杖,,, 各种初始化没有设为1,,,我当时到底在想什么??? 拓展BSGS,这是zky学长讲课的课件截屏: 是不是简单易懂.PS:聪哥说“拓展BSGS是偏题,省选不会考,信我没错”,那 ...
- 【codevs 1565】【SDOI 2011】计算器 快速幂+拓展欧几里得+BSGS算法
BSGS算法是meet in the middle思想的一种应用,参考Yveh的博客我学会了BSGS的模版和hash表模板,,, 现在才会hash是不是太弱了,,, #include<cmath ...
- 【BZOJ-3122】随机数生成器 BSGS
3122: [Sdoi2013]随机数生成器 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1362 Solved: 531[Submit][Sta ...
- BSGS[bzoj2242][bzoj3122]
数论题. 操作一:直接快速幂就好了. 操作二:我用了exgcd,shy和lyz都喜欢欧拉函数...QAQ最后这块还写错了. 对于ax+by=gcd(a,b)的形式,我们可以把他们变成y'x+p'y=1 ...
- 【BZOJ2242】【SDoi2011】计算器 快速幂+EXGCD+BSGS
Description 你被要求设计一个计算器完成以下三项任务: 1.给定y,z,p,计算Y^Z Mod P 的值: 2.给定y,z,p,计算满足xy≡ Z ( mod P )的最小非负整数: 3.给 ...
- BZOJ-2242 计算器 快速幂+拓展欧几里得+BSGS(数论三合一)
污污污污 2242: [SDOI2011]计算器 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 2312 Solved: 917 [Submit][S ...
- bzoj 1467: Pku3243 clever Y 扩展BSGS
1467: Pku3243 clever Y Time Limit: 4 Sec Memory Limit: 64 MB[Submit][Status][Discuss] Description 小 ...
- BSGS模版 a^x=b ( mod c)
kuangbin的BSGS: c为素数: #define MOD 76543 int hs[MOD],head[MOD],next[MOD],id[MOD],top; void insert(int ...
- bzoj 2242: [SDOI2011]计算器 BSGS+快速幂+扩展欧几里德
2242: [SDOI2011]计算器 Time Limit: 10 Sec Memory Limit: 512 MB[Submit][Status][Discuss] Description 你被 ...
- 【poj3358】消因子+BSGS 或 消因子+欧拉定理 两种方法
题意:给你一个分数,求它在二进制下的循环节的长度,还有第一个循环节从哪一位开始. For example, x = 1/10 = 0.0001100110011(00110011)w and 0001 ...
随机推荐
- Message高级特性
一.延迟和定时消息投递 1 xml中配置schedulerSupport属性为true 2 每小时都会发生消息被投递10次,延迟1秒开始,每次间隔1秒: TextMessage message = s ...
- Linux下如何查看tomcat是否启动、查看tomcat启动日志
在Linux系统下,重启Tomcat使用命令的操作! 1.首先,进入Tomcat下的bin目录 cd /usr/local/tomcat/bin 使用Tomcat关闭命令 ./shutdown.sh ...
- 改善C#程序的方法
写在开头: http://www.cnblogs.com/luminji 157个建议_勘误表 一:属性 属性和方法一样.也可以是virtual和abstract. 条款2:运行时常量(read ...
- python打印表格式数据-星号或注释
python打印表格式数据,留出正确的空格,格式化打出 代码如下: def printPicnic(itemsDict,leftWidth,rightWidth): print('PICNIC ITE ...
- 异常信息:javax.net.ssl.SSLHandshakeException: sun.security.validator.ValidatorException: PKIX path building failed
上周五遇到一个问题,工程本地编译运行正常,打包本地tomcat运行也正常.部署到测试环境报错: 2017-05-05 09:38:11.645 ERROR [HttpPoolClientsUtil.j ...
- Javascript--HTML DOM基础知识
1.HTML DOM是什么,以及它的作用: w3c对DOM有一系列的解释和定义,用自己理解的话来说就是:HTML DOM是html的标准对象模型,可以使JavaScript去操作(获取,修改,删除,添 ...
- element之 el-scrollbar组件滚动条的使用
在使用vue + element-ui 搭建后台管理页面的时候,做了一个头部.侧栏.面包屑固定的布局,导航栏和主要内容区域当内容超出时自动滚动.
- C# 中 Linq 操作 DataTable
方法一:更简洁 Console.WriteLine(dt.Rows.OfType<DataRow>().First(x => x.Field<string>(" ...
- python点击短信验证码
代码如下 : import requestsimport time# 手机号码tel=1381380000# 请求地址url="http://192.168.100.101:8080/api ...
- 超详细的Java面试题总结之JavaWeb基础知识总结
ervlet总结: 在Java Web程序中,Servlet主要负责接收用户请求HttpServletRequest,在doGet(),doPost()中做相应的处理,并将回应HttpServletR ...