2427: [HAOI2010]軟體安裝

                                                                    Time Limit: 10 Sec  Memory Limit: 128 MB                                                                                  Submit: 2330  Solved: 978

Description

現在我們的手頭有N個軟體，對於一個軟體

Input

第1行：N, M

Output

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

Sample Input

3 10 5 5 6 2 3 4 0 1 1

Sample Output

5

解析：

       注意圖中是存在環的！所以不能直接進行樹形DP，而是先用Tarjan縮點重建圖再進行DP。

程式碼：

#include <bits/stdc++.h>
using namespace std;

const int Max=505;
int n,m,size,cnt,Index,tot,s;
int f[Max][Max],first[Max],First[Max],c[Max],w[Max],to[Max];
int low[Max],num[Max],father[Max],vis[Max],p[Max],sumw[Max],sumc[Max];
struct shu{int to,next;};
shu edge[Max],Edge[Max];

inline int get_int()
{
	int x=0,f=1;
	char c;
	for(c=getchar();(!isdigit(c))&&(c!='-');c=getchar());
	if(c=='-') f=-1,c=getchar();
	for(;isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+c-'0';
	return x*f;
}

inline void Build(int x,int y)
{
	Edge[++size].next=First[x];
	First[x]=size;
	Edge[size].to=y;
}

inline void build(int x,int y)
{
	edge[++size].next=first[x];
	first[x]=size;
	edge[size].to=y;
}

inline void tarjan(int point)
{
	low[point]=num[point]=++Index;
	vis[point]=1,p[++tot]=point;
	for(int u=First[point];u;u=Edge[u].next)
	{
	  int to=Edge[u].to;
	  if(!num[to]) tarjan(to),low[point]=min(low[point],low[to]);
	  else if(vis[to]) low[point]=min(low[point],num[to]);
	}
	if(low[point]==num[point])
	{
	  cnt++;
	  while(1)
	  {
	  	int x=p[tot--];
	  	vis[x]=0,father[x]=cnt;
	  	sumw[cnt]+=w[x];
	  	sumc[cnt]+=c[x];
	  	if(x==point) break;
	  }
	}
}

inline void dfs(int point)
{
	vis[point]=1,f[point][0]=0;
	for(int u=first[point];u;u=edge[u].next)
	{
	  int to=edge[u].to;
	  if(vis[to]) continue;
	  dfs(to);
	  for(int i=m;i>=0;i--)
	    for(int j=0;j<=i;j++)
	      f[point][i]=max(f[point][i],f[point][i-j]+f[to][j]);
	}
	for(int i=m;i>=0;i--)   //根據題意強制選
	{
	  if(i>=sumw[point]) f[point][i]=f[point][i-sumw[point]]+sumc[point];
	  else f[point][i]=-1e9;  //注意！
	} f[point][sumw[point]]=sumc[point];
}

int main()
{
	n=get_int(),m=get_int();
	for(int i=1;i<=n;i++) w[i]=get_int();
	for(int i=1;i<=n;i++) c[i]=get_int();
	for(int i=1;i<=n;i++)
	{
	  int x=get_int();
	  if(i==x) continue;
	  if(x) Build(x,i);
	}
	size=0;
	for(int i=1;i<=n;i++) if(!num[i]) tarjan(i);
	for(int i=1;i<=n;i++)
	  for(int u=First[i];u;u=Edge[u].next)
        if(father[i]!=father[Edge[u].to]) build(father[i],father[Edge[u].to]),to[father[Edge[u].to]]++;
    for(int i=1;i<=cnt;i++) if(!to[i]) build(0,i);
	memset(vis,0,sizeof(vis));
	dfs(0);
	cout<<f[0][m];
	return 0;
}