【DP~最大子立方体】吃西瓜

 吃西瓜[matrix.pas/c/cpp]

[说明]此题中出现的所有数全为整数

[背景]SubRaY有一天得到一块西瓜,是长方体形的....

[题目描述]SubRaY发现这块西瓜长m厘米,宽n厘米,高h厘米.他发现如果把这块西瓜平均地分成m*n*h块1立方厘米的小正方体,那么每一小块都会有一个营养值(可能为负,因为西瓜是有可能坏掉的,但是绝对值不超过200).

现在SubRaY决定从这m*n*h立方厘米的西瓜中切出mm*nn*hh立方厘米的一块小西瓜(一定是立方体形,长宽高均为整数),然后吃掉它.他想知道他最多能获得多少营养值.(0<=mm<=m,0<=nn<=n,0<=hh<=h.mm,nn,hh的值由您来决定).

换句话说,我们希望从一个m*n*h的三维矩阵中,找出一个三维子矩阵,这个子矩阵的权和最大.

一个2*3*4的例子,最优方案为切红色2*3*1部分

[输入][matrix.in]

首行三个数h,m,n(注意顺序),分别表示西瓜的高,长,宽.

以下h部分,每部分是一个m*n的矩阵,第i部分第j行的第k个数表示西瓜第i层,第j行第k列的那块1立方厘米的小正方体的营养值.

 

[输出][matrix.out]

SubRaY所能得到的最大营养值

 

[样例输入]

2 3 4

4 1 2 8

0 5 -48 4

3 0 1 9

2 1 4 9

1 0 1 7

3 1 2 8

 

[样例输出]

45

 

[数据范围]

对于30%的数据,h=1,1<=m,n<=10

对于全部的数据,1<=h<=32,1<=m,n<=50,保证h<=m,n

 =====================================

============================

{ID:jie19952PROG:LANG:PASCAL}var  h,m,n:longint;  sum,map:array[0..50,0..50,0..50]of longint;  sum1:array[0..50,0..50]of longint;  sum2:array[0..50]of longint;procedure init;begin  assign(input,'matrix.in');  assign(output,'matrix.out');  reset(input); rewrite(output);end;procedure terminate;begin  close(input); close(output);  halt;end;function max(a,b:longint):longint;begin  if a>b then exit(a);  exit(b);end;procedure main;var  h1,m1,m2,n1,n2:longint;  ans:longint;begin  readln(h,m,n);  fillchar(map,sizeof(map),0);  fillchar(sum,sizeof(sum),0);    for h1:=1 to h do    for m1:=1 to m do      for n1:=1 to n do        begin          read(map[h1,m1,n1]);          sum[h1,m1,n1]:=sum[h1,m1,n1-1]+map[h1,m1,n1];        end;          ans:=-maxlongint;  For n1:=1 to n do    for n2:=n1 to n do      begin        //fillchar(sum1,sizeof(sum1),0);        for h1:=1 to h do          begin            sum1[h1,0]:=0;            for m1:=1 to m do              begin                sum1[h1,m1]:=sum1[h1,m1-1]+sum[h1,m1,n2]-sum[h1,m1,n1-1];              end;          end;        sum2[0]:=0;        for m1:=1 to m do          for m2:=m1 to m do            begin              for h1:=1 to h do                begin                  sum2[h1]:=max(sum2[h1-1]+sum1[h1,m2]-sum1[h1,m1-1],sum1[h1,m2]-sum1[h1,m1-1]);                  if sum2[h1]>ans then ans:=sum2[h1];                end;            end;      end;  writeln(ans);end;begin  init;  main;  terminate;end.