【模板】exBSGS/Spoj3105 Mod
【模板】exBSGS/Spoj3105 Mod
题目描述
已知数\(a,p,b\),求满足\(a^x\equiv b \pmod p\)的最小自然数\(x\)。
输入输出格式
输入格式:
每个测试文件中最多包含\(100\)组测试数据。
每组数据中,每行包含\(3\)个正整数\(a,p,b\)。
当\(a=p=b=0\)时,表示测试数据读入完全。
输出格式:
对于每组数据,输出一行。
如果无解,输出No Solution
(不含引号),否则输出最小自然数解。
BSGS
若\(A \perp p\),那么\(\{A^x,x\le \varphi(p)\}\)遍历的剩余系\(\{A^{kx},x\le \varphi(p)\}\)一定也遍历,于是考虑枚举答案
\]
采用分块的思想,设\(t=\sqrt p,x=kt-b\),式子就变成了
\]
\]
我们枚举\(x=0 \sim t\),然后把得到的\(A^xB\)插到\(\tt{Hash}\)表中去。
然后枚举\((A^t)^k\)的\(k\),查询\(\tt{Hash}\)表中有没有\(A^{kt}\)
exBSGS
如果\(p\)不是质数,存在无解的判定\((\gcd(A,p)\nmid B)\)且\(B\not=1\)(\(B=1\)特判\(x=0\))
然后考虑操作一波式子
\]
把\(d\)除掉
\]
设\(C=\frac{A}{d},B'=\frac{B}{d},p'=\frac{p}{d}\)
原方程变为
\]
然后重复是否无解的判断并向下递归,直到\(A\perp p\)或者无解
然后\(BSGS\)即可,而常数\(C\)并不影响我们进行\(BSGS\)
复杂度?显然递归的深度是\(\log\)的,带上BSGS的就可以了。
Code:
#include <cstdio>
#include <cmath>
#include <unordered_map>
std::unordered_map <int,int> Hash;
int gcd(int a,int b){return b?gcd(b,a%b):a;}
#define mul(a,b,p) (1ll*(a)*(b)%p)
int exbsgs(int A,int B,int p)
{
if(B==1) return 0;
int ct=0,d,k=1;
while((d=gcd(A,p))^1)
{
if(B%d) return -1;
B/=d,p/=d,++ct;
k=mul(k,A/d,p);
if(k==B) return ct;
}
int t=sqrt(p)+1,kt=1;
Hash.clear();
for(int i=0;i<t;i++)
{
Hash[mul(kt,B,p)]=i;
kt=mul(kt,A,p);
}
k=mul(k,kt,p);
for(int i=1;i<=t;i++)
{
if(Hash.find(k)!=Hash.end()) return i*t-Hash[k]+ct;
k=mul(k,kt,p);
}
return -1;
}
int main()
{
int a,p,b;
scanf("%d%d%d",&a,&p,&b);
while(a&&p&&b)
{
int ans=exbsgs(a,b,p);
if(~ans) printf("%d\n",ans);
else puts("No Solution");
scanf("%d%d%d",&a,&p,&b);
}
return 0;
}
2018.12.19
【模板】exBSGS/Spoj3105 Mod的更多相关文章
- P4195 【模板】exBSGS/Spoj3105 Mod
传送门 首先要懂得 $BSGS$,$BSGS$ 可以求出关于 $Y$ 的方程 $X^Y \equiv Z (mod\ mo)$ 的最小解,其中 $gcd(X,Z)=1$ $exBSGS$ 算是 $BS ...
- 【BZOJ1467/2480】Pku3243 clever Y/Spoj3105 Mod EXBSGS
[BZOJ1467/2480]Pku3243 clever Y/Spoj3105 Mod Description 已知数a,p,b,求满足a^x≡b(mod p)的最小自然数x. Input ...
- 【bzoj2480】Spoj3105 Mod
2480: Spoj3105 Mod Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 557 Solved: 210[Submit][Status][ ...
- BSGS 扩展大步小步法解决离散对数问题 (BZOJ 3239: Discrete Logging// 2480: Spoj3105 Mod)
我先转为敬? orz% miskcoo 贴板子 BZOJ 3239: Discrete Logging//2480: Spoj3105 Mod(两道题输入不同,我这里只贴了3239的代码) CODE ...
- spoj3105 MOD - Power Modulo Inverted(exbsgs)
传送门 关于exbsgs是个什么东东可以去看看yyb大佬的博客->这里 //minamoto #include<iostream> #include<cstdio> #i ...
- BZOJ2480 Spoj3105 Mod 数论 扩展BSGS
原文链接https://www.cnblogs.com/zhouzhendong/p/BZOJ2480.html 题目传送门 - BZOJ2480 题意 已知数 $a,p,b$ ,求满足 $a^x≡b ...
- 模板BSGS(SDOI2011计算器) 模板EXBSGS
BSGS和EXBSGS是OI中用于解决A^xΞB(mod C)的常用算法. 1.BSGS BSGS用于A,C互质的情况. 令m=sqrt(C),此时x可表示为i*m+j. 式中i和j都<=sqr ...
- [luogu4195 Spoj3105] Mod (大步小步)
传送门 题目描述 已知数a,p,b,求满足a^x≡b(mod p)的最小自然数x. 输入输出格式 输入格式: 每个测试文件中最多包含100组测试数据. 每组数据中,每行包含3个正整数a,p,b. 当a ...
- BZOJ2480 Spoj3105 Mod
乍一看题面:$$a^x \equiv b \ (mod \ m)$$ 是一道BSGS,但是很可惜$m$不是质数,而且$(m, a) \not= 1$,这个叫扩展BSGS[额...... 于是我们需要通 ...
随机推荐
- java事务 深入Java事务的原理与应用
一.什么是JAVA事务 通常的观念认为,事务仅与数据库相关. 事务必须服从ISO/IEC所制定的ACID原则.ACID是原子性(atomicity).一致性(consistency).隔离性 (iso ...
- python3实现合并两个有序数组
很早就听同学和师兄经常说刷题很重要,然而编程能力一直都很渣的我最近才开始从leetcode的初级算法开始.今天遇到的这道题虽然很简单,因为是头一次用自己的方法速度还不错,特此记录一下,还大神们请不要嘲 ...
- 关于linux-centos7 安装完成git后npm突然无法使用问题处理
报错: 解决方法: 查看一下nodejs是否安装,如果没有安装的话安装完成就能解决了
- 6.openldap客户端安装
作者:yaoyao 1.账号登录系统流程讲解 当在客户端输入账号登录系统时.系统根据/etc/nsswitch.conf配置文件获取账号查找顺序,然后在根据pam配置文件调用相关模块,对账号/etc/ ...
- pssh命令详解
基础命令学习目录首页 原文链接:https://www.cnblogs.com/kevingrace/p/6378719.html pssh提供OpenSSH和相关工具的并行版本.包括pssh,psc ...
- Python字典 (dict)
作者博文地址:http://www.cnblogs.com/spiritman/ 字典是Python语言中唯一的映射类型.字典对象是可变的,它是一个容器类型,支持异构.任意嵌套. 创建字典 语法:{k ...
- (转)一个简单的rest_framework demo
转发:https://www.cnblogs.com/fu-yong/p/9100559.html models.py 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 from ...
- Vue实现双向绑定的原理以及响应式数据
一.vue中的响应式属性 Vue中的数据实现响应式绑定 1.对象实现响应式: 是在初始化的时候利用definePrototype的定义set和get过滤器,在进行组件模板编译时实现water的监听搜集 ...
- Daily Scrum (2015/10/31)
这几天我们组的进度有点慢,剩下这一周的我们必须要加油认真对待. 周末这两天我们是这样安排的: 成员 今日任务 时间 明日任务 符美潇 数据库部分代码的编写 1h 每周小组例会 潘礼鹏 团队博客作业 ...
- 第二阶段Sprint冲刺会议6
进展:将“录制”及“保存”整合到一起,修复出现的Bug,使之能够正常运行.