整理的一些模版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+1);
    for(int i=0;i<=lena;i++)
      mat[i].resize(lenb+1,0);
    for(i=1;i<=lena;i++){
      for(j=1;j<=lenb;j++){
        if(a[i-1]==b[j-1])
          mat[i][j]=mat[i-1][j-1]+1;
        else
          mat[i][j]=max(mat[i-1][j],mat[i][j-1]);
        }
    }
    i=lena-1;
    j=lenb-1;
    int k=0;
 while(i!=-1&&j!=-1) {
    if(a[i]==b[j]){
       sav[k++]=a[i];
        i--;
        j--;
  }
  else{
    if(mat[i+1][j+1]==mat[i][j]){
        i--;
        j--;
   }
   else{
    if(mat[i][j+1]>=mat[i+1][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=0,x=0;
    vector<vector<int> >lcs(lena+1);
    for(int i=0;i<=lena;i++)
      lcs[i].resize(lenb+1,0);
    for(i=1;i<=lena;i++){
      for(j=1;j<=lenb;j++)
        if(a[i-1]==b[j-1]){
           lcs[i][j]=lcs[i-1][j-1]+1;
           if(k<lcs[i][j]){
               x=i;
               k=lcs[i][j];
           }
        }
    }
    strncpy(sav,a+(x-k),k);
    return k;
}