[数论][构造][离散对数] Codeforces 913 G. Power Substring

$Solution$

设$a$长度为$n$。
可以尝试构造一个数$b$，使得$a\cdot10^m+b$成为其后缀。
就等价于$x=a\cdot10^m+b\equiv2^k\pmod {10^{n+m}}$
构造$k\ge n+m$，那要满足这样的$k$就等价于$2^{n+m}\mid x\wedge5\nmid x$。
假如说现在有$x'a\cdot10^m$，要怎么构造满足上述两个条件的数$x$呢？
若$2^{n+m}\nmid x'$，则$x=2^{n+m}-x'\bmod2^{n+m}$。
若$5\mid x'$，则$x=x'+2^{n+m}$。
如何解出$k$？
$2^k\equiv x\pmod {10^{n+m}}$，有$2^{n+m}\mid2^k,2^{n+m}\mid x,2^{n+m}\mid10^{n+m}$。
所以$2^{k-n-m}\equiv {x\over{2^{n+m}}}\pmod{5^{n+m}}$。设$y={x\over{2^{n+m}}}$。
$k=n+m+\log_2y$。
这里有一个结论$2$是$5^\alpha$的原根。
一个更广泛的结论：$g$是奇质数$p$的原根，有$g$是$p^\alpha$的原根。
证明：

• 设$d$是$\varphi(p)$的质因子，有$g^{{p^{\alpha-1}(p-1)}\over d}={\big(g^{p-1\over d}\big)}^{p^{\alpha-1}}\equiv g^{p-1\over d}\not\equiv1\pmod{p^\alpha}$
• 只需要证$g^{p^{\alpha - 2}}\not\equiv1\pmod{p^\alpha}$，发现$x^p\equiv1\pmod{p^a}$，$x^p\equiv1\pmod{p^b}是b\le a$的充分条件。所以就等价于$g\not\equiv1\pmod{p^2}$。

但这个模数很大，BSGS是跑不过去的。
设$d_{y,i}$为模$p^i$意义下的$log_gy$，有$d_{i+1}=d_i+j\varphi(p^i),j\in[0,p)$，枚举就好了。

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;

int test, n;
ll a, b, t, l, x, y, MOD, tp, p5;

inline int diglen(ll x) {
    static int cnt;
    for (cnt = 0; x; cnt++) x /= 10;
    return cnt;
}
inline ll Mul(ll a, ll b, ll MOD) {
    ll c = 0;
    while (b > 0) {
        if (b & 1) c = (c + a) % MOD;
        a = (a + a) % MOD; b >>= 1;
    }
    return c;
}
inline ll Pow(ll a, ll b, ll MOD) {
    ll c = 1;
    while (b > 0) {
        if (b & 1) c = Mul(c, a, MOD);
        a = Mul(a, a, MOD); b >>= 1;
    }
    return c;
}

int main(void) {
    freopen("1.in", "r", stdin);
    freopen("1.out", "w", stdout);
    cin >> test;
    while (test--) {
        cin >> a; n = diglen(a);
        tp = 1;
        for (int m = 0; ; m++) {
            x = a; MOD = 1ll << (n + m);
            if (x % MOD != 0) x += MOD - x % MOD;
            if (x % 5 == 0) x += MOD;
            if (x % MOD == 0 && x % 5 != 0 && x - a < tp) {
                y = x >> (n + m);
                t = 0;
                if (y % 5 == 0) t = -1;
                if (y % 5 == 1) t = 0;
                if (y % 5 == 2) t = 1;
                if (y % 5 == 3) t = 3;
                if (y % 5 == 4) t = 2;
                p5 = 5;
                for (int i = 1; i < n + m; i++){
                    while (Pow(2, t, p5 * 5) != y % (p5 * 5))
                        t += p5 / 5 * 4;
                    p5 *= 5;
                }
                cout << t + n + m << endl;
                break;
            }
            a *= 10; tp *= 10;
        }
    }
}