輸出兩個字符串的最長公共子串、最長公共子序列

輸出兩個字符串的最長公共子串和最長公共子序列（不單單是求長度）。

求解兩個字符串的最長公共子串和最長公共子序列在方法上很接近，都是動態規劃。只不過在遞推方程上有一些不同。

 

輸出兩個字符串的最長公共子串

#include <bits/stdc++.h>
using namespace std;
 
string LCS(const string& str1, const string& str2){
    if(str1.empty() || str2.empty()){
        return "";
    }
     
    int indexMax=0, maxn=0;
    vector<vector<int> > L(str1.size(), vector<int>(str2.size(),0) );   //所有初始化爲0
    //L[i][j]表明 str1[0~i-1]和str2[0~j-1] 的最長公共子串的長度  
 
    for(int i=0; i<str1.length(); i++){
        for(int j=0; j<str2.length(); j++){
            if(str1[i] == str2[j] ){
                if(i==0 || j==0){
                    L[i][j]=1;
                }
                else{
                    L[i][j]=L[i-1][j-1]+1;
                }
            }
            //else 是str1[i]!=str2[j]的狀況,這種狀況下L[i][j]=0,因爲初始化已經將其設置爲0,因此這裏再也不寫。
 
            //處理完L[i][j]以後,查看一下是否須要記錄下來
            if(L[i][j] > maxn){
                maxn=L[i][j];   //記錄下最長公共子串的長度
                indexMax=i;     //記錄下出現「最長公共子串」時的末尾字符的位置
            }
        }
    }
    return str1.substr(indexMax+1-maxn, maxn);  
    //字符串截取的長度有(end-start+1) = maxn, 那麼start = end+1-maxn 
    // indexMax - (indexMax+1-maxn) + 1 = maxn, maxn即爲所截取的字符串的長度。
} 
 
int main(){
    string str1, str2;
    string cur;
    while( getline(cin,cur) ){
        if(str1.empty()) str1 = cur;
        else if(str2.empty()) str2 = cur;
         
        if(!str1.empty() && !str2.empty() ){ 
            cout<< LCS(str1, str2) <<endl;
             
            str1.clear();
            str2.clear();
        }
        cur.clear();
    }  
}

　　

 

輸出兩個字符串的最長公共子序列。公共子序列和公共子串的區別在於，序列不要求連續，子串要求是其在原串中是連續的。

參考連接：https://blog.csdn.net/u013074465/article/details/45392687

 

#include <bits/stdc++.h>
using namespace std;

//輸出最長公共子串（LCS）
//二維數組veca記錄的是兩個字符串Xi和Yj的LCS長度
int LCS_length(const string &str1, const string &str2, 
	            vector<vector<int> > &veca, vector<vector<int> > &vecb) {
	int i, j; 
	if (str1 == "" || str2 == "")
		return 0;
 
	for (i = 0; i <= str1.length(); i++) {
		veca[i][0] = 0;
	}
	for (j = 0; j <= str2.length(); j++) {
		veca[0][j] = 0;
	}
	for (i = 1; i <= str1.length(); i++) {
		for (j = 1; j <= str2.length(); j++) {
			//若是Xi-1 == Yj-1，那麼最長子序列爲veca[i - 1][j - 1] + 1
			//此時將vecb[i][j] = 1代表str1[i-1]是子問題LCS的一個元素
			if (str1[i - 1] == str2[j - 1]) {
				veca[i][j] = veca[i - 1][j - 1] + 1;
				vecb[i][j] = 1;
			}
			else {
				if (veca[i - 1][j] >= veca[i][j - 1]) {
					veca[i][j] = veca[i - 1][j];
					vecb[i][j] = 2; //記錄方向,若是爲2,那就往左走
				}
				else {
					veca[i][j] = veca[i][j-1];
					vecb[i][j] = 3; //記錄方向,若是爲3,那就往上走 
				}
			}
		}
	}
	return veca[str1.length()][str2.length()];
}
 
//該函數用於輸出一個LCS的序列
//這裏輸出的順序是先向上尋找，再向左尋找
void PrintOneLCS(vector<vector<int> > &vecb, string &str1, string &str2, int i, int j) {
	if (i == 0 || j == 0)
		return;
	if (vecb[i][j] == 1) {
		PrintOneLCS(vecb, str1,str2, i - 1, j - 1);
		cout << str1[i - 1];  //這裏輸出str1[i-1]和str2[j-1]都是同樣的 
//		cout << str2[j - 1]; 
	}
	else if (vecb[i][j] == 2)
		PrintOneLCS(vecb, str1, str2, i -1, j);
	else
		PrintOneLCS(vecb, str1, str2, i, j - 1);
}

int main() {
	string input;
	getline(cin, input);
	stringstream ss(input);
	string str1, str2;
	ss >> str1;
	ss >> str2;
	//將veca初始化爲一個二維數組,其行列值分別爲str1和str2的長度加1
	//二維數組veca記錄的是兩個字符串Xi和Yj的LCS長度
	//二維數組vecb[i][j]記錄veca[i][j]時所選擇的子問題的最優解
	vector<vector<int> > veca(str1.length() + 1, vector<int>(str2.length() + 1));
	vector<vector<int> > vecb(str1.length() + 1, vector<int>(str2.length() + 1));
	cout << LCS_length(str1, str2, veca, vecb) << endl;
	PrintOneLCS(vecb, str1, str2, str1.length(), str2.length());
	return 0;
}