[bzoj1047][ST表]理想的正方形

Description

　　有一个a*b的整数组成的矩阵，现请你从中找出一个n*n的正方形区域，使得该区域所有数中的最大值和最小值 的差最小。

Input

　　第一行为3个整数，分别表示a,b,n的值第二行至第a+1行每行为b个非负整数，表示矩阵中相应位置上的数。每 行相邻两数之间用一空格分隔。
100%的数据2<=a,b<=1000,n<=a,n<=b,n<=1000

Output

　　仅一个整数，为a*b矩阵中所有“n*n正方形区域中的最大整数和最小整数的差值”的最小值。

Sample Input

5 4 2

1 2 5 6

0 17 16 0

16 17 2 1

2 10 2 1

1 2 2 2

Sample Output

1

题解

原本想写个单调队列
机房有人说这题可以写2维ST表
于是我就写了个2维的ST表
mn[i][j][k],mx[i][j][k]表示(i,j)向左向上2^k的正方形的最小最大值
然后就A了。。
效率感人

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
int mn[1100][1100][11],mx[1100][1100][11];
int bin[11],Log[1100];
int a,b,n;
int mp[1100][1100];
int main()
{
//  freopen("square.in9","r",stdin);
    bin[0]=1;for(int i=1;i<=10;i++)bin[i]=bin[i-1]*2;
    scanf("%d%d%d",&a,&b,&n);
    Log[1]=0;for(int i=2;i<=max(a,b);i++)Log[i]=Log[i/2]+1; 
    for(int i=1;i<=a;i++)for(int j=1;j<=b;j++)scanf("%d",&mp[i][j]),mn[i][j][0]=mx[i][j][0]=mp[i][j];
    for(int i=1;i<=a;i++)
        for(int j=1;j<=b;j++)
            for(int k=1;bin[k]<=min(i,j);k++)
            {
                mn[i][j][k]=min(min(mn[i][j][k-1],mn[i-bin[k-1]][j][k-1]),min(mn[i][j-bin[k-1]][k-1],mn[i-bin[k-1]][j-bin[k-1]][k-1]));
                mx[i][j][k]=max(max(mx[i][j][k-1],mx[i-bin[k-1]][j][k-1]),max(mx[i][j-bin[k-1]][k-1],mx[i-bin[k-1]][j-bin[k-1]][k-1]));
            }
    int ans=999999999;
    for(int i=n;i<=a;i++)
        for(int j=n;j<=b;j++)
        {
            int k=Log[n];
        //  int ax=min(min(mn[i][j][k],mn[i-bin[k]+1][j][k]),min(mn[i][j-bin[k]+1][k],mn[i-bin[k]+1][j-bin[k]+1][k]));
        //  int ay=max(max(mx[i][j][k],mx[i-bin[k]+1][j][k]),max(mx[i][j-bin[k]+1][k],mx[i-bin[k]+1][j-bin[k]+1][k]));
            int ax=min(min(mn[i][j][k],mn[i-(n-bin[k])][j][k]),min(mn[i][j-(n-bin[k])][k],mn[i-(n-bin[k])][j-(n-bin[k])][k]));
            int ay=max(max(mx[i][j][k],mx[i-(n-bin[k])][j][k]),max(mx[i][j-(n-bin[k])][k],mx[i-(n-bin[k])][j-(n-bin[k])][k]));
            ans=min(ans,ay-ax);
        }
    printf("%d\n",ans);
    return 0;
}