整理的一些模版LCS(连续和非连续)

对于连续的最大串，我们称之为子串....非连续的称之为公共序列..

代码：

非连续连续

 int LCS(char a[],char b[],char sav[]){
     int lena=strlen(a);
     int lenb=strlen(b);
     int i,j;
      vector<vector<int> >mat(lena+);
     for(int i=;i<=lena;i++)
       mat[i].resize(lenb+,);
     for(i=;i<=lena;i++){
       for(j=;j<=lenb;j++){
         if(a[i-]==b[j-])
           mat[i][j]=mat[i-][j-]+;
         else
           mat[i][j]=max(mat[i-][j],mat[i][j-]);
         }
     }
     i=lena-;
     j=lenb-;
     int k=;
  while(i!=-&&j!=-) {
     if(a[i]==b[j]){
        sav[k++]=a[i];
         i--;
         j--;
   }
   else{
     if(mat[i+][j+]==mat[i][j]){
         i--;
         j--;
    }
    else{
     if(mat[i][j+]>=mat[i+][j]) i--;
      else j--;
    }
   }
  }
  return k;
 }

连续子串

代码：

 int LCS(char a[],char b[],char sav[]){
     int lena=strlen(a);
     int lenb=strlen(b);
     int i,j,k=,x=;
     vector<vector<int> >lcs(lena+);
     for(int i=;i<=lena;i++)
       lcs[i].resize(lenb+,);
     for(i=;i<=lena;i++){
       for(j=;j<=lenb;j++)
         if(a[i-]==b[j-]){
            lcs[i][j]=lcs[i-][j-]+;
            if(k<lcs[i][j]){
                x=i;
                k=lcs[i][j];
            }
         }
     }
     strncpy(sav,a+(x-k),k);
     return k;
 }