首先根据题意,我们可以一维取1,另外两维取到最大,那这样的方块就是一个截面,而且这样消毒一定最优,那么可以把题意理解为一个立方体,有些点需要被染色,每次可以对一个截面染色,最少要染多少次
如果这道题是二维那就是一道非常经典的二分图匹配问题,对于一个点(x,y),第x行和第y列连边,要覆盖所有的边,最小覆盖=最大匹配数,现在是三维,观察数据范围,因为
a∗b∗c≤5000
,所以最小的数
≤17
,所以
2k
枚举最小的一维切哪些,另外两维的切面就可以转化成二维的二分图匹配问题
换坐标换了好久…各种细节,各种漏洞
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;
}