【文文殿下】ExBSGS

无需逆元版本：

#include<cstdio>
#include<cassert>
#include<cmath>
#include<map>
typedef long long ll;
ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;
}
ll qpow(ll a,ll b,ll p) {ll ret = 1;while(b) {if(b&1) {ret=ret*a%p;}a=a*a%p;b>>=1;}return ret;
}
ll exgcd(ll a,ll b,ll &x,ll &y) {if(b==0) {x=1;y=0;return a;}ll d = exgcd(b,a%b,x,y);ll tmp = x;x = y;y=tmp-a/b*y;return d;
}
inline ll bsgs(ll a,ll b,ll p) {std::map<ll,ll> M;if(p==1) return 0;a%=p;b%=p;ll d,cnt=0,q=1;while((d=gcd(a,p))!=1) {if(b==1) return cnt;if(b%d) return -1;b/=d;p/=d;++cnt;q=a/d*q%p;}ll lmt = ceil(sqrt(p));ll tmp = b%p;for(int i = 0;i<lmt;++i,tmp=tmp*a%p) {M[tmp]=i;}tmp = qpow(a,lmt,p);for(int i = 1;i<=lmt+1;++i) {q=q*tmp%p;if(M.count(q)) {return i*lmt-M[q]+cnt;}}return -1;
}
int main() {ll a,b,p;scanf("%lld%lld%lld",&a,&p,&b);while(a||b||p) {ll ans = bsgs(a,b,p);if(ans==-1) puts("No Solution");else printf("%lld\n",ans);scanf("%lld%lld%lld",&a,&p,&b);}return 0;
}

转载于:https://www.cnblogs.com/Syameimaru/p/10221311.html

【文文殿下】ExBSGS相关推荐

【文文殿下】数论一些经典结论
\(x^2+y^2=n\)的整数解的个数,是n的所有素数中,形如\(4n+1\)的素数的指数+1的乘积,如果有形如\(4n+3\)的素数,指数不为偶数,则无解,证明方法:高斯素数. \(d(ij) = ...
【文文殿下】快速傅里叶变换(FFT)学习笔记
多项式定义形如\(A(x)=\sum_{i=0}^{n-1} a_i x^i\)的式子称为多项式. 我们把\(n\)称为该多项式的次数界. 显然,一个\(n-1\)次多项式的次数界为\(n\). ...
【文文殿下】网络流学习笔记
最大流算法 Dinic 割一个网络的割:存在一个边集,删去集合里的边时,S-T不再连通最小割所有割中边权之和最小的最小割-最大流定理最小割等于最大流二分图匹配 \(M\)为边集\(E\)的 ...
【文文殿下】Manache算法-学习笔记
Manache算法定义:是一个判断回文子串的算法,我们结合例题解释: 题目:给定一个长度为 n 的字符串 S,求其最长回文子串一个字符串是回文的,当且仅当反转后的串与原串完全相等分析:对于这个题 ...
【文文殿下】 [USACO08MAR]土地征用题解
题解斜率优化裸题. 有个很玄学的事情,就是我用\(f[i]=min\{f[j-1]+p[j].y*p[i].x\}\) 会很奇怪的Wa . 明明和\(f[i]=min\{f[j]+p[j+1].y* ...
【文文殿下】浅谈KMP算法next数组与循环节的关系
KMP算法 KMP算法是一种字符串匹配算法,他可以在O(n+m)的时间内求出一个模式串在另一个模式串下出现的次数. KMP算法是利用next数组进行自匹配,然后来进行匹配的. Next数组 Next数 ...
【文文殿下】[BZOJ4008] [HNOI2015] 亚瑟王
题解这是一个经典的概率DP模型设\(f_{i,j}\)表示考虑到前\(i\)张牌,有\(j\)轮没打出牌的可能性,那么显然\(f_{0,r} = 1\). 考虑第\(i+1\)张牌,他可能在剩下的 ...
【文文殿下】Win7打开无线热点
下载程序 : https://files.cnblogs.com/files/Syameimaru/wifi.rar 打开config.ini 修改相关信息,然后运行程序. 运行成功后,打开网络和共享 ...
ad中pcb双面板怎么设置_Altium Designer PCB双面板制作打印操作步骤
Altium Designer PCB双面板制作打印操作步骤百度知道: http://jingyan.baidu.com/article/335530da83441c19cb41c3db.html?s ...

【文文殿下】ExBSGS

无需逆元版本：

【文文殿下】ExBSGS相关推荐

最新文章

热门文章