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;
}