字符串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串长度相同的子串，返回开始位置 
    不保证绝对正确，发生冲突概率为O(sqrt(n)), 
    n为哈希函数的最大值 
\*===================================================*/  
  
#define ull unsigned long long  
const ull B = 1e8+7;    /*according to the book*/  
const int MAXN = 1001;  
char a[MAXN],b[MAXN];  
  
int HashFind()/*返回匹配位置的开始处的下标,-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;  
}  
  
/*===================================================*\ 
    滚动哈希--滚动通过ull的溢出实现,相当于mod 2^64 
 
\*===================================================*/  
  
#define ull unsigned long long  
const int MAXN = 1001;  
const ull B = 1e8+7;    /*according to the book*/  
char a[MAXN],b[MAXN];  
  
int OverLap()/*a的后缀和b的前缀相等的最大长度*/  
{  
    int al=strlen(a),bl=strlen(b),ans=0;  
    int minlen=min(al,bl);  
    ull t=1,ah=0,bh=0;  
    for(int i=1;i<minlen;i++)  
    {  
        ah=ah+a[al-i]*t;  
        bh=bh*B+b[i-1];  
        if(ah == bh)ans=i;  
        t*=B;  
    }  
    return ans;  
}