[BZOJ]1565: [NOI2009]植物大战僵尸 Tarjan+最小割（最大权闭合子图）

Description

题解：

经典的最大权闭合子图模型，不会的可以去网上查找相关资料，挺简单的。这道题呢，对于植物i可以保护植物j，我们可以建一条j->i的有向边（注意同一行相邻2个植物也有保护关系），然后跑一次Tarjan，因为在一个环中的植物我们是不能选到的，对于这些不能选到的植物，我的处理方法是将他们的权值看做负无穷，然后其他的按照最大权闭合子图搞就行了。

代码：

#include<cstdio>#include<cstring>#include<queue>#include<iostream>#include<algorithm>using namespace std;const int inf=2147483647;const int Maxn=35;const int Max=910;struct Edge1{int y,next;}E[400000];struct Edge2{int y,d,next;}e[2110000];int score1[Maxn][Maxn],score2[Max],num[Maxn][Maxn],z=0;int n,m;int last[Max],len=1;int Last[Max],Len=0;int dfn[Max],low[Max],id=0,bel[Max],cnt=0,sta[Max],top=0,p[Max];bool in[Max],vis[Max],hh[Max][Max];void ins(int x,int y,int d){    int t=++len;    e[t].y=y;e[t].d=d;    e[t].next=last[x];last[x]=t;}void addedge(int x,int y,int d){ins(x,y,d);ins(y,x,0);}void Ins(int x,int y){    int t=++Len;    E[t].y=y;E[t].next=Last[x];Last[x]=t;}void Tarjan(int x){    low[x]=dfn[x]=++id;    sta[++top]=x;in[x]=true;    for(int i=Last[x];i;i=E[i].next)    {        int y=E[i].y;        if(!dfn[y])Tarjan(y),low[x]=min(low[x],low[y]);        else if(in[y])low[x]=min(low[x],dfn[y]);    }    if(low[x]==dfn[x])    {        int i;cnt++;        do        {            i=sta[top--];            bel[i]=cnt;            p[cnt]++;            in[i]=false;        }while(i!=x);    }}queue<int>q;int h[Max],st,ed;bool bfs(){    memset(h,0,sizeof(h));h[st]=1;    q.push(st);    while(!q.empty())    {        int x=q.front();q.pop();        for(int i=last[x];i;i=e[i].next)        if(e[i].d&&!h[e[i].y])h[e[i].y]=h[x]+1,q.push(e[i].y);    }    return h[ed];}int dfs(int x,int f){    if(x==ed)return f;    int s=0,t;    for(int i=last[x];i;i=e[i].next)    {        int y=e[i].y;        if(h[y]==h[x]+1&&e[i].d&&s<f)        {            t=dfs(y,min(f-s,e[i].d));            s+=t;e[i^1].d+=t;e[i].d-=t;        }    }    if(s==0)h[x]=0;    return s;}int main(){    int ans=0;    scanf("%d%d",&n,&m);    for(int i=1;i<=n;i++)    for(int j=1;j<=m;j++)    num[i][j]=++z;    for(int i=1;i<=n;i++)    for(int j=1;j<=m;j++)    {        int k;        scanf("%d%d",&score1[i][j],&k);        for(int l=1;l<=k;l++)        {            int x,y;            scanf("%d%d",&x,&y);x++;y++;            Ins(num[x][y],num[i][j]);            hh[num[x][y]][num[i][j]]=true;        }        if(j>1&&!hh[num[i][j-1]][num[i][j]])Ins(num[i][j-1],num[i][j]);    }    for(int i=1;i<=z;i++)if(!dfn[i])Tarjan(i);    st=z+1;ed=z+2;    for(int i=1;i<=n;i++)    for(int j=1;j<=m;j++)    {        int x=num[i][j];        if(p[bel[x]]>1)addedge(x,ed,inf);        else if(score1[i][j]<0)addedge(x,ed,-score1[i][j]);        else ans+=score1[i][j],addedge(st,x,score1[i][j]);        memset(vis,false,sizeof(vis));vis[x]=true;        for(int l=Last[x];l;l=E[l].next)        {            int y=E[l].y;            if(vis[y])break;vis[y]=true;            addedge(x,y,inf);        }    }    while(bfs())ans-=dfs(st,inf);    printf("%d",ans);}