每日一题之 hdu1159 Common Subsequence

Problem Description
A subsequence of a given sequence is the given sequence with some elements (possible none) left out. Given a sequence $X = <x1, x2, ..., xm>$another sequence $Z = <z1, z2, ..., zk>$is a subsequence of X if there exists a strictly increasing sequence <i1,i2,...,ik> < i 1 , i 2 , . . . , i k > <script type="math/tex" id="MathJax-Element-52"> </script> of indices of X such that for all $j = 1,2,...,k, xij = zj$. For example, $Z = <a, b, f, c>$ is a subsequence of $X = <a, b, c, f, b, c>$ with index sequence <1,2,4,6> < 1 , 2 , 4 , 6 > <script type="math/tex" id="MathJax-Element-56"><1, 2, 4, 6></script>. Given two sequences X and Y the problem is to find the length of the maximum-length common subsequence of X and Y.

The program input is from a text file. Each data set in the file contains two strings representing the given sequences. The sequences are separated by any number of white spaces. The input data are correct. For each set of data the program prints on the standard output the length of the maximum-length common subsequence from the beginning of a separate line.

Sample Input
abcfbc abfcab
programming contest
abcd mnp

Sample Output
4
2
0

题意：

给两个字符串，求最长公共子序列的长度。

思路：

令 $X=<x1,x2,...,xm> 和 Y=<y1,y2,...,yn>$ 为两个序列，$Z=<z1,z2,z3,...,zk>$为X和Y的任意LCS。则
如果xm=yn，则zk=xm=yn且Zk−1是Xm−1和Yn−1的一个LCS。
如果xm≠yn，那么zk≠xm，意味着Z是Xm−1和Y的一个LCS。
如果xm≠yn，那么zk≠yn，意味着Z是X和Yn−1的一个LCS。

从上述的结论可以看出，两个序列的LCS问题包含两个序列的前缀的LCS，因此，LCS问题具有最优子结构性质。在设计递归算法时，不难看出递归算法具有子问题重叠的性质。
  设C[i,j]表示Xi和Yj的最长公共子序列LCS的长度。如果i=0或j=0,即一个序列长度为0时，那么LCS的长度为0。根据LCS问题的最优子结构性质，可得如下公式：
  

$$C[i,j]= \begin{cases} 0， & 当 i=0或j=0\\ C[i-1,j-1]+1，&当i,j>0 且x_i=y_j \\ MAX(C[i,j-1],C[i-1,j])&当i,j>0且x_i≠y_j \end{cases}$$
 
 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>

using namespace std;

const int maxn = 1e3+5;

int main()
{
    string str1,str2;

    while(cin >> str1 >> str2) {
        int len1 = str1.length();
        int len2 = str2.length();
        int dp[len1+5][len2+5];
        memset(dp,0,sizeof(dp));


        for (int i = 0; i <= len1; ++i) {
            for (int j = 0; j <= len2; ++j) {
                if (i == 0 || j == 0) 
                    dp[i][j] = 0;
                else if (str1[i-1] == str2[j-1]) {
                    dp[i][j] = dp[i-1][j-1] + 1;
                }
                else 
                    dp[i][j] = max(dp[i][j-1],dp[i-1][j]);
            }
        }
        cout << dp[len1][len2] << endl;
    }

    return 0;
}