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

前言：

好像几天都懒得写要补很多。

题意：

求最大相同正方形。

题解：

二分+hash。
一开始想到要二维hash$O(n^4)$判断。
什么，不会啊。
后来算了算假如用一维的hash满打满算好像也就$O(25^4)$乘点常数。
好像卡一卡能过。
于是竟然真的过了。
orzTYB知道二维怎么写。
code：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#define LL long long
using namespace std;
int n;
int map[55][55][5];
LL ha[55][55][5];
LL bas=10007;//行hash 
const LL mod=1000000007;
LL quick(LL a,int b)
{
    LL ans=1;
    while(b)
    {
        if(b&1) (ans*=a)%=mod;
        (a*=a)%=mod;b>>=1;
    }
    return ans;
}
LL get_hash(int x,int l,int r,int k)
{
    return ((ha[x][r][k]-ha[x][l-1][k]*quick(bas,(r-l+1))%mod)%mod+mod)%mod;
}
bool solve(int x,int y,int x1,int y1,int len)
{
    for(int i=1;i<=len;i++)
        if(get_hash(x+i-1,y-len+1,y,1)!=get_hash(x1+i-1,y1-len+1,y1,2)) return false;
    return true;
}
bool check(int x)
{
    if(x==0) return true;
    for(int i=1;i<=n-x+1;i++)
    for(int j=x;j<=n;j++)
    for(int k=1;k<=n-x+1;k++)
    for(int u=x;u<=n;u++)
        if(solve(i,j,k,u,x)) return true;
    return false;
}
int main()
{
    scanf("%d",&n);
    for(int k=1;k<=2;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++) scanf("%d",&map[i][j][k]);
    for(int k=1;k<=2;k++)
        for(int i=1;i<=n;i++)
        {
            ha[i][0][k]=0;
            for(int j=1;j<=n;j++) ha[i][j][k]=(map[i][j][k]+ha[i][j-1][k]*bas)%mod;
        }
    int l=0,r=n,ans=0;
    while(l<=r) 
    {
        int mid=(l+r)/2;
        if(check(mid)) ans=mid,l=mid+1;
        else r=mid-1;
    }
    //printf("%d\n",solve(2,3,1,2,2));
    printf("%d",ans);
}