POJ Constellations

思路（主要参考：传送门）

代码

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<set>
#define ull unsigned long long
using namespace std;
int N,M,T,P,Q; 
char f[1000][1000];//匹配对象
char p[1000][1000];//匹配模式 
ull Hash[1000][1000],t[1000][1000];
void compute_hash(char a[1000][1000],int n,int m)//a表示要求哈希值的矩阵，并且是n*m的
{
	const ull B1=1e9+7;
	const ull B2=1e9+9;
	
	ull t1=1;//B1的Q次方
	for(int i=0;i<Q;i++) t1*=B1;
	
	//按行方向计算哈希值
	for(int i=0;i<n;i++)
	{
		ull e=0;
		for(int j=0;j<Q;j++) e=e*B1+a[i][j];//求出起始位置为0，长度为Q的字符串的哈希值
		
		for(int j=0;j+Q<=m;j++)//枚举一下起始位置，求出第i行所有起始位置为j，长度为Q的字符串的哈希值
		{
			t[i][j]=e;
			if(j+Q<m) e=e*B1+a[i][j+Q]-a[i][j]*t1;//每次更新一下e，表示当前的哈希值，这里是列数低的是B1进制的高位	
		}	
	} 
	
	ull t2=1;//B2的P次方
	for(int i=0;i<P;i++) t2*=B2;
	
	//按列方向计算哈希值
	for(int j=0;j+Q<=m;j++)
	{
		ull e=0;
		for(int i=0;i<P;i++) e=e*B2+t[i][j];//这里就把第一次的哈希值作为了新矩阵，来再求一次哈希值
		
		for(int i=0;i+P<=n;i++)//同理去枚举一下每一个起始位置
		{
			Hash[i][j]=e;
			if(i+P<n) e=e*B2+t[i+P][j]-t[i][j]*t2;
		}	
	} 
	
}
int main()
{
	multiset<ull> s;
	int count=0;
	while(cin>>N>>M>>T>>P>>Q)
	{
		s.clear();//注意注意
		if(N==M&&N==T&&N==T&&N==P&&N==Q) break;
		for(int i=0;i<N;i++)
			cin>>f[i];
		//将所有匹配模式的哈希值放入一个multiset中，允许元素重复，方便统计、删除元素 
		for(int t=1;t<=T;t++)
		{
			for(int i=0;i<P;i++)
				cin>>p[i];
			compute_hash(p,P,Q);
			s.insert(Hash[0][0]);//因为匹配矩阵哈希值只有一个，就是从0，0开始的哈希值，所以存入multiset中
		}	
		//求匹配对象中相应的哈希值，若multiset中有则删去，说明有匹配模式在匹配对象中出现了
		compute_hash(f,N,M);//求一下对象矩阵的哈希值
		for(int i=0;i+P<=N;i++)//枚举一下起点
			for(int j=0;j+Q<=M;j++)
				s.erase(Hash[i][j]);//删除在multiset中有的哈希值，也就是匹配到的哈希值
				
		int ans=T-s.size();//删除了所有匹配到的哈希值的multiset中哈希值个数被总矩阵数量减去，就是匹配到的矩阵数量
		cout<<"Case "<<++count<<": "<<ans<<endl;
		
	}
}