HDU 5148 Cities (树形DP)*

题目链接：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;
}