树上dp——软件安装

[HAOI2010] 软件安装

题目描述

现在我们的手头有 N N N 个软件,对于一个软件 i i i,它要占用 W i W_i Wi 的磁盘空间,它的价值为 V i V_i Vi。我们希望从中选择一些软件安装到一台磁盘容量为 M M M 计算机上,使得这些软件的价值尽可能大(即 V i V_i Vi 的和最大)。

但是现在有个问题:软件之间存在依赖关系,即软件 i i i 只有在安装了软件 j j j(包括软件 j j j 的直接或间接依赖)的情况下才能正确工作(软件 i i i 依赖软件 j j j)。幸运的是,一个软件最多依赖另外一个软件。如果一个软件不能正常工作,那么它能够发挥的作用为 0 0 0

我们现在知道了软件之间的依赖关系:软件 i i i 依赖软件 D i D_i Di。现在请你设计出一种方案,安装价值尽量大的软件。一个软件只能被安装一次,如果一个软件没有依赖则 D i = 0 D_i=0 Di=0,这时只要这个软件安装了,它就能正常工作。

输入格式

第 1 行: N , M ( 0 ≤ N ≤ 100 , 0 ≤ M ≤ 500 ) N,M(0\leq N\leq 100, 0\leq M\leq 500) N,M(0N100,0M500)

第 2 行: W 1 , W 2 , . . . W i , . . . , W n ( 0 ≤ W i ≤ M ) W_1,W_2, ... W_i, ..., W_n (0\leq W_i\leq M) W1,W2,...Wi,...,Wn(0WiM)

第 3 行: V 1 , V 2 , . . . , V i , . . . , V n ( 0 ≤ V i ≤ 1000 ) V_1, V_2, ..., V_i, ..., V_n (0\leq V_i\leq 1000) V1,V2,...,Vi,...,Vn(0Vi1000)

第 4 行: D 1 , D 2 , . . . , D i , . . . , D n ( 0 ≤ D i ≤ N , D i ≠ i ) D_1, D_2, ..., D_i, ..., D_n (0\leq D_i\leq N, D_i≠i) D1,D2,...,Di,...,Dn(0DiN,Di=i)

输出格式

一个整数,代表最大价值。

样例 #1

样例输入 #1

3 10
5 5 6
2 3 4
0 1 1

样例输出 #1

5

思路

  • 首先,对于本道题,题目描述中并没有说这个东西一定是一棵树,它有可能是一张带环图,也可能是一个森林,因此我们得缩点。
  • 我们发现:这个图建完后是一个带环图。1,2,3号节点互为依赖条件。一旦有其中任何一个点没有被选择,整个环里面的所有内容都无法被使用。我们重新观察这张图,我们发现1,2,3号节点的共性:我们可以对它整个选择,也可以对它整个不选择,对其中的任何一个真子集都没有意义。对于环中点有共性的图,我们考虑用缩点来解决问题。
  • 也就是说我如果要选联通分量内的一个,那么整个联通分量都得选(核心)

代码

#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 220,M = 1010;

int e[N],h[N],hs[N],w[N],ne[N],idx;
int dfn[N],low[N];
int v[N];
bool in_stk[N];
int stk[N],top,scc_size[N];
int scc_cnt,id[N],timecuocuo;
int din[N];
int scc_v[N],scc_w[N];
int f[N][M];//以i为根,体积为j,值为f[i][j]
int n,m;

void add(int h[],int a,int b){
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}

void targin(int u){
    low[u]=dfn[u]=++timecuocuo;
    stk[++top]=u,in_stk[u]=true;
    
    for(int i=h[u];~i;i=ne[i]){
        int j=e[i];
        
        if(!dfn[j]){
            targin(j);
            low[u]=min(low[u],low[j]);
        }else if(in_stk[j]){
            low[u]=min(low[u],dfn[j]);
        }
    }
    
    if(low[u]==dfn[u]){
        scc_cnt++;
        int y;
        do{
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            scc_size[scc_cnt]++;
        }while(y!=u);
    }
}

void dfs(int u){
    for(int i=m;i>=scc_w[u];i--){
        f[u][i]=scc_v[u];
    }
    
    for(int i=hs[u];~i;i=ne[i]){
        int p=e[i];
        
        dfs(p);
        
        for(int j=m;j>=scc_w[u];j--){
            for(int k=0;k<=j-scc_w[u];k++){//选根节点必须滴
                f[u][j]=max(f[u][j],f[p][k]+f[u][j-k]);
            }
        }
    }
}

int main(){
    cin>>n>>m;
    memset(h,-1,sizeof h);
    memset(hs,-1,sizeof hs);
    for(int i=1;i<=n;i++)cin>>w[i];
    
    for(int i=1;i<=n;i++)cin>>v[i];
    
    for(int i=1;i<=n;i++){
        int b;
        cin>>b;
        if(b!=0)
        add(h,b,i);
    }
    
    for(int i=1;i<=n;i++){
        if(!dfn[i]){
            targin(i);
        }
    }
    
    for(int i=1;i<=n;i++){
        for(int j=h[i];~j;j=ne[j]){
            int k=e[j];
            int a=id[i],b=id[k];
            if(a!=b){
                add(hs,a,b);
                din[b]++;
            }
        }
        scc_v[id[i]]+=v[i];
        scc_w[id[i]]+=w[i];
    }
    //超级原点
    for(int i=1;i<=scc_cnt;i++){
        if(!din[i]){
            add(hs,0,i);
        }
    }
    
    dfs(0);
    
    cout<<f[0][m];
    
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

green qwq

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值