Bzoj 4950 (二分图最大匹配)

Description

那是春日里一个天气晴朗的好日子,你准备去见见你的老朋友Patrick,也是你之前的犯罪同伙。Patrick在编程竞赛

上豪赌输掉了一大笔钱,所以他需要再干一票。为此他需要你的帮助,虽然你已经金盆洗手了。你刚开始很不情愿,

因为你一点也不想再回到那条老路上了,但是你觉得听一下他的计划也无伤大雅。在附近的一个仓库里有一批货物,

包含一些贵重的消费性部件,Patrick企图从中尽可能多地偷些东西出来。这意味着要找一条进去的路,弄晕安保人

员,穿过各种各样的激光射线,你懂的,都是常见的抢劫技术。然而,仓库的核心装备了一套Patrick搞不定的安保系

统。这也是他需要你帮助他的地方。这批货物被放置在一些巨大的立方体箱里,每个箱子的尺寸都是相同的。这些

箱子堆放成许多整齐的堆,每个箱子可以表示成一个三维的网格。安保系统每个小时会用三台相机对这堆货物进行

一次拍照,相机分别为:前置相机(front camera),侧置相机(side camera)和顶置相机(top camera)。前置相机的照

片显示了每一行最高的那堆箱子的高度,侧置相机显示了每一列最高的那堆箱子的高度,顶置相机显示了每个位置是

否存在一堆箱子。如果安保系统发现任何一张照片出现了变化,它会立即拉响警报。一旦 Patrick 进去了,他会确

定每堆箱子的高度并且发给你。图1显示了一种网格可能的放置,以及每台相机会得到的视图。

图 1. 网格的高度值与对应的相机视图。

图 2. 洗劫后网格可能的高度值。

Patrick想尽可能多偷走一些箱子。由于他不能弄坏安保系统,他准备重新安排剩余每堆箱子的放置,使得下一次相

机取像时会得到相同的照片,从而骗过安保系统。在上面的例子中,他可以偷走九个箱子。图2显示了一种可能的剩

余箱子的安置方案能使得安保系统认为与原安置情况相同。Patrick想请你帮他确定在保证能骗过安保系统的情况

下他最多能偷走多少个箱子。你会帮他干完这最后一票么?

Input

第一行包含两个整数r(1≤r≤100)和c(1≤n≤100),分别表示网格的行数与列数。

接下来r行,每行包含c个整数,表示对应行上每堆立方体箱的高度(箱子的数量)。

所有的高度在0到10^9之间 (含边界) 。

Output

输出在不被发现的情况下最多能偷走多少箱子。

Sample Input

样例1

5 5

1 4 0 5 2

2 1 2 0 1

0 2 3 4 4

0 3 0 3 1

1 2 2 1 1

样例2

2 3

50 20 3

20 10 3

Sample Output

样例1

9

样例2

30

---------------------------------------------------------------------

题解:

  对于俯视图,我们使有箱子的格子不被拿完即可。对于侧视图和正视图,使行和列的最大值保持不变即可。

  注意到某一行和某一列的最大值可能相同,对此,我们将尽量保留这样的最大值。

  将最大值相等的行和列连边,跑一遍二分图最大匹配即可。

 1 #include<stdio.h>
 2 #include<string.h>
 3 #include<algorithm>
 4 using namespace std;
 5 #define MAXN 100000+10
 6 typedef long long LL;
 7 struct ed{LL v,next;}edge[MAXN];
 8 LL match[220],head[220],vis[220],n,m;
 9 LL map[220][220],sum=0,mh[220],ml[220];
10 void add(LL u,LL v){
11     static LL tot=0;
12     edge[++tot].v=v;
13     edge[tot].next=head[u];
14     head[u]=tot;
15 }
16 bool dfs(LL u){
17     for(LL i=head[u];i;i=edge[i].next){
18         LL v=edge[i].v;
19         if(vis[v])continue;
20         vis[v]=1;
21         if(!match[v]||dfs(match[v])){
22             match[v]=u;
23             return true;
24         }
25     }
26     return false;
27 }
28 int main(){
29     scanf("%lld%lld",&n,&m);
30     for(LL i=1;i<=n;i++)
31         for(LL j=1;j<=m;j++){
32             scanf("%lld",&map[i][j]);
33             mh[i]=max(mh[i],map[i][j]);
34             ml[j]=max(ml[j],map[i][j]);
35             if(map[i][j])sum+=map[i][j]-1;
36         }
37     for(LL i=1;i<=n;i++)
38         for(LL j=1;j<=m;j++)
39             if(mh[i]&&map[i][j]&&mh[i]==ml[j])add(i,n+j);
40     for(LL i=1;i<=n;i++)if(mh[i])sum-=mh[i]-1;
41     for(LL i=1;i<=m;i++)if(ml[i])sum-=ml[i]-1;
42     for(LL i=1;i<=n;i++){
43         memset(vis,0,sizeof(vis));
44         if(mh[i]&&dfs(i))sum+=mh[i]-1;
45     }
46     printf("%lld",sum);
47     return 0;
48 }

  

时间: 2024-08-06 20:02:10

Bzoj 4950 (二分图最大匹配)的相关文章

BZOJ 1143 CTSC2008 祭祀river 二分图最大匹配

题目大意:给定一个拓扑图,求一个最大的点集,点集中的点两两不可达 这实际上就是让你求传递闭包后图的最大点独立集- - 利用二分图最大匹配就能搞- - #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define M 110 using namespace std; int n,m,ans; bool map[M][M]; int state[M],r

BZOJ 1191 HNOI2006 超级英雄Hero 二分图最大匹配

题目大意:给定n个锦囊和m个问题,每个问题可以使用给定的两个锦囊之一,必须连续答题,求最多答上多少题 二分图最大匹配,每出现一个问题就向给定的两个锦囊连边,然后匈牙利算法寻找增广路,如果找不到就break 我这SB居然一开始在两个锦囊之间练了条边,然后二分答案--伤不起啊-- #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define M 1010 usi

BZOJ 1854 SCOI2010 游戏 二分图最大匹配/并查集

题目大意:给定n个武器,每个武器有两个属性,只能使用其中一个,要求选择一些武器 可以造成形如1 2 3 4的伤害 求最大伤害 题目大意我没写明白还是去看原题把QAQ 做法1: 同 1191 每个武器向两个属性连边 然后从1~10000枚举属性 跑二分图最大匹配 无法匹配则输出答案 #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define M 100100

BZOJ 1191 HNOI 2006 超级英雄Hero 二分图最大匹配

题目大意:闯关答题,每一个题可以用两个trick中的一个来解决,一个trick最多只能使用一次,问最多可以连续答对多少题. 思路:我一开始就想到了二分图最大匹配,但是思路完全想歪了.我看每个题有两个trick可以用,就用这个来拆点建图,显然是错的.. 正确的是用每个题和每个trick来建边,来一个问题就建两条边,然后看能不能找到增广路,如果不能就无法答对这个题,输出. CODE: #include <cstdio> #include <cstring> #include <i

bzoj 1854: [Scoi2010]游戏 (并查集||二分图最大匹配)

链接: https://www.lydsy.com/JudgeOnline/problem.php?id=1854 写法1: 二分图最大匹配 思路:  将武器的属性对武器编号建边,因为只有10000种属性,我们直接对1-10000跑二分图匹配,同时用时间戳优化匹配. 实现代码: #include<bits/stdc++.h> using namespace std; const int M = 1e6+10; vector<int>g[M]; int vis[M],pre[M],t

匈牙利算法dfs模板 [二分图][二分图最大匹配]

最近学了二分图最大匹配,bfs模板却死活打不出来?我可能学了假的bfs 于是用到了dfs模板 寻找二分图最大匹配的算法是匈牙利算法 匈牙利算法的主要程序是寻找增广路 寻找增光路是过程是:从一个未经配对的点出发,历经未配边.匹配边.未配边.匹配边.未配边....最终到达一个未配点的过程,只要把路径中的未配边和匹配边的“身份”对调,匹配就加一了.这就是一个寻找增广路的过程,通过不断寻找增广路,可以找到最大的匹配. 1 #include<cstdio> 2 #include<cstring&g

图论——LCA、强联通分量、桥、割顶、二分图最大匹配、网络流

A: 交通运输线 时间限制: 5 Sec  内存限制: 128 MB 题目描述 战后有很多城市被严重破坏,我们需要重建城市.然而,有些建设材料只能在某些地方产生.因此,我们必须通过城市交通,来运送这些材料的城市.由于大部分道路已经在战争期间完全遭到破坏,可能有两个城市之间没有道路.当然在运输线中,更不可能存在圈. 现在,你的任务来了.给你战后的道路情况,我们想知道,两个城市之间是否存在道路,如果存在,输出这两个城市之间的最短路径长度. 输入 第一行一个整数Case(Case<=10)表示测试数据

POJ 2226二分图最大匹配

匈牙利算法是由匈牙利数学家Edmonds于1965年提出,因而得名.匈牙利算法是基于Hall定理中充分性证明的思想,它是二部图匹配最常见的算法,该算法的核心就是寻找增广路径,它是一种用增广路径求二分图最大匹配的算法. #include<stdio.h> #include<string.h> #include<stdlib.h> int n1,n2; char map[1005][1005]; //数组开大点 int mapx[1005][1005],mapy[1005]

【Codevs1922】骑士共存问题(最小割,二分图最大匹配)

题意: 在一个n*n个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示.棋盘上某些方格设置了障碍,骑士不得进入. 对于给定的n*n个方格的国际象棋棋盘和障碍标志,计算棋盘上最多可以放置多少个骑士,使得它们彼此互不攻击. n<=200,m<=n^2 思路:经典的二分图最大匹配问题,采用黑白点染色的思想. 如果按照相邻点黑白不同染色,可以发现每次跳到的点必定与现在所在点不同色,二分图最大匹配即可. 这里用最小割来解决,因为不能允许任何黑白点之间的任何一条边有流量,符合最小割的思想. 1