ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

P4362 [NOI2002] 贪吃的九头龙 题解

2021-11-15 22:33:20  阅读:166  来源: 互联网

标签:大头 min int 题解 ch edge NOI2002 P4362 dp


Description

Luogu传送门

Solution

非常巧妙地树形 \(dp\),刚开始完全不知道如何设状态 QwQ。

言归正传,我们对于每个点可以随便分配个某一个头,只有大头需要必须有 \(k\) 个果子,所以我们的状态里只需要大头就够了。

那么设 \(dp_{x, i}\) 表示以 \(x\) 为根的子树中给大头 \(i\) 个果子最小的难受值。

但是我们发现根本没法转移,所以考虑再多记录一点东西。

设:\(dp_{x, i, 0/1}\),0/1 表示 \(x\) 节点是否分给大头,那么就有转移方程:

\[dp_{x, i, 0} = min\{ dp_{y, j, 0} + dp_{x, i - j, 0} + (m == 2) \times edge_{x, y},\ \ dp_{y, j, 1} + dp_{x, i - j, 0}\} \]

解释一下:

  • 第一项:\(x\) 和它的一个儿子 \(y\) 都不给大头,那么就得把它们分给其他不同的头,但如果此时 \(m = 2\),即除了大头只有一个头,那就只能加上边的难受值了。
  • 第二项:\(y\) 给大头,\(x\) 不给大头,\(edge_{x, y}\) 没有贡献。

再来看 \(dp_{x, i, 1}\) 的转移:

\[dp_{x, i, 1} = min\{dp_{y, j, 1} + dp_{x, i - j, 1} + edge_{x, y},\ \ dp_{y, j, 0} + dp_{x, i - j, 1} \} \]

这个就比较好理解了吧,跟上面的差不多,不多解释了。

但是直接转移会有问题,因为在同一轮转移过程中 \(dp\) 数组的值是不能变的,不然就有问题,所以转移前先给它转存一下,\(f_{i, 0/1} = dp_{x, i, 0/i}\),然后用 \(f_{i, 0/1}\) 去转移。

再放一遍方程吧:

\[dp_{x, i, 0} = min\{ dp_{y, j, 0} + f_{i - j, 0} + (m == 2) \times edge_{x, y},\ \ dp_{y, j, 1} + f_{i - j, 0}\} \]

\[dp_{x, i, 1} = min\{dp_{y, j, 1} + f_{i - j, 1} + edge_{x, y},\ \ dp_{y, j, 0} + f_{i - j, 1} \} \]

Code

#include <bits/stdc++.h>

using namespace std;

namespace IO{
    inline int read(){
        int x = 0;
        char ch = getchar();
        while(!isdigit(ch)) ch = getchar();
        while(isdigit(ch)) x = (x << 3) + (x << 1) + ch - '0', ch = getchar();
        return x;
    }

    template <typename T> inline void write(T x){
        if(x > 9) write(x / 10);
        putchar(x % 10 + '0');
    }
}
using namespace IO;

const int N = 310;
int n, m, k;
struct node{
    int v, w;
};
vector <node> g[N];
int head[N], tot;
int dp[N][N][2], f[N][2];

inline void dfs(int x, int fa){
    dp[x][0][0] = dp[x][1][1] = 0;
    for(auto nxt : g[x]){
        int y = nxt.v;
        if(y == fa) continue;
        dfs(y, x);
        memcpy(f, dp[x], sizeof(dp[x]));
        memset(dp[x], 0x3f, sizeof(dp[x]));
        for(int i = 0; i <= k; ++i)
            for(int j = 0; j <= i; ++j){
                dp[x][i][0] = min(dp[x][i][0], min(dp[y][j][0] + f[i - j][0] + (m == 2) * nxt.w, dp[y][j][1] + f[i - j][0]));
                dp[x][i][1] = min(dp[x][i][1], min(dp[y][j][1] + f[i - j][1] + nxt.w, dp[y][j][0] + f[i - j][1]));
            }
    }
}

int main(){
    n = read(), m = read(), k = read();
    if(n - k < m - 1) return puts("-1"), 0;
    for(int i = 1; i < n; ++i){
        int u = read(), v = read(), w = read();
        g[u].push_back((node){v, w});
        g[v].push_back((node){u, w});
    }
    memset(dp, 0x3f, sizeof(dp));
    dfs(1, 1);
    printf("%d\n", dp[1][k][1]);
    return 0;
}

\[\_EOF\_ \]

标签:大头,min,int,题解,ch,edge,NOI2002,P4362,dp
来源: https://www.cnblogs.com/xixike/p/15558882.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有