首页 > 代码库 > poj 2417
poj 2417
http://poj.org/problem?id=2417
求b^l = n%p;
已知b,n,p求l的一个算法
参考:http://
1 #include <stdio.h> 2 #include <string.h> 3 #include <math.h> 4 #include <map> 5 using namespace std; 6 #define ll long long 7 8 ll b,n,p,l,x; 9 bool flag; 10 11 12 map<ll,ll>s; 13 14 15 long long qsm(long long x) //快速幂 16 { 17 long long sum=1; long long aa=b; 18 while (x>0) 19 { 20 if (x&1) 21 sum=(sum*aa)%p; 22 x=x>>1; 23 aa=(aa*aa)%p; 24 } 25 return sum; 26 } 27 28 void bs() 29 { 30 long long ans; 31 for (ll i=0;i<=x;i++) 32 { 33 if (i==0) 34 { 35 ans=n%p; s[ans]=i; continue; 36 } 37 ans=(ans*b)%p; 38 s[ans]=i; 39 } 40 } 41 42 void gs() 43 { 44 long long t=qsm(x), ans=1; 45 for (ll i=1;i<=x;i++) 46 { 47 ans=(ans*t)%p; 48 if (s[ans]) 49 { 50 int t=i*x-s[ans]; 51 printf("%d\n",(t%p+p)%p); 52 flag=true; 53 break; 54 } 55 } 56 } 57 58 int main() 59 { 60 //freopen("in.txt","r",stdin); 61 while(~scanf("%lld%lld%lld",&p,&b,&n)) 62 { 63 64 s.clear(); 65 flag = false; 66 x = sqrt(p)+1; 67 bs(); 68 gs(); 69 if(!flag) 70 printf("no solution\n"); 71 } 72 return 0; 73 }
blog.csdn.net/clover_hxy/article/details/50683832
poj 2417
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。