我想用C ++解决三个int序列的最长公共子序列。这个问题很经典：
任务即可。给定三个序列A =（a1，a2，...，an），B =（b1，b2，...，bm），并且C =（c1，c2，...，cl），找到长度其 最长的公共子序列，即最大的非负整数p 这样就存在指数1≤i1<1。 i2＆lt; ···＆lt; ip≤n，1≤j1< j2＆lt; ···＆lt; jp≤m，1≤k1<1。 k2＆lt; ···＆lt; kp≤l使得ai1 = bj1 = ck1 ,. 。 。 ，aip = bjp = ckp
输入格式。第一行：n。第二行：a1，a2 ,. 。 。 ，一个。第三行：m。第四行：b1，b2 ,. 。 。 ，bm。第五行：l。第六 line：c1，c2 ,. 。 。 ，cl。
约束即可。 1≤n，m，l≤100; -109＆lt; ai，bi，ci＆lt; 109.
输出格式。输出p。
对于超过29个案例的28个我的解决方案是可以的，但是，我找不到错误
这是我的代码：
#include <iostream>
#include <vector>
#include<algorithm>
using namespace std;
using std::vector;
int get_position(int alpha,vector<pair<int,int>> vect, int A,int B ) {
if (alpha == 8) {
int gamma = 0;
}
if (A != B) {
if (vect[A].first == alpha) {
return vect[A].second;
}
if (vect[B].first == alpha) {
int petitb = B;
while ((petitb > A) && (vect[petitb].first == alpha))
petitb--;
if (vect[petitb].first == alpha) return vect[petitb].second;
else {
petitb++;
if (vect[petitb].first == alpha) return vect[petitb].second;
}
return vect[B].second;
}
int mid = (A + B) / 2;
if (vect[mid].first == alpha) {
return vect[mid].second;
}
if (vect[mid].first > alpha) {
if (mid == B) return -1;
return get_position(alpha, vect, A, mid);
}
else {
if (mid == A) return -1;
return get_position(alpha, vect, mid, B);
}
}
else {
if (vect[A].first == alpha) {
return vect[A].second;
}
}
return -1;
}
int get_relative_pos(int alpha, vector<pair<int, int>> vect,vector<int> V, int A, int B) {
int absolute = get_position(alpha, vect, 0, vect.size()-1);
if ((absolute >= A) && (absolute <= B)) {
return absolute;
}
if (absolute > B)
return -1;
if (absolute < A) {
if (V[A] == alpha) { return A; }
else { return -1; }
}
return -1;
}
bool comp1(pair<int, int>A, pair<int, int>B) {
if (A.first == B.first)
return A.second < B.second;
return A.first < B.first;
}
int explore(vector<int> &a, vector<int> &b, vector<int> &c, int PA, int PB, int PC) {
int length = 0;
while ((PA<a.size()) && (PB<b.size()) && (PC<c.size())) {
if ((a[PA] == b[PB]) && (a[PA] == c[PC])) {
length++;
PA++;
PB++;
PC++;
}
else {
return length;
}
}
return length;
}
int explore2(vector<int> &a, vector<int> &b, vector<int> &c, int PA, int PB, int PC, vector<pair<int, int>> SA, vector<pair<int, int>> SB, vector<pair<int, int>> SC) {
int length = 0;
while ((PA < a.size()) && (PB < b.size()) && (PC < c.size())) {
if ((a[PA] == b[PB]) && (a[PA] == c[PC])) {
length++;
PA++;
PB++;
PC++;
}
else {
//return length;
int Bprim = get_relative_pos(a[PA], SB,b, PB, SB.size() - 1);//get_position(a[PA], SB, PB, SB.size() - 1);
int Cprim = get_relative_pos(a[PA], SC,c, PC, SC.size() - 1);//get_position(a[PA], SC, PC, SC.size() - 1);
if ((Bprim != (-1)) && (Cprim != (-1))) {
PB = Bprim;
PC = Cprim;
}
else { PA++; }
}
}
return length;
}
int cleanData(vector<int> &a, vector<int> &b, vector<int> &c, vector<pair<int, int>> SA, vector<pair<int, int>> SB, vector<pair<int, int>> SC)
{
int indexA, indexB, indexC;
for (int i = 0; i < a.size(); i++) {
indexB = get_position(a[i],SB,0,SB.size()-1);
indexC = get_position(a[i], SC, 0, SC.size() - 1);
if ((indexB == (-1)) || (indexC == (-1)))
a.erase(a.begin()+i);
}
for (int i = 0; i < b.size(); i++) {
indexA = get_position(b[i], SA, 0, SA.size() - 1);
indexC = get_position(b[i], SC, 0, SC.size() - 1);
if ((indexA == (-1)) || (indexC == (-1)))
b.erase(b.begin() + i);
}
for (int i = 0; i < c.size(); i++) {
indexB = get_position(c[i], SB, 0, SB.size() - 1);
indexC = get_position(c[i], SC, 0, SC.size() - 1);
if ((indexB == (-1)) || (indexC == (-1)))
c.erase(c.begin() + i);
}
return 0;
}
int lcs3(vector<int> &a, vector<int> &b, vector<int> &c) {
vector<pair<int, int>> SA;
vector<pair<int, int>> SB;
vector<pair<int, int>> SC;
for (int i = 0; i < a.size(); i++) SA.push_back(make_pair(a[i],i));
for (int i = 0; i < b.size(); i++) SB.push_back(make_pair(b[i], i));
for (int i = 0; i < c.size(); i++) SC.push_back(make_pair(c[i], i));
sort(SA.begin(),SA.end());
sort(SB.begin(), SB.end());
sort(SC.begin(), SC.end());
/*cleanData(a,b,c,SA,SB,SC);
SC.clear();
SA.clear();
SB.clear();
for (int i = 0; i < a.size(); i++) SA.push_back(make_pair(a[i], i));
for (int i = 0; i < b.size(); i++) SB.push_back(make_pair(b[i], i));
for (int i = 0; i < c.size(); i++) SC.push_back(make_pair(c[i], i));
sort(SA.begin(), SA.end());
sort(SB.begin(), SB.end());
sort(SC.begin(), SC.end());*/
//write your code here
int ptrA, ptrB, ptrC;
int currentMax = 0;
int lcs = 0;
ptrA = 0;
ptrB = 0;
ptrC = 0;
int val = 0;
for (int i = 0; i < a.size(); i++) {
val = a[i];
ptrB = get_position(val,SB,0,SB.size()-1);
ptrC = get_position(val, SC, 0, SC.size() - 1);
if ((ptrB != (-1)) && (ptrC!=(-1))) {
//currentMax = explore(a, b, c, i, ptrB, ptrC);
currentMax = explore2(a, b, c, i, ptrB, ptrC, SA, SB, SC);
lcs = max(lcs, currentMax);
}
if ((a.size() - i) < lcs)return lcs;
}
return lcs;
}
int main() {
size_t an;
std::cin >> an;
vector<int> a(an);
for (size_t i = 0; i < an; i++) {
std::cin >> a[i];
}
size_t bn;
std::cin >> bn;
vector<int> b(bn);
for (size_t i = 0; i < bn; i++) {
std::cin >> b[i];
}
size_t cn;
std::cin >> cn;
vector<int> c(cn);
for (size_t i = 0; i < cn; i++) {
std::cin >> c[i];
}
std::cout << lcs3(a, b, c) << std::endl;
return 0;
}
答案 0 :(得分：0)
你的代码对于这样的问题看起来真的很复杂，我猜这个问题不仅仅是一个实现，而是解决这个问题的整个想法。
三序列问题的LCS类似于两个序列的LCS。在这两种情况下，我们使用动态编程方法 我们有两个序列：
for i in 1 .. A.size()
for j in 1 .. B.size()
if A[i] = B[j]
lcs_len[i, j] = lcs_len[i - 1, j - 1] + 1
else
lcs_len[i, j] = max(lcs_len[i, j - 1], lcs_len[i - 1, j])
稍后，如果你想得到lcs（不仅仅是它的长度），你只需要从位置(A.size(), B.size())开始回溯，然后转向之前的值。您可以阅读更多相关信息here。
现在有三个序列问题，你只需要再添加一个循环和if语句，你应该能够想出这个！如果没有，这是一个spoiler。