[POJ2152]消防站解题报告

最新推荐文章于 2019-05-21 09:07:00 发布

原创最新推荐文章于 2019-05-21 09:07:00 发布 · 1.1k 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#dp

DP 专栏收录该内容

43 篇文章

订阅专栏

本文详细介绍了解决POJ2152问题的过程，从最初设计的状态转移导致的时间复杂度过高，到阅读论文后优化状态为f(i,j)，表示j控制i时的最小代价。文章强调了动态规划中状态设计的重要性，并列举了编程过程中的几个常见错误，如树的边数计算、重复节点处理以及多组数据处理的疏忽，这些都是导致WA的原因。" 2626200,214930,中国IT业低迷原因探讨,"['行业分析', '教育体制', '技术发展', '市场竞争', '创新环境']

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

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

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

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

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

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

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

#include<cstdio>
#include<cstring>
#include<iostream>
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]]<dfnr[node])
				if(dfnx[q[h]]>=dfnl[succ[i]]&&dfnx[q[h]]<dfnr[succ[i]])f[node][q[h]]=min(f[node][q[h]]+f[succ[i]][q[h]],inf);
				else f[node][q[h]]=min(f[node][q[h]]+min(f[succ[i]][q[h]],g[succ[i]]),inf);
		if(dfnx[q[h]]>=dfnl[node]&&dfnx[q[h]]<dfnr[node])g[node]=min(g[node],f[node][q[h]]+w[q[h]]);
		//printf(":%d(%d)\n",q[h],f[node][q[h]]);
	}
	while(t--)p[q[t]]=0;
	//printf("---:%d\n",g[node]);
}
inline void in(int &x){
	char c=getchar();
	while(c<'0'||c>'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);
	}
}