[ACNOI2022]进制转换

题目

题目背景
我昨天还以为，$\mathsf{H}\mathsf{a}\mathsf{n}\mathsf{d}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{e}\mathsf{v}\mathsf{i}\mathsf{l}$ 是 日象 家的天才，不可能有人战胜它的 马桶白眼；今天的 $\mathsf{D}\mathsf{D}(\mathsf{X}\mathsf{Y}\mathsf{X})$ 的洞见就超越了它。他的体内，有着比人类更高明、更强大之物，炎翼鸟 也只是祂的宠物——祂是，普林斯普！

“你的那双眼睛，比我的 马桶白眼 更厉害，” 质朴的 $\mathsf{H}\mathsf{a}\mathsf{n}\mathsf{d}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{e}\mathsf{v}\mathsf{i}\mathsf{l}$ 浮现在 $\mathsf{D}\mathsf{D}(\mathsf{X}\mathsf{Y}\mathsf{X})$ 眼前，“现在 $\mathsf{O}\mathsf{n}\mathsf{e}\mathsf{I}\mathsf{n}\mathsf{D}\mathsf{a}\mathsf{r}\mathsf{k}$ 正身处光明之中，而你能将他打入黑暗。”

题目描述
给出十进制数 $y$ 和 $b$ 进制数 $l$，保证 $l$ 的数码都在 $[0,9]$ 中。你需要选择 $b$ 使得 $y_{(10)}⩾l_{(b)}$ 且 $y_{(b)}$ 的每一位都在 $[0,9]$ 之中。求最大的 $b$ 。这里 $x_{(m)}$ 表示将 $x$ 的十进制表示 $x={\sum }_{i⩾0}{a}_{i}10^i(a_i\in [0,9])$ 重写为 $m$ 进制下 $x_{(m)}={\sum }_{i⩾0}{a}_i{m}^i$ 。

数据范围与约定
为了避免某些问题，保证 $10⩽l⩽y⩽10^{18}$ 。数据组数 $T⩽7\times 10^4$ 。

思路

乍看之下，好像 $y⩾l_{(b)}$ 是很简单的，直接二分就得到了 $b$ 的范围。然而下一个问题是：每个数位都在 $[0,9]$ 该怎么翻译？

难以数位 $\mathtt{d}\mathtt{p}$ 的时候，不妨从数学的角度考虑。第 $k$ 位上的数其实就是
$$\lfloor \frac{y}{b^k}\rfloor -b\lfloor \frac{y}{b^{k+1}}\rfloor ⩽9$$

只能想到整除分块，对于 $y⩾b^4$ 就用上式找到合法的 $b$，而 $y<b^4$ 就直接枚举 $y_{(b)}$ 数位。时间复杂度较高。

优化点在于，可行的 $b$ 其实非常非常少。上式中，很多整除分块的结果都是解不出一个合法的 $b$ 的。我们只要能够充分剪枝就成功了。

第一个想法是，最高位的限制最强，所以我们要求 $b^k$ 是严格最高位。但是这还不够强。问题还是在整除分块上——所以要大胆弃用原方法，转而 枚举该数位值。设为 $a$，则
$$a{b}^k⩽y⩽a{b}^k+9\sum _{i=0}^{k-1}{b}^i$$

显然 $b⩾10$ 且 $k⩾1$ 。我们会惊讶地发现，可行的 $b$ 的数量 少得可怜。因为后面简单放缩（或者直接用 $b$ 进制数思考）可以发现
$$y⩽a{b}^k+10b^{k-1}$$

也就是让它再进一位。那么，若找到最小的 $b$ 使得原不等式成立，则 $a\cdot (b+\delta {)}^k⩾a{b}^k+ak\delta b^{k-1}$，在 $\delta ⩾\frac{10}{ak}$ 时就会超过 $a{b}^k+10b^{k-1}$ 即 $y$ 的上界，所以可行的 $b$ 的数量只有 $\mathcal{O}(\delta )=\mathcal{O}(\frac{10}{ak})$ 个！

于是只需要找到 $b_{max}=\sqrt[k]{\frac{y}{a}}$，然后暴力判断 $b\in [b_{max}-\delta ,b_{max}]$ 是否为解。时间复杂度几乎是 $\mathcal{O}(10{\log }_{10}l\ln y)$ 的。

代码

即使分析出来很正确，还是需要剪枝卡常，这就很不厚道了……

#include <cstdio> // JZM yydJUNK!!!
#include <iostream> // XJX yyds!!!
#include <algorithm> // XYX yydLONELY!!!
#include <cstring> // (the STRONG long for LONELINESS)
#include <cctype> // ZXY yydSISTER!!!
#include <cmath>
using namespace std;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
typedef long long llong;
# define double (long double)
inline llong readint(){
	llong a = 0; int c = getchar(), f = 1;
	for(; !isdigit(c); c=getchar())
		if(c == '-') f = -f;
	for(; isdigit(c); c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(const llong &x){
	if(x > 9) writeint(x/10);
	putchar(int((x%10)^48));
}

bool check(llong y, const llong &l, const llong &b){
	llong all = 0, v = 1;
	for(int d; y; y/=b){
		if((d = int(y%b)) > 9) return false;
		all += d*v, v *= 10;
	}
	return all >= l; // satisfying
}
inline llong qkpow(llong b, int q){
	llong a = 1; // no modules
	for(; q; q>>=1,b*=b) if(q&1) a *= b;
	return a;
}

llong powten[19];
int main(){
	rep(i,int(powten[0]=1),18) powten[i] = 10*powten[i-1];
	for(int T=int(readint()); T; --T){
		llong y = readint(), l = readint();
		for(int k=1; powten[k<<1]<y; ++k) rep(a,1,9){
			if(a+1 <= (l-1)/powten[k]) continue;
			llong b = llong(round(pow(double(y)/a,double(1)/k)+1));
			llong lb = llong(round(pow((y-(qkpow(b,k)
				-1)/(b-1)*9)/a,double(1)/k)-1));
			for(lb=max(lb,9ll); b!=lb; --b) // bottom line
				if(check(y,l,b)){ writeint(b); goto SCHEME; }
		}
		for(int b=100; true; --b) // bruteforce
			if(check(y,l,b)){ writeint(b); break; }
		SCHEME: putchar('\n');
	}
	return 0;
}

后记

据称，想到这个做法的原因是，一个数字的最高位对数字大小的影响最大。或许与本题有些相似，或许没有。