Luogu P4398 [JSOI2008]Blue Mary的战役地图 矩阵哈希

其实可以二分矩阵边长但是我太懒了$qwq$。

把每个子矩阵扔到$map$里，然后就没了

#include<cstdio>
#include<map>
#include<iostream>
#define ull unsigned long long
#define R register int
using namespace std;
inline int g() {
    R ret=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
    do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
} const int B=,B2=,N=;
map<ull,bool> mp;
int n;
ull p[N],p2[N],a[N][N];
inline ull calc(int x,int y,int l) {return a[x][y]-a[x-l][y]*p2[l]-a[x][y-l]*p[l]+a[x-l][y-l]*p[l]*p2[l];}
signed main() {
    n=g(); p[]=p2[]=; for(R i=;i<=n;++i) p[i]=p[i-]*B; for(R i=;i<=n;++i) p2[i]=p2[i-]*B2;
    for(R i=;i<=n;++i) for(R j=;j<=n;++j) a[i][j]=g();
    for(R i=;i<=n;++i) for(R j=;j<=n;++j) a[i][j]=a[i][j-]*B+a[i][j];
    for(R i=;i<=n;++i) for(R j=;j<=n;++j) a[i][j]=a[i-][j]*B2+a[i][j];
    for(R l=;l<n;++l) for(R i=l;i<=n;++i) for(R j=l;j<=n;++j) mp[calc(i,j,l)]=true;
    for(R i=;i<=n;++i) for(R j=;j<=n;++j) a[i][j]=g();
    for(R i=;i<=n;++i) for(R j=;j<=n;++j) a[i][j]=a[i][j-]*B+a[i][j];
    for(R i=;i<=n;++i) for(R j=;j<=n;++j) a[i][j]=a[i-][j]*B2+a[i][j];
    for(R l=n;l;--l) for(R i=l;i<=n;++i) for(R j=l;j<=n;++j) if(mp.find(calc(i,j,l))!=mp.end()) {printf("%d\n",l); return ;}
}

---

2019.07.11