模式匹配中的kmp算法

//:kmp
//kmp算法是效率最高的模式匹配算法
//其主要思想是：（假设i，j分别指向主串mStr和模式sStr的当前试匹配位置）当mStr[i]!=mStr[j]时并不
//按传统的回溯i,j ，而是i不变j重一个算好的数组回溯，该数组可使得模式串向右
//移动尽可能多的距离。
//主要原理：假设当i，j失配的时候应该用模式的k字符与主串的j字符进行比较
//则：  sStr[1~k-1] = mStr[i-k+1 ~ i-1]
//又已知：sStr[j-k+1 ~ j-1]=mStr[i-k+1 ~ i-1]
//所以： sStr[1~k-1] = sStr[j-k+1 ~ j-1];
//所以可以根据字串计算出每个字符对应的最大k,储存在next[]中,简便方法为：
//1,初始化 next[0] = -1 ,i,j分别代表sStr[0~i-1] = sStr[j-i-1,j-1]情况下的i,j且j+1是当前失配位置
//2，如果sStr[i] ==  sStr[j] ; 则sStr[0~i] = sStr[j-i-1,j]，如果此时sStr[0~i+1] = sStr[j-i-1,j+1]
//则当j+1位置失配时，推到i+1位置也会失配,所以应使next[j+1] = next[i+1]如果sStr[0~i+1] != sStr[j-i-1,j+1]
//则next[j+1]=i+1;
//3,如果j<strlen(sStr)重复2
//4，退出


#include<iostream>
#include<cstdio>
int * getIndex(const char *pstr,int *next)//及时一个串的模式匹配数组
{
    next[0]=-1,next[1]=0;
   
    int i = 0,j=1;
   
    while(j < strlen(pstr))
    {
        if(i == -1 || pstr[i]==pstr[j])
        {
            //如果i==-1则将模式的第一个字符和主串的下一个字符比较
            //如果pstr[i]==pstr[j]比较后续字符
            ++i,++j;
            if(pstr[i] == pstr[j])
            {//sStr[0~i] = sStr[j-i-1,j]
                next[j] = next[i];
            }
            else
            {//sStr[0~i] != sStr[j-i-1,j]
                next[j]=i;
            }
        }
        else
        {
            i = next[i];
        }
    }
    return next;
}

int kmp(const char * mainSour,const char * subStr,int pos)
{
    int res = -1;
    int * next = new int[strlen(subStr)];
    //计算数组
    getIndex(subStr,next);
   
    int i=pos,j=0,end= strlen(subStr);
   
    while(i<end)//未到达主串的尾部
    {
        if(mainSour[i] == subStr[j] || j == -1)
        {
            ++i,++j;
            if(!subStr[j])
            {
                res = i-strlen(subStr);
                break;
            }
        }
        else
        {//失配则回溯j
            j = next[j];
        }
    }
    return res;
}

int main()
{
    std::cout<<kmp("ababcabcacbab","abcac",0)<<std::endl;
    return 0;
}