HDU 2870 Largest Submatrix

最新推荐文章于 2020-09-12 20:37:46 发布

原创最新推荐文章于 2020-09-12 20:37:46 发布 · 320 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#动态规划

动态规划专栏收录该内容

50 篇文章

订阅专栏

本文介绍了一种解决最大子矩阵问题的算法实现，通过三次遍历分别针对三种字符组合，利用动态规划思想预处理得到每个位置的最大高度、左右扩展范围，并最终计算出最大子矩阵的面积。

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

找最大子矩阵。

先将整个图转化为全部用a或b或c表示的形式，再对此图进行操作。

定义数组h[i][j]表示以点i,j处的最大高度，r[i][j]表示从点i,j开始向左扩展最多扩展多少位，l[i][j]表示向左扩展多少位。

最后结果为max(h[i][j]*(l[i][j]+r[i][j]-1))。

代码：

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int a[1001][1001];
int map[1001][1001];
int h[1001][1001];
int l[1001][1001];
int r[1001][1001];
int main()
{
	ios::sync_with_stdio(false);
	int n,m;
	while(cin>>n>>m)
	{
		memset(a,0,sizeof(a));
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=m;j++)
			{
				char c;cin>>c;
				a[i][j]=c;
			}
		}
		int sum=0;
		for(int k=1;k<=3;k++)
		{
			memset(map,0,sizeof(map));
			memset(h,0,sizeof(h));
			memset(l,0,sizeof(l));
			memset(r,0,sizeof(r));
			for(int i=1;i<=n;i++)
			{
				for(int j=1;j<=m;j++)
				{
					if(k==1&&(a[i][j]=='a'||a[i][j]=='w'||a[i][j]=='y'||a[i][j]=='z'))map[i][j]=1;
					else if(k==2&&(a[i][j]=='b'||a[i][j]=='w'||a[i][j]=='y'||a[i][j]=='z'))map[i][j]=1;
					else if(k==3&&(a[i][j]=='c'||a[i][j]=='x'||a[i][j]=='y'||a[i][j]=='z'))map[i][j]=1;
					else map[i][j]=0;
				}
			}
			for(int j=1;j<=m;j++)
			{
				for(int i=1;i<=n;i++)
				{
					if(map[i][j]==1)h[i][j]=h[i-1][j]+1;
				}
			}
			for(int i=1;i<=n;i++)
			{
				for(int j=1;j<=m;j++)
				{
					if(map[i][j]==1)
					{
						l[i][j]++;
						int t=j-1;
						while(h[i][t]>=h[i][j])
						{
							l[i][j]+=l[i][t];
							t-=l[i][t];
						}
					}
				}
			}
			for(int i=1;i<=n;i++)
			{
				for(int j=m;j>=1;j--)
				{
					if(map[i][j]==1)
					{
						r[i][j]++;
						int t=j+1;
						while(h[i][t]>=h[i][j])
						{
							r[i][j]+=r[i][t];
							t+=r[i][t];
						}
					}
				}
			}
			int tot=0;
			for(int i=1;i<=n;i++)
			{
				for(int j=1;j<=m;j++)
				{
					tot=max(tot,(l[i][j]+r[i][j]-1)*h[i][j]);
				}
			}
			sum=max(sum,tot);
		}
		cout<<sum<<endl;
	}
	return 0;
}