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;
}