1 /** 这道题甚是巧妙啊，也是看了别人的代码才A的……，我相信如果你看了以后也会茅塞顿开的 */
 2 #include <iostream>
 3 #include <cstdio>
 4 using namespace std;
 5 #define N 1005
 6
 7 int head1, tail1, head2, tail2, qMax[N][N], qMin[N][N], matrix[N][N];
 8
 9 struct node{
10     int v, id;
11 }q1[N], q2[N];
12
13 void addMax( int v, int id ) {  //维护队列最大值
14     while( head1 < tail1 && q1[tail1-1].v <= v ) tail1--;
15     q1[tail1].v = v; q1[tail1++].id = id;
16 }
17
18 int getMaxNum( int id ) {
19     while( head1 < tail1 && q1[head1].id < id ) head1++;
20     return q1[head1].v;
21 }
22
23 void addMin( int v, int id ) {
24     while( head2 < tail2 && q2[tail2-1].v >= v ) tail2--;
25     q2[tail2].v = v; q2[tail2++].id = id;
26 }
27
28 int getMinNum( int id ) {
29     while( head2 < tail2 && q2[head2].id < id ) head2++;
30     return q2[head2].v;
31 }
32
33 int main() {
34     int n, m, r, c, i, j;
35     while(scanf("%d%d%d%d", &n, &m, &r, &c) != EOF) {
36         for( i = 0; i < n; ++i )
37         for( j = 0; j < m; ++j ) scanf("%d", &matrix[i][j]);
38
39         for( j = 0; j < m; ++j ) {
40             head1 = tail1 = head2 = tail2 = 0;
41             for( i = 0; i < r - 1; ++i ) {
42                 addMax( matrix[i][j], i );
43                 addMin( matrix[i][j], i );
44             }
45             /** 这里qMax[i][j]表示第j列, max( matrix[ i->(i+r) ][j] )中的最大值
46             *   qMin[i][j],同上。
47             */
48             for( i = r - 1; i < n; ++i ) {
49                 addMax( matrix[i][j], i );
50                 qMax[i-r+1][j] = getMaxNum( i-r+1 );
51                 addMin( matrix[i][j], i );
52                 qMin[i-r+1][j] = getMinNum( i-r+1 );
53             }
54         }
55         /** 这里qMax[i][j]表示第j列, max( qMax[i][ j->(j+c) ] )中的最大值
56         *   即qMax[i][j]表示以(i,j)为左上角的一个r*c子矩阵中的最大值
57         *   qMin[i][j],同上。
58         */
59         for( i = 0; i < n - c + 1; ++i ) {
60             head1 = tail1 = head2 = tail2 = 0;
61             for( j = 0; j < c - 1; ++j ) {
62                 addMax( qMax[i][j], j );
63                 addMin( qMin[i][j], j );
64             }
65             for( j = c - 1; j < m; ++j ) {
66                 addMax( qMax[i][j], j );
67                 qMax[i][j-c+1] = getMaxNum(j-c+1);
68                 addMin( qMin[i][j], j );
69                 qMin[i][j-c+1] = getMinNum(j-c+1);
70             }
71         }
72         int res = qMax[0][0] - qMin[0][0];
73         for( i = 0; i < n - r + 1; ++i )
74         for( j = 0; j < m - c + 1; ++j )
75         res = res > qMax[i][j] - qMin[i][j] ? res : qMax[i][j] - qMin[i][j];
76         printf("%d\n", res);
77     }
78     return 0;
79 }