HDU - 2859 Phalanx(动态规划/哈希表)

题目链接:点击查看

题目大意:给定一个整数n,然后给出一个n*n的方阵,求方阵中最大的对称子方阵(对称指的是以右上角至左下角为对角线)

分析:这个题网上有个一般做法,相当于一个动态规划的小模拟,时间复杂度为n^{3},这里不再赘述,一会直接上代码。

然后还有一个方法,是那天讲题的时候,鑫爷以非同常人的脑回路提出的,可以先用n^{2}打一个哈希表,然后再遍历每一个点进行

判断,如果判断相等则直接更新,如果不相等则可以二分查找最长的相等子序列,这样时间复杂度就下降到了n^{2}logn了,提交的

时候比一般方法快了十倍还多,不得感叹真的tql。

直接上代码吧:

一般做法:

#include<iostream>
#include<cstdio> 
#include<string>
#include<cstring>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<sstream>
#include<cmath>
using namespace std;

typedef long long LL;

const int inf=0x3f3f3f3f;

const int N=1e3+100;

char maze[N][N];

int dp[N][N];

int main()
{
//	freopen("input.txt","r",stdin);
	int n;
	while(scanf("%d",&n)!=EOF&&n)
	{
		for(int i=0;i<n;i++)
			scanf("%s",maze+i);
		fill(dp[0],dp[0]+N*N,1);
		int ans=1;
		for(int i=1;i<n;i++)
			for(int j=0;j<n;j++)
			{
				int x=i;
				int y=j;
				while(maze[x][j]==maze[i][y])
				{
					x--;
					y++;
					if(x<0||y>=n)
						break;
				}
				x=i-x;//所记录的最长相等子序列的长度
				if(x>dp[i-1][j+1])
				{
					dp[i][j]=dp[i-1][j+1]+1;
				}
				else
					dp[i][j]=x;
				ans=max(dp[i][j],ans);
			}
		cout<<ans<<endl;
	}
	
	
	
	
	
	
	
	
	
	
	
	return 0;
} 

哈希表+二分:

#include<iostream>
#include<cstdio> 
#include<string>
#include<cstring>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<sstream>
#include<cmath>
using namespace std;

typedef long long LL;

typedef unsigned long long ULL;

const int inf=0x3f3f3f3f;

const int N=1e3+100;

ULL p[N],hash_1[N][N],hash_2[N][N];

char maze[N][N];

int dp[N][N];

void hash_(int n)
{
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			hash_1[i][j]=hash_1[i][j-1]*131+(maze[i][j]-'a'+1);
		}
	for(int j=1;j<=n;j++)
		for(int i=n;i>=1;i--)
		{
			hash_2[i][j]=hash_2[i+1][j]*131+(maze[i][j]-'a'+1);
		}
}

void initp()
{
	p[0]=1;
	for(int i=1;i<N;i++)
		p[i]=p[i-1]*131;
}

void init(int n)
{
	memset(maze,0,sizeof(maze));
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			dp[i][j]=1;
	memset(hash_1,0,sizeof(hash_1));
	memset(hash_2,0,sizeof(hash_2));
}

int main()
{
//	freopen("input.txt","r",stdin);
	int n;
	initp();
	while(scanf("%d",&n)!=EOF&&n)
	{
		init(n);
		for(int i=1;i<=n;i++)
			scanf("%s",maze[i]+1);
		hash_(n);
		int mmax=1;
		for(int i=2;i<=n;i++)
			for(int j=n-1;j>=1;j--)
			{
				int k=dp[i-1][j+1];
				ULL a=hash_1[i][j+k]-hash_1[i][j-1]*p[k+1];
				ULL b=hash_2[i-k][j]-hash_2[i+1][j]*p[k+1];
				if(a==b)
					dp[i][j]=k+1;
				else
				{
					int l=1,r=k;
					int ans=0;
					while(l<=r)
					{
						int mid=(l+r)/2;
						ULL a=hash_1[i][j+mid]-hash_1[i][j-1]*p[mid+1];
						ULL b=hash_2[i-mid][j]-hash_2[i+1][j]*p[mid+1];
						if(a==b)
						{
							ans=mid;
							l=mid+1;
						}
						else
						{
							r=mid-1;
						}	
					}
					dp[i][j]=ans+1;
				}
				mmax=max(mmax,dp[i][j]);
			}
		cout<<mmax<<endl;
	}
	
	
	
	
	
	
	
	
	
	return 0;
} 

 

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Frozen_Guardian

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值