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;  
}   

[cpp] view plain copy
/*===================================================*\  
    从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;    
}