最长公共子序列LCS

Posted on | In |

二维数组c[][]记录最长公共子序列的长度，b[][]记录最长子序列的来源
c[][]右下角的值即为最长子序列的长度

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
#include<iostream>
#define MAX 100
using namespace std;
int c[MAX][MAX];
int b[MAX][MAX];

void lcs(string s1,string s2){
	int m=s1.length();
	int n=s2.length();
	for(int i=0;i<=m;i++)
		{
			c[i][0]=b[i][0]=0;
			c[0][i]=b[0][i]=0;
		}
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			if(s1[i-1]==s2[j-1]){
				c[i][j]=c[i-1][j-1]+1;
				b[i][j]=1;
			}
			else{
				if(c[i-1][j]>c[i][j-1]){
					c[i][j]=c[i-1][j];
					b[i][j]=3;
				}
				else{
					
					c[i][j]=c[i][j-1];
					b[i][j]=2;
				}
			}		
		}
	}
}
void findPath(int x,int y,string s1){
	if(x==0||y==0) return ;
	if(b[x][y]==1){
		findPath(x-1,y-1,s1);
		cout<<s1[x-1]<<" ";
	}
	else if(b[x][y]==2){
		findPath(x,y-1,s1);
	}
	else{
		findPath(x-1,y,s1);
	}
	
}
int main(){
	string s1,s2;
	freopen("12-2.txt","r",stdin);
	cin>>s1>>s2;
	
	lcs(s1,s2);
		int m=s1.length();
	int n=s2.length();
	cout<<"最长公共子序列的长度:"<<c[m][n]<<endl;
	cout<<"最长公共子序列是:"<<endl;
	findPath(m,n,s1);
	return 0;	
}

image
image