Yukicoder No.417 チューリップバブル

,

http://yukicoder.me/problems/no/417

solution

よくある感じの木DP。ある頂点$i$からの部分木で時間$t$かけて得られる金額の最大値でDP。$\mathrm{dp}_i : M \to \mathbb{N}$を$N$回mergeするので全体で$O(NM^2)$。 使った辺は必ず$2$度使うことを考えると、DP表の大きさは正確には$\lfloor \frac{M}{2}\rfloor + 1$で済む。

implementation

#include <iostream>
#include <vector>
#include <algorithm>
#include <tuple>
#include <functional>
#define repeat(i,n) for (int i = 0; (i) < (n); ++(i))
#define whole(f,x,...) ([&](decltype((x)) y) { return (f)(begin(y), end(y), ## __VA_ARGS__); })(x)
typedef long long ll;
using namespace std;
template <class T> void setmax(T & a, T const & b) { if (a < b) a = b; }
const ll inf = ll(1e18)+9;
int main() {
    // input
    int n, m; cin >> n >> m;
    vector<int> u(n); repeat (i,n) cin >> u[i];
    vector<vector<pair<int,int> > > g(n);
    repeat (i,n-1) {
        int a, b, c; cin >> a >> b >> c;
        g[a].emplace_back(b, c);
        g[b].emplace_back(a, c);
    }
    // compute
    function<vector<ll> (int, int)> dfs = [&](int i, int p) {
        vector<ll> cur(m/2+1, - inf);
        cur[0] = u[i];
        for (auto it : g[i]) {
            int j, c; tie(j, c) = it;
            if (j == p) continue;
            vector<ll> nxt(m/2+1, - inf);
            vector<ll> prv = dfs(j, i);
            repeat (x,m/2+1) {
                setmax(nxt[x], cur[x]);
                repeat (y,m/2+1) {
                    if (x+y+c >= m/2+1) break;
                    setmax(nxt[x+y+c], cur[x] + prv[y]);
                }
            }
            cur.swap(nxt);
        }
        return cur;
    };
    vector<ll> dp = dfs(0, -1);
    // output
    cout << *whole(max_element, dp) << endl;
    return 0;
}