题目描述

最近在生物实验室工作的小 T 遇到了大麻烦。由于实验室最近升级的缘故，他的分格实验皿是一个长方体，其尺寸为 a∗b∗ca*b*ca∗b∗c。为了实验的方便，它被划分为 a∗b∗ca*b*ca∗b∗c 个单位立方体区域，每个单位立方体尺寸为 1∗1∗11*1*11∗1∗1，并用 (i,j,k)(i,j,k)(i,j,k) 标识一个单位立方体。这个实验皿已经很久没有人用了。现在，小 T 被导师要求将其中一些单位立方体区域进行消毒操作（每个区域可以被重复消毒）。

而由于严格的实验要求，他被要求使用一种特定的 F 试剂来进行消毒。这种 F 试剂特别奇怪，每次对尺寸为 x∗y∗zx*y*zx∗y∗z 的长方体区域（它由 x∗y∗zx*y*zx∗y∗z 个单位立方体组成）进行消毒时，只需要使用 min(x,y,z)min(x,y,z)min(x,y,z) 单位的 F 试剂。F 试剂的价格不菲，这可难倒了小 T。

现在请你告诉他，最少要用多少单位的 F 试剂。

解析

非暴力，不合作

首先可以有一个结论：**每次使min(x,y,z)=1min(x,y,z)=1min(x,y,z)=1,一定是不劣的
所以我们就每次一面一面的涂
看一个《弱化版》的题目

在一N∗NN*NN∗N个的矩阵中，有 KKK个格子中有杂物，现在你有一种能力，一次可以消除一行或一列格子中的杂物，问你至少需要几次可以将这些杂物全部消完。

这题应该是二分图的入门题了，把每个点的x坐标与y坐标相连，跑二分图最大匹配即可
不难发现，消毒这题应该就是消除杂物的升级版，从二维变成了三维
然鹅很快我们就发现并不能推广到k维。。。
当然如果您能发明三分图匹配本题就和喝水一样
那么我们怎么办呢？
然后就点开了题解
还是暴力的思想了
考虑到数据范围：
a∗b∗c<=5000a*b*c<=5000a∗b∗c<=5000

那么a、b、c中的最小值应该不超过17
所以我们考虑暴力枚举最小的一维的选取状态，然后每次跑一遍匈牙利取答案最小值即可

另外本题还有一个很巧妙的实现技巧
先把坐标存到三个一维数组里
把最小的一维swap到a的位置同时把对应的那一位的数组swap到第一位
代码实现就变得很简单了

代码

#include<bits/stdc++.h>
using namespace std;
const int N=5020;
#define ll long long
ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();};while(isdigit(c)){x=x*10+c-'0';c=getchar();};return x*f;
}
int n,m,l;
struct node{int from,to,nxt;
}p[N];
int fi[N],cnt=-1;
void addline(int x,int y){p[++cnt]=(node){x,y,fi[x]};fi[x]=cnt;
}
int vis[N],mat[N];
bool dfs(int x,int tim){if(vis[x]==tim) return false;vis[x]=tim;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(!mat[to]||dfs(mat[to],tim)){mat[to]=x;return true;}}return false;
}
int a,b,c;
int hungary(){int res=0;memset(vis,0,sizeof(vis));memset(mat,0,sizeof(mat));for(int i=1;i<=b;i++){if(dfs(i,i)){res++;//printf("  ok:%d\n",i);}}return res;
}int q[4][N],num;
bool ok[35];
int ans;
int calc(){memset(fi,-1,sizeof(fi));cnt=-1;for(int i=1;i<=num;i++){if(ok[q[1][i]]) continue;int x=q[2][i],y=q[3][i];addline(x,y+b);addline(y+b,x);//    printf("x=%d y=%d\n",x,y);}return hungary();
}
void find(int k,int val){if(k>a){//printf("ok");//printf("-------------val=%d\n",val);//for(int i=1;i<=a;i++) printf("%d ",ok[i]);//printf("\n");ans=min(ans,calc()+val);//printf("---ans=%d\n\n",ans);return;}find(k+1,val);ok[k]=1;find(k+1,val+1);ok[k]=0;
}
int main(){int T=read();while(T--){a=read();b=read();c=read();num=0;ans=5000;int mn=min(a,min(b,c));for(int i=1;i<=a;i++){for(int j=1;j<=b;j++){for(int k=1;k<=c;k++){int x=read();if(!x) continue;q[1][++num]=i;q[2][num]=j;q[3][num]=k;}}}if(mn==b){swap(a,b);swap(q[1],q[2]);}else if(mn==c){swap(a,c);swap(q[1],q[3]);}find(1,0);printf("%d\n",ans);}return 0;
}
/**/

YBTOJ洛谷P3231：消毒（二分图匹配）相关推荐

YbtOJ#573-后缀表达【二分图匹配】
正题题目链接:https://www.ybtoj.com.cn/contest/115/problem/2 题目大意给出一个包含字母变量和若干种同级操作符的后缀表达式.求一个等价的表达式满足该表达 ...
YBTOJ洛谷P4074：糖果公园（树上莫队）
文章目录解析 update: 代码所谓树上莫队,就是在树上的莫队 (逃) 传送门解析似乎就是树上的这道题考虑如何转化为序列问题呢? 考虑dfs序但是又一个问题... 似乎这条链的dfs序不 ...
YBTOJ洛谷P4551：最长异或路径（trie树）
洛谷传送门文章目录题目描述解析代码题目描述解析本题关键就在于一点: 若把每个点的深度dep[i]定义为从根到节点边权的异或和那么i到j的路径异或和可以表示为: dep[i] ^ dep ...
BZOJ.3140.[HNOI2013]消毒(二分图匹配匈牙利)
题目链接不难想到每次一定是切一片. 如果是平面,很容易想到直接做二分图匹配.对于3维的? 可以发现min(a,b,c)的最大值只有\(\sqrt[3]{n}≈17\),我们暴力枚举这一最小值代表的是 ...
YBTOJ洛谷P4298：祭祀（二分图匹配）
文章目录题目描述解析解析题目描述在遥远的东方,有一个神秘的民族,自称Y族.他们世代居住在水面上,奉龙王为神.每逢重大庆典, Y族都会在水面上举办盛大的祭祀活动.我们可以把Y族居住地水系看成一 ...
YBTOJ：消除格子（二分图匹配）
文章目录题目描述解析代码题目描述在一个n*n的矩阵中,有 k个格子中有杂物,现在你有一种能力,一次可以消除一行或一列格子中的杂物,问你至少需要几次可以将这些杂物全部消完. 解析看起来像是网 ...
YBTOJ洛谷P4068：数字配对（网络流）
文章目录题目描述解析代码题目描述有n种数字,第 i 种数字是 ai,有 bi个,权值是 ci. 若两个数字 ai,aj 满足, ai 是 aj 的倍数,且 ai/aj 是一个质数,那么这两个 ...
YBTOJ洛谷P2223：软件开发（费用流）
文章目录题目描述解析代码题目描述某软件公司正在规划一项nnn天的软件开发计划,根据开发计划第iii天需要nin_ini个软件开发人员,为了提高软件开发人员的效率,公司给软件人员提供了很多的 ...
YBTOJ：伞兵空降（二分图匹配）
文章目录题目描述解析代码题目描述有n个点和m条边的有向无环图,在这张图上的某些点上空投伞兵,使伞兵可以走到图上所有的点. 且每个点只能被一个伞兵走一次.问至少需要放多少伞兵. 解析考虑一开 ...

YBTOJ洛谷P3231：消毒（二分图匹配）

文章目录

题目描述

解析

代码

YBTOJ洛谷P3231：消毒（二分图匹配）相关推荐

最新文章

热门文章