字符串hash(前缀，区间)。kmp求min循环节长度。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int N =2e3+10,P=1e9+7,seed=31;
//除预处理hash外(O(n^2))，hash调用O(1)
//分别对行列字符串hash压缩成一维。再对一维的数值跑kmp
ll n,q,hr[N][N],hc[N][N],p[N],a[N],nex[N];
string s[N];
int kmp(int n){//1~n,n<2e3+7（自己和自己匹配求nex数值）
    memset(nex,0,sizeof(nex));
    int j=0;
    for(int i=1;i<=n;i++){
        while(j&&a[i+1]!=a[j+1])j=nex[j];//下一个不匹配就跳
        if(a[i+1]==a[j+1])j++;
        nex[i+1]=j;
    }
    return n-nex[n];//返回该序列的最小循环节长度
}
int main() {
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    p[0]=1;//p[i]:seed的i次方
    for(int i=1;i<=2007;i++){p[i]=p[i-1]*seed%P;}
    cin>>n>>q;
    for(int i=1;i<=n;i++){cin>>s[i];s[i]=" "+s[i];}
    for(int i=1;i<=n;i++){//前缀hash(一维)
        for(int j=1;j<=n;j++){
            hr[i][j]=(hr[i][j-1]*seed%P+(s[i][j]-'a'+1))%P;
            hc[i][j]=(hc[i-1][j]*seed%P+(s[i][j]-'a'+1))%P;
            //hash一般值很大，都要取余使他的值落在P区间内，但P太小或冲突
        }
    }
    int x1,x2,y1,y2;
    while(q--){
        cin>>x1>>y1>>x2>>y2;
        int m=0;
        for(int i=1;i<=n;i++)a[i]=0;//存每行（y1~y2区间子串的hash值）
        for(int i=x1;i<=x2;i++){//区间hash(一维)
            a[++m]=(hr[i][y2]-hr[i][y1-1]*p[y2-y1+1]%P+P)%P;
        }
        int tx=kmp(m);//求矩阵(各行子串)的最短循环节长度
        m=0;for(int i=1;i<=n;i++)a[i]=0;//对列
        for(int i=y1;i<=y2;i++){
            a[++m]=(hc[x2][i]-hc[x1-1][i]*p[x2-x1+1]%P+P)%P;
        }
        int ty=kmp(m);
        cout<<tx*ty<<endl;
    }
    return 0;
}