最长公共子序列（LCSubstring）和最长公共子串（LCSunsequence）问题

最长公共子序列（Longest Common Subsequence）与最长公共子串（Longest Common Substring）的区别： 子序列问题只需保持字符相对顺序一致，并不要求连续，而子串问题要求字符在原字符串中是连续的。

一、最长公共子序列（Longest Common Subsequence）（动态规划)

根据下列公式来编写程序：

程序步骤如下：
1、构造dp数组（二维数组多出一行一列来保存dp初始状态的0值）；
2、正序构造dp数组：根据公式给dp数组赋值（注意下标值(str1[i-1]与str2[j-1]）与dp[i][j]对应)）；
3、结果输出：倒序遍历dp数组，如果当前位置的str1[i-1]==str2[j-1]， 则将该点对应的值插入的结果字符串的头部并且i--; j--;，如果当前位置的dp值与上方或左方的dp值相同，则向上或向左移动。如此循环，直到i < 0 || j < 0。 （因为最长公共子序列不一定唯一，这样只是输出一种方案，所有的LCS输出见PrintAllLCS()）。
4、输出LCSubsequence；

#include <iostream>
#include <string>
#include <vector>
#include <set>
using namespace std;
set<string> allRes;

int Max(int a, int b) {
    return a > b ? a : b;
}

void PrintAllLCSubsequence(string &str1, string &str2, int i, int j,
                 vector<vector<int>> &dpvec, string res) {
    //print all LCS
    if (i == 0 || j == 0) {
        if (allRes.find(res) == allRes.end()) {
            cout << res << endl;
            allRes.insert(res);
        }
        return ;
    }
    if (i > 0 && j > 0) {
        if (str1[i-1] == str2[j-1]) {
            string tmp = str1[i-1] + res;
            PrintAllLCSubsequence(str1, str2, i-1, j-1, dpvec, tmp);
        } else {
            if (dpvec[i][j] == dpvec[i-1][j])
                PrintAllLCSubsequence(str1, str2, i-1, j, dpvec, res);

            if (dpvec[i][j] == dpvec[i][j-1])
                PrintAllLCSubsequence(str1, str2, i, j-1, dpvec, res);
        }
    }
}

void PrintOneLCSubsequence(string &str1, string &str2, vector<vector<int>> &dpvec) {
    //print one LCS
    string res = "";
    int i = str1.size();
    int j = str2.size();

    while (i > 0 && j > 0) {
        if (str1[i-1] == str2[j-1]) {
            res = str1[i-1] + res;
            i--;
            j--;
        } else {
            if (dpvec[i][j] == dpvec[i-1][j])
                i--;
            else if (dpvec[i][j] == dpvec[i][j-1])
                j--;
        }
    }
    cout << res << endl;
}

int LCSubsequence(string &str1, string &str2) {
    int len1 = str1.size();
    int len2 = str2.size();
    if (len1 == 0 || len2 == 0)
        return 0;

    vector<vector<int>> dp(len1+1, vector<int>(len2+1));
    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-1][j], dp[i][j-1]);
        }
    }
    PrintOneLCSubsequence(str1, str2, dp);
    //PrintAllLCSubsequence(str1, str2, len1, len2, dp, "");
    return dp[len1][len2];
}

二、最长公共子串（Longest Common Substring）（动态规划）

这里的最大公共字串要求的字串是连续的。
求字串的方法和求子序列方法类似：

当str1[i] == str2[j]时，子序列长度dp[i][j] = dp[i - 1][j - 1] + 1；只是当str1[i] != str2[j]时，dp[i][j] = 0，而不是max{dp[i-1][j], dp[i][j-1]}。

#include <iostream>
#include <string>
#include <vector>
using namespace std;

void PrintOneLCSubstring(string &str1, string &str2, vector<vector<int>> &dpvec, int maxLen) {
     //print LCS
    if (maxLen != 0) {
        for (int i = str1.size(); i > 0; i--) {
            for (int j = str2.size(); j > 0; j--) {
                if (dpvec[i][j] == maxLen) {
                    string res = "";
                    int a = i;
                    int b = j;
                    while (a > 0 && b > 0) {
                        res = str1[a-1] + res;
                        --a;
                        --b;
                    }
                    cout << res << endl;
                    return ;
                }
            }
        }
    }

    return ;

}
int LCSubstring(string &str1, string &str2) {
    int len1 = str1.size();
    int len2 = str2.size();
    if (len1 == 0 || len2 == 0)
        return 0;
    int maxLen = 0;
    vector<vector<int>> dp(len1+1, vector<int>(len2+1));
    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;
                if (dp[i][j] > maxLen)
                    maxLen = dp[i][j];
            } else
                dp[i][j] = 0;
        }
    }
    PrintOneLCSubstring(str1, str2, dp, maxLen);
    return maxLen;
}

int main() {
	string X = "asdfas";
	string Y = "werasdfaswer";
	cout << "The length of LCS is " << LCSubstring(X, Y);
	cout << endl;
	return 0;
}