HihoCoder - 1398 - 网络流五·最大权闭合子图

题目：HihoCoder - 1398

题意：题目说的很清楚了，题解也给了，就是最大权闭合子图，也有证明（虽然我看不懂）

题解：

就讲一下最大权闭合子图的做法吧：

引入闭合子图的概念 : 通俗点说就是选出一个图的子图，使得子图中的所有点出度指向的点依旧在这个子图内，则说明此子图是闭合子图。

最大权闭合子图 : 假设每个点具有点权值，在一个图的所有闭合子图中，点权之和最大的即是最大权闭合子图。

求取最大权闭合子图的权值之和是有一个结论的

一、先抽象出一个超级源、汇点

二、将权值为正的点和超级源点连接、容量为权值

三、将权值为负的点和超级汇点连接、容量为权值的绝对值

四、然后除了源、汇之外的点原本怎么连就怎么连、且容量为无穷大

五、最大权闭合子图权值  =  所有权值为正的权值之和  -  最大流

按照以上方法建图，然后求出最大流（最小割）就差不多好啦！

代码：

#include<bits/stdc++.h>
#define N 100005
#define INF 0x7f7f7f7f
using namespace std;
int n,m,s,t,c,sum;
int head[N],cur[N],deep[N];
struct ljh
{
	int to,next,w;
}e[2*N];
inline void add(int x,int y,int z)
{
	e[c].next=head[x];
	e[c].w=z;
	e[c].to=y;
	head[x]=c++;
}
bool bfs(int s,int t)
{
	memset(deep,0,sizeof(deep));
	deep[s]=1;
	queue<int>q;
	q.push(s);
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		for(int i=head[x];i!=-1;i=e[i].next)
		{
			int nex=e[i].to;
			if(!deep[nex]&&e[i].w>0)
			{
				deep[nex]=deep[x]+1;
				q.push(nex);			
			}
		}
	}
	return deep[t];
}
int dfs(int x,int t,int maxflow)
{
	if(x==t)return maxflow;
	int ans=0;
	for(int i=cur[x];i!=-1;i=e[i].next)
	{
		int nex=e[i].to;
		if(e[i].w<=0||ans>=maxflow||deep[nex]!=deep[x]+1)continue;
		cur[x]=i;
		int k=dfs(nex,t,min(e[i].w,maxflow-ans));
		e[i].w-=k;
		e[i^1].w+=k;
		ans+=k;
	}
	if(ans==0)deep[x]=-2;
	return ans;
}
int Dinic(int s,int t)
{
	int ans=0;
	while(bfs(s,t))
	{
		memcpy(cur,head,sizeof(head));
		ans+=dfs(s,t,INF);
	}
	return ans;
}
int main()
{
	while(~scanf("%d%d",&n,&m))
	{
		int x,y,z;
		s=0;t=n+m+1;
		c=0;
		sum=0;
		memset(head,-1,sizeof(head));
		for(int i=1;i<=m;i++)
		{
			scanf("%d",&x);
			add(n+i,t,x);
			add(t,n+i,0);
		}
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&x,&y);
			sum+=x;
			add(s,i,x);
			add(i,s,0);
			for(int j=1;j<=y;j++)
			{
				scanf("%d",&z);
				add(i,n+z,INF);
				add(n+z,i,0);
			}
		}
		printf("%d\n",sum-Dinic(s,t));
	}
}