【NOI2009】植物大战僵尸【最大权闭合子图】【网络流】

传送门

不晓得为什么最后一个点T了，，

首先学习最大权闭合子图。

最大权闭合子图

在有点权的有向图里面，选择一些点出来，这些点的出边都指向这些点。不指向别人。这是一个闭合子图。

点权和最大，最大权闭合子图。

求法：

我们设S连向所有点权为正的点，边权为点权。负的连向T，边权为绝对值。其它边为inf。所以最小割一定在那些特殊边上。

正点权和-最小割（即最大流）即最大权闭合子图的权值。

反正割掉正的，就是减小贡献。割掉负的，就是选上负的，同样减小贡献。（瞎扯，不会）

总之就是正和-最小割。

本题

然后这道题很符合模型的是，所有能保护这个植物的植物都要被消灭，然后这个植物右边的植物也要被消灭，这个植物才能被消灭。

如果不是闭合子图，那伸出去的那个植物也要被消灭。所有又是一个闭合子图。

将保护关系反向连边，走模板。

注意可能有环，先正向拓扑查一下，有标记的才能用。

#include<bits/stdc++.h>
using namespace std;
#define in read()
int in{
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){
		ch=getchar();if(ch=='-')f=-1;
	}
	while(isdigit(ch)){
		cnt=cnt*10+ch-48;
		ch=getchar();
	}return cnt*f;
}int sum;
const int inf=0x3f3f3f3f;
int first[1003],nxt[500003],to[500003],tot,rd[1003],w[500003];
void add(int a,int b){
	//cout<<"add "<<a<<" "<<b<<endl;
	rd[b]++;nxt[++tot]=first[a];first[a]=tot;to[tot]=b;
}
void add(int a,int b,int c){
	//cout<<" "<<a<<" "<<b<<endl;
	nxt[++tot]=first[a];first[a]=tot;to[tot]=b;w[tot]=c;
	nxt[++tot]=first[b];first[b]=tot;to[tot]=a;w[tot]=0;
}
struct node{
	int u,v;
}edge[200003];int cnt;
int n,m;
int get(int x,int y){
	return x*31+y;
}int flag[1003];
int score[1003];queue<int> q;
void topo(){
	for(register int i=1;i<=n;++i)for(register int j=1;j<=m;++j)if(!rd[get(i,j)])q.push(get(i,j));
	while(!q.empty()){
		int u=q.front();q.pop();flag[u]=1;
		for(register int i=first[u];i;i=nxt[i]){
			int v=to[i];--rd[v];if(!rd[v])q.push(v);
		}
	}
}
int s=0,t;
int maxflow,dep[1003],vis[1003];
bool bfs(int S,int T){
	memset(dep,0,sizeof(dep));
	while(!q.empty())q.pop();
	q.push(S);dep[S]=1;
	while(!q.empty()){
		int u=q.front();q.pop();
		for(register int i=first[u];i;i=nxt[i]){
			int v=to[i];
			if(w[i]&&(!dep[v])){
				dep[v]=dep[u]+1;q.push(v);
			}
		}
	}return dep[T]!=0;
}
int dfs(int now,int T,int limit){
	//cout<<now<<" "<<T<<endl;
	if(now==T||!limit)return limit;
	int flow=0,f;
	for(register int i=first[now];i;i=nxt[i]){
		int v=to[i];
		if(dep[v]==dep[now]+1&&(f=dfs(v,T,min(limit,w[i])))){
			flow+=f;w[i]-=f;w[i^1]+=f;limit-=f;if(!limit)break;
		}
	}return flow;
}
void dinic(int S,int T){
	while(bfs(S,T)){
		maxflow+=dfs(S,T,0x3f3f3f3f);
	}
}
signed main(){
	//freopen("1.in","r",stdin);
	n=in;m=in;t=get(n,m)+33;
	for(register int i=1;i<=n;++i){
		for(register int j=1;j<=m;++j){
			score[get(i,j)]=in;int x=in;
			for(register int k=1;k<=x;k++){
				int xx=in+1;int yy=in+1;edge[++cnt]={get(i,j),get(xx,yy)};
			}
		}
	}
	for(register int i=1;i<=n;++i)
		for(register int j=1;j<m;++j)
			edge[++cnt]={get(i,j+1),get(i,j)};
	for(register int i=1;i<=cnt;++i)add(edge[i].u,edge[i].v);
	//for(register int i=1;i<=n;++i)for(register int j=1;j<=m;++j)cout<<rd[get(i,j)]<<" ";cout<<endl;
	topo(); //for(register int i=1;i<=n;++i)for(register int j=1;j<=m;++j)cout<<flag[get(i,j)]<<" ";cout<<endl;
	memset(first,0,sizeof(first));tot=1;
	for(register int i=1;i<=cnt;++i)
		if(flag[edge[i].u]&&flag[edge[i].v])add(edge[i].v,edge[i].u,inf);
//	for(register int i=1;i<=n;++i)
//		for(register int j=1;j<m;++j)
//			if(flag[get(i,j)]&&flag[get(i,j+1)])
//				add(get(i,j),get(i,j+1),inf);
	for(register int i=1;i<=n;++i)
		for(register int j=1;j<=m;++j)
			if(flag[get(i,j)])
				if(score[get(i,j)]>=0)add(s,get(i,j),score[get(i,j)]),sum+=score[get(i,j)];
				else add(get(i,j),t,-score[get(i,j)]);
	dinic(s,t);//cout<<sum<<" "<<maxflow<<endl;
	cout<<sum-maxflow;
	return 0;
}