[SPOJ MINSUB Largest Submatrix]单调栈

[SPOJ MINSUB Largest Submatrix]单调栈

知识点：data structure 单调栈

1. 题目链接

[SPOJ MINSUB Largest Submatrix]

2. 题意描述

给定一个$R*C$的矩阵，矩阵元素都是非负整数。要你找一个满足下列条件的一个子矩阵。
- 面积大于$K$
- 保证子矩阵的最小值最大的情况下，让面积也最大。

输出这个子矩阵的最小值和面积。
$(R, C\le 1000, K\le R*C, 矩阵元素\le 10^9)$

3. 解题思路

可以对子矩阵的最小值进行二分。然后令矩阵中所有小于当前二分的值的元素为0，否则为1。
然后就是一个单调栈的简单应用——求元素全为1的子矩阵的最大面积。
详细可以参考《 [51nod 1158 全是1的最大子矩阵] 单调栈》

4. 实现代码

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;
typedef long double LB;
typedef pair<int, int> PII;
typedef pair<LL, LL> PLL;

const int INF = 0x3f3f3f3f;
const LL INFL = 0x3f3f3f3f3f3f3f3fLL;
const LB eps = 1e-8;
const int MX = 1000 + 5;

int row, col, K;
int G[MX][MX], H[MX], L[MX], R[MX];
struct SNode {
    int val, pos;
};
SNode stk[MX]; int top;

int calc(int base) {
    int ret = 0;
    memset(H, 0, sizeof(H));
    for(int i = 1; i <= row; ++i) {
        for(int j = 1; j <= col; ++j) H[j] = (G[i][j] >= base) ? H[j] + 1 : 0;
        top = 0;
        for(int j = 1; j <= col; ++j) {
            L[j] = j;
            while(top && H[j] <= (stk[top - 1]).val) L[j] = stk[-- top].pos;
            stk[top ++] = SNode{H[j], L[j]};
        }

        top = 0;
        for(int j = col; j >= 1; --j) {
            R[j] = j;
            while(top && H[j] <= (stk[top - 1]).val) R[j] = stk[-- top].pos;
            stk[top ++] = SNode{H[j], R[j]};
        }

        for(int j = 1; j <= col; ++j) ret = max(ret, H[j] * (R[j] - L[j] + 1));
    }
    return ret;
}

int main() {
#ifdef ___LOCAL_WONZY___
    freopen("input.txt", "r", stdin);
#endif // ___LOCAL_WONZY___
    int _; scanf("%d", &_);
    while(_ --) {
        scanf("%d %d %d", &row, &col, &K);
        for(int i = 1; i <= row; ++i) {
            for(int j = 1; j <= col; ++j) {
                scanf("%d", &G[i][j]);
            }
        }
        int lb = 0, ub = INF, mid, ans = 0;
        while(lb <= ub) {
            mid = (lb + ub) >> 1;
            LL ret = calc(mid);
            if(ret >= K) lb = mid + 1, ans = mid;
            else ub = mid - 1;
        }
        printf("%d %d\n", ans, calc(ans));
    }
    return 0;
}