1567: [JSOI2008]Blue Mary的战役地图
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
1 2 3
4 5 6
7 8 9
5 6 7
8 9 1
2 3 4
Sample Output
2
题解:
二维hash
#include<bits/stdc++.h> using namespace std; #pragma comment(linker, "/STACK:102400000,102400000") #define ls i<<1 #define rs ls | 1 #define mid ((ll+rr)>>1) #define pii pair<int,int> #define MP make_pair typedef long long LL; typedef unsigned long long ULL; const long long INF = 1e18+1LL; const double pi = acos(-1.0); const int N = 2e5+10, M = 1e3+20,inf = 2e9; /*inline LL read() { LL x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; }*/ const ULL mod1 = 133700613937LL, mod2 = 10031LL; ULL sqr1[60],sqr2[60],ahas[60][60],bhas[60][60]; LL a[60][60],b[60][60]; int n; map<ULL, int > mp; int ok(int w) { mp.clear(); for(int i = w; i <= n; ++i) { for(int j = w; j <= n; ++j) { mp[sqr1[n - i + 1] * sqr2[n - j + 1] * (ahas[i][j] + ahas[i-w][j-w] - ahas[i-w][j] - ahas[i][j-w])] = 1; } } for(int i = w; i <= n; ++i) { for(int j = w; j <= n; ++j) { if(mp[sqr1[n - i + 1] * sqr2[n - j + 1] * (bhas[i][j] + bhas[i-w][j-w] - bhas[i-w][j] - bhas[i][j-w])]) return 1; } } return 0; } int main() { sqr1[0] = 1LL;sqr2[0] = 1LL; for(int i = 1; i <= 59; ++i) sqr1[i] = sqr1[i-1] * mod1, sqr2[i] = sqr2[i-1] * mod2; scanf("%d",&n); for(int i = 1; i <= n; ++i) for(int j = 1; j <= n; ++j) scanf("%lld",&a[i][j]), ahas[i][j] = sqr1[i] * sqr2[j] * a[i][j], ahas[i][j] += ahas[i][j-1] + ahas[i-1][j] - ahas[i-1][j-1]; for(int i = 1; i <= n; ++i) for(int j = 1; j <= n; ++j) scanf("%lld",&b[i][j]), bhas[i][j] = sqr1[i] * sqr2[j] * b[i][j], bhas[i][j] += bhas[i][j-1] + bhas[i-1][j] - bhas[i-1][j-1]; int ans = 0; for(int w = n; w >= 1; w--) { if(ok(w)) { printf("%d\n",w); return 0; } } printf("0\n"); return 0; }