Yukicoder 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;
}