最长公共子序列问题

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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
#include<bits/stdc++.h>
using namespace std;

/*
最长公共子序列问题
2019-12-15 12:48:14

*/

char x[100];
char y[100];
int dp[100][100];
int b[100][100];

int LCS(char x[], char y[],int m, int n)
{
// 初始化
for (int i = 0; i <= m; i++) // 初始化第一列
{
dp[i][0] = 0;
}
for (int i = 0; i <= n; i++)
{
dp[0][i] = 0; // 初始化第一行
}

for (int i = 1; i <= m; i++)
{
for (int j = 1; j <= n; j++)
{
if (x[i] == y[j])
{
dp[i][j] = dp[i - 1][j - 1] + 1;
b[i][j] = 1;
}
else if (dp[i - 1][j] >= dp[i][j - 1])
{
dp[i][j] = dp[i - 1][j];
b[i][j] = 2;
}
else
{
dp[i][j] = dp[i][j - 1];
b[i][j] = 3;

}
}
}

return dp[m][n];

}


// 根据b矩阵得出最长公共子序列
stack<char> s;
void getLCS(int i, int j)
{
// basecase
if (i > 0 && j > 0)
{
if (b[i][j] == 1) {
s.push(x[i]);
getLCS(i - 1, j - 1);
}
if (b[i][j] == 2)
{
getLCS(i - 1, j);
}
if (b[i][j] == 3)
{
getLCS(i, j - 1);
}

}

}



int main()
{
int m, n = 0;
cin >> m >> n;

for (int i = 1; i <= m; i++)
{
cin >> x[i];
}
for (int i = 1; i <= n; i++)
{
cin >> y[i];
}

cout << LCS(x, y, m, n) << endl;



//dp矩阵
for (int i = 0; i <= m; i++)
{
for (int j = 0; j <= n; j++)
{
cout << dp[i][j] << "\t";
}
cout << endl;
}
cout << endl;


getLCS(m, n);

//打印最长公共子序列
while (!s.empty())
{
cout << s.top() << endl;
s.pop();
}

system("pause");

return 0;

}