bzoj1567 [JSOI2008]Blue Mary的战役地图

http://www.elijahqi.win/archives/2985
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
二分一下最大的长度是多少 然后用二维hash判断是否相同 每次用二分查找判断hash值是否相同

#include<cstdio>
#include<algorithm>
#define ll unsigned long long
#define N 55
#define g1 1333
#define g2 23333
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
ll hs1[N][N],hs2[N][N],q[N*N],p1[N],p2[N];int n; 
inline bool check(int md){int cnt=0;
    for (int i=md;i<=n;++i){
        for (int j=md;j<=n;++j){
            ll tmp=hs1[i][j];tmp-=hs1[i-md][j]*p2[md];tmp-=hs1[i][j-md]*p1[md];
            tmp+=hs1[i-md][j-md]*p2[md]*p1[md];q[++cnt]=tmp;
        }
    }sort(q+1,q+cnt+1);
    for (int i=md;i<=n;++i){
        for (int j=md;j<=n;++j){
            ll tmp=hs2[i][j];tmp-=hs2[i-md][j]*p2[md];tmp-=hs2[i][j-md]*p1[md];
            tmp+=hs2[i-md][j-md]*p2[md]*p1[md];if (*lower_bound(q+1,q+cnt+1,tmp)==tmp) return 1;
        }
    }return 0;
}
int main(){
    freopen("bzoj1567.in","r",stdin);
    n=read();p1[0]=p2[0]=1;for (int i=1;i<=n;++i) p1[i]=p1[i-1]*g1,p2[i]=p2[i-1]*g2;
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) hs1[i][j]=hs1[i][j-1]*g1+read();
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) hs2[i][j]=hs2[i][j-1]*g1+read();
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) hs1[i][j]+=hs1[i-1][j]*g2;
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) hs2[i][j]+=hs2[i-1][j]*g2;
    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\n",r);
    return 0;
}