BZOJ 3140 [HNOI2013]消毒
2016-03-08 23:20
246 查看
title: ‘BZOJ 3140 [HNOI2013]消毒’
categories: BZOJ
date: 2016-3-8 21:51:00
tags: [二分图匹配]
由于实验室最近升级的缘故,他的分格实验皿是一个长方体,其尺寸为a*b*c,a、b、c 均为正整数。为了实验的方便,它被划分为a*b*c个单位立方体区域,每个单位立方体尺寸为1*1*1。用(i,j,k)标识一个单位立方体,1 ≤i≤a,1≤j≤b,1≤k≤c。这个实验皿已经很久没有人用了,现在,小T被导师要求将其中一些单位立方体区域进 行消毒操作(每个区域可以被重复消毒)。而由于严格的实验要求,他被要求使用一种特定 的F试剂来进行消毒。 这种F试剂特别奇怪,每次对尺寸为x*y*z的长方体区域(它由x*y*z个单位立方体组 成)进行消毒时,只需要使用min{x,y,z}单位的F试剂。F试剂的价格不菲,这可难倒了小T。现在请你告诉他,最少要用多少单位的F试剂。(注:min{x,y,z}表示x、y、z中的最小者。)
1
4 4 4
1 0 1 1
0 0 1 1
0 0 0 0
0 0 0 0
0 0 1 1
1 0 1 1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
outut.txt
3
首先考虑二维的情况怎么做,相当于是要选出一些行和列涂色,覆盖所有的标记点。我们可以将标记点的横坐标与纵坐标连边,然后求二分图最小点覆盖,这样做的含义是用最少的点(行或列)来覆盖所有的边(标记点),这样直接二分图最大匹配用匈牙利跑就好了。
那三维的该怎么做呢?我们发现a*b*c≤5000,也就是说最小的那一维最多只有17,我们直接考虑2^k枚举哪一层选或不选,之后的跟二维一样。
categories: BZOJ
date: 2016-3-8 21:51:00
tags: [二分图匹配]
Description
最近在生物实验室工作的小T遇到了大麻烦。由于实验室最近升级的缘故,他的分格实验皿是一个长方体,其尺寸为a*b*c,a、b、c 均为正整数。为了实验的方便,它被划分为a*b*c个单位立方体区域,每个单位立方体尺寸为1*1*1。用(i,j,k)标识一个单位立方体,1 ≤i≤a,1≤j≤b,1≤k≤c。这个实验皿已经很久没有人用了,现在,小T被导师要求将其中一些单位立方体区域进 行消毒操作(每个区域可以被重复消毒)。而由于严格的实验要求,他被要求使用一种特定 的F试剂来进行消毒。 这种F试剂特别奇怪,每次对尺寸为x*y*z的长方体区域(它由x*y*z个单位立方体组 成)进行消毒时,只需要使用min{x,y,z}单位的F试剂。F试剂的价格不菲,这可难倒了小T。现在请你告诉他,最少要用多少单位的F试剂。(注:min{x,y,z}表示x、y、z中的最小者。)
Input
第一行是一个正整数D,表示数据组数。接下来是D组数据,每组数据开头是三个数a,b,c表示实验皿的尺寸。接下来会出现a个b 行c列的用空格隔开的01矩阵,0表示对应的单位立方体不要求消毒,1表示对应的单位立方体需要消毒;例如,如果第1个01矩阵的第2行第3列为1,则表示单位立方体(1,2,3)需要被消毒。输入保证满足a*b*c≤5000,T≤3。Output
仅包含D行,每行一个整数,表示对应实验皿最少要用多少单位 的F试剂。Sample
input.txt1
4 4 4
1 0 1 1
0 0 1 1
0 0 0 0
0 0 0 0
0 0 1 1
1 0 1 1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
outut.txt
3
Solution
贪心地考虑,肯定是当一维选最小的时候,剩下两维都取到最大值。首先考虑二维的情况怎么做,相当于是要选出一些行和列涂色,覆盖所有的标记点。我们可以将标记点的横坐标与纵坐标连边,然后求二分图最小点覆盖,这样做的含义是用最少的点(行或列)来覆盖所有的边(标记点),这样直接二分图最大匹配用匈牙利跑就好了。
那三维的该怎么做呢?我们发现a*b*c≤5000,也就是说最小的那一维最多只有17,我们直接考虑2^k枚举哪一层选或不选,之后的跟二维一样。
Code
#include<cstdio> #include<cstdlib> #include<cstring> #include<algorithm> #define maxn 5000+5 using namespace std; struct Point{ int x,y,z; }p[maxn]; struct Edge{ int u,v,nxt; }e[maxn<<2]; int head[maxn],link[maxn],vis[maxn],bin[maxn],ok[maxn]; int a,b,c,A,B,C,cnt,ind,ans=0x3f3f3f3f; void init(){ bin[0]=1; for(int i=1;i<=18;i++) bin[i]=bin[i-1]<<1; } void addedge(int x,int y){ e[ind]=(Edge){x,y,head[x]},head[x]=ind++; } bool match(int x){ for(int i=head[x];i!=-1;i=e[i].nxt) if(!vis[e[i].v]){ vis[e[i].v]=1; if(link[e[i].v]==-1 || match(link[e[i].v])){ link[e[i].v]=x; return true; } } return false; } void work(int s){ int res=0; ind=0; for(int i=1;i<=c;i++) ok[i]=0,head[i]=-1,link[i]=-1; for(int i=1;i<=a;i++) if(s&(1<<(i-1))) res++,ok[i]=1; for(int i=1;i<=cnt;i++) if(!ok[p[i].x]) addedge(p[i].y,p[i].z); for(int i=1;i<=b;i++){ for(int j=1;j<=c;j++) vis[j]=0; match(i); } for(int i=1;i<=c;i++) if(link[i]>0) res++; ans=min(ans,res); } int main(){ int T; init(); scanf("%d",&T); while(T--){ cnt=0; ans=0x3f3f3f3f; scanf("%d%d%d",&a,&b,&c); for(int i=1;i<=a;i++) for(int j=1;j<=b;j++) for(int k=1;k<=c;k++){ A=a; B=b; C=c; int opt; scanf("%d",&opt); if(opt) p[++cnt]=(Point){i,j,k}; if(A>B) swap(A,B),swap(p[cnt].x,p[cnt].y); if(A>C) swap(A,C),swap(p[cnt].x,p[cnt].z); } a=A; b=B; c=C; for(int s=0;s<bin[a];s++) work(s); printf("%d\n",ans); } return 0; }
相关文章推荐
- 冒泡排序--BubbleSort
- 【bzoj4421】【cerc2015】【Digit Division】
- 【杭电oj】1716 - 排列2(字典序输出)
- Linux普通用户使用1024以下端口的问题
- 构建可枚举类型(IEnumerable和IEnumerator接口)
- Gulp实现web服务器
- 蓝桥杯:剪格子
- 高性能MySQL主从复制
- 杭电2089~不要62
- 【Eclipse】1、eclipse软件下载
- 正则--match
- 正则表达式----search
- 一种更清晰的 Android 架构
- [转]考虑 PHP 5.0~5.6 各版本兼容性的 cURL 文件上传
- 5-3 逆序的三位数
- BZOJ 1008 越狱
- Jquery插件
- tomcat启动后访问应用超时问题
- kali windows双系统安装!
- 冒泡排序