题目大意
给你一个矩阵,在里面找两个最大的相同的子正方形矩阵(可重叠,输出正方形边长和两个正方形的左上角坐标)
题目分析
首先二分正方形矩阵的边长,然后哈希每一个这样的正方形矩阵即可。
哈希的方法如下(x是矩阵中的元素)
然后容斥解决每一个子矩阵。(具体看代码,但是…看的时候请冷静点,因为很长)
⎡⎣⎢x1∗p2∗q3x5∗p∗q2x9∗q3x2∗p2∗q2x6∗p∗q2x10∗q2x3∗p2∗qx7∗p∗qx11∗qx4∗p2x8∗px12⎤⎦⎥
注意以下几点:
1.哈希的种子要设大一点(WA*2+调不出来的惨痛教训)
2.在更新了ans的情况下才要更新x1,y1,x2,y2(WA*1的惨痛教训)
代码
#include<algorithm>
#include<iostream>
#include<cstring>
#include<climits>
#include<cstdio>
using namespace std;
#define ull unsigned long long
const int N=505;
int n,m,x1,y1,x2,y2,mod=999997,tot,ans=0;
char s[N][N];
ull ha[N][N],p=1313,q=3131;
int h[1000000],ne[N*N],xx[N*N],yy[N*N];ull num[N*N];
ull sp[N],sq[N];
int add(int x,int y,int len){
ull t=ha[x+len-1][y+len-1]-ha[x-1][y+len-1]*sp[len]-ha[x+len-1][y-1]*sq[len]+ha[x-1][y-1]*sp[len]*sq[len];//通过容斥获得哈希值
int kl=t%mod;
if(kl<0)kl=-kl;
for(int i=h[kl];i!=-1;i=ne[i])//在链表hash里查询
if(num[i]==t){
if(len>ans)x1=xx[i],y1=yy[i];//更新的前提是len>ans
return 1;}
++tot,num[tot]=t,xx[tot]=x,yy[tot]=y;
ne[tot]=h[kl],h[kl]=tot;
return 0;
}
int ok(int len){
int i,j;
memset(h,-1,sizeof(h));tot=0;
for(i=1;i<=n-len+1;++i)//枚举左上角坐标
for(j=1;j<=m-len+1;++j)
if(add(i,j,len)){
if(len>ans)x2=i,y2=j;//更新的前提是len>ans
return 1;
}
return 0;
}
int main()
{
int i,j,l=1,r,mid;
scanf("%d%d",&n,&m);r=min(n,m);
for(i=0;i<n;++i)scanf("%s",s[i]);
sp[1]=p,sq[1]=q;
for(i=2;i<=500;++i)sp[i]=sp[i-1]*p,sq[i]=sq[i-1]*q;
for(i=1;i<=n;++i)
for(j=1;j<=m;++j)
ha[i][j]=ha[i-1][j]*p+ha[i][j-1]*q-ha[i-1][j-1]*p*q+s[i-1][j-1];//通过容斥获得哈希值
while(l<=r){//二分
mid=(l+r)>>1;
if(ok(mid))ans=mid,l=mid+1;
else r=mid-1;
}
printf("%d\n",ans);
if(!ans)return 0;
printf("%d %d\n",x1,y1);
printf("%d %d\n",x2,y2);
return 0;
}