1567: [JSOI2008]Blue Mary的战役地图

被VFK的Hash Killer吓出心理阴影了(不要拿辛普森公式求阴影面积谢谢)。

然而还是硬着头皮YY了一个很脑残的Hash,大概就是求LCP时用的那个。

不过改成二维了,还纠结了一下二维前缀和的构造什么的。

然后有惊无险地A了,顺便膜Rank1的0MS,真是给跪了。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int x=123;
const int y=321;
typedef long long ll;
ll hash[2][55][55],ex[55],ey[55];
int n;
ll gethash(int x1,int y1,int len,int num){
	return hash[num][x1][y1]-hash[num][x1-len][y1]*ex[len]-hash[num][x1][y1-len]*ey[len]+hash[num][x1-len][y1-len]*ex[len]*ey[len];
}
bool check(int len){
	for(int x1=len;x1<=n;x1++)
	for(int y1=len;y1<=n;y1++)
	for(int x2=len;x2<=n;x2++)
	for(int y2=len;y2<=n;y2++)
	if(gethash(x1,y1,len,0)==gethash(x2,y2,len,1))return true;
	return false;
}
int main(){
	scanf("%d",&n);
	ex[0]=ey[0]=1;
	for(int i=1;i<=n;i++)ex[i]=ex[i-1]*x,ey[i]=ey[i-1]*y;
	ll tmp;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++){
		scanf("%lld",&tmp);
		hash[0][i][j]=hash[0][i][j-1]*y+hash[0][i-1][j]*x-hash[0][i-1][j-1]*x*y+tmp;
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++){
		scanf("%lld",&tmp);
		hash[1][i][j]=hash[1][i][j-1]*y+hash[1][i-1][j]*x-hash[1][i-1][j-1]*x*y+tmp;
	}
	if(!check(1)){printf("0");return 0;}
	int l=0,r=n;
	while(l<=r){
		int mid=l+r>>1;
		if(check(mid))l=mid+1;
		else r=mid-1;
	}
	printf("%d",l-1);
	return 0;
}


你可能感兴趣的:(1567: [JSOI2008]Blue Mary的战役地图)