BZOJ3140: [Hnoi2013]消毒

首先根据题意，我们可以一维取1，另外两维取到最大，那这样的方块就是一个截面，而且这样消毒一定最优，那么可以把题意理解为一个立方体，有些点需要被染色，每次可以对一个截面染色，最少要染多少次
如果这道题是二维那就是一道非常经典的二分图匹配问题，对于一个点(x,y)，第x行和第y列连边，要覆盖所有的边，最小覆盖=最大匹配数，现在是三维，观察数据范围，因为$a*b*c≤5000$，所以最小的数$≤17$，所以$2^k$枚举最小的一维切哪些，另外两维的切面就可以转化成二维的二分图匹配问题

---

换坐标换了好久…各种细节，各种漏洞

---

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

void up(int &x,int y){if(x<y)x=y;}
void down(int &x,int y){if(x>y)x=y;}
const int maxn = 60000;
const int maxm = 600000;

struct edge
{
    int y,nex;
    edge(){}
    edge(int _y,int _nex){y=_y;nex=_nex;}
}a[maxm]; int len,fir[maxn];
int n,m,h;

void ins(int x,int y)
{
    a[++len]=edge(y,fir[x]); fir[x]=len;
}
int match[maxn];
bool v[maxn];
bool find_(int x)
{
    if(v[x]) return false;
    v[x]=true;
    for(int k=fir[x];k;k=a[k].nex)
    {
        int y=a[k].y;
        if(!match[y]||find_(match[y]))
        {
            match[y]=x;
            return true;
        }
    }
    v[x]=false;
    return false;
}
struct node
{
    int x,y,z;
}p[maxn]; int N;
void build_(int k)
{
    len=0; for(int i=1;i<=n;i++) fir[i]=0;
    for(int i=1;i<=N;i++)
        if(!(k&(1<<p[i].z-1))) ins(p[i].x,n+p[i].y);
    for(int i=1;i<=m;i++) match[n+i]=0;
}
int solve()
{
    int ret=INT_MAX;
    for(int i=0;i<(1<<h);i++)
    {
        build_(i);
        int k=0;
        for(int j=0;j<h;j++) if(i&(1<<j)) k++;
        for(int j=1;j<=n;j++)
        {
            for(int l=1;l<=n;l++) v[l]=false;
            if(find_(j)) k++;
        }
        down(ret,k);
    }
    return ret;
}
struct abc
{
    int x,c;
}a1[6];
int p1[maxn][6];
bool cmp(abc x,abc y){return x.c>y.c;}

int main()
{
    int t; scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d",&n,&m,&h);
        a1[1].x=1; a1[1].c=n;
        a1[2].x=2; a1[2].c=m;
        a1[3].x=3; a1[3].c=h;
        sort(a1+1,a1+3+1,cmp);

        N=0;
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++) for(int l=1;l<=h;l++)
            {
                int c; scanf("%d",&c);
                if(c)
                {
                    N++;
                    p1[N][1]=i; p1[N][2]=j; p1[N][3]=l;
                }
            }
        }

        n=a1[1].c; m=a1[2].c; h=a1[3].c;
        for(int i=1;i<=N;i++)
        {
            p[i].x=p1[i][a1[1].x];
            p[i].y=p1[i][a1[2].x];
            p[i].z=p1[i][a1[3].x];
        }
        printf("%d\n",solve());
    }

    return 0;
}