最长公共子串与最长公共子序列问题

题目：

1、最长公共子串问题（子串：连续的字符串）str ="abcde"，“ab”"cd"都为子串

2、最长公共子序列问题（子序列：保证相对顺序即可）str="abcde"，"ac" ,"ae"都为子序列。

在概念上可以说，子串是特殊的子序列。

假定有字符串str A , lenA = n, str B , lenB =m;

解法：

两个题目都是动态规划问题。

最长子串的问题，因为是子串是连续的。dp[i][j]的更新条件是

A[i]==B[j]时候，dp[i][j] = A[I-1][B-1]+1。

同时用一个flag变量记录dp中的最大值即可。

代码如下：

    int findLongest(string A, int n, string B, int m) {
        // write code here
        if(m==0||n==0)
            return 0;
        int dp[n][m];
        int ans =0;
        for(int i =0;i<m;i++)//第一行
        {
            if(A[0]==B[i])
            {
                dp[0][i] = 1;
                ans = 1;
            }
            else
                dp[0][i] = 0;
        }
        for(int i =0;i<n;i++)//第一列
        {
            if(A[i]==B[0])
            {
                dp[i][0] = 1;
                ans = 1;
            }
            else
                dp[i][0] = 0;
        }
        for(int i =1;i<n;i++)
        {
            for(int j =1;j<m;j++)
            {
                if(A[i]==B[j])
                {
                    dp[i][j] = dp[i-1][j-1]+1;
                    ans = max(ans,dp[i][j]);
                }
                else
                    dp[i][j] =0;
            }
        }
        return ans;
    }

针对公共子序列:

dp[i][j]可以取两种情况，A[i]!=B[i]的时候，dp[i][j] =max(dp[i-1][j],dp[i][j-1])；

A[i] ==B[j]的时候，dp[i][j] = dp[i-1][j-1] +1;又因为dp[i][j]存放的数字表示：

以A[0]~A[i]的子串 与 B[0]~B[j]的子串的最长子序列。所以

dp[i][j] = max(dp[i-1][j-1] +1 , max(dp[i-1][j],dp[i][j-1]) );

dp[i-1][j-1] +1 表示当子序列的最后一个字母为A[i](或者B[j]，两者相等)的最长子序列。

但是有可能该子序列的长度没有 以A[i-1]或B[j-1]结尾的子序列长度长。

代码：

    int findLCS(string A, int n, string B, int m) {
    // write code here
        if(n==0||m==0)
            return 0;
        int dp[n][m];
        int ans =0;
        dp[0][0] = (A[0]==B[0]?1:0);
        for(int i =1;i<n;i++)
        {
            dp[i][0] = max(dp[i-1][0],A[i]==B[0]?1:0);
            ans = max(ans,dp[i][0]);
        }
        for(int i=1;i<m;i++)
        {
            dp[0][i] = max(dp[0][i-1],A[0]==B[i]?1:0);
            ans = max(ans,dp[0][i]);
        }
        for(int i =1;i<n;i++)
        {
            for(int j =1;j<m;j++)
            {
                dp[i][j] = max(dp[i-1][j],dp[i][j-1]);
                ans = max(ans,dp[i][j]);
                if(A[i]==B[j])
                {
                    dp[i][j] = max(dp[i][j],dp[i-1][j-1]+1);
                    ans = max(ans,dp[i][j]);
                }
            }
        }
        return ans;
    }