转载c++实现kmp

• 1. #include <iostream>  
  2. #include <vector>  
  3. #include <functional>  
  4. #include <algorithm>  
  6. using namespace std;
  7. template<typename FType, typename SType>
  8. struct equal_cmp
  9. {
  10.   bool operator ()(const FType& _F, const SType& _S)const
  11.   {return _F == _S;}  
  12. };  
  15. template<typename RandAccessConcept,
  16.          typename EqualCmp>
  17. inline void     
  18. CalFailFunction(
  19.                 const RandAccessConcept& in,
  20.                 int len,
  21.                 vector<int>& ff,
  22.                 EqualCmp cmp
  23.                 )
  24. { 
  25.     ff.clear();
  26.     ff.reserve(len);
  27.     ff.push_back(-1);
  28.     for (int curr = 1, k = -1; curr < len; ++curr)
  29.     {
  30.         for (;k >= 0 && !cmp(in[k+1], in[curr]);
  31.                 k = ff[k]);
  32.         if (cmp(in[k+1], in[curr])) ++k;
  33.         ff.push_back(k);
  34.     }   
  35. }  
  37. template<typename FRandAccessConcept,
  38.          typename SRandAccessConcept,
  39.          typename  EqualCmp>
  40. inline int      
  41. KMPSearchFirst(
  42.                 const FRandAccessConcept& d,
  43.                 int   dlength,
  44.                 const SRandAccessConcept& p,
  45.                 int   plength,
  46.                 EqualCmp cmp,
  47.                 int start = 0
  48.                 )
  49. {
  50.     vector<int> ff;
  51.     CalFailFunction(p, plength, ff, cmp);
  52.     for (int i = start, j = -1; i < dlength; ++i)
  53.     {
  54.         while (j >= 0 && !cmp(d[i], p[j+1])) j = ff[j];
  55.         if (cmp(d[i], p[j+1])) ++j;
  56.         if (j == plength - 1) return i - plength + 1;
  57.     }
  58.     return -1; 
  59. }    
  61. template<typename FRandAccessConcept,
  62.          typename SRandAccessConcept,
  63.          typename EqualCmp>
  64. inline int      
  65. KMPSearchALL(
  66.                 const FRandAccessConcept& d,
  67.                 int   dlength,
  68.                 const SRandAccessConcept& p,
  69.                 int   plength,
  70.                 vector<int>& result,
  71.                 EqualCmp cmp
  72.                 )
  73. {
  74.     vector<int> ff;
  75.     CalFailFunction(p, plength, ff, cmp);
  76.     result.clear();
  77.     for (int i = 0, j = -1; i < dlength; ++i)
  78.     {
  79.         while (j >= 0 && !cmp(d[i], p[j+1])) j = ff[j];
  80.         if (cmp(d[i], p[j+1])) ++j;
  81.         if (j == plength - 1)
  82.         {
  83.              result.push_back(i - plength + 1);
  84.              j = ff[j];
  85.         }    
  86.     }
  87.     return result.size();
  88. }      
  90. int main()
  91. {
  92.     int d[] = {1,2,3,1,2,3,1,2};
  93.     char p[] = "/1/2/3/1/2";
  94.     vector<int> result;
  95.     int found = KMPSearchALL(d, 8, p, 5, result, equal_cmp<int, char>());
  96.     
  97.     for (int i = 0; i < found; ++i)
  98.         cout << result[i] << " ";
  99.     cout << endl;
  100.     
  101.     return 0;
  102. }    

   

  转自 http://blog.csdn.net/baihacker/archive/2008/08/20/2803362.aspx

   

  我以前用c写过，有空的时候可以贴出来