[BZOJ 1567] Blue Mary的战役地图

Link:

BZOJ 1567 传送门

Solution:

矩阵Hash/二维$Hash$模板题

涉及到需要快速查询、匹配的题目，考虑直接上$Hash$

矩阵$Hash$其实就是每行先各$Hash$一次，对于$x~x+l-1$行取出$y~y+l-1$列的$Hash$值再进行一次$Hash$即可

取出$l~r$的$Hash$值的方法：$Hash[j]-Hash[i-1]*base^{j-i+1}$

注意横竖两次$Hash$乘数要不同！否则碰撞概率极高

Tip: 判断一个数是否在有序序列中的方式：$*lower\_ bound(begin,end,now)==now$

Code:

#include <bits/stdc++.h>

using namespace std;
typedef unsigned long long ull;
const int MAXN=;
const ull hsn1=,hsn2=;
int n,dat[][MAXN][MAXN];
ull pre[MAXN],hs[][MAXN][MAXN],t[MAXN*MAXN];

ull cal(int k,int x,int y,int l)
{
    ull ret=;
    for(int i=x;i<=x+l-;i++)
        ret=ret*hsn2+hs[k][i][y+l-]-hs[k][i][y-]*pre[l];
    return ret;
}

bool check(int x)
{
    int cnt=;
    for(int i=;i<=n-x+;i++) for(int j=;j<=n-x+;j++)
        t[++cnt]=cal(,i,j,x);
    sort(t+,t+cnt+);
    for(int i=;i<=n-x+;i++) for(int j=;j<=n-x+;j++)
    {
        ull cur=cal(,i,j,x);
        if(*lower_bound(t+,t+cnt+,cur)==cur)
            return true;
    }
    return false;
}

int main()
{
    scanf("%d",&n);
    for (int k=;k<;k++) for(int i=;i<=n;i++) for(int j=;j<=n;j++)
        scanf("%d",&dat[k][i][j]);
    for(int k=;k<;k++)
        for(int i=;i<=n;i++)
        {
            hs[k][i][]=;
            for(int j=;j<=n;j++)
                hs[k][i][j]=hs[k][i][j-]*hsn1+dat[k][i][j];
        }

    pre[]=;for(int i=;i<=;i++) pre[i]=pre[i-]*hsn1;
    int l=,r=n;
    while(l<=r)
    {
        int mid=(l+r)>>;
        if(check(mid)) l=mid+;
        else r=mid-;
    }
    printf("%d",r);
    return ;
}