最长公共子序列（LCS）问题－－动态规划

一、问题描述

一、什么是最长公共子序列

    什么是最长公共子序列呢？举个简单的例子吧，一个数列S，若分别是两个或多个已知序列的子序列，且是所有符合条件序列中最长的，则S称为已知序列的最长公共子序列。

  举例如下，如：有两个随机数列，1 2 3 4 5 6 和 3 4 5 8 9，则它们的最长公共子序列便是：3 4 5。

  最长公共子串（Longest Common Substirng）和最长公共子序列（Longest Common Subsequence，LCS）的区别为：子串是串的一个连续的部分，子序列则是从不改变序列的顺序，而从序列中去掉任意的元素而获得新的序列；也就是说，子串中字符的位置必须是连续的，子序列则可以不必连续。

二、蛮力法

   蛮力法是解决最长公共子序列问题最容易想到的方法，即对S的每一个子序列，检查是否为T的子序列，从而确定它是否为S和T的公共子序列，并且选出最长的公共子序列。

   S和T的所有子序列都检查过后即可求出S和T的最长公共子序列。S的一个子序列相应于下标序列1,2，...，n的一个子序列。因此，S共有2^n个子序列。当然，T也有2^m个子序列。

   因此，蛮力法的时间复杂度为O(2^n * 2^m)，这可是指数级别的啊。

三、动态规划方法

   1、序列str1和序列str2

  ·长度分别为m和n；

  ·创建1个二维数组L[m.n]；

    ·初始化L数组内容为0

    ·m和n分别从0开始，m++，n++循环：

       - 如果str1[m] == str2[n]，则L[m,n] = L[m - 1, n -1] + 1；

       - 如果str1[m] != str2[n]，则L[m,n] = max{L[m,n - 1]，L[m - 1, n]}

    ·最后从L[m,n]中的数字一定是最大的，且这个数字就是最长公共子序列的长度

    ·从数组L中找出一个最长的公共子序列

   2、从数组L中查找一个最长的公共子序列

   i和j分别从m，n开始，递减循环直到i = 0，j = 0。其中，m和n分别为两个串的长度。

  ·如果str1[i] == str2[j]，则将str[i]字符插入到子序列内，i--，j--；

  ·如果str1[i] != str[j]，则比较L[i,j-1]与L[i-1,j]，L[i,j-1]大，则j--，否则i--；（如果相等，则任选一个）

 

   根据上图，我们可以得到其中公共子串：B C B A 和 B D A B。

 

二、代码实现

  

#include <iostream>

#include <string>

using namespace std;

int main(int argc, char **argv)

{

string str1 = "ABCBDAB";

string str2 = "BDCABA";

int x_len = str1.length();

int y_len = str2.length();

int arr[50][50] = {{0,0}};

int i = 0;

int j = 0;

for(i = 1; i<= x_len; i++)

{

for(j = 1; j<= y_len; j++)

{

if(str1[i- 1] == str2[j- 1])

{

arr[i][j] = arr[i- 1][j- 1] + 1;

}

else

{

if(arr[i][j- 1] >= arr[i- 1][j])

{

arr[i][j] = arr[i][j- 1];

}

else

{

arr[i][j] = arr[i-1][j];

}

}

}

}

for(i = 0 ; i<= x_len; i++)

{

for( j = 0; j<= y_len; j++)

 {

cout << arr[i][j] << "  ";

 }

cout << endl;

}

for(i = x_len, j= y_len; i>= 1 && j>= 1;)

 {

 if(str1[i- 1] == str2[j- 1])

  {

cout << str1[i- 1] << " ";//倒序打印的

i--;

j--;

}

else

{

 //  if(arr[i][j -1] >= arr[i - 1][j])//打印：B A D B

if(arr[i][j-1] > arr[i- 1][j]) //打印:A B C B

 {

 j--;

}

 else

 {

  i--;

 }

  }

}

cout << endl;

 return 0;

 }