二维哈希+二分
说是二维,其实就是先把列hash了,然后再用列的hash值hash行,这样可以O(n)的计算一个正方形的hash值,然后二分边长,枚举左上角点的坐标然后hash判断即可
只要base选的妙(闭眼随便敲一个数),单模数unsigned long long自然溢出也能过
#include<iostream>
#include<cstdio>
#include<map>
using namespace std;
const int N=55;
int n,a[2][N][N];
unsigned long long b[N],f[2][N][N],h[5005];
map<unsigned long long,bool>mp;
int clc(int t,int x,int y,int l)
{
unsigned long long r=1;
for(int i=x;i<=x+l-1;i++)
r=(r*3454+f[t][i][y+l-1]-f[t][i][y-1]*b[l]);
return r;
}
bool ok(int x)
{
mp.clear();
for(int i=1;i<=n-x+1;i++)
for(int j=1;j<=n-x+1;j++)
mp[clc(0,i,j,x)]=1;
for(int i=1;i<=n-x+1;i++)
for(int j=1;j<=n-x+1;j++)
if(mp[clc(1,i,j,x)])
return 1;
return 0;
}
int main()
{
b[0]=1;
for(int i=1;i<=50;i++)
b[i]=b[i-1]*4322;
scanf("%d",&n);
for(int k=0;k<=1;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
scanf("%d",&a[k][i][j]);
for(int k=0;k<=1;k++)
for(int i=1;i<=n;i++)
{
f[k][i][0]=1;
for(int j=1;j<=n;j++)
f[k][i][j]=f[k][i][j-1]*4322+a[k][i][j];
}
int l=1,r=n,ans=0;
while(l<=r)
{
int mid=(l+r)>>1;
if(ok(mid))
l=mid+1,ans=mid;
else
r=mid-1;
}
printf("%d",ans);
return 0;
}