[BZOJ1047][HAOI2007]理想的正方形

Time Limit: 10 Sec
Memory Limit: 162 MB

Description

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

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

Output
仅一个整数，为ab矩阵中所有“nn正方形区域中的最大整数和最小整数的差值”的最小值。
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

题解：
ST表RMQ，但是由于查询范围是正方形那么我们就可以把ST表压掉一维,时间卡的挺紧的

#include
#define LiangJiaJun main
#define INF 2147483647
#define max(a,b) (a>b?a:b)
#define min(a,b) (a
using namespace std;
inline int read(){
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while('0'<=ch&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x;
}
int a,b,n;
int mm[1004];
int maxn[1004][1004][14],mint[1004][1004][14];
void pre(){
     int k,i,j;
     for(k=1;k<=mm[min(a,b)];k++){
         for(i=1;i+(1<<k)-1<=a;i++){
             for(j=1;j+(1<<k)-1<=b;j++){
                 maxn[i][j][k]=max(max(maxn[i][j][k-1],maxn[i+(1<<(k-1))][j][k-1]),
                                   max(maxn[i][j+(1<<(k-1))][k-1],maxn[i+(1<<(k-1))][j+(1<<(k-1))][k-1]));
                 mint[i][j][k]=min(min(mint[i][j][k-1],mint[i+(1<<(k-1))][j][k-1]),
                                   min(mint[i][j+(1<<(k-1))][k-1],mint[i+(1<<(k-1))][j+(1<<(k-1))][k-1]));
             }
         }
     }
}
int queryMax(int x1,int y1,int x2,int y2){
    int k=mm[n];
    x2=x2-(1<<k)+1;
    y2=y2-(1<<k)+1;
    return max(max(maxn[x1][y1][k],maxn[x1][y2][k]),max(maxn[x2][y1][k],maxn[x2][y2][k]));
}
int queryMin(int x1,int y1,int x2,int y2){
    int k=mm[n];
    x2=x2-(1<<k)+1;
    y2=y2-(1<<k)+1;
    return min(min(mint[x1][y1][k],mint[x1][y2][k]),min(mint[x2][y1][k],mint[x2][y2][k]));
}
int w33ha(){
    for(int i=1;i<=a;i++){
        for(int j=1;j<=b;j++){
            mint[i][j][0]=maxn[i][j][0]=read();
        }
    }
    pre();
    int ans=INF;
    for(int i=1;i+n-1<=a;i++){
        for(int j=1;j+n-1<=b;j++){
            ans=min(ans,queryMax(i,j,i+n-1,j+n-1)-queryMin(i,j,i+n-1,j+n-1));
        }
    }
    printf("%d\n",ans);
    return 0;
}
int LiangJiaJun(){
    mm[0]=-1;
    for(int i=1;i<=1001;i++){
        mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
    }
    while(scanf("%d%d%d",&a,&b,&n)!=EOF)w33ha();
    return 0;
}