bzoj2427 [HAOI2010]軟體安裝

標籤：des   blog   io   ar   os   sp   for   strong   on   

Description

現在我們的手頭有N個軟體，對於一個軟體i，它要佔用Wi的磁碟空間，它的價值為Vi。我們希望從中選擇一些軟體安裝到一台磁碟容量為M電腦上，使得這些軟體的價值儘可能大（即Vi的和最大）。

但是現在有個問題：軟體之間存在依賴關係，即軟體i只有在安裝了軟體j（包括軟體j的直接或間接依賴）的情況下才能正確工作（軟體i依賴軟體j)。幸運的是，一個軟體最多依賴另外一個軟體。如果一個軟體不能正常工作，那麼它能夠發揮的作用為0。

我們現在知道了軟體之間的依賴關係：軟體i依賴軟體Di。現在請你設計出一種方案，安裝價值盡量大的軟體。一個軟體只能被安裝一次，如果一個軟體沒有依賴則Di=0，這時只要這個軟體安裝了，它就能正常工作。

 

Input

第1行：N, M  （0<=N<=100, 0<=M<=500）
      第2行：W1, W2, ... Wi, ..., Wn （0<=Wi<=M ）
      第3行：V1, V2, ..., Vi, ..., Vn  （0<=Vi<=1000 ）
      第4行：D1, D2, ..., Di, ..., Dn （0<=Di<=N, Di≠i ）

 

Output

一個整數，代表最大價值。

 

Sample Input3 10
5 5 6
2 3 4
0 1 1 Sample Output5

 

tarjan縮完點樹形dp

黃巨大說這種題很多但是我是第一次做啊……

#include<cstdio>#include<iostream>#include<cstring>#include<cstdlib>#include<algorithm>#include<cmath>#include<queue>#include<deque>#include<set>#include<map>#include<ctime>#define LL long long#define inf 0x7ffffff#define pa pair<int,int>#define pi 3.1415926535897932384626433832795028841971#define N 110#define M 510using namespace std;inline LL read(){    LL x=0,f=1;char ch=getchar();    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}    return x*f;}inline void write(LL a){    if (a<0){printf("-");a=-a;}    if (a>=10)write(a/10);    putchar(a%10+‘0‘);}inline void writeln(LL a){write(a);printf("\n");}struct edge{int to,next;}e1[2*N],e2[2*N];int head1[N],head2[N];int c[N],w[N],fa[N],belong[N];int C[N],W[N],f[N][M];int dfn[N],low[N];int zhan[N],top;bool inset[N],mrk[N];int size[N];int n,m,cnt1,cnt2,cnt3,tt,ans;inline void ins1(int u,int v){    e1[++cnt1].to=v;    e1[cnt1].next=head1[u];    head1[u]=cnt1;}inline void ins2(int u,int v){    e2[++cnt2].to=v;    e2[cnt2].next=head2[u];    head2[u]=cnt2;}inline void dfs(int x){    low[x]=dfn[x]=++tt;    zhan[++top]=x;inset[x]=1;    for (int i=head1[x];i;i=e1[i].next)    {        if (!dfn[e1[i].to])        {            dfs(e1[i].to);            low[x]=min(low[x],low[e1[i].to]);        }else if (inset[e1[i].to])            low[x]=min(low[x],dfn[e1[i].to]);    }    if (dfn[x]==low[x])    {        cnt3++;        int p=-1;        while (p!=x)        {            p=zhan[top--];            inset[p]=0;            belong[p]=cnt3;            C[cnt3]+=c[p];            W[cnt3]+=w[p];        }    }}inline void tarjan(){    for (int i=1;i<=n;i++)      if (!dfn[i])dfs(i);}inline void treedp(int x){    if(mrk[x])return;mrk[x]=1;    for (int i=head2[x];i;i=e2[i].next)      if (!mrk[e2[i].to])      {        treedp(e2[i].to);        for (int j=m-C[x];j>=0;j--)          for (int k=j;k>=0;k--)            f[x][j]=max(f[x][j],f[x][j-k]+f[e2[i].to][k]);      }    for (int i=m;i>=0;i--)    {        if (i>=C[x]) f[x][i]=f[x][i-C[x]]+W[x];        else f[x][i]=0;    }}int main(){    n=read();m=read();    for (int i=1;i<=n;i++)c[i]=read();    for (int i=1;i<=n;i++)w[i]=read();    for (int i=1;i<=n;i++)    {        fa[i]=read();        if (fa[i])ins1(fa[i],i);    }    tarjan();    for (int i=1;i<=n;i++)    {      for (int j=head1[i];j;j=e1[j].next)    if (belong[e1[j].to]!=belong[i])ins2(belong[i],belong[e1[j].to]);    }    for (int i=1;i<=cnt3;i++)ins2(0,i);    treedp(0);    for (int i=0;i<=m;i++)ans=max(ans,f[0][i]);    printf("%d\n",ans);}

　　

bzoj2427 [HAOI2010]軟體安裝