BZOJ.3140.[HNOI2013]消毒(二分图匹配 匈牙利)

题目链接

不难想到每次一定是切一片。
如果是平面,很容易想到直接做二分图匹配。对于3维的?
可以发现min(a,b,c)的最大值只有\(\sqrt[3]{n}≈17\),我们暴力枚举这一最小值代表的是否选,对于剩下的两维二分图匹配

能用匈牙利用什么网络流啊。。那么麻烦。
在匈牙利求解的过程中记得剪枝!
因为实际点数应该远不如ABC多,所以不要用memset,直接枚举。
woc卡着时过去 那些2000ms怎么过的??网上找不到快点的懒得找

zz的我考场写的每次DFS完重建边+ISAP,成功都T掉(写的方式也很zz。。每次ABC枚举找点,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;
}

原文地址:https://www.cnblogs.com/SovietPower/p/8671252.html

时间: 2024-10-24 21:09:30

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

bzoj 3140: [Hnoi2013]消毒

3140: [Hnoi2013]消毒 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被导师要求将其中一些单位立方体区域进 行消毒操作(每个区域可以被重复消毒).而由于严格的实验要求,

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

The Perfect StallHal Burch Farmer John completed his new barn just last week, complete with all the latest milking technology. Unfortunately, due to engineering problems, all the stalls in the new barn are different. For the first week, Farmer John r

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

Kingdom of Obsession Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 49    Accepted Submission(s): 14 Problem Description There is a kindom of obsession, so people in this kingdom do things very

HDU1507 Uncle Tom&#39;s Inherited Land* 二分图匹配 匈牙利算法 黑白染色

原文链接http://www.cnblogs.com/zhouzhendong/p/8254062.html 题目传送门 - HDU1507 题意概括 有一个n*m的棋盘,有些点是废的. 现在让你用1*2的矩形覆盖所有的不废的点,并且不重叠,问最多可以覆盖多少个1*2的矩形,输出方案,有SPJ. 输入描述: 多组数据,每组首先两个数n,m(如果n和m为0,则结束程序) 然后给出k 然后给出k个二元组(x,y)表示废点的坐标. 题解 按照前两片博文的算法已经不行了,因为方案不对了. 所以我们要进行

矩阵游戏|ZJOI2007|BZOJ1059|codevs1433|luoguP1129|二分图匹配|匈牙利算法|Elena

1059: [ZJOI2007]矩阵游戏 Time Limit: 10 Sec  Memory Limit: 162 MB Description 小Q是一个非常聪明的孩子,除了国际象棋,他还很喜欢玩一个电脑益智游戏--矩阵游戏.矩阵游戏在一个N *N黑白方阵进行(如同国际象棋一般,只是颜色是随意的).每次可以对该矩阵进行两种操作:行交换操作:选择 矩阵的任意两行,交换这两行(即交换对应格子的颜色)列交换操作:选择矩阵的任意行列,交换这两列(即交换 对应格子的颜色)游戏的目标,即通过若干次操作,

网络流24题 第一题 - 洛谷2756 飞行员配对方案 二分图匹配 匈牙利算法

欢迎访问~原文出处--博客园-zhouzhendong 去博客园看该题解 题目传送门 题意概括 裸的二分图匹配 题解 匈牙利算法 上板子 代码 #include <cstring> #include <cstdio> #include <algorithm> #include <cstdlib> #include <cmath> using namespace std; const int N=100+5; int m,n,a,b,match[N

图论-二分图匹配-匈牙利算法

有关概念: 二分图:图G中的点集可以分为两个互不相交的子集,且G中的每条边连接的两个点分别属于这两个子集 二分图匹配:二分图G的子图M中每个结点上只连一条边,则称M为一个匹配 极大匹配:无法再向二分图中加边且满足匹配条件的匹配 最大匹配:所有极大匹配中边数最多的一个 增广路:若P为图G上连接两个未匹配结点的路径,且已匹配边和未匹配边在P上交替出现,则称P为相对于M的一条增广路 匈牙利算法即用来求二分图的极大匹配 思路: 在图G中找出增广路P,对P上每一条边取反(即已匹配边改为未匹配边,未匹配边改

#图# #二分图匹配# #匈牙利算法#

二分图 设G=(V,E)是一个无向图,如果顶点V可分割为两个互不相交的子集(A,B),并且图中的每条边(i,j)所关联的两个顶点i和j分别属于这两个不同的顶点集(i in A,j in B),则称图G为一个二分图. 区别二分图,关键是看点集是否能分成两个独立的点集. 二分图匹配(匈牙利算法) 最大匹配:设 是一个无向图.如顶点集V可分割为两个互不相交的子集 ,选择这样的子集中边数最大的子集称为图的最大匹配问题(maximal matching problem). 完全匹配:如果一个匹配中,  且

【二分图匹配/匈牙利算法】飞行员配对方案问题

P2756 飞行员配对方案问题 确认过眼神, 是二分图匹配板子题啦!!! 跑个匈牙利, 有匹配的输出, 记得先输出外籍飞行员, 因为有spj顺序无所谓啦qwq 最近A的最顺利的题了哈哈哈哈哈哈开心!!!!!!!! 1 #include<cstdio> 2 #include<cstring> 3 #include<iostream> 4 using namespace std; 5 const int sz = 100010; 6 int n, m, num = 0, a