解一类A^x=B(mod C)(C是质数)的方程 （BSGS）

bsgs算法
主要用来解决 A^x=B(mod C)(C是质数)，都是整数，已知A、B、C求x。（poj 2417 Discrete Logging）
具体步骤如下：
先把x=i*m-j，其中m=ceil(sqrt(C))，（ceil是向上取整）。
这样原式就变为A^(i*m-j)=B(mod C)，
再变为A^j×B=A^(m*i) (mod C)。
枚举j(范围0-m),将A^j×B存入hash表
枚举i(范围1-m),从hash表中寻找第一个满足A^j×B=A^(m*i) (mod C)。
此时x=i*m-j即为所求。
在网上题解大多用的是x=i*m+j，其实思路是一样的，只是会牵扯的求逆元，使x=i*m-j就可以轻松避免这个问题了。
那么肯定有人会有疑问为何只计算到m=ceil(sqrt(C))就可以确定答案呢？
x=i*m-j 也就是x 的最大值不会超过p,那超过p的怎么办 ？
有一个公式 a^(k mod p-1)=a^k (mod p) 这个公式的推导需要用到费马小定理
k mod p-1可以看做 k-m（p-1） ,原式可化成 a^k/(a^(p-1))^m=a^k (mod p)
根据费马小定理 a^(p-1)=1 (mod p) 其中p为质数 ，a,p 互质，可得a^k/1^m=a^k (mod p) a^k=a^k (mod p) 得证。
可以简单理解为分块之后的暴力枚举。
下面第一份代码为 i*m-j， 第二份代码为 i*m+j。 存储时Hash，map都行。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<map> 
#define LL long long    
using namespace std; 

LL a,b,c,m,f[10000000];

map<LL,int> mp;    

LL qsm(LL x){
  //快速幂    
    LL sum=1, aa=a;     
    while(x>0){