BZOJ 2427: [HAOI2010]软件安装

来源：互联网发布：html12canvas.js 编辑：程序博客网时间：2024/06/05 08:31

本来是想找道树形依赖随便水一下的顺便试试非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;}

阅读全文

1 0