标签:qpow ch frac P3306 ll 生成器 return SDOI2013 xn
思路:\(BSGS\)
提交:\(1\)次
题解:
原式可以化为\[x_{i+1}+\frac{b}{a-1}=a(x_{i}+\frac{b}{a-1})\mod p\]
这不是等比数列吗?
\[x_{n}+\frac{b}{a-1}=a^{n-1}\cdot (x_{1}+\frac{b}{a-1})\mod p\]
所以有
\[a^{n-1}=(x_{1}+\frac{b}{a-1})^{-1}\cdot (x_{n}+\frac{b}{a-1})\mod p\]
这样我们可以\(BSGS\)
注意特判\(a=1,0\)的情况
#include<cstdio>
#include<iostream>
#include<unordered_map>
#include<cmath>
#define ll long long
#define R register int
using namespace std;
namespace Luitaryi {
template<class I> inline I g(I& x) { x=0; register I f=1;
register char ch; while(!isdigit(ch=getchar())) f=ch=='-'?-1:f;
do x=x*10+(ch^48); while(isdigit(ch=getchar())); return x*=f;
}
int T,p,a,b,x1,xn;
inline ll qpow(ll a,ll b) { register ll ret=1;
for(;b;b>>=1,(a*=a)%=p) if(b&1) (ret*=a)%=p; return ret;
}
inline int BSGS() {
unordered_map<int,int> hsh; hsh.clear();
R t=sqrt(p)+1; R c=(xn+1ll*b*qpow(a-1,p-2))%p*qpow((x1+1ll*b*qpow(a-1,p-2))%p,p-2)%p;
for(R i=1;i<=t;++i) {
R vl=1ll*c*qpow(a,i)%p;
hsh[vl]=i;
} a=qpow(a,t);
if(a==0) return c==0?1:-2;
for(R i=1;i<=t;++i) {
R vl=qpow(a,i);
if(hsh.count(vl)&&i*t-hsh[vl]>=0) return i*t-hsh[vl];
} return -2;
}
inline void main() {
g(T); while(T--) {
g(p),g(a),g(b),g(x1),g(xn);
if(x1==xn) {puts("1"); continue;}
if(a==0) {if(xn==b) puts("2"); else puts("-1"); continue;}
if(a==1&&b==0) {puts("-1"); continue;}
if(a==1) {
printf("%d\n",1ll*((xn-x1)%p+p)%p*qpow(b,p-2)%p+1);
continue;
} printf("%d\n",BSGS()+1);
}
}
} signed main() {Luitaryi::main(); return 0;}
2019.08.24
76
标签:qpow,ch,frac,P3306,ll,生成器,return,SDOI2013,xn 来源: https://www.cnblogs.com/Jackpei/p/11404629.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。