解题报告：一道关于进制转换的题

Problem

给定一个十进制数$y$，请选择一个进制$b$，使得$y$转成$b$进制后大于或等于$b$进制下的$l$，且每位均为$0$到$9$之间的数码。求最大的$b$。

数据范围：$10 \le l \le y \le 10^{18}$

Analysis

乍看起来是道很简单的题。$(y)_{10} \ge (l)_b$有两种考虑方向，可以把$l$按$b$进制解析，得到 $\sum_i l_ib^i \le y$。这样左边对$b$是单调的，可以考虑二分。

另一种方向是把$y$化成$b$进制来与$l$比较，设为$m+1$位数$a$。这样就是要找一个$a \ge l$满足 $\sum_{i=0}^m a_ib^i = y$。左边$a_mb^m$比其余项加起来还大很多，这样可以给$b$很强的限制。

但这只考虑了一条主要限制，细想之下，另一条限制更加棘手，即$y$在$b$进制下每位均在$0$到$9$之间。这使得大部分$b$都不合法，并且分布方式难以捉摸，答案可能离满足$(y)_{10} \ge (l)_b$的最大的$b$相去甚远，因此前一种思路告吹。

所以关键在于解这样一个式子

$$\begin{cases} a_0b^0 + a_1b^1 + \cdots + a_mb^m = y\\0\le a_i\le9\\a \ge l \end{cases}$$
这里$10 \le y,l,b,a \le 10^{18}$，因此$1 \le m \le 18$

Solution 1

把$a_0$移到右边，得 $b|(y-a_0)$，因此把$y,y-1,\cdots,y-9$这10个数用Pollard-Rho分解（其实我不会写rho算法），枚举一个约数$b$来check。

check就是暴力把$y$转成$b$进制，单次复杂度$\log_by\le\lg y$

Solution 2

由于$b^m \le 10^{18}$，所以值得考虑的$b$的分布是极不均匀的。可以在$b$小时枚举$b$，在$b$大时$m$很小从而可以枚举$a$。例如枚举完$a<10^4$后，剩下的一定有$m\ge4$，于是$b \le \sqrt[4]{10^{18}} < 31623$

现在给定一个$a$，如何判断是否存在$b$使$y=(a)_b$。我们发现$a_mb^m \le y < (a_m+1)b^m$并不能给出足够强的限制，因为当$m=1,a_m=1$时它给出的区间可以长达$10^{18}/2$。

但是这里用关于$b$的单调性来二分是没问题的。复杂度约为$10000\log\lg+31613\lg$

Solution 3

$a_mb^m \le y < (a_m+1)b^m$不足以给出足够强的限制，其实是因为我们太粗犷了。式中认为$a_i < b$，但我们知道更强的限制是$a_i \le 9$。所以重新列出式子：

$$a_mb^m \le y \le a_mb^m + 9\times\frac{b^m-1}{b-1}$$

看来$b$确实只能在$\sqrt[m]{y/a_m}$很近的地方，可是如何严谨地证明这一点呢。式中$b$和$b^m$同时出现，不大能解。但是我们只需要估算就行了。

不等式的左边很好解，$b\le\sqrt[m]{y/a_m}$，即$b_{max} = \lfloor\sqrt[m]{y/a_m}\rfloor$是$b$的上界。我们可以把右边分式部分中的$b$换成$b_{max}$，这样合法是因为相当于把$a$末尾一串9999的权值放大。这样就可以把$b$的范围严谨地解出来了：

$$\sqrt[m]\frac{y - 9\times\frac{b_{max}^m-1}{b_{max}-1}}{a_m} \le b \le \sqrt[m]\frac{y}{a_m}$$

可见在$b$较大时，这个上下界的距离很近，几乎是常数（如果有人会算这东西到底有多小，跪求大佬指点）。枚举这个范围内的$b$来check

对$b\le100$的暴力，再枚举$1 \le m \le 8, 0 \le a_m \le 9$计算。复杂度大约不超过$300\lg$

Code

#include <cmath>
#include <cstdio>
typedef long long ll;
inline ll and10(ll l) {return l > 10 ? l : 10;}

ll calc(ll y, ll b) {
    ll a = 0, p = 1;
    do{ ll t = y % b;
        if (t > 9) return -1;
        a += t * p;
        p *= 10;
    } while (y /= b);
    return a;
}

int main() {
    ll y, l;
    scanf("%lld%lld",&y,&l);
    for (int m = log10l(l); m < 9; m++)
        for (int a = 1; a <= 9; a++)
            for (ll br = powl(y / a + 1, 1.0L/m),
                bl = and10(powl((y - (powl(br, m) - 1) / (br - 1) * 9) / a, 1.0L/m));
                br >= bl; br--) if (calc(y, br) >= l) {printf("%lld\n",br); return 0;}
    int b = 100;
    while (calc(y, b) < l) b--;
    printf("%d\n",b);
    return 0;
}