BZOJ 1567: [JSOI2008]Blue Mary的战役地图( 二分答案 + hash )

二分答案, 然后用哈希去判断... 

-------------------------------------------------------------------------

#include<cstdio>

#include<algorithm>

#include<cstring>

#include<set>

 

using namespace std;

 

typedef unsigned long long ull;

 

const int maxn = 59;

const int p = 31, _p = 59;

const int P = 10000007;

 

int N, seq[maxn];

ull ha[maxn][maxn], hb[maxn][maxn], Power[maxn];

set<ull> Bst;

 

void init() {

scanf("%d", &N);

for(int i = 0; i < N; i++) {

ha[i][N] = 0;

for(int j = 0; j < N; j++) scanf("%d", seq + j);

for(int j = N; j--; ) ha[i][j] = ha[i][j + 1] * p + seq[j];

}

for(int i = 0; i < N; i++) {

hb[i][N] = 0;

for(int j = 0; j < N; j++) scanf("%d", seq + j);

for(int j = N; j--; ) hb[i][j] = hb[i][j + 1] * p + seq[j];

}

Power[0] = 1;

for(int i = 1; i <= N; i++)

Power[i] = Power[i - 1] * p;

}

 

ull get(int R, int l, int r, ull h[][maxn]) {

return h[R][l] - h[R][r] * Power[r - l];

}

 

ull Get(int Rl, int Rr, int Cl, int Cr, ull h[][maxn]) {

ull ret = 0;

for(int i = Rl; i < Rr; i++)

ret = ret * _p + get(i, Cl, Cr, h);

return ret;

}

 

bool check(int m) {

Bst.clear();

for(int i = N - m; ~i; i--)

for(int j = N - m; ~j; j--)

Bst.insert(Get(i, i + m, j, j + m, ha));

for(int i = N - m; ~i; i--)

for(int j = N - m; ~j; j--)

if(Bst.find(Get(i, i + m, j, j + m, hb)) != Bst.end()) return true;

return false;

}

 

int main() {

init();

int ans = 0, l = 1, r = N;

while(l <= r) {

int m = (l + r) >> 1;

if(check(m))

ans = m, l = m + 1;

else

r = m - 1;

}

printf("%d ", ans);

return 0;

}

------------------------------------------------------------------------- 

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

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 690  Solved: 395
[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

2

HINT

样例解释：

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

约定：
n<=50

Source