LA4015

此题是2007 成都 区域赛的题目  由于刚接触树形dp 耗时较久。且此题是我在想了一下午无果后 去洗了个澡  哎。。。洗澡的时候就胡思乱想 然后就想了一个类似于最短路的思想来做。

我不知道我这样的理解是否正确。。。。 如果大神发现我理解错了 望加以指正！！非常感谢

想了我一下午。 我开始想到的状态方程是包含当前点的最大距离。所以我开始的时候是用dp【x】【y】表示x到y的最小长度。但是想了想不对啊。他要给定查询条件。

也就是说我从当前X点到各个点的最小长度都算出来。那么我前面想的就有问题了。因为少了一维状态来存各个点的距离

所以此题需要三维。同样第三维保存包含不包含当前节点所以有两个状态。然后第1维表示从当前节点出发 第二维j表示以当前节点出发到J个儿子节点的距离。

那么我当树输进去后 预处理一遍到个点距离就可以了

因为题意是以跟节点出发的

那么我们需要考虑的是。从当前节点出发到子树，那么也可以从子树打个来回然后在到其他子树（我是理解为以当前节点为出发既做出发节点又作中转节点）

那么我如果作为了中转节点，那么我距离就需要加上到其中一个子树的距离的两倍（打个来回）如果不做中转节点，那么距离就只要去的距离就可以了。

方程可以这么想：

y：特定儿子节点；

j：特定儿子节点的访问儿子个数

i：当前节点访问儿子节点个数

dp【当前节点】【子树节点】【0】0表示不含当前节点的距离=min（……，dp【当前节点】【i-j】【0】+dp【y】【j】【0】+len*2）

上面这个表示当前节点是中转节点但是不含当前节点，那么既然是中转节点  len就要来回*2

dp【当前节点】【子树节点】【1】1表示涵盖当前节点的距离=min（……，dp【当前节点】【i-j】【1】+dp【y】【j】【1】+len）

上面这个表示当前节点到Y节点去了就不往回走了。就是从其他节点到了X然后从X到Y节点下去。

dp【当前节点】【子树节点】【1】=min（……，dp【当前节点】【i-j】【0】+dp【y】【j】【0】+len*2）

上面这个表示当前节点访问到了X节点但是同样以X节点做中转节点  len就要来回*2

这里我贴出我的代码。因为树形dp我真没把握完全理解，特别对于刚入门的来说。如果希望大神们能对我的理解加以指正  在此感激不尽

我用前向星见图的，如果边数较多的话我的代码不一定好。我是参照LA2038的建图做的

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
const int MAXN= 600;
const int INF= 0x3f3f3f3f;
struct node
{
	int v,len,next;
}edge[MAXN<<1];
int dp[MAXN][MAXN][2],son[MAXN],head[MAXN<<1];
int cnt=0,n;
void addedge(int u,int v,int val)
{
	edge[cnt].v=v;
	edge[cnt].next=head[u];
	edge[cnt].len=val;
	head[u]=cnt++;
	
	edge[cnt].v=u;
	edge[cnt].next=head[v];
	edge[cnt].len=val;
	head[v]=cnt++;	
}
void init()
{
	memset(head,-1,sizeof(head));
	memset(son,0,sizeof(son));
	for(int i=0;i<=n;i++)
	for(int j=0;j<=n;j++)
	dp[i][j][0]=dp[i][j][1]=INF;
	cnt=0;
}
void dfs(int x,int fa)
{
	dp[x][1][0]=0;dp[x][1][1]=0;
	son[x]=1;
	for(int i=head[x];i!=-1;i=edge[i].next)
	{
		int y=edge[i].v;
		int len=edge[i].len;
		if(y!=fa)
		{
			dfs(y,x);
			son[x]+=son[y];
			for(int i=son[x];i>1;i--)
			for(int j=1;j<=son[y];j++)
			{
				dp[x][i][0]=min(dp[x][i][0],dp[x][i-j][0]+dp[y][j][0]+len*2);
				dp[x][i][1]=min(dp[x][i][1],dp[x][i-j][0]+dp[y][j][1]+len);
				dp[x][i][1]=min(dp[x][i][1],dp[x][i-j][1]+dp[y][j][0]+len*2);
			}
		}
	}
}
int main()
{
	int cas=0;
	while(scanf("%d",&n)!=EOF,n)
	{
		init();

		for(int i=1;i<n;i++)
		{
			int u,v,len;
			scanf("%d%d%d",&u,&v,&len);
			addedge(u+1,v+1,len);
		}
		dfs(1,0);
		int Q;
		scanf("%d",&Q);
		printf("Case %d:\n",++cas);
		for(int i=1;i<=Q;i++)
		{
			int dis;
			scanf("%d",&dis);
			int ans=1;
			for(int i=n;i>=1;i--)
			{
				if(dp[1][i][0]<=dis || dp[1][i][1]<=dis)
				ans=max(ans,i);
			}
			printf("%d\n",ans);
		}		
	}
	return 0;
}