hdu2870子串长度子串最小数

//本题大意:给出一个含有w,x,y,z,a,b,c的矩阵.因为其中w,x,y,z可以转换为a或者b或者c.要求找一个子矩阵.里面全是
//         相同的元素.输出这个最大值.
//解题思路:这个题分为a,b,c来枚举以后就变成类似hdu1505,1506的类型.在一维上面来说.就是解决在一串数字(正数)
//         中寻找一个最大的子串.子串的价值=子串长度*子串中最小的数.
//         暴力做法就是枚举出所有i->j的最小值*j-i+1即可..在这里利用了>符号的传递性.
//         即如果A>B&&B>C那么A>C.  所以以每一个元素为基准..即认定a[i]是所在子串的最小值的话.
//         只需找出该子串中所含的元素个数..就是a[i]左边和右边连续的>=a[i]的个数.
//         用l[i]记录左边满足条件的元素的最右边的位置.r[i]记录右边满足条件的元素的最左边的位置
//         为了避免挨个遍历..在确定l[i]的时候.只需去和l[i-1]比较.如果成立.那么l[i] = l[i-1]
//         即是l[i]的位置至少在l[i-1]左边u后者就是l[i-1].这样递推下去即可..
//题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=2870
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<iostream>
#include<math.h>
#include<algorithm>
#include<vector>
#include<queue>
#include<functional>
using namespace std;
const int mn = 2110;
const int inf = 10000;
int ps[mn][mn];
int l[mn],r[mn],a[mn];
int solve(int row,int col,int cas[][mn])
{
    memset(ps,0,sizeof(ps));
    for(int i = 0 ; i < row ; i++)
    {
        ps[i+1][1] = cas[i][0];
        for(int j = 1 ; j < col ; j++)
        {
            if(ps[i+1][j] + cas[i][j] < 0)
                ps[i+1][j+1] = cas[i][j];
            else ps[i+1][j+1] = ps[i+1][j] + cas[i][j];
        }
    }
    int ans = 0;
    for(int i = 1 ; i <= col ; i++)
    {
        int t = 1;
        for(int j = 1 ; j <= row ; j++)
        {
            if(ps[j][i] < 0)
            {

                int p = 0;
                for(int k = t ; k < j ; k++)
                {
                    p++;
                    a[p] = ps[k][i];
                }
                for(int k = 1 ; k <= p ; k++)
                {
                    r[k] = k;
                    l[k] = k;
                }
                a[0] = -1;
                a[p+1] = -1;
                for(int k = 1 ; k <= p ; k++)
                    while(a[l[k]-1] >= a[k])
                        l[k] = l[l[k]-1];
                for(int k = p ; k >= 1 ; k--)
                    while(a[r[k]+1] >= a[k])
                        r[k] = r[r[k]+1];
                for(int k = 1 ; k <= p ; k++)
                    ans = max(ans,(r[k]-l[k]+1)*a[k]);
                t = j+1;
            }
        }
        int p = 0;
        for(int k = t ; k <= row ; k++)
        {
            p++;
            a[p] = ps[k][i];
        }
        for(int k = 1 ; k <= p; k++)
        {
            r[k] = k;
            l[k] = k;
        }
        a[0] = -1;
        a[p+1] = -1;
        for(int k = 1 ; k <= p ; k++)
         while(a[l[k]-1] >= a[k])
          l[k] = l[l[k]-1];
        for(int k = p ; k >= 1 ; k--)
         while(a[r[k]+1] >= a[k])
          r[k] = r[r[k]+1];
        for(int k = 1 ; k <= p ; k++)
          ans = max(ans,(r[k]-l[k]+1)*a[k]);
    }
    return ans;
}
int cas[mn][mn];
char map[mn][mn];
int main()
{
    int row,col;
    while(scanf("%d%d",&row,&col) != EOF)
    {
        for(int i = 0 ; i < row ; i++)scanf("%s",map[i]);
        for(int i = 0 ; i < row ; i++)
            for(int j = 0 ; j < col ; j++)
                if(map[i][j] == 'x' || map[i][j] == 'b' || map[i][j] == 'c')
                    cas[i][j] = -inf;
                else cas[i][j] = 1;
        int res = 0;
        res = max(res,solve(row,col,cas));
        for(int i = 0 ; i < row ; i++)
            for(int j = 0 ; j < col ; j++)
                if(map[i][j] == 'y' || map[i][j] == 'a' || map[i][j] == 'c')
                    cas[i][j] = -inf;
                else cas[i][j] = 1;
        res = max(res,solve(row,col,cas));
        for(int i = 0 ; i < row ; i++)
            for(int j = 0 ; j < col ; j++)
                if(map[i][j] == 'w' || map[i][j] == 'b' || map[i][j] == 'a')
                    cas[i][j] = -inf;
                else cas[i][j] = 1;
        res = max(res,solve(row,col,cas));
        cout<<res<<endl;
    }
    return 0;
}