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=,P1=,P2=;
 int n,m,r,c,Q;
 ull a[N][N],pw1[N],pw2[N],tot,hs[];
 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[]=pw2[]=;
     rep(i,,) pw1[i]=pw1[i-]*P1,pw2[i]=pw2[i-]*P2;
     rep(i,,n) scanf("%s",mp[i]+);
     rep(i,,n) rep(j,,m) a[i][j]=a[i][j-]*P1+mp[i][j]-'';
     rep(i,,n) rep(j,,m) a[i][j]=a[i-][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+,hs+tot+); memset(a,,sizeof(a));
     for (scanf("%d",&Q); Q--; ){
         rep(i,,r) scanf("%s",mp[i]+);
         rep(i,,r) rep(j,,c) a[i][j]=a[i][j-]*P1+mp[i][j]-'';
         rep(i,,r) rep(j,,c) a[i][j]=a[i-][j]*P2+a[i][j];
         int k=lower_bound(hs+,hs+tot+,a[r][c])-hs;
         if (hs[k]==a[r][c]) puts(""); else puts("");
     }
     return ;
 }