[NOI 2009]植物大战僵尸

类型：最大权值闭合子图

建模分析：

每个植物有正价值和负价值 这就像太空飞行计划问题中的实验和器材一样
而每个植物保护和被保护的关系就是器材和实验的关系
那么我们就可以套用最大权值闭合子图的模型
首先 从源点向正价值植物建一条容量为价值的边
从每个负价值植物向汇点建一条容量为价值绝对值的边
然后从被保护的植物向保护的植物连一条容量为inf的边
需要注意的是，我们从右方进攻也是一种保护和被保护的关系
答案即为所有正价值植物价值总和-网络最大流
然而这样建完模型之后，样例都不对？？？？？？？
到底是哪里出错了？
观察一下样例，我发现了一个问题，如果一个植物保护了它右边的植物。
那么这两个植物直接便形成了环.
我们无论如何都吃不掉这些植物的.
我们跑最大流建立的模型必须是DAG.
于是我们得先把图中的环去掉.
拓扑排序以及Tarjan都可以解决这个问题.
拓扑好写而且是队列实现，于是就用了拓扑排序.
去掉环之后重建图即可.

易错点：

1.边的数量比较多 足足有$8*(10^5)$
2.价值总和是在重建图中正价值的价值总和
3.我们可以不进攻，也就是说当我们得到的答案是负数的时候，其实答案是0，不要忘记特判

题目链接：

COGS
Luogu
BZOJ
这次BZOJ好评，前两个OJ给的空间太大了，我MLE没检测出来

Ac Code：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <queue>
#include <map>
#include <cstring>
#define il inline
#define ll long long
using namespace std;
const int inf=0x7fffffff;
const int maxm=220000;
int head[maxm],to[maxm*2],cap[maxm*2],net[maxm*2],deep[800],cnt=1;
il void add(int x,int y,int c){cnt++,to[cnt]=y,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
queue <int> dl;
int vis[maxm],cx[maxm],id[60][60],r[maxm];
int n,m,s,t,sum=0;
il bool BFS(int s,int t)
{
    while(!dl.empty()) dl.pop();
    memset(deep,-1,sizeof(deep));
    dl.push(s),deep[s]=0;
    while(!dl.empty())
    {
        int x=dl.front();dl.pop();
        for(int i=head[x];i;i=net[i])
         if(cap[i]>0&&deep[to[i]]==-1)
          dl.push(to[i]),deep[to[i]]=deep[x]+1;
    }
    return deep[t]==-1?0:1;
}
int dfs(int now,int flow,int t)
{
    if(now==t) return flow;
    int w,used=0;
    for(int i=head[now];i;i=net[i])
    {
        int v=to[i];
        if(deep[v]==deep[now]+1&&cap[i])
        {
            w=dfs(v,min(flow-used,cap[i]),t);
            cap[i]-=w;
            cap[i^1]+=w;
            used+=w;
            if(used==flow) return flow;
        }
    }
    if(!used) deep[now]=-1;
    return used;
}
il int dinic(int s,int t)
{
    int maxflow=0;
    while(BFS(s,t)) maxflow+=dfs(s,inf,t);
    return maxflow;
}
inline void adx(int x,int y,int cax)
{
    add(x,y,cax),add(y,x,0);
}
int headx[maxm],tox[maxm*2],netx[maxm*2],cntx;
il void addx(int x,int y){cntx++,tox[cntx]=y,netx[cntx]=headx[x],headx[x]=cntx;}
void build()
{
    scanf("%d%d",&n,&m);
    int bx=0;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      id[i][j]=++bx;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
     {
        int w;
        scanf("%d%d",&cx[id[i][j]],&w);
        for(int k=1,x,y;k<=w;k++)
        {
            scanf("%d%d",&x,&y);
            x++,y++;
            addx(id[i][j],id[x][y]);
            r[id[x][y]]++;
        }
        if(j!=m) addx(id[i][j+1],id[i][j]),r[id[i][j]]++;
     }
}
void TO()
{
    while(!dl.empty()) dl.pop();
    for(int i=1;i<=n*m;i++)
     if(!r[i]) dl.push(i),vis[i]=1;
    while(!dl.empty())
    {
        int x=dl.front();
        dl.pop();
        for(int i=headx[x];i;i=netx[i])
        {
            int tmp=tox[i];
            if(vis[tmp]) continue;
            r[tmp]--;
            if(r[tmp]==0) dl.push(tmp),vis[tmp]=1;
        }
    }
}
void rebuild()
{
    s=n*m+5,t=n*m+6;
    for(int i=1;i<=n*m;i++)
    if(vis[i])
    {
        cx[i]<0?adx(i,t,-cx[i]):adx(s,i,cx[i]);
        if(cx[i]>0) sum+=cx[i];
        for(int j=headx[i];j;j=netx[j])
        {
            int tmp=tox[j];
            if(!vis[tmp]) continue;
            adx(tmp,i,inf);
        }
    } 
}
int main()
{
   //freopen("pvz.in","r",stdin);
   //freopen("pvz.out","w",stdout);
   build();//建图
   TO();//拓扑排序
   rebuild();//重建
   return printf("%d\n",max(sum-dinic(s,t),0))*0;
}