[HNOI2013]消毒

如果只有二维，那么我们建二分图求最大匹配，O(n^2)的复杂度可以实现。而现在是三维，发现不会三分图，好像也没有三分图这个算法…

再思考一下对于二维，是否有别的暴力实现方法？

我们可以暴力枚举哪些行是被消掉的，然后对于不是剩下的黑点，枚举每一列是否需要再消一次。复杂度：O(2^n*m)

由于a * b * c<=5000，可得min(a,b,c)<=17,所以现在我们也可以暴力枚举每一层是否消掉，然后对于没有消掉的层，进行二分图求最大匹配，因为现在三维就变成二维的了。

我们将最小的那一维作为立方体的高，这样可以保证复杂度为：O(2^17*295)。

#include <bits/stdc++.h>
using namespace std;
int T,a,b,c,x,tot,ans;
int num[4][5005],match[5005];
bool pd[18],vis[5005];
int cnt,head[5005];
struct edge{int next,to;}e[5005];

inline void add(int u,int v)
{
	cnt++;
	e[cnt].next=head[u];
	e[cnt].to=v;
	head[u]=cnt;
}

int dfs(int u)
{
	for (register int i=head[u]; i; i=e[i].next)
	if (!vis[e[i].to])
	{
		vis[e[i].to]=true;
		if (!match[e[i].to] || dfs(match[e[i].to]))
		{
			match[e[i].to]=u;
			return 1;
		}	
	}
	return 0;
}

inline int solve(int p)
{
	int sum=0;
	for (register int i=1; i<=a; ++i) pd[i]=false;
	for (register int i=1; i<=a; ++i) 
	if ((1<<(i-1))&p) pd[i]=true;
	else sum++;
	cnt=0;
	for (register int i=1; i<=b; ++i) head[i]=0;
	for (register int i=1; i<=tot; ++i) if (pd[num[1][i]]) add(num[2][i],num[3][i]);
	for (register int i=1; i<=b; ++i) match[i]=0;
	for (register int i=1; i<=b; ++i)
	{
		for (register int j=1; j<=c; ++j) vis[j]=false;	
		sum+=dfs(i);
	}
	return sum;
}

int main(){
	scanf("%d",&T);
	while (T--)
	{
		ans=2e9;
		tot=0;
		scanf("%d%d%d",&a,&b,&c);
		int minn=min(a,min(b,c));
		for (register int i=1; i<=a; ++i)
		for (register int j=1; j<=b; ++j)
		for (register int k=1; k<=c; ++k)
		{
			scanf("%d",&x);
			if (!x) continue;
			num[1][++tot]=i;
			num[2][tot]=j;
			num[3][tot]=k;
		}
		if (minn==b) swap(a,b),swap(num[1],num[2]);
		else if (minn==c) swap(a,c),swap(num[1],num[3]);
		int max_statue=(1<<a)-1;
		for (register int p=0; p<=max_statue; ++p) ans=min(ans,solve(p));
		printf("%d\n",ans);
	}
return 0;	
}