本来是想找道 树形依赖随便水一下的 顺便试试非dfs序的那种解法
谁知道这题居然有环,靠,那我打多个tarjan咯。。。
这1987B 让我最近提交记录的代码长度继续保持在了2000B以内
QAQ (惨啊 天天刷水题
#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof a)
using namespace std;
const int N=102,inf=1e9+7;
char B[1<<14],*S=B,*T=B;
#define gc (S==T&&(T=(S=B)+fread(B,1,1<<14,stdin),S==T)?-1:*S++)
inline int read(){
int x=0,f=1; char ch=gc;
while(ch<'0' || ch>'9'){if(ch=='-')f=-1; ch=gc;}
while(ch>='0' && ch<='9'){x=(x<<1)+(x<<3)+ch-'0'; ch=gc;}
return x*f;
}
struct node{int x,y,nex;}a[N];
int w[N],v[N],len,fir[N],m,f[N][N*5];
void ins(int x,int y){
a[++len]=(node){x,y,fir[x]},fir[x]=len;
}
void dp(int x){
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]=-inf;
for(int k=fir[x];k;k=a[k].nex){
int y=a[k].y;
for(int i=0;i<=m;++i)f[y][i]=f[x][i];
dp(y);
for(int i=0;i<=m;++i)f[x][i]=max(f[x][i],f[y][i]);
}
}
int dfn[N],low[N],st[N],tp,id,cnt,bl[N];
bool in[N],o[N];
void dfs(int x){
low[x]=dfn[x]=++id,st[++tp]=x,in[x]=1;
for(int k=fir[x];k;k=a[k].nex){
int y=a[k].y;
if(!dfn[y]) dfs(y),low[x]=min(low[x],low[y]);
else if(in[y])low[x]=min(low[x],dfn[y]);
}
if(dfn[x]==low[x]){
++cnt; int s=0,i,nw=0,nv=0;
do{
i=st[tp--],s++,bl[i]=x;
in[i]=0,nw+=w[i],nv+=v[i];
}while(i!=x);
if(s!=1) w[x]=nw,v[x]=nv,o[x]=1;
}
}
int main(){
int n=read(); m=read();
for(int i=1;i<=n;++i) w[i]=read();
for(int i=1;i<=n;++i) v[i]=read();
for(int i=1;i<=n;++i){ int x=read(); ins(x,i); }
for(int i=1;i<=n;++i) if(!dfn[i]) dfs(i);
int e=len; len=0; me(fir,0);
for(int i=1;i<=e;++i){
int x=a[i].x,y=a[i].y;
if(bl[x]!=bl[y])ins(bl[x],bl[y]);
}
for(int i=1;i<=n;++i) if(o[i]) ins(0,i);
dp(0);
printf("%d\n",f[0][m]);
return 0;
}