当前位置: 代码迷 >> 综合 >> HDU 5148 Cities (树形DP)*
  详细解决方案

HDU 5148 Cities (树形DP)*

热度:17   发布时间:2023-11-15 14:47:07.0

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5148

#include<bits/stdc++.h>
using namespace std;#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define read(x,y) scanf("%d%d",&x,&y)
#define ll unsigned long long#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
const int  maxn =2e3+5;
const int mod=1e9+7;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意:从一棵树中选出k个城市,
这k个点构成的权重就是任意点对距离和,问最小值。首先不难发现这k个点一定是要连在一起的,
那么我们就可以设置dp状态了,
令dp[u][p]为以u为根节点在u子树中选择p个的最小值,
那么思路就清晰了,假设更新到当前节点的第i个子节点,
之前的i-1个都更新过了,现在的dp[u]状态数组是之前更新的,
那么我们在更新当前节点v时,就要考虑如何利用v完成的状态去进一步更新u,
首先我们可以枚举从v子树中选择p个,利用之前的u子树的j个去更新j+p状态,
注意j的状态要倒着扫描,根据DP更新时间序。
时间复杂度:O(k^2*n)
*/ll dp[maxn][60],INF;///树形DP结构
///前向星
struct edge{int u,nxt, w;}e[maxn<<1];
int head[maxn],tot=0;
void init(){memset(head,-1,sizeof(head));tot=0;}
void add(int x,int y,int z){e[tot]=edge{y,head[x],z};head[x]=tot++;}
///数据域
int n,k,x,y,z;
ll ans;
///树形DP
void Dp(int u,int pre)
{dp[u][1]=0;for(int i=head[u];~i;i=e[i].nxt){int v=e[i].u;if(v==pre) continue;Dp(v,u);for(int j=k-1;j;j--)///选择已经更新过的j状态{for(int p=1;p+j<=k;p++)///选择当前子树中的p个点{if(dp[u][j]==INF||dp[v][p]==INF) continue;dp[u][j+p]=min(dp[u][j+p],dp[u][j]+dp[v][p]+2LL*p*(k-p)*e[i].w);///j是已经用过的状态}}}ans=min(ans,dp[u][k]);
}
int main()
{int t;scanf("%d",&t);while(t--){init(); scanf("%d%d",&n,&k);for(int i=1;i<n;i++){scanf("%d%d%d",&x,&y,&z);add(x,y,z),add(y,x,z);}memset(dp,0xf,sizeof(dp));INF=dp[0][0];ans=INF;  Dp(1,-1);printf("%I64d\n",ans);}return 0;
}