题目链接

不难想到每次一定是切一片。

如果是平面,很容易想到直接做二分图匹配。对于3维的?

可以发现min(a,b,c)的最大值只有\(\sqrt[3]{n}≈17\),我们暴力枚举这一最小值代表的是否选,对于剩下的两维二分图匹配

能用匈牙利用什么网络流啊。。那么麻烦。

在匈牙利求解的过程中记得剪枝!

因为实际点数应该远不如ABC多,所以不要用memset,直接枚举。

woc卡着时过去 那些2000ms怎么过的??网上找不到快点的懒得找。

zz的我考场写的每次DFS完重建边+ISAP,成功都T掉(写的方式也很zz。。每次A*B*C枚举找点,3次讨论3个DFS)

感觉数据一点梯度没有全都5s过。。

//2120kb	9300ms
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=5005,M=20005; int A,B,C,Time,bit[1<<18],Enum,H[N],nxt[M],to[M],match[N],use[N]/*空间换时间,匈牙利常用套路?*/;
struct Node
{
int x,y,z;
Node() {}
Node(int x,int y,int z):x(x),y(y),z(z) {}
}node[N]; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline void AddEdge(int u,int v){
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
}
bool Can(int x)
{
for(int i=H[x]; i; i=nxt[i])
if(use[to[i]]!=Time)
{
use[to[i]]=Time;
if(!match[to[i]]||Can(match[to[i]]))
{
match[to[i]]=x; return 1;
}
}
return 0;
}
void Work()
{
A=read(),B=read(),C=read();
int cnt=0, f=1, lim=1<<A, Max=B;
if(B<A&&B<C) f=2,lim=1<<B,Max=A;
else if(C<A&&C<B) f=3,lim=1<<C,Max=A;
for(int v,i=1; i<=A; ++i)
for(int j=1; j<=B; ++j)
for(int k=1; k<=C; ++k)
if(v=read()){
if(f==1) node[++cnt]=Node(i,j,k);//直接交换顺序即可
else if(f==2) node[++cnt]=Node(j,i,k);
else node[++cnt]=Node(k,i,j);
}
int Ans=10000;
for(int res,s=0; s<lim; ++s)
{
Enum=0, res=bit[s];//Time并不需要清空
// memset(H,0,sizeof H), memset(match,0,sizeof match);
for(int i=1; i<=Max; ++i) H[i]=match[i]=0;//here!
for(int i=1; i<=cnt; ++i)
if(!(s&(1<<node[i].x-1))) AddEdge(node[i].y,node[i].z);
++Time;
for(int i=1; i<=Max; ++i,++Time)
if(Can(i)){
if(++res>=Ans) break;//剪枝!
}
Ans=std::min(res,Ans);
}
printf("%d\n",Ans);
} int main()
{
for(int i=1; i<(1<<17); ++i) bit[i]=bit[i>>1]+(i&1);
int T=read();
while(T--) Work();
return 0;
}

BZOJ.3140.[HNOI2013]消毒(二分图匹配 匈牙利)的更多相关文章

  1. bzoj 3140: [Hnoi2013]消毒

    3140: [Hnoi2013]消毒 Description 最近在生物实验室工作的小T遇到了大麻烦. 由于实验室最近升级的缘故,他的分格实验皿是一个长方体,其尺寸为a*b*c,a.b.c 均为正整数 ...

  2. [bzoj]1059矩阵游戏<二分图匹配*匈牙利算法>

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1059 初见此题,我觉得这是水题,我认为只要每一行和每一列至少存在一个黑格就可以出现对角线, ...

  3. [HNOI2013] 消毒 - 二分图匹配

    容易发现 \(a,b,c\) 中至少有一个 \(\leq 17\) 不妨将其调剂为 \(a\),那么我们可以暴力枚举哪些 \(x\) 片片要被直接削掉,剩下的拍扁成二维情况 二维时,如果有一个格子是 ...

  4. 【刷题】BZOJ 3140 [Hnoi2013]消毒

    Description 最近在生物实验室工作的小T遇到了大麻烦. 由于实验室最近升级的缘故,他的分格实验皿是一个长方体,其尺寸为abc,a.b.c 均为正整数.为了实验的方便,它被划分为abc个单位立 ...

  5. [BZOJ3140][HNOI2013]消毒(二分图最小点覆盖)

    3140: [Hnoi2013]消毒 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1621  Solved: 676[Submit][Status] ...

  6. HDU 5943 Kingdom of Obsession 【二分图匹配 匈牙利算法】 (2016年中国大学生程序设计竞赛(杭州))

    Kingdom of Obsession Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Oth ...

  7. USACO 4.2 The Perfect Stall(二分图匹配匈牙利算法)

    The Perfect StallHal Burch Farmer John completed his new barn just last week, complete with all the ...

  8. [ZJOI2009]假期的宿舍 二分图匹配匈牙利

    [ZJOI2009]假期的宿舍 二分图匹配匈牙利 一个人对应一张床,每个人对床可能不止一种选择,可以猜出是二分图匹配. 床只能由本校的学生提供,而需要床的有住校并且本校和外校两种人.最后统计二分图匹配 ...

  9. 3140:[HNOI2013]消毒 - BZOJ

    题目描述 Description 最近在生物实验室工作的小 T 遇到了大麻烦. 由于实验室最近升级的缘故,他的分格实验皿是一个长方体,其尺寸为 a*b*c,a.b.c均为正整数.为了实验的方便,它被划 ...

随机推荐

  1. LOJ 3093: 洛谷 P5323: 「BJOI2019」光线

    题目传送门:LOJ #3093. 题意简述: 有 \(n\) 面玻璃,第 \(i\) 面的透光率为 \(a\),反射率为 \(b\). 问把这 \(n\) 面玻璃按顺序叠在一起后,\(n\) 层玻璃的 ...

  2. Django配置图片上传

    本文首先实现django中上传图片的过程,然后解决富文本编辑器文件上传的问题. 一. 上传图片 1.在 settings.py 中配置MEDIA_URL  和 MEDIA_ROOT 在 D:\blog ...

  3. v4l2文档之——color and format【转】

    转自:https://blog.csdn.net/zoe6553/article/details/17715407 v4l2文档第五A--颜色与格式    颜色与格式这是不定期发布的关于写视频驱动程序 ...

  4. Nodejs+定时截图+发送邮件

    功能 每天定时截图,并把截到的图片自动通过邮件发送. 说明 代码注释已经非常详细,就不多做说明,需要的朋友自己查看代码即可,主文件Mail.js,截图文件capturePart1.js,capture ...

  5. HTML学习笔记06-连接

    HTML超链接 HTML使用标签<a>来设置文本超链接. 超链接可以是文字,也可以是图片,点击这些内容跳转到新的文档或当前文档的某个部分 代码类似这样: <a href=" ...

  6. vw 解决方案

    vw 解决方案 1. 安装并配置PostCss插件 复制代码代码如下: npm i postcss-aspect-ratio-mini postcss-px-to-viewport postcss-w ...

  7. Python-元类 单例

    2.元类 用于创建类的类 叫元类 默认元类都是type 主要控制类的创建__init__ 类的实例化__call__ 3.单例 一个类只有一个实例 什么时候用,大家公用一个打印机,没有必要每个人创建新 ...

  8. 中文多分类 BERT

    直接把自己的工作文档导入的,由于是在外企工作,所以都是英文写的 Steps: git clone https://github.com/google-research/bert prepare dat ...

  9. keras2.0的一些变化

    keras 变化太快了https://github.com/fchollet/keras/wiki/Keras-2.0-release-notes

  10. linux 下 eclipse 安装

    下载: 官网选择相应安装包下载,我这里下了tar.gz包 安装: tar xzvf eclipse-inst-linux64.tar.gz 设置环境变量 export JAVA_HOME=/usr/l ...