BZOJ 2427: [HAOI2010]軟件安裝

本來是想找道 樹形依賴隨便水一下的 順便試試非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;
}