BZOJ 1565 [NOI2009]植物大战僵尸

title: ‘BZOJ 1565 [NOI2009]植物大战僵尸’
categories: BZOJ
date: 2016-3-20 9:30:00
tags: [最小割]

---

Description

Input

Output

仅包含一个整数，表示可以获得的最大能源收入。注意，你也可以选择不进行任何攻击，这样能源收入为0。

Sample

input.txt
3 2
10 0
20 0
-10 0
-5 1 0 0
100 1 2 1
100 0

outut.txt
25

Solution

没刷过几道网络流的题= = 这题应该算比较好想的，原题貌似要输出方案，OJ比较良心只要输出最大收益就好了。可以看得出，当被保护的点连向保护它的点，左边的点连向右边的点后，明显是要求一个最大权闭合图。这个问题直接将源点连向正权点，负权点连向汇点，原图中的边连INF的边，求最小割就好了，割掉正权点的边就相当于不选，割掉负权点的边就相当于选了，并且一个割也对应一个闭合图，具体证明可以参照论文胡伯涛《最小割模型在信息学竞赛中的应用》。
然而这题还有一个要注意的地方，环上的点是不能选的，同时能走到环的点也不能选。这个可以把边反向建然后跑一遍拓扑，拓扑不到的点判掉就好了。

Code

（代码教丑，谨慎观看）

#include<queue>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>

#define maxs 1000+5
#define maxm 1000000+5 
#define INF 0x3f3f3f3f

using namespace std;

struct Graph{
    int u[maxm],v[maxm],cap[maxm],nxt[maxm],head[maxs];
    int ind;
    Graph(){ ind=0; memset(head,-1,sizeof(head)); }
    void addedge(int x,int y,int c){
        u[ind]=x; v[ind]=y; cap[ind]=c; nxt[ind]=head[x]; head[x]=ind++;
        u[ind]=y; v[ind]=x; cap[ind]=0; nxt[ind]=head[y]; head[y]=ind++;
    }
}G,Gr;

bool ok[maxs];
int sc[maxs],du[maxs],bow[maxs],stack[maxs],lv[maxs];
int n,m,s,t,ind,top,ans;

inline int pos(int i,int j){
    return (i-1)*m+j;
}


void Toposort(){
    for(int i=1;i<=n*m;i++)
        if(!du[i]) stack[++top]=i;
    while(top){
        int x=stack[top--]; ok[x]=1;
        for(int i=Gr.head[x];i!=-1;i=Gr.nxt[i])
            if(Gr.cap[i])
                if((--du[Gr.v[i]])==0) stack[++top]=Gr.v[i];          
    }
    for(int i=1;i<=n*m;i++)
        if(ok[i] && sc[i]>0) ans+=sc[i];
}

bool bfs(){
    queue <int> q;
    memset(lv,0,sizeof(lv));
    lv[s]=1; q.push(s);
    while(!q.empty()){
        int x=q.front(); q.pop();
        for(int i=G.head[x];i!=-1;i=G.nxt[i])
            if(ok[G.v[i]] && !lv[G.v[i]] && G.cap[i]){
                lv[G.v[i]]=lv[x]+1;
                q.push(G.v[i]);
            }
    }
    if(!lv[t]) return false;
    return true;

}

int dfs(int x,int t,int f){
    if(x==t) return f;
    for(int &i=bow[x];i!=-1;i=G.nxt[i])
        if(ok[G.v[i]] && lv[G.v[i]]>lv[x] && G.cap[i]){
            int fff=dfs(G.v[i],t,min(G.cap[i],f));
            if(fff){
                G.cap[i]-=fff;
                G.cap[i^1]+=fff;
                return fff;
            }
        }
    return 0;
}

void Dinic(){
    while(bfs()){
        memcpy(bow,G.head,sizeof(G.head));
        while(true){
            int fff=dfs(s,t,INF);
            if(fff) ans-=fff;
            else break;
        }
    }
}

int main(){
    scanf("%d%d",&n,&m); s=0; t=n*m+1; ok[s]=ok[t]=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            int sum;
            scanf("%d%d",&sc[pos(i,j)],&sum);
            if(sc[pos(i,j)]>0) G.addedge(s,pos(i,j),sc[pos(i,j)]);
            else if(sc[pos(i,j)]<0) G.addedge(pos(i,j),t,-sc[pos(i,j)]);    
            for(int k=1;k<=sum;k++){
                int x,y;
                scanf("%d%d",&x,&y); x++; y++;
                G.addedge(pos(x,y),pos(i,j),INF);
                Gr.addedge(pos(i,j),pos(x,y),INF);
                du[pos(x,y)]++;
            }
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<m;j++){
            G.addedge(pos(i,j),pos(i,j+1),INF);
            Gr.addedge(pos(i,j+1),pos(i,j),INF),du[pos(i,j)]++;
        }
    Toposort();
    Dinic();
    printf("%d",ans);
    return 0;
}