hdu 2870 动态规划

最新推荐文章于 2019-12-18 20:09:00 发布

onthewaytotop

最新推荐文章于 2019-12-18 20:09:00 发布

阅读量246

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/onthewaytotop/article/details/50811534

本文介绍了一种算法，该算法通过分析字符矩阵中特定字符的分布来寻找最大的子矩阵，具体关注了'a'、'b'、'c'等字符，并采用动态规划的方法进行高效计算。

//1505 更加强版，分为a，b，c三种情况~
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int a[1005][1005],l[1005],r[1005];
char map[1005][1005];
int main()
{
	int m,n;
	while(~scanf("%d%d",&m,&n))
	{
		int i,j;
		int ans=0,res;
		for(i=1;i<=m;i++)
		{
			scanf("%s",map[i]+1);
		}
		for(j=1;j<=n;j++) a[0][j]=0;
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++)
			if(map[i][j]=='a'||map[i][j]=='w'||map[i][j]=='y'||map[i][j]=='z')
			a[i][j]=a[i-1][j]+1;
			else
			a[i][j]=0;
		}
		int t;
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++) l[j]=r[j]=1;
			for(j=2;j<=n;j++)
			{
				t=j-1;
				while(a[i][t]>=a[i][j]&&t>=1)
				{
					l[j]+=l[t];
					t=t-l[t];
				}
			}
			for(j=n-1;j>=1;j--)
			{
				t=j+1;
				while(a[i][t]>=a[i][j]&&t<=n)
				{
					r[j]+=r[t];
					t=t+r[t];
				}
			}
			for(j=1;j<=n;j++)
			{
				res=a[i][j]*(l[j]+r[j]-1);
				if(ans<res) ans=res;
			}
		}
		for(j=1;j<=m;j++) a[0][j]=0;
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++)
			if(map[i][j]=='b'||map[i][j]=='w'||map[i][j]=='x'||map[i][j]=='z')
			a[i][j]=a[i-1][j]+1;
			else
			a[i][j]=0;
		}
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++) l[j]=r[j]=1;
			for(j=2;j<=n;j++)
			{
				t=j-1;
				while(a[i][t]>=a[i][j]&&t>=1)
				{
					l[j]+=l[t];
					t=t-l[t];
				}
			}
			for(j=n-1;j>=1;j--)
			{
				t=j+1;
				while(a[i][t]>=a[i][j]&&t<=n)
				{
					r[j]+=r[t];
					t=t+r[t];
				}
			}
			for(j=1;j<=n;j++)
			{
				res=a[i][j]*(l[j]+r[j]-1);
				if(ans<res) ans=res;
			}
		}
		for(j=1;j<=m;j++) a[0][j]=0;
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++)
			if(map[i][j]=='c'||map[i][j]=='x'||map[i][j]=='y'||map[i][j]=='z')
			a[i][j]=a[i-1][j]+1;
			else
			a[i][j]=0;
		}
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++) l[j]=r[j]=1;
			for(j=2;j<=n;j++)
			{
				t=j-1;
				while(a[i][t]>=a[i][j]&&t>=1)
				{
					l[j]+=l[t];
					t=t-l[t];
				}
			}
			for(j=n-1;j>=1;j--)
			{
				t=j+1;
				while(a[i][t]>=a[i][j]&&t<=n)
				{
					r[j]+=r[t];
					t=t+r[t];
				}
			}
			for(j=1;j<=n;j++)
			{
				res=a[i][j]*(l[j]+r[j]-1);
				if(ans<res) ans=res;
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}