POJ 2019: Cornfields

题目链接：http://poj.org/problem?id=2019

题目大意：

给定一个N*N的矩阵，

查询若干个B*B的矩阵中最大值和最小值的差值。

算法：

由于对于每个矩阵，查询的子矩阵的大小是固定的。

所以一看就是基础数据结构的练习。

不过我暂时没空管这个，我做这道题是为了写一个二维RMQ的模板

二维RMQ的原理也很简单。

首先像一维RMQ一样，处理出以每个点为左上角的大小为(2^i)*(2^j)的矩阵中的最值。

查询的时候，找出使得(2^k<b)的最大的k。

然后就可以把b*b的矩阵划分成4个边长为2^k的小矩阵(有重叠)

代码如下：

#include<cstdio>
#include<iostream>
#include<sstream>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
#include<set>
#include<map>
#define INF 0x3f3f3f3f
#define eps 1e-8
using namespace std;

const int MAXN = 250;
const int MAXD = 8;

int mx[MAXN][MAXN];
int rmq1[MAXN][MAXN][MAXD][MAXD];
int rmq2[MAXN][MAXN][MAXD][MAXD];

inline int cal(int x)
{
    int ret = 0;
    while((1 << ret) <= x)
    {
        ret ++;
    }
    return ret - 1;
}

int main()
{
    int n, b, q;
    scanf("%d %d %d", &n, &b, &q);
        for(int i = 0; i < n; i ++)
        {
            for (int j = 0; j < n; j ++)
            {
                scanf("%d", &mx[i][j]);
                rmq1[i][j][0][0] = rmq2[i][j][0][0] = mx[i][j];
            }
        }
        for (int i = n - 1; i >= 0; i --)
        {
            for (int j = n - 1; j >= 0; j --)
            {
                for (int k = 0; i + (1 << k) <= n; k ++)
                {
                    for (int l = 0; j + (1 << l) <= n; l ++)
                    {
                        if (k)
                        {
                            rmq1[i][j][k][l] = min(rmq1[i][j][k - 1][l], rmq1[i + (1 << (k - 1))][j][k - 1][l]);
                            rmq2[i][j][k][l] = max(rmq2[i][j][k - 1][l], rmq2[i + (1 << (k - 1))][j][k - 1][l]);
                        }
                        else if (l)
                        {
                            rmq1[i][j][k][l] = min(rmq1[i][j][k][l - 1], rmq1[i][j + (1 << (l - 1))][k][l - 1]);
                            rmq2[i][j][k][l] = max(rmq2[i][j][k][l - 1], rmq2[i][j + (1 << (l - 1))][k][l - 1]);
                        }
                    }
                }
            }
        }
        int k, l;
        k = l = cal(b);
        while(q --)
        {
            int i, j;
            scanf("%d%d", &i, &j);
            i --;
            j --;
            int ans1 = INT_MAX, ans2 = INT_MIN;
            ans1 = min(ans1, rmq1[i][j][k][l]);
            ans1 = min(ans1, rmq1[i + b - (1 << k)][j][k][l]);
            ans1 = min(ans1, rmq1[i][j + b - (1 << l)][k][l]);
            ans1 = min(ans1, rmq1[i + b - (1 << k)][j + b - (1 << l)][k][l]);
            ans2 = max(ans2, rmq2[i][j][k][l]);
            ans2 = max(ans2, rmq2[i + b - (1 << k)][j][k][l]);
            ans2 = max(ans2, rmq2[i][j + b - (1 << l)][k][l]);
            ans2 = max(ans2, rmq2[i + b - (1 << k)][j + b - (1 << l)][k][l]);
            printf("%d\n", ans2 - ans1);
        }
    return 0;
}