#最大流，拓扑排序，最大权闭合子图#JZOJ 2856 洛谷 2805 植物大战僵尸

题目链接

---

题目大意

给一个图，选择一个点必须选择该点的后继，问可以获得的最大点权和

---

分析

其实之前我就觉得这道题目非常难，首先植物大战僵尸这么名字非常熟悉，而且题目又丑又长，所以说也不敢做，现在？？？
首先这道题目可以转换成最大权闭合子图这个东东，然后源点连正权点，汇点连负权点，原来的边流量inf，那么就可以求出答案。
恭喜你做错了，没错，还要来个拓扑排序把环去掉

---

代码

#include <cstdio>
#include <cctype>
#include <queue>
#define rr register
using namespace std;
struct node{int y,w,next;}e[300001];
queue<int>q,tq[601]; int dis[601],ls[601],a[601],p[601],n,m,tot,s,t,ans,k=1;
inline signed iut(){
    rr int ans=0,f=1; rr char c=getchar();
    while (!isdigit(c)) f=(c=='-')?-f:f,c=getchar();
    while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
    return ans*f;
}
inline void add(int x,int y,int w){
	e[++k]=(node){y,w,ls[x]}; ls[x]=k;
	e[++k]=(node){x,0,ls[y]}; ls[y]=k;
}
inline signed bfs(int s){
    for (rr int i=1;i<=t;++i) dis[i]=0;
    queue<int>q; q.push(s); dis[s]=1;
    while (q.size()){
        rr int x=q.front(); q.pop();
        for (rr int i=ls[x];i;i=e[i].next)
        if (e[i].w>0&&!dis[e[i].y]){
            dis[e[i].y]=dis[x]+1;
            if (e[i].y==t) return 1;
            q.push(e[i].y);
        }
    }
    return 0;
}
inline signed dfs(int x,int now){
    if (x==t||!now) return now;
    rr int rest=0,f;
    for (rr int i=ls[x];i;i=e[i].next)
    if (e[i].w>0&&dis[e[i].y]==dis[x]+1){
        rest+=(f=dfs(e[i].y,min(now-rest,e[i].w)));
        e[i].w-=f; e[i^1].w+=f;
        if (now==rest) return rest;
    }
    if (!rest) dis[x]=0;
    return rest;
}
signed main(){
    n=iut(); m=iut();
    for (rr int i=1;i<=n;++i)
    for (rr int j=1;j<=m;++j){
        a[(i-1)*m+j]=iut();
        for (rr int t=iut();t;--t){
        	rr int x=iut()+1,y=iut()+1;
        	++ls[(x-1)*m+y],tq[(i-1)*m+j].push((x-1)*m+y);
		}
		if (j<m) ++ls[(i-1)*m+j],tq[(i-1)*m+j+1].push((i-1)*m+j);
    }
    for (rr int i=1;i<=n*m;++i) if (!ls[i]) q.push(i);
    while (q.size()){
    	rr int x=q.front(),cnt=tq[x].size(); q.pop();
    	for (rr int i=1;i<=cnt;++i){
    		rr int y=tq[x].front(); tq[x].pop(); tq[x].push(y);
    		if (!(--ls[y])) q.push(y);
		}
	}
	for (rr int i=1;i<=n*m;++i) if (!ls[i]) p[i]=++tot; else ls[i]=0;
	s=tot+1,t=tot+2;
	for (rr int i=1;i<=n*m;++i)
	if (p[i]){
		if (a[i]<0) add(p[i],t,-a[i]);
		    else add(s,p[i],a[i]),ans+=a[i];
		while (tq[i].size()){
			rr int y=tq[i].front(); tq[i].pop();
			if (p[y]) add(p[y],p[i],23333333);
		}
	}
	while (bfs(s)) ans-=dfs(s,23333333);
	return !printf("%d",ans);
}