字串 KMP 演算法

 

#include <iostream>

 

using namespace std;

 

void fail(string & s, int f[]);

int fastFind(string &src,string &pat,int *fail);

 

/*

fail[i] = k (0<=k<i 且滿足P[0]P[1]……P[k] = P[i-k]P[i-k+1]……P[i] 的最大整數)

         = -1  (others)

*/

 

void fail(string & s, int f[])

{

     int lengthP=s.length();

     f[0]=-1;

     for(int j=1;j<lengthP;j++)

     {

         int i=f[j-1]; // previous value

 

         while(s[j] != s[i + 1] && i > -1)

         {

              i=f[i];

         }

 

         if(s[j]== s[i + 1])

              f[j]=i+1; 

         else

              f[j]=-1;

     }

     cout << endl;

}

 

int fastFind(string &src,string &pat,int *fail)

{

     int posP=0,posT=0;

     int lengthP=pat.length(),lengthT=src.length();

     while(posP<lengthP && posT<lengthT)

     {

         if(pat.at(posP)==src.at(posT))

         {

              posP++;

              posT++;//匹配成功一位就同時向後移一位

         }

         else if(posP==0)

         {

              posT++;//第一個字元都不匹配時，被匹配字串向後移

         }

         else

         {

              posP=fail[posP-1]+1;//否則，通過失效函數找到重新匹配的起始位置

         }

     }

 

     if(posP<lengthP)

     {

         return -1; // Not found

     }

     else

     {

         return posT - lengthP;//匹配成功時返回被匹配串中的匹配上的起始位置

     }

}

 

int main( void )

{

     string src = "aabbaabcaabb";

     string t = "aabc";

     int a[20];

     fail(t, a);

     int pos = fastFind(src, t, a);

     cout << "pos = " << pos << endl;

 

     for(int k = 0; k < t.length(); ++k)

     {

         cout << k << " : " << a[k] << endl;

     }

}