poj 3690 Constellations 二维滚动哈希匹配

//poj 3690
//sep9
#include <iostream>
#include <set>
using namespace std;
typedef unsigned long long ull;

const int MAXN=1024;
const int MAXT=128;
const int MAXP=64;
const int B1=9973;
const int B2=100000007;
char sky[MAXN][MAXN];
char patterns[MAXT][MAXN][MAXN];
int n,m,t,p,q;
ull hash[MAXN][MAXN];

void compute_hash(char a[MAXN][MAXN],int row,int col)
{
	ull x=1;
	for(int i=0;i<q;++i)	
		x=x*B1;
	for(int i=0;i<row;++i){
		ull e=0;
		for(int j=0;j<q;++j)
			e=B1*e+a[i][j];
		for(int j=0;j+q<=col;++j){
			hash[i][j]=e;
			if(j+q<m)
		 		e=B1*e+a[i][j+q]-x*a[i][j];
		}	
	}
	ull y=1;
	for(int i=0;i<p;++i)
		y=y*B2;
	for(int i=0;i+q<=col;++i){
		ull e=0,tmp;
		for(int j=0;j<p;++j)
			e=B2*e+hash[j][i];
		for(int j=0;j+p<=row;++j){
			tmp=hash[j][i];
			hash[j][i]=e;
			if(j+p<n)
				e=B2*e+hash[j+p][i]-y*tmp;
		}
	}
}

int main()
{
	int cases=0;
	while(scanf("%d%d%d%d%d",&n,&m,&t,&p,&q)==5){
		if(n==0&&m==0&&t==0&&p==0&&q==0) break;
		for(int i=0;i<n;++i)
			scanf("%s",sky[i]);
		for(int i=0;i<t;++i)
			for(int j=0;j<p;++j)
				scanf("%s",patterns[i][j]);
		multiset<ull> myset;
		for(int i=0;i<t;++i){
			compute_hash(patterns[i],p,q);
			myset.insert(hash[0][0]);	
		}
		compute_hash(sky,n,m);
		for(int i=0;i+p<=n;++i)
			for(int j=0;j+q<=m;++j)
				myset.erase(hash[i][j]);	
		printf("Case %d: %d\n",++cases,t-myset.size());
	}
}