BSGS(进阶篇，解法+题目)

进阶篇

问题：

$$x^a\equiv b(modp)$$

解法一：

设$g$是p的一个原根，根据原根的性质存在$c$，满足$g^c\equiv x$成立，同理存在$t$，满足$g^t\equiv b$成立。
$$\begin{gathered} \therefore (g^c{)}^a\equiv b(modp) \\ (g^a{)}^c\equiv b(modp) \\ {g}^a已知，所以我们就能用基础篇直接求解c，也就是一个特解。 \end{gathered}$$

解法二：

设$g$是p的一个原根，根据原根的性质存在$c$，满足$g^c\equiv x$成立，同理存在$t$，满足$g^t\equiv b$成立。
$$\begin{gathered} \therefore g^{ac}\equiv g^t(modp) \\ 根据阶的性质有 \\ \therefore ac\equiv t(mod\phi (p)) \\ 根据exgcd求出c，也是一个特解。 \end{gathered}$$

求全部的解：

我们在已知一个特解$g^c$的情况下，我们要得到全部解。
$$\begin{gathered} \because g^{\phi (p)}\equiv 1(modp) \\ \therefore \forall t\in \mathbb{Z},x^a\equiv g^{ca+\frac{t\phi (p)}{a}}(modp) \\ \therefore \forall t\in \mathbb{Z}且a|t\phi (p),x\equiv g^{c+\frac{t\phi (p)}{a}}(modp) \\ \because a|t\phi (p) \\ \therefore \frac{a}{gcd(a,\phi (p))}|t \\ \therefore t=\frac{a}{gcd(a,\phi (p))}i \\ \therefore 全部的解为：\forall i\in \mathbb{Z},x\equiv g^{c+\frac{\phi (p)}{gcd(a,\phi (p))}i}(modp) \\ （解的个数为gcd(a,\phi (p))） \end{gathered}$$
​

题目：

P3306 [SDOI2013] 随机数生成器

根据题意将原式化成等比数列的形式。再去推式子。

$$\begin{gathered} x_{i+1}\equiv a{x}_i+b(modp) \\ {x}_{i+1}+\frac{b}{a-1}\equiv a(x_i+\frac{b}{a-1})(modp) \\ {x}_n\equiv t(modp) \\ {a}^{n-1}(x_1+\frac{b}{a-1})-\frac{b}{a-1}\equiv t(modp) \\ {a}^{n-1}\equiv (t+\frac{b}{a-1})\ast inv(x_1+\frac{b}{a-1})(modp) \end{gathered}$$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll qpow(ll x, ll y, ll mod) {
    ll ans = 1;
    while(y) {
        if(y & 1) ans = ans * x % mod;
        x = x * x % mod;
        y >>= 1;
    }
    return ans;
}

unordered_map<ll, ll> mp;
ll bsgs(ll a, ll b, ll p) {
    if(a % p == 0) return -1;
    mp.clear();
    ll k = ceil(sqrt(p));
    for(int i=0; i<=k; i++) {
        mp[b] = i;
        b = b * a % p;
    }
    ll aa = qpow(a, k, p), A = aa;
    for(int i=1; i<=k; i++) {
        if(mp[aa]) {
            return 1ll * i * k - mp[aa] + 1;
        }
        aa = aa * A % p;
    }
    return -1;
}
int main() {
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
    int t;
    cin >> t;
    while(t--) {
        ll p, a, b, x, t;
        cin >> p >> a >> b >> x >> t;
        if(x == t) {
            cout << 1 << endl;
            continue;
        }
        if(a == 0) {
            if(b == t) cout << 2 << endl;
            else cout << -1 << endl;
            continue;
        }
        if(a == 1) {
            if(b == 0) cout << -1  << endl;
            else {
                ll k = qpow(b, p-2, p);
                cout << ((t-x+p)%p*k)  % p + 1 << endl;
            }
            continue;
        }

        ll tmp = b*qpow(a-1, p-2, p) % p;
        t = (t + tmp) % p;
        t = t * qpow((x+tmp)%p, p-2, p) % p;
        ll ans = bsgs(a, t, p);
        cout << ans << endl;
    }
    return 0;
}

261. Discrete Roots模板

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll qpow(ll x, ll y, ll mod) {
	ll ans = 1;
	while(y) {
		if(y & 1) ans = ans * x % mod;
		x = x * x % mod;
		y >>= 1;
	}
	return ans;
}
ll G(ll p) {//求原根
	if(p == 2) return 1;
	vector<ll> tmp;
	int phi = p-1, n = phi;
	for(int i=2; 1ll*i*i<=n; i++) {
		if(n % i == 0) {
			tmp.push_back(i);
			while(n % i == 0) n /= i;
		}
	}
	if(n > 1) tmp.push_back(n);
	for(int i=1; i<=p; i++) {
		bool flag = 0;
		for(auto it : tmp) {
			if(qpow(i, phi/it, p) == 1) {
				flag = 1;
				break;
			}
		}
		if(!flag) return i;
	}
	return -1;
}
unordered_map<ll, ll> mp;
ll bsgs(ll a, ll b, ll p) {
	if(a % p == 0 && b != 0) return -1;
	mp.clear();
	int k = ceil(sqrt(p));
	for(int i=0; i<=k; i++) {
		mp[b] = i;
		b = b * a % p;
	}
	ll aa = qpow(a, k, p), A = aa;
	for(int i=1; i<=k; i++) {
		if(mp[A]) return 1ll * i * k - mp[A];
		A = A * aa % p;
	}
	return -1;
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
	ll p, k, a;
	while(cin >> p >> k >> a) {
		if(a == 0) {
			cout << "1\n0\n";
			continue;
		}
		ll g = G(p), gk = qpow(g, k, p);
		ll x0 = bsgs(gk, a, p), t = bsgs(g, a, p);
		if(x0 == -1) {
			cout << "0\n";
			continue;
		}

		vector<ll> ans;
		ll d = __gcd(p-1, k), mod = (p-1)/d;
		x0 = x0 % mod;
		for(int i=0; i<d; i++) {//最多只有d个不同的解。
			ans.push_back(qpow(g, ( x0+i*mod%(p-1) )%(p-1), p));
		}
		sort(ans.begin(), ans.end());
		cout << ans.size() << endl;
		for(auto it : ans) cout << it << ' ';
		cout << endl;
	}
	return 0;
}