【BZOJ5070】危险的迷宫 最小费用最大流

【BZOJ5070】危险的迷宫

Description

JudgeOnline/upload/201710/55.doc

Input

第一行是两个整数A与B(1≤A,B≤10)，中间用空格分隔，表示该迷宫是A行B列的。

第2行至第A+1行，每行有B个1至100以内的整数，表示该迷宫每一格的危险程度。

以下一行是一个整数K。接着K行每行有四个整数X0,Y0,X1,Y1,

(1 ≤X0,X1≤A, 1≤Y0,Y1≤B) ，(X0,Y0),(X1,Y1)为相邻的两格，这两格相通。

接着一行是一个整数N（0≤N≤A*B/2），表示有N个出口与入口。

以下N行，每行有两个整数X0,Y0，表示每个入口的行列位置。

以下还有N行，每行有两个整数X1,Y1，表示每个出口的行列位置。

Output

输出仅一个数，若队员们不能全部到达指定目标位置，则输出-1；

否则输出所有队员所经过的所有单元格的危险程度之和。

Sample Input

3 4
20 30 40 30
30 60 20 20
20 15 20 20
13
1 1 2 1
1 2 1 3
1 2 2 2
1 3 1 4
1 4 2 4
2 1 2 2
2 1 3 1
2 2 2 3
2 3 2 4
2 4 3 4
3 1 3 2
3 2 3 3
3 3 3 4
2
1 1
1 2
2 3
3 4

Sample Output

235

题解：直接连边跑最小费用最大流即可。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#define P(X,Y) ((X-1)*B+Y)
#define Q(X,Y) ((X-1)*B+Y+A*B)
using namespace std;
int A,B,n,m,sum,ans,cnt,S,T;
int cost[100000],flow[100000],next[100000],head[100000],to[100000],dis[100000],inq[100000],pv[100000],pe[100000];
queue<int> q;
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')	f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+gc-'0',gc=getchar();
	return ret*f;
}
inline void add(int a,int b,int c,int d)
{
	to[cnt]=b,cost[cnt]=c,flow[cnt]=d,next[cnt]=head[a],head[a]=cnt++;
	to[cnt]=a,cost[cnt]=-c,flow[cnt]=0,next[cnt]=head[b],head[b]=cnt++;
}
inline int bfs()
{
	memset(dis,0x3f,sizeof(dis));
	dis[S]=0,q.push(S);
	int u,i;
	while(!q.empty())
	{
		u=q.front(),q.pop(),inq[u]=0;
		for(i=head[u];i!=-1;i=next[i])	if(dis[to[i]]>dis[u]+cost[i]&&flow[i])
		{
			dis[to[i]]=dis[u]+cost[i],pe[to[i]]=i,pv[to[i]]=u;
			if(!inq[to[i]])	inq[to[i]]=1,q.push(to[i]);
		}
	}
	return dis[T]<0x3f3f3f3f;
}
int main()
{
	A=rd(),B=rd(),S=0,T=2*A*B+1;
	int i,j,a,b,c,d,mf;
	memset(head,-1,sizeof(head));
	for(i=1;i<=A;i++)	for(j=1;j<=B;j++)	add(P(i,j),Q(i,j),rd(),1);
	m=rd();
	for(i=1;i<=m;i++)
	{
		a=rd(),b=rd(),c=rd(),d=rd();
		add(Q(a,b),P(c,d),0,1),add(Q(c,d),P(a,b),0,1);
	}
	n=rd();
	for(i=1;i<=n;i++)	a=rd(),b=rd(),add(S,P(a,b),0,1);
	for(i=1;i<=n;i++)	a=rd(),b=rd(),add(Q(a,b),T,0,1);
	while(bfs())
	{
		mf=1<<30;
		for(i=T;i!=S;i=pv[i])	mf=min(mf,flow[pe[i]]);
		sum+=mf,ans+=dis[T]*mf;
		for(i=T;i!=S;i=pv[i])	flow[pe[i]]-=mf,flow[pe[i]^1]+=mf;
	}
	if(sum<n)	printf("-1");
	else	printf("%d",ans);
	return 0;
}