zoj 2788 网络流sap+最小割

建图:(另设源点s)

如果某房间i有入侵者,则建立一条<s,i>,权值为无限大。

如果某房间i有控制面板,控制进入j房间的门,则建立两条边<i,j>,<j,i>,权值分别为无限大和1

建图完毕就可以求网络流解决这个最小割问题了。

(为什么又是最小割- - 。题目的本质为将源点和汇点分开- - 。所以属于最小割问题)

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
#include <vector>
#include <cstdio>
#include <string>
#include<string.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
#define INF 999999999;
#define set(x,i) memset(x,i,sizeof(x));
const int maxn=10000;
struct edge
{
	int to,cost,next;
}e[maxn];
int head[maxn],cnt;
int pre[maxn];
int gap[maxn];
int cur[maxn];
int level[maxn];
int s,t,tot;
template<class T>
inline void read(T &m)
{
	T x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	m=x*f;
}
void addedge(int a,int b,int c)
{
	e[cnt].to=b;
	e[cnt].cost=c;
	e[cnt].next=head[a];
	head[a]=cnt++;

	e[cnt].to=a;
	e[cnt].cost=0;
	e[cnt].next=head[b];
	head[b]=cnt++;
}
void build_graph(int n,int m)  //题目的边序号是从0-(n-1)
{
	s=0,t=m+1,tot=n+1,cnt=0;
	set(head,-1);
	set(gap,0);
	set(level,0);
	gap[0]=tot;
	for(int i=1;i<=n;i++)
	{
		char ch[5];
		scanf("%s",ch);
		if(ch[0]=='I')
			addedge(s,i,999999);//与源点相连的边<s,i>
		int k,kk;
		read(k);
		while(k--)
		{
			read(kk);
			addedge(i,kk+1,99999);  //边<i,j>
			addedge(kk+1,i,1);       //边<j,i>
		}
	}
}
int sap()
{
	for(int i=0;i<tot;i++)
		cur[i]=head[i];
	int u=s,ans=0,inf=INF;
	while(level[s]<tot)
	{
		bool flag=0;
		for(int &i=cur[u];i!=-1;i=e[i].next)
		{
			
			int v=e[i].to;
			if(e[i].cost>0&&level[u]==level[v]+1)
			{
				flag=1;
				pre[v]=u;
				u=v;
				if(e[i].cost<inf)
					inf=e[i].cost;
				if(u==t)
				{
					while(u!=s)
					{
						u=pre[u];
						e[cur[u]].cost-=inf;
						e[cur[u]^1].cost+=inf;
					}
					ans+=inf;
					inf=INF;
				}
				break;
			}
		}
		if(!flag)
		{
			if(--gap[level[u]]==0)
				break;
			int minn=tot;
			for(int i=head[u];i!=-1;i=e[i].next)
				if(e[i].cost>0&&level[e[i].to]<minn)
				{
					minn=level[e[i].to];
					cur[u]=i;
				}
			level[u]=minn+1;
			gap[minn+1]++;
			if(u!=s)
				u=pre[u];
		}
	}
	return ans;
}
int main()
{
	int cas,n,m;
	read(cas);
	while(cas--)
	{
		read(n),read(m);
		build_graph(n,m);
		int ans=sap();
		if(ans>=99999)//如果求出的最小割为无限大,即始终进入有控制面板的门
			printf("PANIC ROOM BREACH\n");
		else
			printf("%d\n",ans);
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值