省选专练NOI2009植物大战僵尸

#include<iostream>

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
struct Front_star{
	int u,v,w,nxt;
}e[400100];
int cnt=1;
int first[1001]={0};
void addedge(int u,int v,int w){
	cnt++;
	e[cnt].u=u;
	e[cnt].v=v;
	e[cnt].w=w;
	e[cnt].nxt=first[u];
	first[u]=cnt;
} 
int n,m;
int S,T;
int mp[5001]={0};
int rd[5001]={0};
int cur[5001]={0};//当前弧优化，准确说就是不走边权为0的边了。 
void add(int u,int v,int w){
	addedge(u,v,w);
	addedge(v,u,0);
	rd[u]++;
}
int d[5001]={0};
queue<int> q;
int bfs(){
	for(int i=S;i<=T;i++)
		if(d[i]!=-2)
			d[i]=-1;
	 q.push(S);
	 d[S]=1;
	 while(!q.empty()){
	 	int x=q.front();
	 	q.pop();
	 	for(int i=first[x];i>1;i=e[i].nxt){
	 		if(e[i].w>0&&d[e[i].v]==-1){
	 			d[e[i].v]=d[x]+1;
				 q.push(e[i].v); 
	 		}
	 	}
	 }
	 return d[T]>0;
}
int dfs(int st,int ed,int nowdat){
	int dat=0;
	if(st==ed)
		return nowdat;
	for(int i=cur[st];i>1&&nowdat>dat;i=e[i].nxt){
		cur[st]=i;
		int v=e[i].v;
		if(e[i].w>0&&d[st]+1==d[e[i].v]){
			int now=min(nowdat-dat,e[i].w);
			now=dfs(v,ed,now);
			dat+=now;
			e[i].w-=now;
			e[i^1].w+=now;
			if(e[i].w){
				cur[st]=i;
			}
		}
	}
	if(!dat)
		d[st]=-3;
	return dat;
}
int dinic(){
	int ans=0;
	while(bfs()){
		for(int i=S;i<=T;i++){
			cur[i]=first[i];
		}
		ans+=dfs(S,T,0x3f3f3f3f);
	}
	return ans;
}
int main(){
	scanf("%d%d",&n,&m);
	S=0;
	T=n*m+1;
	for(int i=1;i<=n*m;i++){
		scanf("%d",&mp[i]);
		if(mp[i]>0)
			add(S,i,mp[i]);
		else
			add(i,T,-mp[i]);
//		cout<<"$%#$%"<<endl;
		int tot;
		scanf("%d",&tot);
		for(int j=1;j<=tot;j++){
			int dx,dy;
			scanf("%d%d",&dx,&dy);
			add(dx*m+dy+1,i,INF);
		}
		if(i%m)//不换行
			add(i,i+1,INF); 
	}
	int q[5001]={0};//Topsort去除环
	int he=0,ta=0;
	for(int i=S;i<=T;i++){
		if(!rd[i]){
			q[ta++]=i;
		}
		d[i]=-2;
	}
	int ans=0;
	while(he<ta){
		int x=q[he++];
		d[x]=0;
		if(mp[x]>0)//对无前效应的点，贪心得出可获得的值。
			 ans+=mp[x];
		for(int i=first[x];i;i=e[i].nxt){
			if(i%2==1){//如果不是反向边 
				if(!--rd[e[i].v])
					q[ta++]=e[i].v;
			}
		}
	}
//	cout<<ans<<endl;
	ans=ans-dinic();
	cout<<ans;
}

典型的最大权闭合图问题。

注意不能用tarjan缩点

先Topsort 缩一个点，再跑dinic 完