BZOJ 2159 「国家集训队」Crash 的文明世界(第二类斯特林数,换根DP)BZOJ计划

Posted 繁凡さん

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ 2159 「国家集训队」Crash 的文明世界(第二类斯特林数,换根DP)BZOJ计划相关的知识,希望对你有一定的参考价值。

整理的算法模板合集: ACM模板

点我看算法全家桶系列!!!

实际上是一个全新的精炼模板整合计划


题目链接

https://hydro.ac/d/bzoj/p/2159

hydro 的 BZOJ 修复工程 !(我也去领了一点题慢慢修着玩,这题就是我修的嘿嘿嘿)

Problem

Crash 小朋友最近迷上了一款游戏——文明5 (Civilization V)。在这个游戏中,玩家可以建立和发展自己的国家,通过外交和别的国家交流,或是通过战争征服别的国家。

现在 Crash 已经拥有了一个 n n n 个城市的国家,这些城市之间通过道路相连。由于建设道路是有花费的,因此 Crash 只修建了 n − 1 n-1 n1 条道路连接这些城市,不过可以保证任意两个城市都有路径相通。

在游戏中,Crash 需要选择一个城市作为他的国家的首都,选择首都需要考虑很多指标,有一个指标是这样的:

S ( i ) = ∑ j = 1 n d i s t ( i , j ) k S(i) = \\sum_{j = 1}^{n}{\\rm dist}(i, j) ^ k S(i)=j=1ndist(i,j)k

其中 S ( i ) S(i) S(i) 表示第 i i i 个城市的指标值, d i s t ( i , j ) {\\rm dist}(i, j) dist(i,j) 表示第 i i i 个城市到第 j j j 个城市需要经过的道路条数的最小值, k k k 为一个常数且为正整数。

因此 Crash 交给你一个简单的任务:给出城市之间的道路,对于每个城市,输出这个城市的指标值,由于指标值可能会很大,所以你只需要输出这个数   m o d     10007 \\bmod\\ 10007 mod 10007 的值。

1 ≤ n ≤ 5 × 1 0 4 1\\le n\\le 5\\times 10^4 1n5×104 1 ≤ k ≤ 150 1\\le k\\le 150 1k150

Solution

没有边权的无根树^ q ^

k = 1 k=1 k=1 时显然是经典的换根DP。

考虑如何处理 d i s t k \\mathrm{dist}^k distk

我们知道第二类斯特林数 S ( n , m ) = { n m } S(n,m)= \\left \\{ {\\begin{matrix} n\\\\ m \\end{matrix}} \\right \\} S(n,m)={nm} 满足

m n = ∑ i = 0 n { i n } i ! ( m i ) m^{n}=\\sum_{i=0}^{n} \\left \\{ {\\begin{matrix} i\\\\ n \\end{matrix}} \\right \\} i ! {m\\choose i} mn=i=0n{in}i!(im)

则有:

S x = ∑ i = 1 n d i s ( i , x ) k = ∑ i = 1 n ∑ j = 0 k { k j } ( d i s ( i , x ) j ) j ! = ∑ j = 0 k { k j } j ! ∑ i = 1 n ( d i s ( i , x ) j ) = ∑ j = 0 k { k j } j ! ∑ i = 1 n ( ( d i s ( i , x ) − 1 j ) + ( d i s ( i , x − 1 ) j − 1 ) ) \\begin{aligned}S_{x} &=\\sum_{i=1}^{n} \\mathrm{d i s}(i, x)^{k} \\\\&=\\sum_{i=1}^{n} \\sum_{j=0}^{k} \\left \\{ {\\begin{matrix} k\\\\ j \\end{matrix}} \\right \\} {\\mathrm{d i s}(i, x)\\choose j}j ! \\\\&=\\sum_{j=0}^{k}\\left\\{\\begin{array}{c}k \\\\j\\end{array}\\right\\} j ! \\sum_{i=1}^{n} {\\mathrm{d i s}(i, x)\\choose j} \\\\&=\\sum_{j=0}^{k}\\left\\{\\begin{array}{c}k \\\\j\\end{array}\\right\\} j ! \\sum_{i=1}^{n}\\left({\\mathrm{d i s}(i, x)-1\\choose j}+{\\mathrm{d i s}(i, x-1)\\choose j-1}\\right)\\end{aligned} Sx=i=1ndis(i,x)k=i=1nj=0k{kj}(jdis(i,x))j!=j=0k{kj}j!i=1n(jdis(i,x))=j=0k{kj}j!i=1n((jdis(i,x)1)+(j1dis(i,x1)))

考虑继续做换根DP。

f [ x , j ] f[x,j] f[x,j] 表示以 x x x 为根的子树中,上述表达式的后半部分的答案。

显然有转移方程:

f [ x , j ] = ∑ s o n ∈ x _ s o n f [ s o n , j ] + f [ s o n , j − 1 ] f[x,j]=\\sum\\limits_{\\mathrm{son}\\in \\mathrm{x\\_{son}}}f[\\mathrm{son}, j]+f[\\mathrm{son}, j-1] f[x,j]=sonx_sonf[son,j]+f[son,j1]

预处理之后直接进行换根DP即可。

时间复杂度 O ( n k ) O(nk) O(nk)

Code

以上是关于BZOJ 2159 「国家集训队」Crash 的文明世界(第二类斯特林数,换根DP)BZOJ计划的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 2159 「国家集训队」Crash 的文明世界(第二类斯特林数,换根DP)BZOJ计划

bzoj 2159: Crash 的文明世界

bzoj 2159: Crash 的文明世界

BZOJ2159 : Crash 的文明世界

bzoj2159: Crash 的文明世界

bzoj 2159 Crash 的文明世界 —— 第二类斯特林数+树形DP