最长公共子序列

Sample Input

  
  

   
   abcfbc abfcab
programming contest 
abcd mnp
  
  

Sample Output

  
  

   
   4
2
0
  
  

题意：求最长公共子序列。

分析：假设S1...Si，T1...Ti对应两个字符串序列，定义dp[i][j]为LCS的长度，则S1...S(i+1)和T1...T(j+1)对应的公共子序列有2种状 态：1.当S(i+1)=T(j+1)时，相当于在S1...Si和T1...Tj 的最长公共子列末尾追加上S(i+1)（这个思想很重要）。2.如果不等的话，就是相当于求S1...Si和T1...T(j+1)的公共子序列或者求S1...S(i+1)和T1...Tj的公共子序列。

利用dp思想，下面给个复杂度O（m*n）的算法：

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int main()
{
    int dp[500][500];
    char a[500], b[500];

    while(~scanf("%s%s", a,b))
    {
        memset(dp, 0, sizeof(dp));
        int n = strlen(a), m = strlen(b);
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<m; j++)
            {
                if(a[i] == b[j]) dp[i+1][j+1] = dp[i][j]+1;
                else dp[i+1][j+1] = max(dp[i+1][j], dp[i][j+1]);
            }
        }
        printf("%d\n", dp[n][m]);
    }
    return 0;
}

好像还有n*logn算法，不过现在还不会，等碰到了再说吧~

还有具体的思想可以借鉴下面的文章：

动态规划 (Dynamic Programming) 之 最长公共子序列(Longest Common Subsequence)