二维矩阵匹配问题,至今不知道Q的范围是多少,反正是要求做到读入复杂度。
二维Hash:就是一维的等效拓展,注意两维的Base不能相同。
其余就是一维Hash和二维前缀和的结合,可以自然溢出,据说概率很科学。
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 #define rep(i,l,r) for (int i=(l); i<=(r); i++) 5 typedef unsigned int ull; 6 typedef long long ll; 7 using namespace std; 8 9 const int N=1010,P1=2,P2=9191891; 10 int n,m,r,c,Q; 11 ull a[N][N],pw1[N],pw2[N],tot,hs[1000500]; 12 char mp[N][N]; 13 14 int main(){ 15 freopen("bzoj2462.in","r",stdin); 16 freopen("bzoj2462.out","w",stdout); 17 scanf("%d%d%d%d",&n,&m,&r,&c); pw1[0]=pw2[0]=1; 18 rep(i,1,1000) pw1[i]=pw1[i-1]*P1,pw2[i]=pw2[i-1]*P2; 19 rep(i,1,n) scanf("%s",mp[i]+1); 20 rep(i,1,n) rep(j,1,m) a[i][j]=a[i][j-1]*P1+mp[i][j]-'0'; 21 rep(i,1,n) rep(j,1,m) a[i][j]=a[i-1][j]*P2+a[i][j]; 22 rep(i,r,n) rep(j,c,m) 23 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]; 24 sort(hs+1,hs+tot+1); memset(a,0,sizeof(a)); 25 for (scanf("%d",&Q); Q--; ){ 26 rep(i,1,r) scanf("%s",mp[i]+1); 27 rep(i,1,r) rep(j,1,c) a[i][j]=a[i][j-1]*P1+mp[i][j]-'0'; 28 rep(i,1,r) rep(j,1,c) a[i][j]=a[i-1][j]*P2+a[i][j]; 29 int k=lower_bound(hs+1,hs+tot+1,a[r][c])-hs; 30 if (hs[k]==a[r][c]) puts("1"); else puts("0"); 31 } 32 return 0; 33 }