【BZOJ2427】【HAOI2010】软件安装

无力吐槽……

原题：

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

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

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

0<=N<=100, 0<=M<=500

恩这道题刚开始的时候我是有思路的，就是树上有限制的背包嘛

但是为了防止写了数h+然后发现思路歪了还是去看了一下题解，然后才发现题中可以有环

环之间相互依赖，要选一个则其它必选，就可以直接看成一个点

然后就是tarjian缩强连通分量辣（环也是强连通分量

但是晚上调了1h+，知道错误在哪里但是怎么都想不明白为什么，因为很晚了急着走所以就去网上对着AC代码改程序过掉了，然而依旧不能想明白为什么，电脑关机一下错误的代码也没有了，心好累，感觉应该想不明白这个问题了

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int oo=168430090;
int rd(){int z=0,mk=1;  char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')mk=-1;  ch=getchar();}
    while(ch>='0'&&ch<='9'){z=(z<<3)+(z<<1)+ch-'0';  ch=getchar();}
    return z*mk;
}
struct ddd{int nxt,y;}e[5100000];  int lk[1100000],ltp=0,indgr[210];
inline void ist(int x,int y){  e[++ltp].nxt=lk[x],lk[x]=ltp,e[ltp].y=y,++indgr[y];}
int n,m;
int w[210],v[210];
int dfn[210],low[210],dfscnt=0;
int stck[210],tp=0;  bool vstd[210];
int grp[210],grpcnt=0;
int f[210][510];
bool flg[210];
void tj(int x){
    flg[x]=true;
    dfn[x]=low[x]=++dfscnt;
    stck[++tp]=x,vstd[x]=true;
    for(int i=lk[x];i;i=e[i].nxt){
        if(!dfn[e[i].y]){  tj(e[i].y);  low[x]=min(low[x],low[e[i].y]);}
        else if(vstd[e[i].y])  low[x]=min(low[x],dfn[e[i].y]);
    }
    if(dfn[x]==low[x]){
        ++grpcnt;  int tmp;
        do{  tmp=stck[tp--],vstd[tmp]=false,grp[tmp]=grpcnt+n;}while(tmp!=x);
    }
}
void dfs(int x){
    for(int i=lk[x];i;i=e[i].nxt){
        dfs(e[i].y);
        for(int j=m;j>=0;--j)for(int k=j;k>=0;--k)
            f[x][j]=max(f[x][j],f[x][j-k]+f[e[i].y][k]);
    }
    for(int i=m;i>=w[x];--i)  f[x][i]=f[x][i-w[x]]+v[x];
    for(int i=0;i<w[x];++i)  f[x][i]=0;
}
int main(){//freopen("ddd.in","r",stdin);
    memset(vstd,0,sizeof(vstd));
    memset(f,0,sizeof(f));
    memset(flg,0,sizeof(flg));
    cin>>n>>m;
    for(int i=1;i<=n;++i)  w[i]=rd();
    for(int i=1;i<=n;++i)  v[i]=rd();
    for(int i=1;i<=n;++i)  ist(rd(),i);
    for(int i=1;i<=n;++i)if(!dfn[i])  tj(i);
    for(int i=1;i<=n;++i){
        w[grp[i]]+=w[i],v[grp[i]]+=v[i];
        for(int j=lk[i];j;j=e[j].nxt)if(grp[e[j].y]!=grp[i])
            ist(grp[i],grp[e[j].y]);
    }
    for(int i=1;i<=grpcnt;++i)if(!indgr[i+n])  ist(grpcnt+n+1,i+n);
    dfs(grpcnt+n+1);
    int mx=0;
    for(int i=0;i<=m;++i)  mx=max(mx,f[grpcnt+n+1][i]);
    cout<<mx<<endl;
    return 0;
}