FZU 2080 最大差值 二维单调队列（甚是巧妙……）

/** 这道题甚是巧妙啊，也是看了别人的代码才A的……，我相信如果你看了以后也会茅塞顿开的 */
#include <iostream>
#include <cstdio>
using namespace std;
#define N 1005

int head1, tail1, head2, tail2, qMax[N][N], qMin[N][N], matrix[N][N];

struct node{
    int v, id;
}q1[N], q2[N];

void addMax( int v, int id ) {  //维护队列最大值
    while( head1 < tail1 && q1[tail1-1].v <= v ) tail1--;
    q1[tail1].v = v; q1[tail1++].id = id;
}

int getMaxNum( int id ) {
    while( head1 < tail1 && q1[head1].id < id ) head1++;
    return q1[head1].v;
}

void addMin( int v, int id ) {
    while( head2 < tail2 && q2[tail2-1].v >= v ) tail2--;
    q2[tail2].v = v; q2[tail2++].id = id;
}

int getMinNum( int id ) {
    while( head2 < tail2 && q2[head2].id < id ) head2++;
    return q2[head2].v;
}

int main() {
    int n, m, r, c, i, j;
    while(scanf("%d%d%d%d", &n, &m, &r, &c) != EOF) {
        for( i = 0; i < n; ++i )
        for( j = 0; j < m; ++j ) scanf("%d", &matrix[i][j]);

        for( j = 0; j < m; ++j ) {
            head1 = tail1 = head2 = tail2 = 0;
            for( i = 0; i < r - 1; ++i ) {
                addMax( matrix[i][j], i );
                addMin( matrix[i][j], i );
            }
            /** 这里qMax[i][j]表示第j列, max( matrix[ i->(i+r) ][j] )中的最大值 
            *   qMin[i][j],同上。
            */
            for( i = r - 1; i < n; ++i ) {
                addMax( matrix[i][j], i );
                qMax[i-r+1][j] = getMaxNum( i-r+1 );
                addMin( matrix[i][j], i );
                qMin[i-r+1][j] = getMinNum( i-r+1 );
            }
        }
        /** 这里qMax[i][j]表示第j列, max( qMax[i][ j->(j+c) ] )中的最大值
        *   即qMax[i][j]表示以(i,j)为左上角的一个r*c子矩阵中的最大值 
        *   qMin[i][j],同上。
        */
        for( i = 0; i < n - c + 1; ++i ) {
            head1 = tail1 = head2 = tail2 = 0;
            for( j = 0; j < c - 1; ++j ) {
                addMax( qMax[i][j], j );
                addMin( qMin[i][j], j );
            }
            for( j = c - 1; j < m; ++j ) {
                addMax( qMax[i][j], j );
                qMax[i][j-c+1] = getMaxNum(j-c+1);
                addMin( qMin[i][j], j );
                qMin[i][j-c+1] = getMinNum(j-c+1);
            }
        }
        int res = qMax[0][0] - qMin[0][0];
        for( i = 0; i < n - r + 1; ++i )
        for( j = 0; j < m - c + 1; ++j )
        res = res > qMax[i][j] - qMin[i][j] ? res : qMax[i][j] - qMin[i][j];
        printf("%d
", res);
    }
    return 0;
}