POJ2559
題意：給你一個柱狀圖，讓你找出圖中能組成的矩形面積最大為多少。可轉化為對于一個序列，對于任意區間都有一個貢獻，區間最小值乘于區間的長度為該區間的貢獻，求這個序列貢獻最大值為多少。
思路：

對序列中每個元素 i ，以它為最小值向兩邊擴展，當前的貢獻值為擴展長度 len 乘 元素值 val, 求出最大的貢獻值。
如果要一直擴展下去，必須當前元素值不小于擴展元素值。
如果我們從第 1 個元素計算到第 n 個元素，計算到元素 i 時，
如果我們能快速的找出左邊第一個比 i 小的元素，就可以完成向左的擴展，
而右擴展，這個過程本來就是一直在向右走，所以直接擴展即可。
用單調遞減棧來維護這個過程，棧內儲存的為坐標i。
在每個元素入棧的時候，更新入棧元素向左擴展的最小位置
在每個元素出棧的時候，更新棧內元素向右擴展的最大位置

AC代碼:

#include <cstdio>
#include <iostream>
#include <stack>
#include <string.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> iip;
typedef pair<ll, ll> llp;
const int MAXN = 200005;
const int INF = 0x3f3f3f3f;
const int MOD = 998244353;

int i, j, n, a[MAXN], l[MAXN], r[MAXN], top;
ll ans;
stack<int> s;

int main()
{
    while(scanf("%d", &n) != EOF && n) {
        for(i = 0; i < n; i++) {
            scanf("%d", &a[i]);
        }
        ans = 0;
        for(i = 0; i < n; i++) {
            while(!s.empty() && a[i] < a[s.top()]) {
                r[s.top()] = i;
                s.pop();
            }
            if(s.empty())   l[i] = 0;
            else if(a[i] == a[s.top()])  l[i] = l[s.top()];
            else  l[i] = s.top() + 1;
            s.push(i);
        }
        while(!s.empty()) {
            r[s.top()] = i;
            s.pop();
        }
        for(i = 0; i < n; i++) {
            ans = max(ans, 1ll * a[i] * (r[i]-l[i]));
        }
        printf("%lld\n", ans);
//        cout.flush();
    }
    return 0;
}

POJ3494
題意：給你一個n*m的矩陣，矩陣的每個位置值為0或者1，問你在這個矩陣中全部由1組成的最大的矩形面積為多少。
思路：矩陣中的子矩陣等同于序列中的子序列，只不過此題要做一個預處理，將這個矩陣分為以每一行為x軸的n個柱狀圖，對于每一個等于1的點，它的高度都等于上一行同一列的點的高度加一。初始化后，再對n個柱狀圖進行如上題POJ2559的處理。
例如：

例子

初始化后

以第一行為 x 軸

以第二行為 x 軸

以第三行為 x 軸

AC代碼：

#include<cstdio>
#include<stack>
#include<string.h>
#include <iostream>
using namespace std;

typedef long long ll;
typedef pair<int, int> iip;
typedef pair<ll, ll> llp;
const int MAXN = 2005;
int i, j, n, m, a[MAXN][MAXN] = {0}, cnt[MAXN][MAXN] = {0}, ans = 0, l[MAXN], r[MAXN];
stack<int> s;

void init() {
    for(i = 1; i <= n; i++) {
        for(j = 1; j <= m; j++) {
            if(a[i][j] == 1) {
                cnt[i][j] = cnt[i-1][j] + 1;
            }
        }
    }
}

int main() {
    while(scanf("%d%d", &n, &m) != EOF) {
        memset(cnt, 0, sizeof(cnt));
        for(i = 1; i <= n; i++) {
            for(j = 1; j <= m; j++) {
                scanf("%d", &a[i][j]);
                if(a[i][j] == 1) {
                    cnt[i][j] = 1;
                }
            }
        }
        init();
        ans = 0;
        for(i = 1; i <= n; i++) {
            int tmp = 1;
            for(j = 1; j <= m; j++) {
                if(cnt[i][j] == 0) {
                    while(!s.empty()) {
                        r[s.top()] = j;
                        s.pop();
                    }
                    tmp = j+1;
                } else {
                    while(!s.empty() && cnt[i][j] < cnt[i][s.top()]) {
                        r[s.top()] = j;
                        s.pop();
                    }
                    if(s.empty())   l[j] = tmp;
                    else if(cnt[i][j] == cnt[i][s.top()])   l[j] = l[s.top()];
                    else l[j] = s.top() + 1;
                    s.push(j);
                }
            }
            while(!s.empty()) {
                r[s.top()] = j;
                s.pop();
            }
            for(j = 1; j <= m; j++) {
                ans = max(ans, cnt[i][j] * (r[j] - l[j]));
            }
         }
        printf("%d\n", ans);
//        cout.flush();
    }
    return 0;
}