bzoj2427-tarjan+树形dp+背包

传送门
一道练习综合代码能力的好题
先把被依赖的点连向依赖这个点的点，再tarjan缩点，缩点后会形成多棵树，对每棵树进行dp，最后进行一次背包即可，注意要倒着枚举。
代码：

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
struct wp{
    int v,w;
}a[110],b[110];
int low[110],dfn[110];
int cnt,fa[110];
int n,m,z[110];
struct edg{
    int to,next;
}e[110];
int size,head[110],now;
int vis[110];
int ta,st[110];
int tot,rot[110];
int f[110][510];
int du[110];
int ans[110][510];
void add(int x,int y)
{
    size++;
    e[size].to=y;
    e[size].next=head[x];
    head[x]=size;
}
int dfs(int x,int f_a)
{
    for (int i=head[x];i;i=e[i].next)
    {
        int y=e[i].to;
        if (y!=f_a) 
        {
            dfs(y,x);
            for (int j=m;j>=0;j--)
            {
                for (int k=j;k>=0;k--)
                f[x][j]=max(f[x][j],f[x][k]+f[y][j-k]);
            }
        }
    }
    for (int i=m;i>=0;i--)
        if (i>=b[x].w)
            f[x][i]=f[x][i-b[x].w]+b[x].v;
        else f[x][i]=0;
}

void tarjan(int x,int &now)
{
    vis[x]=1;
    st[++ta]=x; 
    dfn[x]=low[x]=++now;
    for (int i=head[x];i;i=e[i].next)
    {
        int y=e[i].to;
        if (vis[y]==0)
        {
            tarjan(y,now);
            low[x]=min(low[x],low[y]);
        }
        else if (vis[y]==1) low[x]=min(low[x],dfn[y]); 
    }
    if (low[x]==dfn[x])
    {
        cnt++;
        while (st[ta]!=x)
        {
            fa[st[ta]]=cnt;
            vis[st[ta]]=2;
            ta--;
        }
        fa[st[ta]]=cnt;
        vis[st[ta]]=2;
        ta--;
    }
}
int main()
{
//  freopen("2427.in","r",stdin);
//  freopen("2427.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        scanf("%d",&a[i].w);
    for (int i=1;i<=n;i++)
        scanf("%d",&a[i].v);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&z[i]);
        if (z[i])add(z[i],i);
    }
    for (int i=1;i<=n;i++)
    {
        if (vis[i]==0) now=0,tarjan(i,now); 
    }
    memset(head,0,sizeof(head));
    memset(e,0,sizeof(e));
    size=0;

    for (int i=1;i<=n;i++)
    {
        if (z[i])
        { 
            if (fa[z[i]]!=fa[i])
                add(fa[z[i]],fa[i]),du[fa[i]]++; 
        }
        b[fa[i]].w+=a[i].w;
        b[fa[i]].v+=a[i].v;
    }
    for (int i=1;i<=cnt;i++)
        if (du[i]==0)
            dfs(i,0),rot[++tot]=i;
    for (int i=1;i<=tot;i++)
    {
        for (int j=m;j>=0;j--)
        {
            ans[i][j]=ans[i-1][j];
            for (int k=j;k>=0;k--)
            {
                ans[i][j]=max(ans[i][j],ans[i-1][j-k]+f[rot[i]][k]);
            }
        }
    }
    int maxn=0;
    for (int i=0;i<=m;i++)
        maxn=max(maxn,ans[tot][i]);
    printf("%d",maxn);

}