bzoj 1565: [NOI2009]植物大战僵尸

题意：

在n*m的地图上放僵尸，在吃一个植物前必须先吃掉若干给定的植物，然后可以获得一个收益（可以使负的），问最大。

题解：

tarjan+网络流。
首先明显是一个最大权闭合子图。
那么就将正权连原点，负权连汇点，若x需要先吃y则x->y。
然后最小割。
但是有一个问题，就是假如一堆正权的点连成环那么就gg了。
所以先强连通，将所有成环的点的得分改为负无穷就好了。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;
const int inf=(1<<28);
struct node{
    int x,y,c,next,other;
}a[2110000];int len=0,last[1010];
int dfn[1010],low[1010],id=0,sta[1010],cnt=0,belong[1010],top=0,pe[1010];
bool v[1010];
bool map[35][35][35][35];
int n,m,sc[35][35],att[35][35];
int qs(int x,int y){return (x-1)*m+y;}
void insert(int x,int y)
{
    a[++len].y=y;
    a[len].next=last[x];last[x]=len;
}
int h[1010],s[1010],st,ed;
void ins(int x,int y,int c)
{
    int k1=++len;
    a[len].x=x;a[len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;
    int k2=++len;
    a[len].x=y;a[len].y=x;a[len].c=0;
    a[len].next=last[y];last[y]=len;
    a[k1].other=k2;a[k2].other=k1;
}
bool bt_h()
{
    memset(h,0,sizeof(h));
    int l=1,r=2;s[l]=st;h[st]=1;
    while(l!=r)
    {
        int x=s[l];
        for(int i=last[x];i;i=a[i].next)
        {
            int y=a[i].y;
            if(h[y]==0&&a[i].c>0) h[y]=h[x]+1,s[r++]=y;
        }
        l++;
    }
    return h[ed]!=0;
}
int findflow(int x,int f)
{
    if(x==ed) return f;
    int t,ans=0;
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        if(h[x]+1==h[y]&&a[i].c>0&&ans<f)
        {
            ans+=(t=findflow(y,min(a[i].c,f-ans)));
            a[i].c-=t;a[a[i].other].c+=t;
        }
    }
    if(ans==0) h[x]=0;
    return ans;
}
void dfs(int x)
{
    dfn[x]=low[x]=++id;
    v[x]=true;sta[++top]=x;
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        if(!dfn[y])
        {
            dfs(y);
            low[x]=min(low[x],low[y]);
        }
        else 
            if(v[y])
                low[x]=min(low[x],dfn[y]);
    }
    if(dfn[x]==low[x])
    {
        int i,num=0;cnt++;
        do
        {
            i=sta[top--];
            belong[i]=cnt;
            v[i]=false;
            num++;
        }while(i!=x);
        pe[cnt]=num;
    }
}
int main()
{
    scanf("%d %d",&n,&m);
    memset(map,false,sizeof(map));
    for(int x=1;x<=n;x++)
        for(int y=1;y<=m;y++)
        {
            scanf("%d %d",&sc[x][y],&att[x][y]);
            if(y!=m) map[x][y][x][y+1]=true;
            for(int i=1;i<=att[x][y];i++)
            {
                int xx,yy;scanf("%d %d",&xx,&yy);
                xx++;yy++;
                map[xx][yy][x][y]=true;
            }
        }
    for(int x=1;x<=n;x++)
    for(int y=1;y<=m;y++)
        for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(map[x][y][i][j]) insert(qs(x,y),qs(i,j));
    memset(v,false,sizeof(v));
    memset(dfn,0,sizeof(dfn));
    for(int i=1;i<=n*m;i++) if(!dfn[i]) dfs(i);
    len=0;memset(last,0,sizeof(last));st=0;ed=n*m+1;
    int ans=0;
    for(int x=1;x<=n;x++)
        for(int y=1;y<=m;y++)
        {
            if(pe[belong[qs(x,y)]]!=1) sc[x][y]=-inf;
            if(sc[x][y]>=0) ins(st,qs(x,y),sc[x][y]),ans+=sc[x][y];
            else ins(qs(x,y),ed,-sc[x][y]);
            for(int i=1;i<=n;i++)
                for(int j=1;j<=m;j++)
                    if(map[x][y][i][j]) ins(qs(x,y),qs(i,j),inf);
        }
    while(bt_h()) ans-=findflow(st,inf);
    printf("%d",ans);
}