欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  IT编程

BZOJ1812: [Ioi2005]riv(树形dp)

程序员文章站 2022-07-12 08:32:55
题意 "题目链接" Sol 首先一个很显然的思路是直接用$f[i][j] / g[i][j]$表示$i$的子树中选了$j$个节点,该节点是否选的最小权值。但是直接这样然后按照树形背包的套路转移的话会有一种情况无法处理,就是说该节点不选,儿子节点也不选,这样我们就不清楚儿子节点的子节点的贡献了 一种暴 ......

题意

题目链接

sol

首先一个很显然的思路是直接用\(f[i][j] / g[i][j]\)表示\(i\)的子树中选了\(j\)个节点,该节点是否选的最小权值。但是直接这样然后按照树形背包的套路转移的话会有一种情况无法处理,就是说该节点不选,儿子节点也不选,这样我们就不清楚儿子节点的子节点的贡献了

一种暴力的做法是钦定该节点选,并重新枚举子树中的所有节点,转移出dp值之后背包合并

最后再把\(0\)号节点的合并一次

#include<bits/stdc++.h>
#define chmin(x, y) (x = x < y ? x : y)
#define chmax(x, y) (x = x > y ? x : y)
using namespace std;
const int maxn = 1001, inf = 2e9 + 10;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int n, k, a[maxn], dis[maxn], siz[maxn], f[maxn][maxn], g[maxn][maxn], ans;
vector<int> v[maxn];
void dfs2(int x, int fa, int root) {
    f[x][0] = dis[root] * a[x];
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; if(to == fa) continue;
        dfs2(to, x, root);
        for(int j = min(k, siz[x]); ~j; j--) 
            for(int k = 0; k <= min(j, siz[to]); k++) 
                chmax(f[x][j], f[to][k] + f[x][j - k]);
    }
    for(int i = siz[x]; i; i--) chmax(f[x][i], g[x][i]);
}
void dfs(int x, int fa) {
    dis[x] += dis[fa];  siz[x] = 1; 
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i];
        if(to == fa) continue;
        dfs(to, x); siz[x] += siz[to];
    }
    g[x][0] = 0;
    memset(f, 0, sizeof(f));
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; if(to == fa) continue;
        dfs2(to, x, x);
        for(int j = min(k, siz[x]); ~j; j--) 
            for(int k = 0; k <= min(j, siz[to]); k++) 
                chmax(g[x][j], f[to][k] + g[x][j - k]);
    }
    for(int i = min(k, siz[x]); i; i--) g[x][i] = g[x][i - 1] + dis[x] * a[x];
}
int main() {
    n = read(); k = read();
    for(int i = 1; i <= n; i++) {
        a[i] = read(); int fa = read(); dis[i] = read();
        v[fa].push_back(i); 
    }
    dfs(0, -1);
    for(int i = 1; i <= n; i++) ans += dis[i] * a[i];
    memset(f, 0, sizeof(f));
    int x = 0;
    for(int i = 0; i < v[x].size(); i++) {
        int to = v[x][i]; 
        dfs2(to, x, x);
        for(int j = min(k, siz[x]); ~j; j--) 
            for(int k = 0; k <= min(j, siz[to]); k++) 
                chmax(f[x][j], f[to][k] + f[x][j - k]);
    }
    int out = inf;
    for(int i = 0; i <= k; i++) out = min(out, ans - f[0][i]);
    printf("%d\n", out);
    return 0;
}