hdu3630

/*
分析：
    二维RMQ。
    对于所求的区域内、以每一个非lake点所在行为底部、先向上使劲儿扩展、然后再努力向
两侧扩展、可以得到一个不含lake的矩形。枚举出所有的这种矩形（必定涵盖ans情况）、最佳
的既为ans。n*m*q次枚举、RMQ询问O(1)、sum询问O(1)，时间复杂度为O(n*m*q)。
    感谢u010228612的提醒，这题矩阵中的每个node的值一直不被修改、所以sum是不变的、
是可以预处理的，我竟然在用二维树状数组维护一个动态矩阵、活生生乘了俩log(n)、更犀利
的是经过一个晚上加一个早上的优化，竟然让我用1400MS给a了， 。。。
   下面是改过的、312MS：。。。


                                                                        2013-07-02
*/

#include"iostream"
#include"cstdio"
#include"cmath"
#include"cstring"
using namespace std;
const int N=155;

int n,m,q,s;
int dp[N][N][9][9],le[N][N],rg[N][N],sum[N][N];
int lowbit[N];
inline int max(int a,int b){
    return a>b?a:b;
}

void RMQ()
{
    int j1,j2,row,col;
    int mx=int(log(1.0*n)/log(2.0));
    int my=int(log(1.0*m)/log(2.0));
    for(j1=0;j1<=mx;j1++)
    {
        for(j2=0;j2<=my;j2++)
        {
            if(!j1 && !j2)    continue;
            for(row=1;row+(1<<j1)-1<=n;row++)
            {
                for(col=1;col+(1<<j2)-1<=m;col++)
                {
                    if(!j1)
                        dp[row][col][j1][j2]=max(dp[row][col][j1][j2-1],dp[row][col+(1<<(j2-1))][j1][j2-1]);
                    else
                        dp[row][col][j1][j2]=max(dp[row][col][j1-1][j2],dp[row+(1<<(j1-1))][col][j1-1][j2]);
                }
            }
        }
    }
}
int power[N];
inline int query(int x1,int y1,int x2,int y2)
{
    int kx=power[x2-x1+1];
    int ky=power[y2-y1+1];
    int m1=dp[x1][y1][kx][ky];
    int m2=dp[x2-(1<<kx)+1][y1][kx][ky];
    int m3=dp[x1][y2-(1<<ky)+1][kx][ky];
    int m4=dp[x2-(1<<kx)+1][y2-(1<<ky)+1][kx][ky];
    return max(max(m1,m2),max(m3,m4));
}
void init()
{
    int i,l,temp;
    for(i=0;i<=150;i++) power[i]=log(double(i))/log(2.0);
    for(i=1;i<=n;i++)
    {
        temp=0;
        for(l=1;l<=m;l++)
        {
            scanf("%d",&dp[i][l][0][0]);
            temp+=dp[i][l][0][0];
            sum[i][l]=sum[i-1][l]+temp;
        }
    }

    for(i=1;i<=n;i++)   {le[i][0]=1;rg[i][m+1]=m;}
    for(i=1;i<=n;i++)
    {
        for(l=1;l<=m;l++)
        {
            if(dp[i][l][0][0]==-1)    le[i][l]=l+1;
            else    le[i][l]=le[i][l-1];
        }
        for(l=m;l>=1;l--)
        {
            if(dp[i][l][0][0]==-1)    rg[i][l]=l-1;
            else    rg[i][l]=rg[i][l+1];
        }
    }
    RMQ();
}
int main()
{
    int T,Case;
    int i,l;
    int up;
    int x1,y1,x2,y2,lmax,rmin;
    int t,temp,ans;
    cin>>T;
    for(i=1;i<=150;i++) lowbit[i]=i&(-i);
    for(Case=1;Case<=T;Case++)
    {
        if(Case>1)  printf("\n");
        printf("Case %d:\n",Case);
        scanf("%d%d%d%d",&n,&m,&q,&s);
        init();
        while(q--)
        {
            ans=0;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            for(l=y1;l<=y2;l++)
            {
                up=x1;
                lmax=y1;
                rmin=y2;
                for(i=x1;i<=x2;i++)
                {
                    if(le[i][l]==l+1)   {up=i+1;lmax=y1;rmin=y2;continue;}
                    if(lmax<le[i][l])   lmax=le[i][l];
                    if(rmin>rg[i][l])   rmin=rg[i][l];
                    temp=query(up,lmax,i,rmin);
                    t=(s-1)*temp+sum[i][rmin]-sum[i][lmax-1]-sum[up-1][rmin]+sum[up-1][lmax-1];
                    if(ans<t)   ans=t;
                }
            }
            printf("%d\n",ans);
        }
    }
    return 0;
}