BSGS算法

BSGS算法介绍

ＢＳＧＳ是用来解决形如 y^x ≡ ｚ mod p 的问题的，它能够在有解的情况下求出来最小值。这里前提是gcd(y, p) == 1，否则就需要ｅｘＢＳＧＳ解决。

算法过程

首先我们让 x = i*m-k
　这样上面的式子就变为了　y^i*m/y^k ≡ z mod p, 通过变换得到 y^i*m ≡ z * y^k mod p。其中ｍ等于ceil(sqrt( p))，也就是向上取整根号ｐ，ｋ的范围是0～m-1。
　通过枚举得到所有y^k ｋ∈ [0, m)，并使用ｈａｓｈ表或者ｍａｐ存储当前ｋ值，然后在枚举左边的ｉ，如果在ｈａｓｈ表或者ｍａｐ里面能够找到，就代表找到了最小的解。输出i * m-map[num]就行了。
　
　１，首先说为什么m等于ceil(sqrt§)呢？
　没有为啥，ｍ等于什么都行，但是当他等于ceil(sqrt§)时是最快的。
　
　２，ｉ的取值范围为什么是[1, m]？
　首先需要保证ｘ为正数，所以ｉ一定大于０。然后通过欧拉定理的应用我们知道 a^{t mod φ( p)} = a^t mod p, 这样就保证了ｘ一定是小于ｐ的，因为ｐ为素数的时候欧拉函数值最大为p-1。所以i * m-k <= p，而ｋ的取值范围为ｋ∈ [0, m)，当ｋ = 0的时候，ｉ * m <= p, 又因为我们知道ｍ取值为sqrt( p ), 所以i的取值范围是[1, m]。
　
　３，为什么得到的一定是最小解？
　这个算法叫做大小步算法，在记录所有y^k k∈[0, m)的时候就是小步走，小步走时，如果遇到取模后相同的值就会用当前的ｋ值替换到hash表里面的，这样最后找到的一定是最大的ｋ，因此 i * m - k一定是最小值也就是最小解。

代码实现

// y^x = z mod p
void slove(ll y, ll z, ll p)
{
    map<ll, int> s;
	//特判一下，如果ｙ是ｐ的倍数一定无解
    if (y%p == 0)
    {
        printf("Orz, I cannot find x!\n");
        return;
    }

    y%=p, z%=p;
	//特判
    if (z == 1)
    {
        printf("0\n");
    }
    else
    {	//求出来ｍ的值，加一效果同ceil一样
        int m = sqrt(p)+1;
        ll t = z;
		//求解ｙ^k，并记录下来，ｋ是从[0, m)的，
		//这里为了让最后一个ｔ值也记录下来，就循环到了ｍ
        for (int i=0; i<=m; i++)
        {
            s[t] = i+1; //记录位置加一，一会在加上来，因为判断需要用到值，没办法从０开始
            t = t*y%p;
        }

        ll q = qul(y, m, p);
        t = 1;
        bool flag = false;

        for (int i=1; i<=m; i++)
        {
            t = t*q%p;

            if (s[t])
            {
                printf("%d\n", i*m-s[t]+1);
                flag = true;
                break;
            }
        }

        if (!flag)
                printf("Orz, I cannot find x!\n");
    }
}