【BZOJ】1565: [NOI2009]植物大战僵尸 -最大权闭合图&最小割

传送门：bzoj1565

感觉很多题解都没说清楚啊(或者是我理解不了)QAQ

---

题解

引用一下将狼踩尽 19891101的题解的图。

最大权闭合图
设一个有向图(不一定是$DAG$)中，点集为$V$，每个点都有权值$a_i$(可正可负)，现在要求选出一个总权值最大的子图(点集的子集)满足：对于所有有向边集$E=(u,v)$如果$u$在子图上$v$就必须在。

一个例子：

我们考虑转成网络流模型：建立源点$S$和汇点$T$，从$S$向所有权值为正的点连一条流量为$a_i$的边，每个权值为负的点向$T$连一条流量为$-a_i$的边。而原图中的每条边保持不变，但流量设为$+\infty$。

得到了这样一张图。可以证明，任意简单割分割开后$S$所在的点集是一个闭合子图。(简单割，割集中任意边都与$S$或$T$相连)

因为一个割集中的边必然不包含原图中的边，而选择了某个点，它所指向的所有点必然也和它在同一个点集中，这也就满足了闭合子图的要求。
设所有权值为正的点权值之和为$sum$，割边边权之和为$cost$，这样的一个闭合子图的总权值即为$sum-cost$。
1.所有被割掉的权值为正的点显然就不在这个子图里了，所以要减去它们的权值。
2.因为割集的存在使得$S,T$不连通，所以所有权值为负的点连向$T$的边若没有被割掉，则说明它不在子图中，反之则在子图中。
所以$sum-cost$相当于除开所有未选的权值为正的点后，减去了所选的权值为负的点值。

求最大权闭合图:$sum$不变，所以使$cost$尽量小，那么就是最小割(最大流)。

而按照这道题的要求建图可能存在环，环上的点都不能选，需要拓扑排序判一下。

---

代码

#include<bits/stdc++.h>
#define trs(r,c) ((r)*m+(c))
using namespace std;
const int N=1000,M=5e5,inf=0x3f3f3f3f;

int n,m,val[N],ind[N],vs[N],ans,dep[N];
int head[N],to[M],w[M],nxt[M],cot,tot,S,T;

struct L{int u,v;}le[M];
queue<int>que;

inline void lk(int u,int v)
{le[++cot]=(L){u,v};nxt[cot]=head[u];head[u]=cot;ind[v]++;}
inline void lk(int u,int v,int vv)
{
   to[++tot]=v;nxt[tot]=head[u];head[u]=tot;w[tot]=vv;
   to[++tot]=u;nxt[tot]=head[v];head[v]=tot;w[tot]=0;
}

inline bool bfs()
{
	memset(dep,0xff,sizeof(dep));dep[S]=1;que.push(S);
	int i,j,x;
	for(;que.size();){
		x=que.front();que.pop();
		for(i=head[x];i;i=nxt[i]){
			j=to[i];if((!w[i]) || (~dep[j])) continue;
			dep[j]=dep[x]+1;que.push(j);
		}
	}
	return (dep[T]!=-1);
}

int dfs(int x,int f)
{
	if(x==T) return f;
	int i,j,res,ss=0;
	for(i=head[x];i;i=nxt[i]){
		j=to[i];if((!w[i])||(dep[j]!=dep[x]+1)) continue;
		res=dfs(j,min(f-ss,w[i]));if(!res) continue;
		w[i]-=res;w[i^1]+=res;ss+=res;if(ss==f) return f;
	}
	if(!ss) dep[x]=-1;return ss;
}

int main(){
	int i,j,num,x,y,bs=0;
	scanf("%d%d",&n,&m);S=n*m;T=S+1;
	for(i=0;i<n;++i)
	  for(j=0;j<m;++j){
	     scanf("%d%d",&val[bs],&num);
	     if(j+1<m) lk(bs+1,bs);
	     for(;num;--num) {
	        scanf("%d%d",&x,&y);lk(bs,trs(x,y));
		 }
	     bs++;
	  }
	for(i=0;i<S;++i) if(!ind[i]) que.push(i);
	for(;que.size();){
		x=que.front();que.pop();vs[x]=1;
		for(i=head[x];i;i=nxt[i]){
			j=le[i].v;ind[j]--;
			if(!ind[j]) que.push(j);
		}
	}
	memset(head,0,sizeof(head));tot=1;
	for(i=1;i<=cot;++i) if(vs[le[i].u]&&vs[le[i].v]) lk(le[i].v,le[i].u,inf);
	for(i=0;i<S;++i) if(vs[i]) {if(val[i]>0) lk(S,i,val[i]),ans+=val[i];else lk(i,T,-val[i]);}
	for(;bfs();) ans-=dfs(S,inf);
	printf("%d",ans);
	return 0;
}