字符创匹配（Rabin-Karp 算法）（hash）

<span style="font-size:18px;">//题目链接：http://poj.org/problem?id=3690

//复杂度:如果一一检查，会退化为O(mn) 不同字符串哈希值冲突的概率非常低，
//通常可以忽视，实际不做朴素的检查，所以是O(m+n)。

//哈希函数：H(C)=(c1*pow(b,m-1)+c2*pow(b,m-2)+c3*pow(b,m-3)+...+cm*pow(b,0))mod h.（b为哈希的基数）
//递推：H(S[k+1...k+m])=( H(S[k...k+m-1])*b-sk*pow(b,m)+sk+m )mod h
//取h=pow(2,64),用自然溢出省去模运算

#include <iostream>
#include <cstdio>
#include <cstring>
#include <set>
using namespace std;
typedef unsigned long long ull;
const int maxn=1005,maxt=105;

const ull b1=9973,b2=100000007;
ull tmp[maxn][maxn],Hash[maxn][maxn];
char field[maxn][maxn],model[maxt][maxn][maxn];
multiset<ull> mul;
int P,Q;

void solve(char M[maxn][maxn],int r,int c)
{
   int i,j;
   ull t=1,e;
   for(i=0;i<Q;i++)t*=b1;
   for(i=0;i<r;i++)
   {
      e=0;
      for(j=0;j<Q;j++)
      {
         e=e*b1+M[i][j];
      }
      for(j=0;j+Q<=c;j++)
      {
         tmp[i][j]=e;
         if(j+Q<c)
         {
            e=e*b1-M[i][j]*t+M[i][j+Q];
         }
      }
   }
   t=1;
   for(i=0;i<P;i++)t*=b2;
   for(i=0;i<c-Q+1;i++)
   {
      e=0;
      for(j=0;j<P;j++)
      {
         e=e*b2+tmp[j][i];
      }
      for(j=0;j+P<=r;j++)
      {
         Hash[j][i]=e;
         if(j+P<r)
         {
            e=e*b2-tmp[j][i]*t+tmp[j+P][i];
         }
      }
   }
}

int main()
{
   int N,M,T,i,j,kase=1;
   while(scanf("%d%d%d%d%d",&N,&M,&T,&P,&Q)!=EOF)
   {
      mul.clear();
      if(N==0 && M==0 && T==0 && P==0 && Q==0)break;
      getchar();
      for(i=0;i<N;i++) scanf("%s",field[i]);
      for(i=0;i<T;i++)
      {
         getchar();
         for(j=0;j<P;j++)
         {
            scanf("%s",model[i][j]);
         }
      }
      for(i=0;i<T;i++)
      {
         solve(model[i],P,Q);
         mul.insert(Hash[0][0]);
      }
      solve(field,N,M);
      for(i=0;i<N-P+1;i++)
      {
         for(j=0;j<M-Q+1;j++)
         {
            mul.erase(Hash[i][j]);
         }
      }
      printf("Case %d: %d\n",kase++,T-mul.size());
   }
   return 0;
}

</span>