[BZOJ2228][Zjoi2011][单调队列]礼物(gift)

题意

---

在一个p*q*r的正方体内找到一个a*a*b的子正方体，满足子正方体内均为N，求最大的4ab

---

考虑以一个面作为底面，可以正方体中每个点为左下角的正方形的边长（正方面平行于底面），然后考虑同一列中的点，可以把问题就转化成求一个序列，序列的最小值和区间长度相乘的结果最小。用单调队列做。
转动下正方体就可以计算所有情况啦

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#define N 210

using namespace std;

char w[N*N][N];
char A[N][N][N],B[N][N][N];
int p,q,r,t,Ans=0;
int bianc[N][N][N];
int Q[N],f[N];

inline void Checking(){
    memset(bianc,0,sizeof(bianc));
    for(int k=1;k<=r;k++)
        for(int i=1;i<=p;i++)
            for(int j=1;j<=q;j++)
                if(A[i][j][k]=='N')
                    bianc[i][j][k]=min(min(bianc[i][j-1][k],bianc[i-1][j][k]),bianc[i-1][j-1][k])+1;
    for(int i=1;i<=p;i++)
        for(int j=1;j<=q;j++){
            Q[t=0]=0; bianc[i][j][0]=-1;
            for(int k=1;k<=r;k++){
                while(t&&bianc[i][j][k]<=bianc[i][j][Q[t]]) t--;
                if(A[i][j][k]=='N') f[k]=k-Q[t]; else f[k]=0;
                Q[++t]=k;
            }
            Q[t=0]=r+1; bianc[i][j][r+1]=-1;
            for(int k=r;k;k--){
                while(t&&bianc[i][j][k]<=bianc[i][j][Q[t]]) t--;
                if(A[i][j][k]=='N') f[k]+=Q[t]-k;
                Q[++t]=k;
            }
            for(int k=1;k<=r;k++) Ans=max(Ans,4*bianc[i][j][k]*(f[k]-1));
        }
}

inline void Changing(){
    for(int i=1;i<=p;i++)
        for(int j=1;j<=q;j++)
            for(int k=1;k<=r;k++)
                B[j][k][i]=A[i][j][k];
    int k=p; p=q; q=r; r=k;
    for(int i=1;i<=p;i++)
        for(int j=1;j<=q;j++)
            for(int k=1;k<=r;k++)
                A[i][j][k]=B[i][j][k];
}

int main(){
    scanf("%d%d%d",&p,&q,&r);
    for(int i=1;i<=p*q;i++) scanf("%s",w[i]+1);
    for(int i=1;i<=p;i++)
        for(int j=1;j<=q;j++)
            for(int k=1;k<=r;k++) A[i][j][k]=w[j*p+i-p][k];
    Checking();
    Changing();
    Checking();
    Changing();
    Checking();
    printf("%d\n",Ans);
    return 0;
}

跑的好慢啊……