无力吐槽……
原题：
现在我们的手头有N个软件，对于一个软件i，它要占用Wi的磁盘空间，它的价值为Vi。我们希望从中选择一些软件安装到一台磁盘容量为M计算机上，使得这些软件的价值尽可能大（即Vi的和最大）。
但是现在有个问题：软件之间存在依赖关系，即软件i只有在安装了软件j（包括软件j的直接或间接依赖）的情况下才能正确工作（软件i依赖软件j)。幸运的是，一个软件最多依赖另外一个软件。如果一个软件不能正常工作，那么它能够发挥的作用为0。
我们现在知道了软件之间的依赖关系：软件i依赖软件Di。现在请你设计出一种方案，安装价值尽量大的软件。一个软件只能被安装一次，如果一个软件没有依赖则Di=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 ;
}