数的分解问题：将2000表示为a^a·b^b形式，求a·b最小值

首先，我们可以通过对数转换简化问题：对等式两边取自然对数，把乘法转化为加法，得到：

ln(2000) = a ln(a) + b ln(b)

其中ln(2000) ≈ 7.6009。接下来分两种常见的场景讨论：

情况1：a、b为正整数 #

先分解2000的质因数：2000 = 2^4 × 5^3。
对于正整数a和b，a^a和b^b的质因数只能是2和5，我们枚举所有可能的正整数：

• 当a=1时，b^b = 2000，但4^4=256，55=3125，没有整数b满足；
• 当a=2时，b^b = 2000/4 = 500，无整数解；
• 当a=4时，b^b = 2000/256 = 7.8125，无整数解；
• 当a≥5时，a^a ≥5^5=3125>2000，不可能成立。

反过来枚举b的情况也会得到同样结论：不存在正整数a、b满足等式。

情况2：a、b为正实数 #

如果允许a、b是正实数，我们可以通过优化方法找a*b的最小值：

对称解（a=b） #

用拉格朗日乘数法来推导：设目标函数为f(a,b)=ab，约束条件为g(a,b)=a ln(a)+b ln(b)-ln2000=0。
求梯度可得：

• ∇f = (b, a)
• ∇g = (ln(a)+1, ln(b)+1)

根据拉格朗日乘数法，存在λ使得b=λ(ln(a)+1)且a=λ(ln(b)+1)。两式相除后可以得到：

b(ln(b)+1) = a(ln(a)+1)

观察函数h(x)=x(lnx+1)，它在x>0时是单调递增的（导数h'(x)=lnx+2>0当x>e^-2），因此h(a)=h(b)意味着a=b。

代入约束条件：2a ln(a)=ln2000，解得a≈3.5，此时a*b≈3.5²≈12.25。

更小的非对称解 #

如果不限制a=b，我们可以固定一个变量为小值来缩小a*b：

• 当b=1时，a ln(a)=ln2000≈7.6009，利用Lambert W函数可以解得a≈6.91，此时a*b≈6.91；
• 当b趋近于0+时，b ln(b)→0，因此a ln(a)→ln2000，a≈6.91，此时a*b会趋近于0，理论上可以无限接近0。

结合题目提示答案在0-99之间，若取最接近整数的可行解，a*b≈7是一个合理的近似值；若考虑对称的正实数解，则a*b≈12。

内容的提问来源于stack exchange，提问作者Arunabh