字符串HASH模板

/* 
    从b串中寻找和a串长度相同的子串，返回开始位置 
    不保证绝对正确，发生冲突概率为O(sqrt(n)), 
    n为哈希函数的最大值 
*/ 
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#define ull unsigned long long  
const ull B = 1e8+7;     
const int MAXN = 5010;  
char a[MAXN],b[MAXN];  
  
int find()/*返回匹配位置的开始处的下标,-1表示没有匹配*/  
{  
    int al=strlen(a),bl=strlen(b);  
    if(al>bl)return -1;  
    ull t=1,ah=0,bh=0;  
    for(int i=0;i<al;i++)  
    {  
        t*=B;  
        ah=ah*B+a[i];  
        bh=bh*B+b[i];  
    }  
    for(int i=0;i+al<=bl;i++)  
    {  
        if(ah == bh)return i;  
        if(i+al < bl)bh=bh*B+b[i+al]-b[i]*t;  
    }  
    return -1;  
}  

int main(){
	freopen("strhash.in","r",stdin);
	freopen("strhash.out","w",stdout);	
	gets(a);
	gets(b);
	printf("%d\n",find());
	return 0;
} /* 
    从b串中寻找和a串长度相同的子串，返回开始位置 
    不保证绝对正确，发生