BZOJ3140: [Hnoi2013]消毒（匈牙利）

最新推荐文章于 2022-01-07 23:05:11 发布

Dreamstar_DS

最新推荐文章于 2022-01-07 23:05:11 发布

阅读量313

点赞数

分类专栏：二分图文章标签： OI

本文链接：https://blog.csdn.net/Dreamstar_DS/article/details/82500109

版权

二分图专栏收录该内容

1 篇文章 0 订阅

订阅专栏

3140: [Hnoi2013]消毒

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 1808 Solved: 711
[Submit][Status][Discuss]
Description

最近在生物实验室工作的小T遇到了大麻烦。
由于实验室最近升级的缘故，他的分格实验皿是一个长方体,其尺寸为abc，a、b、c 均为正整数。为了实验的方便，它被划分为abc个单位立方体区域，每个单位立方体尺寸
为111。用(i,j,k)标识一个单位立方体，1 ≤i≤a，1≤j≤b，1≤k≤c。这个实验皿已经很久没有人用了，现在，小T被导师要求将其中一些单位立方体区域进行消毒操作（每个区域可以被重复消毒）。而由于严格的实验要求，他被要求使用一种特定的F试剂来进行消毒。这种F试剂特别奇怪，每次对尺寸为xyz的长方体区域（它由xyz个单位立方体组成）进行消毒时，只需要使用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)需要被消毒。输入保证满足abc≤5000,T≤3。
Output

仅包含D行，每行一个整数，表示对应实验皿最少要用多少单位的F试剂。
Sample Input

4 4 4

1 0 1 1

0 0 1 1

0 0 0 0

0 0 1 1

1 0 1 1

0 0 0 0

1 0 0 0

0 0 0 0

1 0 0 0

Sample Output

如果有人会“三分图匹配”不妨写下…
好吧肯定我们不需要开发一种新算法，条件min{a , b , c}决定了这道题是二分图匹配，而题目中的a * b *c <= 5000是一个很重要的信息，它表示{a , b , c}中最小的数<=17，因此我们可以枚举最小的数，然后再来跑二分图匹配
进行了一些剪枝，详见Code
注意：本题不要用memset！枚举时用memset显然要T，另外vis直接用++vis0的方式表示，不要去清空

AC Code：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define rg register
#define il inline
#define maxn 5005
#define lid id << 1
#define rid (id << 1) | 1
#define ll long long
using namespace std;
il int read(){rg int x = 0 , w = 1;rg char ch = getchar();while (ch < '0' || ch > '9'){if (ch == '-') w = -1;ch = getchar();}while (ch >= '0' && ch <= '9'){x = (x << 3) + (x << 1) + ch - '0';ch = getchar();}return x * w;}
struct edge{
    int to , next;  
}e[maxn << 5];
struct point{
    int x , y , z;  
}p[maxn];
int head[maxn << 1] , cnt , tot , match[maxn << 1] , ans , mi[maxn];
ll vis[maxn << 1] , vis0;
void add(int u,int v){
    e[++cnt].to = v;
    e[cnt].next = head[u];
    head[u] = cnt;
}
int find(int u){
    for (rg int i = head[u] ; i ; i = e[i].next){
        rg int to = e[i].to;
        if (vis[to] < vis0){
            vis[to] = vis0;
            if (!match[to] || find(match[to])){
                match[to] = u;
                return 1;
            }
        }
    }
    return 0;
}
void work(int x,int a,int b,int c,int tot){
    rg int res = 0;
    cnt = 0;
    for (rg int i = 0;i < a;++i) if (x & mi[i]) ++res;
    if (res >= ans) return;
    for (rg int i = 1;i <= b;++i) head[i] = 0;
    for (rg int i = 1;i <= c;++i) match[i] = 0;
    for (rg int i = 1;i <= tot;++i)
        if (!(x & mi[p[i].x - 1])) add(p[i].y , p[i].z);
    for (rg int i = 1;i <= b;++i){
        ++vis0;
        res += find(i);
        if (res >= ans) return;
    }
    ans = (res < ans)?res:ans;
}
bool cmp(point a , point b){
    if (a.x!=b.x) return a.x<b.x;
    if (a.y!=b.y) return a.y<b.y;
    return a.z<b.z;
}
int main(){
    mi[0] = 1;for(rg int i = 1;i <= 19;++i) mi[i] = (mi[i - 1] << 1); 
    rg int t = read();
    while (t--){
        rg int a = read() , b = read() , c = read();
        int num , x , y , z;
        cnt = tot = 0;
        int id;
        for (rg int i = 1;i <= a;++i) for (rg int j = 1;j <= b;++j) for (rg int k = 1;k <= c;++k){
            if (read()) p[++tot].x = i , p[tot].y = j , p[tot].z = k;
        }
        if (b < a) {swap(a , b);for (rg int i = 1;i <= tot;++i) swap(p[i].x , p[i].y);}
        if (c < a) {swap(a , c);for (rg int i = 1;i <= tot;++i) swap(p[i].x , p[i].z);}
        sort(p+1,p+1+tot,cmp);
        ans = 999999999;
        for (rg int i = 0;i < mi[a];++i) work(i , a , b , c , tot);
        printf("%d\n",ans);
    }
    return 0;
}