[HNOI2013]消毒

题目描述

最近在生物实验室工作的小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中的最小 者。)

输入输出格式

输入格式：

第一行是一个正整数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。

输出格式：

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

思路：看了好多博客终于想明白了。。在二维上直接就是最小点覆盖，三维比二维对了一个操作，枚举哪一层直接用一的费用去消毒，对于我来说难点并不在这，而是剩下的那些层为什么可以压缩在一层上进行匹配，发挥一下空间想象力，对于那些覆盖在同一位置上的点，可以用费用为1去进行消毒(沿枚举方向)。。(感觉说的不太清楚啊

#include<bits/stdc++.h>
#define read(x) scanf("%d",&x)
using namespace std;
const int maxn=1e5+3;
struct node{
    int x,y,z;
}p[maxn];
struct edge{
    int u,v,next;
}e[maxn];
int head[maxn],link[maxn],bt[maxn];
int vis[maxn],ok[maxn];
int a,b,c,A,B,C;
int cnt=0;
int ind,ans=0x3f3f3f3f;
void init() {
    bt[0]=1;
    for(int i=1;i<=19;i++) {
        bt[i]=bt[i-1]<<1;
    }
}
void addedge(int f,int t) {
    e[ind]=(edge){f,t,head[f]};
    head[f]=ind++;
}
bool match(int x) {
    for(int i=head[x];i!=-1;i=e[i].next) {
        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 solve(int x) {
    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(x&(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);
//    cout<<res<<endl;
}
int main() {
    int T;
    init();
    read(T);
    while(T--) {
        cnt=0;
        ans=0x3f3f3f3f;
        read(a);
        read(b);
        read(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 pos;
                    read(pos);
                    if(pos)p[++cnt]=(node){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 i=0;i<bt[a];i++) {
            solve(i);
//            if(i==0)cout<<ans<<endl;
        }
        printf("%d\n",ans);
    }
}