bzoj1565:[NOI2009]植物大战僵尸

传送门

对于一些植物存在保护它的植物，比如能攻击到它的和它右边的植物都能保护它
这就符合最大权闭合子图的定义
按照最大权闭合子图的建模方式，正权连源点，负权连汇点，中间就对于每个植物连向保护它的植物就好了
可惜有些保护是成环的，环内的植物以及被这个环保护的植物你一个也碰不了
所以先反向建图，跑一边拓扑排序，将环处理出来，建模的时候忽略掉这些点就好了
代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
#include<vector>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
void print(int x){
    if(x<0)putchar('-'),x=-x;
    if(x>9)print(x/10);
    putchar(x%10+'0');
}
#define rg register
const int maxn=1250,inf=1e9;queue<int>q;bool vis[maxn];vector<int>g[maxn];
int dis[maxn],in[maxn],cnt=1,n,m,ans,pre[maxn*maxn],sum,x[30][40],r[30][40][700],c[30][40][700];
int cur[maxn],s,t,a[21][31],nxt[maxn*maxn],h[maxn],v[maxn*maxn];
int id(int x,int y){return (x-1)*m+y;}
void add(int x,int y,int z)
{
	pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt,v[cnt]=z;
	pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt,v[cnt]=0;
}
bool bfs()
{
	memset(dis,0,sizeof dis);
	q.push(s),dis[s]=1;
	while(!q.empty())
	{
		int x=q.front();q.pop();
		for(rg int i=h[x];i;i=nxt[i])
			if(v[i]&&!dis[pre[i]]&&!in[pre[i]])dis[pre[i]]=dis[x]+1,q.push(pre[i]);
	}
	return dis[t];
}
int dfs(int x,int flow)
{
	if(x==t||!flow)return flow;
	int f=flow;
	for(rg int &i=cur[x];i;i=nxt[i])
		if(dis[pre[i]]>dis[x]&&v[i])
		{
			int y=dfs(pre[i],min(f,v[i]));
			f-=y,v[i]-=y,v[i^1]+=y;
			if(!f)return flow;
		}
	if(f==flow)dis[x]=-1;
	return flow-f;
}
int main()
{
	read(n),read(m);t=n*m+1;
	for(rg int i=1;i<=n;i++)
		for(rg int j=1;j<=m;j++)
		{
			read(a[i][j]),read(x[i][j]);
			for(rg int k=1;k<=x[i][j];k++)read(r[i][j][k]),read(c[i][j][k]),r[i][j][k]++,c[i][j][k]++,in[id(r[i][j][k],c[i][j][k])]++,g[id(i,j)].push_back(id(r[i][j][k],c[i][j][k]));
			if(j!=1)in[id(i,j-1)]++,g[id(i,j)].push_back(id(i,j-1));
		}
	queue<int>qq;
	for(rg int i=1;i<=n*m;i++)if(!in[i])qq.push(i),vis[i]=1;
	while(!qq.empty())
	{
		int x=qq.front(),w=g[x].size();qq.pop();
		for(rg int i=0;i<w;i++)if(!(--in[g[x][i]]))qq.push(g[x][i]),vis[g[x][i]]=1;
	}
	for(rg int i=1;i<=n;i++)
		for(rg int j=1;j<=m;j++)
			if(vis[id(i,j)])
			{
				if(a[i][j]>0)add(s,id(i,j),a[i][j]),sum+=a[i][j];if(j!=m)add(id(i,j),id(i,j+1),inf);
				for(rg int k=1;k<=x[i][j];k++)add(id(r[i][j][k],c[i][j][k]),id(i,j),inf);
				if(a[i][j]<0)add(id(i,j),t,-a[i][j]);
			}
	for(;bfs();ans+=dfs(s,inf))memcpy(cur,h,sizeof h);
	printf("%d\n",sum-ans);
}