参考：http://www.bubuko.com/infodetail-3544495.html
ac代码+解析：

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
typedef unsigned long long ull;
#define base1 3
#define base2 5//这两个数尽量为奇数
int n,m,p,q;
ull hash_code1;//小矩阵的哈希值
char ch1[2001][2001];
char ch2[2001][2001];
ull hash1[2001][2001];
ull hash2[2001][2001];

int main()
{
    cin>>n>>m>>p>>q;
    for(int i=0;i<n;i++)
    {
        ull sum=0;
        for(int j=0;j<m;j++)
        {
            cin>>ch1[i][j];
            sum=sum*base1+ch1[i][j];
        }
        hash_code1=hash_code1*base2+sum;
    }
    for(int i=0;i<p;i++)
    {
        for(int j=0;j<q;j++)
        {
            cin>>ch2[i][j];
        }
    }
    ull now=1;
    for(int i=0;i<m;i++)
    {
        now*=base1;
    }
    for(int i=0;i<p;i++)
    {
        ull sum=0;
        for(int j=0;j<m;j++)
        {
            sum=sum*base1+ch2[i][j];
        }
        hash1[i][m-1]=sum;//每行长度为m的字符串的哈希值,长度小于m没必要求
        for(int j=m;j<q;j++)
        {
            hash1[i][j]=hash1[i][j-1]*base1+ch2[i][j]-ch2[i][j-m]*now;//这里求的是以ch2[i][j]为结尾的长度为m的字符串的哈希值
        }
    }
    now=1;
    int ans=0;
    for(int i=0;i<n;i++)
    {
        now*=base2;
    }
    for(int i=m-1;i<q;i++)
    {
        ull sum=0;
        for(int j=0;j<n;j++)
        {
            sum=sum*base2+hash1[j][i];//二维哈希
        }
        hash2[n-1][i]=sum;//右下角为ch2[n-1][i]行数为n列数为m的小矩阵的哈希值
        if(hash2[n-1][i]==hash_code1)
        {
            ans++;
        }
        for(int j=n;j<p;j++)
        {
            hash2[j][i]=hash2[j-1][i]*base2+hash1[j][i]-hash1[j-n][i]*now;
            if(hash2[j][i]==hash_code1)
            {
                ans++;
            }
        }
    }
    printf("%d\n",ans);
    return 0;
}