【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=;
 int rd(){int z=,mk=;  char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')mk=-;  ch=getchar();}
     while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-'';  ch=getchar();}
     return z*mk;
 }
 struct ddd{int nxt,y;}e[];  int lk[],ltp=,indgr[];
 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[],v[];
 int dfn[],low[],dfscnt=;
 int stck[],tp=;  bool vstd[];
 int grp[],grpcnt=;
 int f[][];
 bool flg[];
 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>=;--j)for(int k=j;k>=;--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=;i<w[x];++i)  f[x][i]=;
 }
 int main(){//freopen("ddd.in","r",stdin);
     memset(vstd,,sizeof(vstd));
     memset(f,,sizeof(f));
     memset(flg,,sizeof(flg));
     cin>>n>>m;
     for(int i=;i<=n;++i)  w[i]=rd();
     for(int i=;i<=n;++i)  v[i]=rd();
     for(int i=;i<=n;++i)  ist(rd(),i);
     for(int i=;i<=n;++i)if(!dfn[i])  tj(i);
     for(int i=;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=;i<=grpcnt;++i)if(!indgr[i+n])  ist(grpcnt+n+,i+n);
     dfs(grpcnt+n+);
     int mx=;
     for(int i=;i<=m;++i)  mx=max(mx,f[grpcnt+n+][i]);
     cout<<mx<<endl;
     return ;
 }