HDU 2870 Largest Submatrix DP求最大子矩阵

题目:http://acm.hdu.edu.cn/showproblem.php?pid=2870

题意:给出一个字符矩阵Matritx[][], 求其最大子矩阵, 详细如题.

解题思路:

　　将全部字符依次转化a, b, c, 再分别求出这三个矩阵的最大子矩阵即可.

　　于是, 问题转化为求矩阵中最大的子矩阵了.

　　设置一个变量Num[][]记录位置的最大高度, Num[i][j]表示Matritx[i][j]位置上的最大高度

　　这样, 只要枚举以各个Num[i][j]为矩阵最小高度, 分别向前后推进扩展矩阵, 如果Num[i][j + 1] >= Num[i][j]则可以向前扩展, 同理Num[i][j - 1] >= Num[i][j]则可以向后扩展, 用H[j], T[j]分别表示当前位置 j 向前向后扩展的最大位置

　　这样, 可以得到一个DP式, 用以优化算法计算过程

　　　　H[j] = H[H[j] - 1], if Num[i][j] <= Num[i][H[j] - 1];

　　　　T[j] = T[T[j] + 1], if Num[i][j] <= Num[i][T[j] + 1];

　　那么, 以Num[i][j]为最小高度的最大矩阵为 Num[i][j] * (T[j] - H[j] + 1), 枚举所有的i, j, 取最大解.

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#include <vector>
#include <map>

using namespace std;

struct TT
{
    int a;
    int b;
    int c;
}M[300];
const int MAX = 1000 + 10;
int n, m;
int Num[3][MAX][MAX];
int H[MAX], T[MAX];

void init()
{
    for(int i = 0; i <= n; ++i)
    {
        for(int j = 1; j <= m; ++j)
        {
            Num[0][i][j] = 0;
            Num[1][i][j] = 0;
            Num[2][i][j] = 0;
        }
    }
}

int solve()
{
    int Ans = 0;
    for(int i = 0; i < 3; ++i)
    {
        for(int j = 1; j <= n; ++j)
        {
            Num[i][j][0] = Num[i][j][m + 1] = -1;
            for(int k = 1; k <= m; ++k)
                H[k] = T[k] = k;
            for(int k = 2; k <= m; ++k)
            {
                while(Num[i][j][k] <= Num[i][j][H[k] - 1])
                    H[k] = H[H[k] - 1];
            }
            for(int k = m - 1; k >= 1; --k)
            {
                while(Num[i][j][k] <= Num[i][j][T[k] + 1])
                    T[k] = T[T[k] + 1];
            }
            for(int k = 1; k <= m; ++k)
                if(Num[i][j][k] * (T[k] - H[k] + 1) > Ans)
                    Ans = Num[i][j][k] * (T[k] - H[k] + 1);
        }
    }
    return Ans;
}

bool isValid(char c)
{
    if(c == 'a' || c == 'b' || c == 'c' ||
       c == 'w' || c == 'x' || c == 'y' ||
       c == 'z')
        return true;
    return false;
}
int main()
{
    freopen("in.txt","r",stdin);
    char c[MAX];
    M['a'].a = M['a'].b = M['a'].c = 0;
    M['b'].a = M['b'].b = M['b'].c = 1;
    M['c'].a = M['c'].b = M['c'].c = 2;
    M['w'].a = M['w'].c = 0;
    M['w'].b = 1;
    M['x'].a = M['x'].c = 1;
    M['x'].b = 2;
    M['y'].a = M['y'].c = 0;
    M['y'].b = 2;
    M['z'].a = 0;
    M['z'].b = 1;
    M['z'].c = 2;

    while(scanf("%d%d", &n, &m) != EOF)
    {
        init();
        for(int i = 1; i <= n; ++i)
        {
            scanf("%s", c);
            for(int j = 0; j < m; ++j)
            {
                if(!isValid(c[j]))
                    continue;
                Num[M[c[j]].a][i][j + 1] = Num[M[c[j]].a][i - 1][j + 1] + 1;
                Num[M[c[j]].b][i][j + 1] = Num[M[c[j]].b][i - 1][j + 1] + 1;
                Num[M[c[j]].c][i][j + 1] = Num[M[c[j]].c][i - 1][j + 1] + 1;
            }
        }
        printf("%d\n", solve());
    }
    return 0;
}

转载自： http://www.cnblogs.com/Kenfly/archive/2011/04/18/2020110.html