3140:[HNOI2013]消毒 - BZOJ

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

输入描述 Input Description
输入文件第一行是一个正整数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。
输出描述 Output Description
仅包含 D 行，每行一个整数，表示对应实验皿最少要用多少单位的 F 试剂。
样例输入 Sample Input
1
4 4 4
1 0 1 1
0 0 1 1
0 0 0 0
0 0 0 0
0 0 1 1
1 0 1 1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
样例输出 Sample Output
3
数据范围及提示 Data Size & Hint
对于区域(1,1,3)-(2,2,4)和(1,1,1)-(4,4,1)消毒，分别花费2个单位和1个单位的F试剂。

终于过了，fillchar没有一点优势啊
还是用各种for来把赋初值快一些

二维的怎么做大家都知道，就是二分图匹配，但是这是三维
怎么办呢？
！！！！强行转化成二维，枚举第三维的情况，那这需要2^k所以我们找最小的k
因为a*b*c<=5000，所以min{a,b,c}最大是17左右
然后开始枚举，但是还是很吃力所以对于把一个数组重新赋值的事最好避免

 var
     f:array[..,..]of longint;
     first,next,last,cover,bool:array[..]of longint;
     flag,sum:array[..]of boolean;
     a,b,c,tot,num,min,aug,ans,match,t,time,time2,maxa,maxb,maxc:longint;

 procedure swap(var x,y:longint);
 var
     t:longint;
 begin
     t:=x;
     x:=y;
     y:=t;
 end;

 procedure init;
 var
     i,j,k,s,ai,bi,ci:longint;
 begin
     read(ai,bi,ci);
     a:=;
     b:=;
     c:=;
     num:=;
     maxa:=;
     maxb:=;
     maxc:=;
     for i:= to ai do
       for j:= to bi do
         for k:= to ci do
           begin
             read(s);
             if s= then
             begin
               inc(num);
               f[num,]:=i;
               f[num,]:=j;
               f[num,]:=k;
               if i>maxa then maxa:=i;
               if j>maxb then maxb:=j;
               if k>maxc then maxc:=k;
             end;
           end;
     min:=ai;
     if bi<min then min:=bi;
     if ci<min then min:=ci;
     if (bi<=ai)and(bi<=ci) then
       begin
         swap(a,b);
         swap(maxa,maxb);
       end
     else
       if (ci<=ai)and(ci<=bi) then
       begin
         swap(a,c);
         swap(maxa,maxc);
       end;
     fillchar(sum,sizeof(sum),false);
     for i:= to num do
       sum[f[i,a]]:=true;
 end;

 function path(x:integer):boolean;
 var
     i:integer;
 begin
     i:=first[x];
     while i<> do
       begin
         if bool[last[i]]<>time2 then
           begin
             bool[last[i]]:=time2;
             if (cover[last[i]]=)or(path(cover[last[i]])) then
             begin
               cover[last[i]]:=x;
               exit(true);
             end;
           end;
         i:=next[i];
       end;
     exit(false);
 end;

 procedure insert(x,y:integer);
 begin
     inc(tot);
     last[tot]:=y;
     next[tot]:=first[x];
     first[x]:=tot;
 end;

 procedure work;
 var
     i:integer;
 begin
     if aug>ans then exit;
     tot:=;
     match:=;
     for i:= to maxb do
       first[i]:=;
     for i:= to maxc do
       cover[i]:=;
     for i:= to num do
       if flag[f[i,a]] then insert(f[i,b],f[i,c]);
     for i:= to maxb do
       begin
         inc(time2);
         if path(i) then inc(match);
       end;
     if ans>aug+match then ans:=aug+match;
 end;

 procedure try(x:integer);
 begin
     if x>min then work
     else
       begin
         flag[x]:=true;
         try(x+);
         if sum[x] then
         begin
           flag[x]:=false;
           inc(aug);
           try(x+);
           dec(aug);
         end;
       end;
 end;

 begin
     read(t);
     for time:= to t do
       begin
         init;
         ans:=;
         try();
         writeln(ans);
       end;
 end.