luogu P2598 [ZJOI2009]狼和羊的故事 |网络流最小割

题目描述

“狼爱上羊啊爱的疯狂,谁让他们真爱了一场;狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想:狼和羊如此和谐,为什么不尝试羊狼合养呢?说干就干! Orez的羊狼圈可以看作一个n*m个矩阵格子,这个矩阵的边缘已经装上了篱笆。可是Drake很快发现狼再怎么也是狼,它们总是对羊垂涎三尺,那首歌只不过是一个动人的传说而已。所以Orez决定在羊狼圈中再加入一些篱笆,还是要将羊狼分开来养。 通过仔细观察,Orez发现狼和羊都有属于自己领地,若狼和羊们不能呆在自己的领地,那它们就会变得非常暴躁,不利于他们的成长。 Orez想要添加篱笆的尽可能的短。当然这个篱笆首先得保证不能改变狼羊的所属领地,再就是篱笆必须修筑完整,也就是说必须修建在单位格子的边界上并且不能只修建一部分。

输入格式

文件的第一行包含两个整数n和m。接下来n行每行m个整数,1表示该格子属于狼的领地,2表示属于羊的领地,0表示该格子不是任何一只动物的领地。

输出格式

文件中仅包含一个整数ans,代表篱笆的最短长度。


#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=4e4+5,M=6*N,inf=1<<30;
inline int read(){
	int x=0; char c=getchar();
	while(c<'0'||c>'9')c=getchar();
	while(c>='0'&&c<='9'){ x=(x<<1)+(x<<3)+(c^48); c=getchar(); }
	return x;
}
int nxt[M],head[N],go[M],edge[M],tot=1;
inline void add(int u,int v,int w){
	nxt[++tot]=head[u],head[u]=tot,go[tot]=v,edge[tot]=w;
	nxt[++tot]=head[v],head[v]=tot,go[tot]=u,edge[tot]=0;
}
int n,m,vis[105][105],s,t;
int d[N];
inline bool bfs(){
	memset(d,0,sizeof(d));
	queue<int>q; q.push(s); d[s]=1;
	while(q.size()){
		int u=q.front(); q.pop();
		for(int i=head[u];i;i=nxt[i]){
			int v=go[i];
			if(edge[i]&&!d[v]){
				d[v]=d[u]+1;
				q.push(v);
				if(v==t)return 1;	
			}
		}
	}
	return 0;
}
int dinic(int u,int flow){
	if(u==t)return flow;
	int rest=flow;
	for(int i=head[u];i&&rest;i=nxt[i]){
		int v=go[i];
		if(edge[i]&&d[v]==d[u]+1){
			int k=dinic(v,min(rest,edge[i]));
			if(!k)d[v]=-1;
			edge[i]-=k;
			edge[i^1]+=k;
			rest-=k;
		}
	}
	return flow-rest;
}

#define id(x,y) (x-1)*m+y
signed main(){
	n=read(),m=read();
	for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)vis[i][j]=read();	
	s=0,t=n*m*2+2;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++){
		if(vis[i][j]==1)add(s,id(i,j),inf);
		else if(vis[i][j]==2)add(id(i,j),t,inf);
		if(j>1)add(id(i,j),id(i,j-1),1);
		if(i>1)add(id(i,j),id(i-1,j),1);
		if(i<n)add(id(i,j),id(i+1,j),1);
		if(j<m)add(id(i,j),id(i,j+1),1);
	}
	int flow=0,maxflow=0;
	while(bfs())while(flow=dinic(s,inf))maxflow+=flow;
	cout<<maxflow<<endl;
}
原文地址:https://www.cnblogs.com/naruto-mzx/p/13061754.html