[POJ2152] Fire [树形dp]

[$\mathfrak{L}\mathfrak{i}\mathfrak{n}\mathfrak{k}$]

---

题意：在树里面选一些点建firehouse，满足所有的点在D_i范围内都有firehouse。
求mincost。

---

这个题目很容易让人联想到某道消防局的设立
当然实际上差了很多（

each vertice都要被control，朴素地想就是对于每一个点枚举一个firehouse
然后就搞出来一个n!做法当然是不行的啦。

所以得有更加确定的状态。尝试约束更新方向。
这道题从上往下很不好处理约束，考虑从下往上。

考虑x被某个点i覆盖。
要明确的是，这个时候复杂度达到了上限，不能再加状态了。
那么，把这个设置为状态f(x,i)
表示以x为根的子树被完全control，x被i control的最小花费。（最优子结构）
同时，f(x,i)需要由f(child,…)更新而来（以覆盖子树）
现在的问题就是，w(i)有没有被重复计算。

那么对于每个子树里的点考虑它是不是被i更新,
当然不能真的考虑子树，会爆炸的。
先只考虑儿子看看，说不定可以提前消除影响？
是的话在累加的时候消除重复部分（w(i)）和不是的情况比较
这样的话初始化f(x,i)为w(i)，枚举每个子结点，
f(x,i)加上min{f(child,i)-w(i),min{f(child)}}
（其中min{f(child)}可以另外保存；而且它就是每个点的答案）
发现这种情况下并不用对于x考虑子树内每个点，因为非direct son的重复早就被消除了。

---

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define add_edge(a,b,c) nxt[++tot]=head[a],head[a]=tot,to[tot]=b,val[tot]=c
int t, n, tot=0;
int head[1005];
int nxt[2005];
int to[2005];
int val[2005];
int w[1005], d[1005];
int f[1005][1005];
int ans[1005];
int dis[1005][1005];
void workDis(const int &root, int x, int fa) {
	for (int i = head[x]; i; i = nxt[i]) {
		if (to[i] == fa) continue;
		dis[root][to[i]] = dis[root][x] + val[i];
		workDis(root, to[i], x);
	}
}
void dfs(int x, int fa) {
	ans[x] = 0x3f3f3f3f;
	for (int i = head[x]; i; i = nxt[i]) {
		if (to[i] == fa) continue;
		dfs(to[i], x);
	}
	for (int i = 1; i <= n; ++i) {
		if (dis[x][i] > d[x]) continue;
		f[x][i] = w[i];
		for (int j = head[x]; j; j = nxt[j]) {
			if (to[j] == fa) continue;
			f[x][i] += min(f[to[j]][i] - w[i], ans[to[j]]);
		}
		ans[x] = min(ans[x], f[x][i]);
	}
}
int main() {
	scanf("%d", &t);
	while (t--) {
		tot = 0;
		memset(f, 0x3f, sizeof(f));
		memset(head, 0, sizeof(head));
		scanf("%d", &n);
		for (int i = 1; i <= n; ++i) scanf("%d", &w[i]);
		for (int i = 1; i <= n; ++i) scanf("%d", &d[i]);
		for (int a, b, c, i = 1; i < n; ++i) {
			scanf("%d%d%d", &a, &b, &c);
			add_edge(a, b, c);
			add_edge(b, a, c);
		}
		for (int i = 1; i <= n; ++i) {
			workDis(i, i, 0);
		}
		dfs(1, 0);
		printf("%d\n", ans[1]);
	}
	return 0;
}