Rabin-Karp字符串匹配算法c源代码

/*Rabin-Karp字符串匹配算法*/
/*Author:Milo
 *Date:2011/12/22
 */
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
/*参数说明
*********************
T:要查找的字符串
P:模式字符串
d:查找字符串的基数
q:一个素数
*********************
*/
static void rabin_karp_macher(const char *T,const char *P,int d,int q)
{
	int n=strlen(T);
	int m=strlen(P);
	int p=0;/*存放将模式字符串转换为一个大整数后对q求模的结果*/
	int ts=0;/*存放T中与P等长度的一段字符串转换为大整数后对q求模的结果*/
	int i;
	for(i=0;i<m;i++){
		p=(d%q*p+P[i]-'a')%q;/*求p*/
		ts=(d%q*ts+T[i]-'a')%q;/*求以T[0]开始的长为m的字符串转换为大整数后对q求模的结果*/
	}
	int s=0;
	int h=((int)pow((double)d,m-1))%q;/*求d^(m-1)%q的值*/
	while(s<=n-m){
		if(p==ts){/*匹配字符串的模值与待匹配字符串的子串模值相等，当不一定两个字符串相等*/
			if(!strncmp(T+s,P,m))/*此时需要看两个字符串是否完全相等*/
				printf("match at index: %d\n",s+1);
		}
		if(s<n-m){/*调整待匹配字符串的字串，前进一位，求新字串的大整数模值*/
			ts=(d*(ts-h*(T[s]-'a'))+(T[s+m]-'a'))%q;
			if(ts<0)/*模值有可能为0，需要调整为整数*/
				ts=ts+q;
		}
		s++;
	}
}

int main(int argc,char *argv[])
{
	char T[]="lqmwjjszbhellowjj";
	char P[]="wjj";
	rabin_karp_macher(T,P,26,29);
	return 0;
}