[POJ2152]消防站解题报告

做这道题的时候，我一开始状态设计的太紧了，我设f(i)表示以i为根的子树自己搞定自己需要的最小代价，然后我枚举控制根的点。。然后发现转移是O(N^3)的。

看了论文之后发现其实改变一下状态就可以了，其实显然的一件事就是我们现在是在枚举两个点之间的关系了，那么我们不妨设f(i,j)为j控制i，且以i为根的子树都被完全搞定的最小代价；那么状态转移方程就呼之欲出了。

①这道题告诉我最大的一点就是当我们发现转移代价原大于状态代价的时候，我们不妨将加强状态以试图削弱转移！

我写代码的时候还犯了一些SB错误：

②树的边有N-1条。。而不是N条！

③为了避免一个点在队列中重复出现，我们要用bool数组判重或保留父亲什么的。。我本来记得判来着，后来不知道为啥sb忘判了。

④然后交到poj上的时候。。没看见多组数据直接交又WA了一次。

#include
#include
#include
using namespace std;
int f[1005][1005],g[1005],w[1005],d[1005];
int next[2005],succ[2005],ptr[1005],etot,l[2005];
inline void addedge(int u,int v,int weight){
	next[etot]=ptr[u],ptr[u]=etot,succ[etot]=v,l[etot++]=weight;
}
int dfna[1005],dfnl[1005],dfnr[1005],dfntot,dfnx[1005];
int q[1005],dis[1005];
bool p[1005];
#define inf 100000000
inline void dfs(int node){
	dfna[dfntot]=node,dfnx[node]=dfntot,dfnl[node]=dfntot++;
	int i;
	for(i=ptr[node];i;i=next[i])
		if(dfnx[succ[i]]==0)
			dfs(succ[i]);
	dfnr[node]=dfntot;
	int h=0,t=1;
	q[0]=node,p[node]=1;
	//printf("------%d:%d-[%d,%d)-------\n",node,dfnx[node],dfnl[node],dfnr[node]);
	for(;h!=t;++h){
		for(i=ptr[q[h]];i;i=next[i])
			if(dis[h]+l[i]<=d[node]&&!p[succ[i]]){
				dis[t]=dis[h]+l[i];
				q[t++]=succ[i];
				p[succ[i]]=1;
			}
		f[node][q[h]]=0;
		for(i=ptr[node];i;i=next[i])
			if(dfnx[succ[i]]>=dfnl[node]&&dfnx[succ[i]]=dfnl[succ[i]]&&dfnx[q[h]]=dfnl[node]&&dfnx[q[h]]'9')c=getchar();
	x=0;
	for(;c>='0'&&c<='9';c=getchar())x=x*10+c-'0';
}
int main(){
	freopen("poj_2152.in","r",stdin);
	freopen("poj_2152.out","w",stdout);
	int T;
	in(T);
	int N,i;
	int u,v,l;
	int ans;
	while(T--){
		etot=1,dfntot=1;
		memset(ptr,0,sizeof(ptr));
		memset(f,60,sizeof(f));
		memset(g,60,sizeof(g));
		memset(dfnx,0,sizeof(dfnx));
		in(N);
		for(i=1;i<=N;++i)in(w[i]);
		for(i=1;i<=N;++i)in(d[i]);
		for(i=N;--i;){
			in(u),in(v),in(l);
			addedge(u,v,l),addedge(v,u,l);
		}
		dfs(1);
		ans=0x7fffffff;
		for(i=N;i;--i)ans=min(ans,f[1][i]+w[i]);
		printf("%d\n",ans);
	}
}