BZOJ2462[Beijing2011]矩阵模板(二维Hash)

二维矩阵匹配问题，至今不知道Q的范围是多少，反正是要求做到读入复杂度。

二维Hash：就是一维的等效拓展，注意两维的Base不能相同。

其余就是一维Hash和二维前缀和的结合，可以自然溢出，据说概率很科学。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
typedef unsigned int ull;
typedef long long ll;
using namespace std;

const int N=1010,P1=2,P2=9191891;
int n,m,r,c,Q;
ull a[N][N],pw1[N],pw2[N],tot,hs[1000500];
char mp[N][N];

int main(){
    freopen("bzoj2462.in","r",stdin);
    freopen("bzoj2462.out","w",stdout);
    scanf("%d%d%d%d",&n,&m,&r,&c); pw1[0]=pw2[0]=1;
    rep(i,1,1000) pw1[i]=pw1[i-1]*P1,pw2[i]=pw2[i-1]*P2;
    rep(i,1,n) scanf("%s",mp[i]+1);
    rep(i,1,n) rep(j,1,m) a[i][j]=a[i][j-1]*P1+mp[i][j]-'0';
    rep(i,1,n) rep(j,1,m) a[i][j]=a[i-1][j]*P2+a[i][j];
    rep(i,r,n) rep(j,c,m)
        hs[++tot]=a[i][j]-a[i][j-c]*pw1[c]-a[i-r][j]*pw2[r]+a[i-r][j-c]*pw1[c]*pw2[r];
    sort(hs+1,hs+tot+1); memset(a,0,sizeof(a));
    for (scanf("%d",&Q); Q--; ){
        rep(i,1,r) scanf("%s",mp[i]+1);
        rep(i,1,r) rep(j,1,c) a[i][j]=a[i][j-1]*P1+mp[i][j]-'0';
        rep(i,1,r) rep(j,1,c) a[i][j]=a[i-1][j]*P2+a[i][j];
        int k=lower_bound(hs+1,hs+tot+1,a[r][c])-hs;
        if (hs[k]==a[r][c]) puts("1"); else puts("0");
    }
    return 0;
}