KMP算法

题目来源】https://vjudge.net/problem/HDU-1711 
【题意】 
给出两列数字,问第二个序列是否存在于第一个序列中,若存在,输出最小的左端点,不存在,输出-1 。 
【思路】 
字符串匹配问题。这几天会把这个专题做一下。。这是第一题。。。 
起初学习KMP的时候,是在大一过完年刚来的时候,那时候感觉可以学的,但是学长的经验告诉我们,不适合我们学。 
但是,现在的我觉得都可以学,并没有什么神奇的地方,只是通过一个next数组记录相应的编号而已。 
在这里,分享一下哔哩哔哩上的KMP视频,感觉不错的。 
总结两点: 
1、KMP算法里利用前缀和后缀来缩短匹配次数。 
2、通过next数组以O(n)的时间复杂度高效的查询前缀和后缀。

【代码】

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<limits.h>
#include<algorithm>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int mod=1e9+7;
const double esp=1e-5;
typedef unsigned long long ll;
typedef long long LL;
int ma[1000000+10],mo[10000+10];
int naxt[10000+10];
int n,m;
void get_next()
{
    int i=0,j;
    j=naxt[0]=-1;
    while(i<m)
    {
        while(-1!=j&&mo[i]!=mo[j]) j=naxt[j];
        naxt[++i]=++j;
    }
}
int kmp()
{
    int i=0,j=0;
    int ans=0;
    get_next();
    while(i<n)
    {
        while(-1!=j&&ma[i]!=mo[j]) j=naxt[j];
        i++;
        j++;
        if(j>=m)
        {
            return i;
        }
    }
    return -1;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++)
            scanf("%d",&ma[i]);
        for(int i=0;i<m;i++)
            scanf("%d",&mo[i]);
        int ans=kmp();
        if(ans>=0)
            ans-=m-1;//知道右端点,求左端点。
        printf("%d\n",ans);

    }

//自己写的一些代码,不喜勿喷
//kmp模板、
#include <cstdio>
#include <iostream>
#include  <string>
using namespace std;
#define INF 0x3f3f3f3f//无穷大
#define MOD 0xc0c0c0c0//无穷小
const int maxn=100000+10;
int next[maxn];
string patten,text;
void get_next()
{
    //得到一个临时数组,对于patten来说
    int p_len=patten.size();
    int j=-1,i=0;
    next[0]=-1;//important 很重要
    while(i<p_len){
        if(j==-1||patten[i]==patten[j])
            next[++i]=++j;
        else j=next[j];
    }
}
int kmp()//在 text中patten第一次出现的位置
{
    get_next();
    int i=0;//text的下标
    int j=0;//patten的小标
    int t_len=text.size();
    int p_len=patten.size();
    while(i<t_len&&j<p_len){
        if(j==-1||text[i]==patten[j]){//patten的第一个字母不匹配或text[i]==patten[j]
            i++;j++;
        }
        else j=next[j];//不符合进行跳转
    }
    if(j==p_len)//能够遍历到patten的末尾代表着匹配成功
        return i-j;
    return -1;//没有匹配成功

}
int kmp_count()
{
    int ans=0;
    int i=0,j=0;
    if(text.size()==1&&patten.size()==1){
        if(text[0]==patten[0])  return 1;
        else return 0;
    }
    get_next();//得到一个临时数组
    for(i=0;i<text.size();i++){
        while(j>0&&text[i]!=patten[j])
            j=next[j];
        if(text[i]==patten[j])  j++;
        if(j==patten.size()){
            ans++; j=next[j];
        }
    }
    return ans;
}
int main()
{
    text="catcatcat";
    patten="cat";
    printf("%d\n",kmp());
    printf("%d\n",kmp_count());
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值