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

将矩阵hash。s[0]忘了弄成0，输出中间过程发现了。

hash。sort。判重。大概这样子的步骤吧。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
#define ll unsigned long long
#define rep(i,n) for(int i=1;i<=n;i++)
#define clr(x,c) memset(x,c,sizeof(x))
#define REP(i,s,t) for(int i=s;i<=t;i++)
int read(){
	int x=0;char c=getchar();bool f=true;
	while(!isdigit(c)) {
		if(c=='-') f=false;c=getchar();
	}
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
	return f?x:-x;
} 
ll n,f[2][55][55],mm[55],s[3000];
void init(){
	n=read();
	clr(f,0);clr(s,0);
	rep(i,n) rep(j,n) f[0][i][j]=f[0][i][j-1]*103+read();
	rep(i,n) rep(j,n) f[1][i][j]=f[1][i][j-1]*103+read();
	mm[1]=1;REP(i,2,54) mm[i]=mm[i-1]*103;
}
ll get(int x,int y,int t,int op){
	ll ans=1;
	REP(i,x,x+t-1) ans=(ans*127+f[op][i][y+t-1]-f[op][i][y-1]*mm[t+1]);
	return ans;
}
bool check(int x){
	s[0]=0;
	rep(i,n-x+1) rep(j,n-x+1) s[++s[0]]=get(i,j,x,0);
	sort(s+1,s+s[0]+1);
	//rep(i,s[0]) printf("%d ",s[i]);printf("
");
	rep(i,n-x+1) rep(j,n-x+1) {
		ll tmp=get(i,j,x,1);
		//printf("=>%d
",tmp);
		if(*lower_bound(s+1,s+s[0]+1,tmp)==tmp) return true;
	}
	return false;
}
void work(){
	int l=1,r=n,ans=0,mid;
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid)) ans=mid,l=mid+1;
		else r=mid-1;
	}
	printf("%d
",ans);
}
int main(){
	init();work();return 0;
}

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

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 800 Solved: 462
[Submit][Status][Discuss]

Description

Blue Mary最近迷上了玩Starcraft(星际争霸) 的RPG游戏。她正在设法寻找更多的战役地图以进一步提高自己的水平。由于Blue Mary的技术已经达到了一定的高度,因此，对于用同一种打法能够通过的战役地图，她只需要玩一张，她就能了解这一类战役的打法，然后她就没有兴趣再玩儿这一类地图了。而网上流传的地图有很多都是属于同一种打法，因此Blue Mary需要你写一个程序，来帮助她判断哪些地图是属于同一类的。具体来说，Blue Mary已经将战役地图编码为n*n的矩阵，矩阵的每个格子里面是一个32位（有符号）正整数。对于两个矩阵，他们的相似程度定义为他们的最大公共正方形矩阵的边长。两个矩阵的相似程度越大，这两张战役地图就越有可能是属于同一类的。

Input

第一行包含一个正整数n。以下n行，每行包含n个正整数，表示第一张战役地图的代表矩阵。再以下n行，每行包含n个正整数，表示第二张战役地图的代表矩阵。

Output

仅包含一行。这一行仅有一个正整数，表示这两个矩阵的相似程度。

Sample Input

3
1 2 3
4 5 6
7 8 9
5 6 7
8 9 1
2 3 4

Sample Output

HINT

样例解释：

子矩阵：
5 6
8 9
为两个地图的最大公共矩阵

约定：
n<=50

Source

[Submit][Status][Discuss]