标签:blog http os io for 2014 ar amp
http://poj.org/problem?id=3690
UVA还有一道也是这样的题,LRJ给的算法是AC自动机----我还没写过,今天用HASH搞了这道题
思路很清晰,但是处理起来还因为HASH函数写混WA了几次。。。
文本矩阵n*m T个匹配矩阵p*q
思路:
1、把每一行处理出长为q的hash值
2、对于1中得到的p个哈希值在算一次哈希,这样就把一个矩阵用一个hash值替代了
3、把所有的匹配矩阵压入multiset,然后对于文本矩阵的每一个p*q的子矩阵,算出矩阵哈希值,从multiset中删去该值
4、ans=T-multiset.size()
本题可以作为字符矩阵匹配模板
#include <cstdio> #include <cstring> #include <algorithm> #include <string> #include <iostream> #include <cmath> #include <map> #include <set> #include <queue> using namespace std; #define ls(rt) rt*2 #define rs(rt) rt*2+1 #define ll long long #define ull unsigned long long #define rep(i,s,e) for(int i=s;i<e;i++) #define repe(i,s,e) for(int i=s;i<=e;i++) #define CL(a,b) memset(a,b,sizeof(a)) #define IN(s) freopen(s,"r",stdin) #define OUT(s) freopen(s,"w",stdin) const ull ba1=9973; const ull ba2=1e8+7; const int MAXT=1000+10; const int MAXP=MAXT;//55; char text[MAXT][MAXT]; char pat[100+5][MAXP][MAXP]; ull ha[MAXT][MAXT],tmp[MAXT][MAXT]; ull base1[MAXT],base2[MAXT]; int N,M,T,P,Q; void init() { base1[0]=1; base2[0]=1; for(int i=1;i<MAXT;i++) { base1[i]=base1[i-1]*ba1; base2[i]=base2[i-1]*ba2; } } void calhash(char mat[][MAXT],int n, int m) { ull e; //行 for(int i=0;i<n;i++) { for(int j=0;j<Q;j++) tmp[i][j]=(j==0)?mat[i][j]:(tmp[i][j-1]*ba1+mat[i][j]); for(int j=Q;j<m;j++) tmp[i][j]=tmp[i][j-1]*ba1+mat[i][j]-mat[i][j-Q]*base1[Q]; } //列 for(int j=Q-1;j<m;j++) { for(int i=0;i<P;i++) ha[i][j]=(i==0)?tmp[i][j]:(ha[i-1][j]*ba2+tmp[i][j]); for(int i=P;i<n;i++) ha[i][j]=ha[i-1][j]*ba2+tmp[i][j]-tmp[i-P][j]*base2[P]; } } int solve() { multiset<ull>app; for(int i=0;i<T;i++) { calhash(pat[i],P,Q); app.insert(ha[P-1][Q-1]); /*printf("#T--%d# ",i); for(int i=0;i<P;i++) { for(int j=0;j<Q;j++) cout << tmp[i][j] << "/" << ha[i][j] << " "; putchar('\n'); }*/ } calhash(text,N,M); for(int i=P-1;i<N;i++) for(int j=Q-1;j<M;j++) app.erase(ha[i][j]); return T-app.size(); } int main() { //IN("poj3690.txt"); init(); int ic=0,ans; while(~scanf("%d%d%d%d%d",&N,&M,&T,&P,&Q) && N+M+T+P+Q) { for(int i=0;i<N;i++) scanf("%s",text[i]); for(int i=0;i<T;i++) { for(int j=0;j<P;j++) scanf("%s",pat[i][j]); } printf("Case %d: %d\n",++ic,solve()); } return 0; }
poj 3690 字符矩阵匹配----HASH算法,布布扣,bubuko.com
标签:blog http os io for 2014 ar amp
原文地址:http://blog.csdn.net/u011026968/article/details/38495821