【SCOI2019】超矩形 （min25筛）（狄利克雷卷积）（数数）

传送门

• $T\le 100,n\le 1e6,\sum k\le 1e6$
• 简要题解：
  首先可以暴力 $dp$，$f_{i,j}$ 表示到了第 $i$ 维，乘积为 $j$ 的所有方案数的和
  $$f_{i,j}=\sum _{l|j}{f}_{i-1,l}(\frac{j}{l}{)}^{b_i}$$
  写成前缀和的形式就是
  $$S_{i,n}=\sum _{t=1}^n{t}^{b_i}{S}_{i-1,\lfloor \frac{n}{t}\rfloor }$$
  于是可以 $O(T\sum k{n}^{3/4})$ 做
  然后你发现把这个写成个类似生成函数的东西就是
  $$f=\prod _i{g}_{b_i}^{r_i}$$
  其中 $g_{b_i}(t)={\sum }_{x=1}^n{x}^{b_i}{t}^i$
  乘法及幂为狄利克雷卷积，于是可以对每一种 $b_i$ 处理， $O(Tnlog(n)log(k))$ 做
• 注意到 $f$ 其实是一个积性函数
  我们考虑用 $min\_25$ 筛来求出 $f$ 的前缀和
  只需要搞定 $p$ 及 $p^k$ 处的取值
• $p$ 只可能在一种 $b_i$ 中出现，对于某一种 $b_i$，$f(p)=r_i\ast p^{b_i}$
  但 $p^k$ 可以在多种 $b_i$ 中出现，我们先求出一种的，然后将多种的卷积起来（指数做加法卷积）$f(p^k)=\left(\genfrac{}{}{0ex}{}{k+r_i-1}{k}\right)\ast p^{k\ast b_i}$，对于这个预处理，复杂度为 $\sum lo{g}_p(n{)}^2$ 可以忽略不计
  所以我们可以在 $O(T\frac{n^{3/4}}{log(n)})$ 的时间内解决这个问题