【BZOJ5070】危险的迷宫 最小费用最大流
最大 最小 费用 危险 迷宫
2023-09-11 14:15:24 时间
【BZOJ5070】危险的迷宫
Description
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
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; }
相关文章
- hdu1358 最小循环节,最大循环次数 KMP
- 【网络流24题】魔术球问题(最大流)
- 【网络流24题】餐巾计划问题(最小费用最大流)
- 【BZOJ1877】[SDOI2009]晨跑 最小费用最大流
- 《图论资料》最大权闭合子图
- k8s【资源管理(resources)】4--LimitRange为配置命名空间内存最小和最大约束
- 二叉树——求最大节点和最小节点
- 《模式识别》学习笔记(十七)含拒绝判断的最小损失准则,最小最大损失准则
- [算法]最大连续子数组和,最长重复子串,最长无重复字符子串
- UVa 1161 Objective: Berlin (最大流)
- 聚类算法-最大最小距离算法(实例+代码)
- 最大匹配数,最小路径覆盖数,最大独立数,最小点覆盖数 定理总结
- Flexbox 的最小宽度和最大宽度声明在 Safari 上不起作用?为什么?
- 力扣解法汇总1599. 经营摩天轮的最大利润
- 求方阵的鞍点(即在行最小列最大的那个点)
- 【leetcode 53】 最大子序和
- 【bzoj3261】最大异或和 可持久化Trie树
- 【bzoj1565】[NOI2009]植物大战僵尸 拓扑排序+最大权闭合图
- 【bzoj1738】[Usaco2005 mar]Ombrophobic Bovines 发抖的牛 Floyd+二分+网络流最大流
- 滴滴今年发力智慧交通、国际化 “最大裁员”不攻自破