hdu1505 City game(最大子矩阵)

题意：

Bob打算在城市中寻找一块长方形的区域，已知城市中有些地已经被占用，要求我们帮助他选出最大的长方形区域。R代表已被占用，F代表未被占用。

分析：

这是hdu1506的升级版，即从一维变成了二维。但我们依旧可以把它转换成一维，即以每一列为一个底，算出此列中长方形的高度，然后依次按行(即按高度)dp，找出最大值。

—————>

在每一行找出当前行中第j个点向左向右延伸的最大长度，然后乘以高度即是面积。
注意最终答案要乘以3！！！！
状态转移方程：
if(h[j-1] >= h[j]) L[j] = L[j-1];
if(h[j+1] >= h[j]) R[j] = R[j-1];

具体代码：

#include <set>
#include <numeric>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cctype>
#include <string>
#include <sstream>
#include <map>
#include <functional>

using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
#define REP(idx1,num1) for(int idx1=0;idx1<(num1);idx1++)
#define pb push_back
#define empb emplace_back
#define mp make_pair
#define mem(s) memset(s,0,sizeof(s));
const double EPS = 1e-6;
const int maxn = 1000 + 10;
int d[maxn][maxn];
int l[maxn],r[maxn];
int r1,c1;
void debug()
{
    for(int i = 0; i < r1; ++i){
        for(int j = 0; j < c1; ++j)
            cout <<d[i][j]<<" ";
        cout<<endl;
    }
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&r1,&c1);
        mem(d);
        for(int i = 0; i < r1; ++i)
            for(int j = 0; j < c1; ++j){
                char ch;cin >> ch;
                //cout << ch << endl;
                if(ch == 'F') d[i][j] = (i == 0) ? 1 : d[i-1][j]+1;
            }
        //debug();
        int max_area = 0;
        //从每一行开始
        for(int i = 0; i < r1; ++i){
            for(int j = 0; j < c1; ++j){
                int t = j;
                while(t >= 1 && d[i][t-1] >= d[i][j])
                    t = l[t-1];
                l[j] = t;
            }
            for(int j = c1-1; j >= 0; --j){
                int t = j;
                while(t < c1-1 && d[i][t+1] >= d[i][j])
                    t = r[t+1];
                r[j] = t;
            }
            /*
            for(int j = 0; j < c1; ++j)
                cout << l[j] << " ";
            cout << endl;
            */
            for(int j = 0; j < c1; ++j){
                int area = (r[j] - l[j] + 1) * d[i][j];
                if(max_area < area) max_area = area;
            }
        }
        printf("%d\n",max_area*3);
    }
    return 0;
}