最长公共子序列打印不正确

时间:2013-11-12 21:24:15

标签: c++ string algorithm

我一直在尝试实施最长公共子序列几个小时。我已经检查过LCSLength函数返回正确的长度,但它没有正确打印序列。

int max(int a , int b)
{
    if(a>b) 
        return a;
    else 
        return b;
}
/* The printing function that is not functioning Properly */
string backtrack(vector< vector<int> > C, string X,string Y,int i,int j)
{
    if( i==0  || j==0)
        return "";
    else if( X[i] == Y[j])
        return backtrack(C,X,Y,i-1,j-1) + X[i];
    else
        if( C[i][j-1] >= C[i-1][j] )
            return backtrack(C,X,Y,i,j-1);
        else
            return backtrack(C,X,Y,i-1,j);

}



   /* It correctly returns the subsequence length. */
   int LCSLength(string s1,string s2)
   {   
        vector< vector <int> >  C (s1.size()+1 , vector<int> (s2.size() +1 ) );
        int i,j;
        int m = s1.size();
        int n = s2.size();
        for(i =0; i<=m; i++){
        for( j =0; j<=n; j++){
            if( i==0 || j==0)
                C[i][j] = 0;
            else if( s1[i-1] == s2[j-1] )
                C[i][j] = C[i-1][j-1] +1;
            else
                C[i][j] = max(C[i][j-1],C[i-1][j]);
        }
    }




    cout << backtrack(C,s1,s2,m,n);


    return C[m][n];
}

我遵循此处给出的伪代码:http://en.wikipedia.org/wiki/Longest_common_subsequence_problem

你会感激的任何帮助。

经测试的案例:

 cout << LCSLength("HUMAN", "CHIMPANZEE") << endl;

返回4但字符串序列不正确。

  cout << LCSLength("AAACCGTGAGTTATTCGTTCTAGAA","CACCCCTAAGGTACCTTTGGTTC") << endl;

返回14

感谢。

1 个答案:

答案 0 :(得分:2)

我猜你的停止条件是:

if( i==0  || j==0)
     return "";

您永远无法到达X [0]或Y [0]处的角色。您打印输出错误的示例(“MAN”而不是“HMAN”)与此匹配,因为H是第一个字符。

请注意,您链接的wiki值将字符串定义为X[1..m] and Y[1..n],这可能不是在c ++中执行此操作的直观方式。

尝试切换为-1作为停止条件,或者在开头填充字符串。