ural1486 二分+二维哈希

题目大意

给你一个矩阵，在里面找两个最大的相同的子正方形矩阵（可重叠，输出正方形边长和两个正方形的左上角坐标）

题目分析

首先二分正方形矩阵的边长，然后哈希每一个这样的正方形矩阵即可。
哈希的方法如下（x是矩阵中的元素）
然后容斥解决每一个子矩阵。（具体看代码，但是…看的时候请冷静点，因为很长）

$$\begin{bmatrix} x_1*p^2*q^3 & x_2*p^2*q^2 & x_3*p^2*q & x_4*p^2 \\ x_5*p*q^2 & x_6*p*q^2 & x_7*p*q & x_8*p \\ x_9*q^3 & x_{10}*q^2 & x_{11}*q & x_{12} \end{bmatrix}$$
注意以下几点：
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;
}