bzoj 1567: [JSOI2008]Blue Mary的战役地图【二分+hash】

二维哈希+二分
说是二维,其实就是先把列hash了,然后再用列的hash值hash行,这样可以O(n)的计算一个正方形的hash值,然后二分边长,枚举左上角点的坐标然后hash判断即可
只要base选的妙(闭眼随便敲一个数),单模数unsigned long long自然溢出也能过

#include<iostream>
#include<cstdio>
#include<map>
using namespace std;
const int N=55;
int n,a[2][N][N];
unsigned long long b[N],f[2][N][N],h[5005];
map<unsigned long long,bool>mp;
int clc(int t,int x,int y,int l)
{
	unsigned long long r=1;
	for(int i=x;i<=x+l-1;i++)
		r=(r*3454+f[t][i][y+l-1]-f[t][i][y-1]*b[l]);
	return r;
}
bool ok(int x)
{
	mp.clear();
	for(int i=1;i<=n-x+1;i++)
		for(int j=1;j<=n-x+1;j++)
			mp[clc(0,i,j,x)]=1;
	for(int i=1;i<=n-x+1;i++)
		for(int j=1;j<=n-x+1;j++)
			if(mp[clc(1,i,j,x)])
				return 1;
	return 0;
}
int main()
{
	b[0]=1;
	for(int i=1;i<=50;i++)
		b[i]=b[i-1]*4322;
	scanf("%d",&n);
	for(int k=0;k<=1;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				scanf("%d",&a[k][i][j]);
    for(int k=0;k<=1;k++)
		for(int i=1;i<=n;i++)
		{
			f[k][i][0]=1;
			for(int j=1;j<=n;j++)
				f[k][i][j]=f[k][i][j-1]*4322+a[k][i][j];
		}
	int l=1,r=n,ans=0;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(ok(mid))
			l=mid+1,ans=mid;
		else 
			r=mid-1;
	}
	printf("%d",ans);
	return 0;
}
posted @ 2018-09-07 17:14  lokiii  阅读(111)  评论(0编辑  收藏  举报