題面

bsgs問(wèn)題。因?yàn)閜可能不為質(zhì)數(shù)，所以我們將原先解題的式子變形
每次除以p與a的最大公約數(shù)，直到最大公約數(shù)為1或b不能整除為止

代碼

#include<iostream> #include<algorithm> #include<cstdio> #include<cstring> #include<cmath> #include<map> #define LL long longusing namespace std;LL a,b,m,p,now,ans; bool flag; map<LL,int> mp;inline LL fast_pow(LL a,LL b){LL ret=1;LL aa=a;for(;b;b>>=1){if(b&1) ret=ret*aa%p;aa=aa*aa%p;}return ret; }int main(){while(~scanf("%lld%lld%lld",&a,&p,&b)){if(a==0 && b==0 && p==0) break;if(a%p==0){puts("No Solution");continue;}if(b==1) {puts("0");continue;}flag=false;a%=p;b%=p;LL t=1,cnt=0;for(register int i=__gcd(a,p);i!=1;i=__gcd(a,p)){if(b%i){puts("No Solution");flag=1;break;}p/=i;t=t*a/i%p;b/=i;cnt++;if(b==t) {printf("%lld\n",cnt);flag=1;break;}}if(flag) continue;mp.clear();now=b;mp[now]=0;m=ceil(sqrt(p));for(register int i=1;i<=m;i++){now=now*a%p;mp[now]=i;}now=t;LL k=fast_pow(a,m);for(register int i=1;i<=m;i++){now=now*k%p;if(mp[now]){flag=true;ans=i*m-mp[now]+cnt;printf("%lld\n",ans);break;}} if(!flag) puts("No Solution"); }return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/sdfzsyq/p/9677003.html

總結(jié)

以上是生活随笔為你收集整理的LUOGU P4195 Spoj3105 Mod的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。