最长公共子序列的问题

谈谈最长公共子序列

   最长公共子序列（LCS）问题是比较常见的问题，它是这么定义的：一个数列 ，如果分别是两个或多个已知数列的子序列，且是所有符合此条件序列中最长的，则 称为已知序列的最长公共子序列。当然注意它与最长公共子串的区别，LCS可以是不连续的。

   举个例子：比如说两个串：L1={a,b,f,c,a,b}，L2={a,b,c,f,b,c}。那么他们的最长公共子序列的长度就是4。{a,b,f,b}就是一个最长公共子序列。

   但是怎么求出最长公共子序列的长度呢？容易想到的是暴力解决，但这也是最容易排除的方法，因为时间复杂度太高，不是高效的算法。那在这里我们用动态规划来解决这个问题。

   还是以上文中的两个序列作为依据，令F[i][j]表示扫描到L1[i],L2[j]的最长最序列的长度。

   写出状态转移方程：F[i][j]=  F[i-1][j-1]+1  (L1[i]==L2[j])

                               max(F[i-1][j],F[i][j-1]) (L1[i]!=L2[j])

   有了状态转移方程就可以写出代码了（具体的说明示例代码中有注释）

   

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
using namespace std;
#define max(a,b) (a)>(b)?(a):(b)

const int MAX=1000+5;

char L1[MAX],L2[MAX];
int  g[MAX],g1[MAX];//滚动数组代替二维数组，减少了空间的开销

void LSC()
{
   int l1=strlen(L1),l2=strlen(L2);

   for(int i=1;i<=l1;i++)
   {
       for(int j=1;j<=l2;j++)
       {
           if(L1[i-1]==L2[j-1])//状态转移方程
           {
               g1[j]=g[j-1]+1;
           }else
           {
               g1[j]=max(g1[j-1],g[j]);
           }
       }
       for(int j=1;j<=l2;j++)//滚动数组，每次更新一次状态
           {
//               cout<<g1[j];
               g[j]=g1[j];
           }
//           cout<<endl;
   }

   cout<<g1[l2]<<endl;
}

int main()
{
//    freopen("s","r",stdin);
    memset(g,0,sizeof(g));//数组清零

    cin>>L1>>L2;
//    cout<<L1<<endl<<L2<<endl;
    LSC();

    return 0;
}