- T ≤ 100 , n ≤ 1 e 6 , ∑ k ≤ 1 e 6 T\le 100,n\le 1e6,\sum k\le 1e6 T≤100,n≤1e6,∑k≤1e6
- 简要题解:
首先可以暴力 d p dp dp, f i , j f_{i,j} fi,j 表示到了第 i i i 维,乘积为 j j j 的所有方案数的和
f i , j = ∑ l ∣ j f i − 1 , l ( j l ) b i f_{i,j}=\sum_{l|j}f_{i-1,l}(\frac{j}{l})^{b_i} fi,j=l∣j∑fi−1,l(lj)bi
写成前缀和的形式就是
S i , n = ∑ t = 1 n t b i S i − 1 , ⌊ n t ⌋ S_{i,n}=\sum_{t=1}^nt^{b_i}S_{i-1,\lfloor\frac{n}{t}\rfloor} Si,n=t=1∑ntbiSi−1,⌊tn⌋
于是可以 O ( T ∑ k n 3 / 4 ) O(T\sum kn^{3/4}) O(T∑kn3/4) 做
然后你发现把这个写成个类似生成函数的东西就是
f = ∏ i g b i r i f=\prod_i g_{b_i}^{r_i} f=i∏gbiri
其中 g b i ( t ) = ∑ x = 1 n x b i t i g_{b_i}(t)=\sum_{x=1}^nx^{b_i}t^i gbi(t)=∑x=1nxbiti
乘法及幂为狄利克雷卷积,于是可以对每一种 b i b_i bi 处理, O ( T n l o g ( n ) l o g ( k ) ) O(Tnlog(n)log(k)) O(Tnlog(n)log(k)) 做 - 注意到
f
f
f 其实是一个积性函数
我们考虑用 m i n _ 25 min\_25 min_25 筛来求出 f f f 的前缀和
只需要搞定 p p p 及 p k p^k pk 处的取值 -
p
p
p 只可能在一种
b
i
b_i
bi 中出现,对于某一种
b
i
b_i
bi,
f
(
p
)
=
r
i
∗
p
b
i
f(p)=r_i*p^{b_i}
f(p)=ri∗pbi
但 p k p^k pk 可以在多种 b i b_i bi 中出现,我们先求出一种的,然后将多种的卷积起来(指数做加法卷积) f ( p k ) = ( k + r i − 1 k ) ∗ p k ∗ b i f(p^k)=\binom{k+r_i-1}{k}*p^{k*b_i} f(pk)=(kk+ri−1)∗pk∗bi,对于这个预处理,复杂度为 ∑ l o g p ( n ) 2 \sum log_p(n)^2 ∑logp(n)2 可以忽略不计
所以我们可以在 O ( T n 3 / 4 l o g ( n ) ) O(T\frac{n^{3/4}}{log(n)}) O(Tlog(n)n3/4) 的时间内解决这个问题
【SCOI2019】超矩形 (min25筛)(狄利克雷卷积)(数数)
最新推荐文章于 2024-07-14 21:17:23 发布