【字符串二维hs】BZOJ 2462 BeiJing 2011 矩阵模板

2462: [BeiJing2011]矩阵模板

Time Limit: 2 Sec  Memory Limit: 128 MB
Submit: 1929  Solved: 820
[Submit][Status][Discuss]
Description

给定一个M行N列的01矩阵，以及Q个A行B列的01矩阵，你需要求出这Q个矩阵哪些在
原矩阵中出现过。 
   所谓01矩阵，就是矩阵中所有元素不是0就是1。 
 

Input

输入文件的第一行为M、N、A、B，参见题目描述。 
接下来M行，每行N个字符，非0即1，描述原矩阵。 
接下来一行为你要处理的询问数Q。 
接下来Q个矩阵，一共Q*A行，每行B个字符，描述Q个01矩阵。 
 

Output

你需要输出Q行，每行为0或者1，表示这个矩阵是否出现过，0表示没有出现过，1表
示出现过。

Sample Input

3 3 2 2 

111 

000 

111 

3 

11 

00 

11 

11 

00 

11 
Sample Output

1 

0 

1 

 
HINT

对于100%的数据，N,M<=1000 A,B<=100

Source

Day4

二维hs就是横着一个seed，纵着一个seed，然后合并查询和矩阵前缀和没有区别，十分easy

注意n和m不要弄混了

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1111
#define p 13131
#define q 1313131
#define ull unsigned long long
using namespace std;
int n,m,a,b,qq;
ull hs[N][N],tt[N][N],sep[N],seq[N];
char in[N][N],ques[150][150];
inline ull geths(int l1,int l2,int r1,int r2)
    {return hs[r2][l2]-hs[r2][l1-1]*seq[b]-hs[r1-1][l2]*sep[a]+hs[r1-1][l1-1]*sep[a]*seq[b];}
int main()
{
    scanf("%d%d%d%d",&n,&m,&a,&b);
    sep[0]=seq[0]=1;
    for(int i=1;i<N;i++)sep[i]=sep[i-1]*p,seq[i]=seq[i-1]*q;
    for(int i=1;i<=m;i++)
        scanf("%s",in[i]+1);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            hs[i][j]=hs[i-1][j]*p+hs[i][j-1]*q-hs[i-1][j-1]*p*q+in[i][j];
    scanf("%d",&qq);
    for(int que=1;que<=qq;que++)
    {
        for(int i=1;i<=a;i++)
            scanf("%s",ques[i]+1);
        memset(tt,0,sizeof(tt));
        for(int i=1;i<=a;i++)
            for(int j=1;j<=b;j++)
                tt[i][j]=tt[i-1][j]*p+tt[i][j-1]*q-tt[i-1][j-1]*p*q+ques[i][j];
        ull val=tt[a][b];
        bool fl=0;
        for(int l1=1;l1+b-1<=n;l1++)
        {
            int l2=l1+b-1;
            for(int r1=1;r1+a-1<=m;r1++)
            {
                int r2=r1+a-1;
                if(geths(l1,l2,r1,r2)==val)
                    {fl=1;break;}
            }
            if(fl)break;
        }
        if(fl)printf("1
");
        else printf("0
");
    }
    return 0;
}