[Zjoi2011]礼物(gift)

2228: [Zjoi2011]礼物(gift)
Time Limit: 10 Sec  Memory Limit: 259 MB
Submit: 32  Solved: 19
[Submit][Status][Discuss]
Description
小白的生日就要到了，小蓝决定送一件自己亲手做的手工艺品使自己的礼物与众不同。具体来说，小蓝已经通过某种方式制作出了一个p×q×r的木块（由pqr个单位小木块组成）。但由于小蓝手艺不精，现在这个木块中的有些单位小木块是有问题的(有裂缝、里面是空心等等)，这样的礼物小蓝是不可能直接送出去的。于是小蓝决定在这个木块中再挖出一个a×a×b的子木块(即要求挖出的长方体木块存在两条长度相等的相邻边)，当然这个子木块中是不能包含有问题的单位小木块的。为了使这个木块上能包含更多的图案，小蓝希望从所有可行的方案中挑取4ab的值最大的方案。但小蓝光检测木块中哪些地方有问题就已经耗尽了体力，作为小蓝的好友，你能帮帮小蓝吗？
Input
每个输入文件中仅包含一个测试数据。第一行包含三个由空格隔开的正整数，p,q,r。 接下来有pq行，每行包含r个字符，每个字符只可能是’P’(Poor)或者’N’(Nice)，表示该单位小木块有问题或者没问题。具体的说，第1+(yp+x-p)行的第z个字符描述的是坐标为(x,y,z)的小木块情况。(1<=x<=p,1<=y<=q,1<=z<=r)
Output
输出文件仅包含一个整数，表示最佳方案的4ab的值。
Sample Input
3 2 5
PNNNN
PNNNN
NPPNP
PNNNP
NNNNP
PPNNP

Sample Output
24

【数据范围】
对于100%的数据，0
HINT

Source
day1

 

题意:

这道题并不难,综合了 家的范围 http://hi.baidu.com/hccz95/blog/item/10f7fff9533ba30609244df9.html 和 zoj1985两题

详见程序

/************************************************************** Problem: 2228 User: hccz95 Language: Pascal Result: Accepted Time:6560 ms Memory:24028 kb ****************************************************************/ uses math; {$r-,s-,q-} const maxn = 150 ; type garr = array[ 0..maxn,0..maxn,0..maxn ] of Char ; var wid , hei , len , ans : Longint ; c , L , R , stk : array[ 0..maxn ] of Longint ; g : array[ 0..2 ] of Garr ; f : array[ 0..maxn,0..maxn,0..maxn ] of Longint ; function get( n : Longint ) : Longint ; var i , top : Longint ; begin top := 0 ; stk[ top ] := 0 ; for i := 1 to n do begin while ( c[ i ] <= c[ stk[ top ] ] ) and ( top > 0 ) do dec( top ); L[ i ] := i - stk[ top ] ; inc( top ); stk[ top ] := i ; end; top := 0 ; stk[ top ] := n + 1 ; for i := n downto 1 do begin while ( c[ i ] <= c[ stk[ top ] ] ) and ( top > 0 ) do dec( top ); R[ i ] := stk[ top ] - i ; inc( top ); stk[ top ] := i ; end; top := 0 ; for i := 1 to n do top := max( top , c[ i ] * ( L[ i ] + R[ i ] - 1 ) ); get := top ; end; procedure work( p , q , r : Longint ; var g : Garr ); var i , j , k : Longint ; begin fillchar( f , sizeof( f ) , 0 ); for i := 1 to p do for j := 1 to q do for k := 1 to r do if g[ i ][ j ][ k ] = 'N' then f[ i ][ j ][ k ] := min( f[ i ][ j - 1 ][ k - 1 ] , min( f[ i ][ j - 1 ][ k ] , f[ i ][ j ][ k - 1 ] ) ) + 1 ; for j := 1 to q do for k := 1 to r do begin for i := 1 to p do c[ i ] := f[ i ][ j ][ k ] ; ans := max( ans , get( p ) ); end; end; procedure main; var i , j , k : Longint ; begin readln( wid , hei , len ); for i := 1 to hei do for j := 1 to wid do begin for k := 1 to len do begin read( g[ 0 ][ i ][ j ][ k ] ); g[ 1 ][ j ][ i ][ k ] := g[ 0 ][ i ][ j ][ k ] ; g[ 2 ][ k ][ i ][ j ] := g[ 0 ][ i ][ j ][ k ] ; end; readln; end; ans := 0 ; work( hei , wid , len , g[ 0 ] ); work( wid , hei , len , g[ 1 ] ); work( len , hei , wid , g[ 2 ] ); writeln( ans * 4 ); end; begin main; end.