《編程之美》：最長公共子序列(串)

求最長公共子串有兩種情況：1, 不要求子串連續，兩個字符串BDCABA和ABCBDAB，字符串BCBA和BDAB都是它們的最長公共子串。2, 要求子串連續。對於上面兩個串來說，BD和AB是它們的最長公共子串。
一，有關連續的已經寫過。設定一個矩陣p，p[i][j]表示以str1[i]和str2[j]結尾的子串的最大公共子串的長度。所以有遞推式：
如果 str1[i] != str2[j], p[i][j] = 0.
如果 str1[i]  = str2[j], p[i][j] = p[i-1][j-1] + 1.
二，如果不要求子串連續，處理有一些小的區別。還是設定一個矩陣p，p[i][j]表示str1[i]和str2[j]之前的子串最大公共子串的長度(不需要一定包含str1[i]或str2[j])，也就是說即使str1[i]和str2[j]不相等，也不能把p[i][j]清零。因爲題目要求的子串不一定是連續的。所以這種情況的遞推式爲：
如果 str1[i] != str2[j], p[i][j] = max(p[i-1][j],  p[i][j-1, p[i-1][j-1]).
如果 str1[i]  = str2[j], p[i][j] = p[i-1][j-1] + 1.
而p[i-1][j-1]肯定小於前兩者，所以第一個條件可改爲：p[i][j] = max(p[i-1][j],  p[i][j-1])。
爲了避免重複計算或者浪費不必要的空間，可以從開始往後遞推。#include <iostream> #include <string.h> using namespace std; enum DRECT { LEFT = 1, UP, LEFT_UP }; int max(int i, int j) { return i > j ? i : j; } void print(int** direct, const char* str1, int row, int col) { if (row < 0 || col < 0) { return; } if (direct[row][col] == LEFT_UP) { print(direct, str1, row - 1, col - 1); cout << str1[row]; } else if (direct[row][col] == LEFT) { print(direct, str1, row, col - 1); } else if (direct[row][col] == UP) { print(direct, str1, row - 1, col); } } int maxCommStr(char* str1, char* str2) { int len1, len2, i, j; len1 = strlen(str1); len2 = strlen(str2); int** pf = new int*[len1]; //根據字符串長度分配一個二維數組 for (i = 0; i<len1; i++) { pf[i] = new int[len2]; } for (i = 0; i<len1; i++) { for (j = 0; j<len2; j++) { pf[i][j] = 0; //數組初始化 } } int** direction = new int*[len1]; //根據字符串長度分配一個二維數組,用來存放方向 for (i = 0; i<len1; i++) { direction[i] = new int[len2]; } for (i = 0; i<len1; i++) { for (j = 0; j<len2; j++) { direction[i][j] = 0; //數組初始化 } } for (i = 0; i<len1; i++) { for (j = 0; j<len2; j++) { if (i == 0 || j == 0) //處理以第一個字符結尾(長度爲1)的情況 { if (str1[i] == str2[j]) { pf[i][j] = 1; direction[i][j] = LEFT_UP; } else { pf[i][j] = 0; } } else { if (str1[i] == str2[j]) { pf[i][j] = pf[i-1][j-1] + 1; direction[i][j] = LEFT_UP; } else if (pf[i-1][j] > pf[i][j-1]) { pf[i][j] = pf[i-1][j]; direction[i][j] = UP; } else { pf[i][j] = pf[i][j-1]; direction[i][j] = LEFT; } } } } print(direction, str1, len1 - 1, len2 - 1); cout << endl; return pf[len1-1][len2-1]; } int main() { char* str1 = "BDCABA"; char* str2 = "ABCBDAB"; cout << maxCommStr(str1, str2) << endl; return 0; }