背包问题(dsu on tree)

背包问题

解：

感觉就是一个暴力乱搞，我们直接暴力插入一个点进背包的复杂度是$m$的。先考虑一个暴力做法：每次清空背包，然后把这个子树的点全部加进背包求最值。这样的复杂度是$n^2m$的。但是考虑这个方法并卡不满，因为每棵子树的大小不能达到$n$。在二叉树的情况下，这个做法是$nlog^nm$的。就可以通过这道题。所以我们暴力可以拿到60分。
考虑优化这个做法，我们把前面的背包记录下来，每次用最大的子树继续做背包，我们就会去掉一个最大的子树的复杂度。
经过重链剖分的证明，我们可以发现每个点最多被重新放进背包$log_n$次，莫名复杂度就对了。而且还贼快。
code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
struct lxy{
    int to,next;
}eg[30005];
int n,m,cnt,a1;
int v[30005],head[30005],size[30005],wson[30005],root[30005];
int w[30005],ans[30005];
int f[105][505];
stack <int> d;

int readit(){
    int a=0;char s=getchar();
    while(s<'0'||s>'9') s=getchar();
    while(s>='0'&&s<='9')
      a=a*10+s-'0',s=getchar();
    return a;
}

void add(int op,int ed){
    eg[++cnt].next=head[op];
    eg[cnt].to=ed;
    head[op]=cnt;
}

void insert(int aim,int u){
    for(int i=m;i>=v[u];i--)
      if(f[aim][i]<f[aim][i-v[u]]+w[u])
        f[aim][i]=f[aim][i-v[u]]+w[u];
}

void dfs(int u){
    size[u]=1;int p=0;
    for(int i=head[u];i!=-1;i=eg[i].next){
        dfs(eg[i].to);
        size[u]+=size[eg[i].to];
        if(size[eg[i].to]>p) p=size[eg[i].to],wson[u]=eg[i].to;
    }
}

void baoli(int aim,int u){
    insert(aim,u);
    for(int i=head[u];i!=-1;i=eg[i].next)
      baoli(aim,eg[i].to);
}

void solve(int u)
{
    if(wson[u]!=0)solve(wson[u]);
    for(int i=head[u];i!=-1;i=eg[i].next)if(eg[i].to!=wson[u]){
        solve(eg[i].to);
        memset(f[root[eg[i].to]],0,sizeof(f[root[eg[i].to]]));d.push(root[eg[i].to]);
        baoli(root[wson[u]],eg[i].to);
    }
    if(wson[u]==0) root[u]=d.top(),d.pop();
    else root[u]=root[wson[u]];root[wson[u]]=0;
    insert(root[u],u);
    ans[u]=f[root[u]][m];
}

int main()
{
    memset(head,-1,sizeof(head));
    for(int i=100;i>=1;i--) d.push(i);
    n=readit();m=readit();
    for(int i=1;i<=n;i++)
      v[i]=readit(),w[i]=readit();
    for(int i=2;i<=n;i++)
      a1=readit(),add(a1,i);
    dfs(1);solve(1);
    for(int i=1;i<=n;i++) printf("%d ",ans[i]);
    return 0;
}