poj 2417

http://poj.org/problem?id=2417

求b^l = n%p；

已知b,n,p求l的一个算法

参考：http://

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <map>
using namespace std;
#define ll long long

ll b,n,p,l,x;
bool flag;


map<ll,ll>s;


long long  qsm(long long x)  //快速幂
{
  long long sum=1; long long aa=b;
  while (x>0)
   {
     if (x&1)
      sum=(sum*aa)%p;
     x=x>>1;
     aa=(aa*aa)%p;
   }
  return sum;
}

void bs()
{
     long long ans;
     for (ll i=0;i<=x;i++)
      {
         if (i==0)
          {
            ans=n%p; s[ans]=i; continue;
          }
         ans=(ans*b)%p;
         s[ans]=i;
      }
}

void gs()
{
     long long t=qsm(x), ans=1;
     for (ll i=1;i<=x;i++)
      {
        ans=(ans*t)%p;
        if (s[ans])
         {
            int t=i*x-s[ans];
            printf("%d
",(t%p+p)%p);
            flag=true;
            break;
         }
      }
}

int main()
{
    //freopen("in.txt","r",stdin);
    while(~scanf("%lld%lld%lld",&p,&b,&n))
    {

        s.clear();
        flag  = false;
        x = sqrt(p)+1;
        bs();
        gs();
        if(!flag)
            printf("no solution
");
    }
    return 0;
}

blog.csdn.net/clover_hxy/article/details/50683832