KMP算法,C语言实现

研究了两天,终于把KMP算法整明白了,贴个代码
参考链接:https://www.bilibili.com/video/BV1iJ411a7Kb?p=1&vd_source=b5fa7d9756fa49e4552b8f222264fc66

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

/* 构建 d 数组*/
static void get_next(int* d,const char* p) {
	int i = 1;
	int j = 0;
	
	d[0] = 0;					// 第一个字符不参与匹配过程

	int len_p = strlen(p);
	while (i < len_p) {	
		if (p[i] == p[j]) {		
			d[i++] = ++j;		// 字符匹配成功,此时的d[i]就是匹配成功的前后缀长度(j+1)
		} else {
			if (j > 0) {
				// 字符匹配失败,但此时并不能直接将d[i]置0,因为在p[0] ~ p[j-1]之间可能还存在公共前后缀,所以还要继续匹配
				j = d[j - 1];	
			} else {
				i++;
			}
		}
	}
}

/* kmp模式匹配 */
static int kmp(const char* t, const char* p) {
	int* d;						// d数组(表示p[0] ~ p[i] 范围之间的最大公共前后缀) 需要初始化为 0
	int i = 0;					// 主串下标索引
	int j = 0;					// 模式串下标索引

	int len_t = strlen(t);		// 主串长度
	int len_p = strlen(p);		// 模式串长度

	/*********************  打印区 方便测试  *************************/
	for ( i = 0; i < len_t; i++) {
		printf("%2d ", i);
	}
	printf("\n");
	for (i = 0; i < len_t; i++) {
		printf("%2c ", t[i]);
	}
	printf("\n");
	for (i = 0; i < len_p; i++) {
		printf("%2c ", p[i]);
	}
	printf("\n");
	/*********************  打印区 *************************/

	d = (int*)malloc(sizeof(int) * len_p);
	if (!d) {
		printf("分配失败!");
		exit(0);
	}
	memset(d, 0, len_p * sizeof(int));		

	get_next(d, p);

	i = 0;
	while (i < len_t) {
		if (t[i] == p[j]) {			// 字符匹配
			i++;
			j++;
			if (j == len_p) {		// 找到了
				//return i - j;
				return i - len_p;
			}
		} else {					// 失配
			if (j > 0) {
				j = d[j - 1];		// 保持i不变,将 j 指向 d[j - 1] (p[0] ~ p[j-1] 范围之间的最大公共前后缀的长度)
			} else {
				i++;				// j == 0,说明没有公共前后缀,将i往后面移动一位
			}
		}
	}

	free(d);
	return -1;						// 没找到,返回-1
}

int main() {
	printf("%d\n**************\n", kmp("ABABABABC", "ABABC"));
	printf("%d\n**************\n", kmp("ABABABABC", "ABAB"));
	printf("%d\n**************\n", kmp("ABABCABAB", "ABAB"));
	printf("%d\n**************\n", kmp("AAAAAAA", "AAA"));
	printf("%d\n**************\n", kmp("ABABABC", "ABABC"));
	printf("%d\n**************\n", kmp("XYXZdeOXZZKWXYZ", "WXYZ"));
	printf("%d\n**************\n", kmp("GCAATGCCTATGTGACCTATGTG", "TATGTG"));
	printf("%d\n**************\n", kmp("AGATACGATATATAC", "ATATA"));
	printf("%d\n**************\n", kmp("CATCGCGGAGAGTATAGCAGAGAG", "GCAGAGAG"));
	return 0;
}

---------------------------------------------------------------------------

测试结果:
在这里插入图片描述
在这里插入图片描述

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值