最长公共子序列(longest common subsequence) 动态规划

问题

题目：[Common Subsequence]

思路

1. 状态定义 ：$dp[i][j]$表示$X_{i}$,$Y_{j}$的最长公共子序列长度。
2. 转移函数
   $$dp[i][j] = \left\{ \begin{aligned} 0, && i = 0, j = 0 \\ dp[i-1][j-1] + 1, && x_{i} = y_{j} \\ max(dp[i-1][j], dp[i][j-1]), && x_{i} \ne y_{j} \end{aligned} \right. \tag{1}$$

代码

#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
int cal_lcs( const std::string& X, const std::string Y );

int main()
{
    std::string s1;
    std::string s2;
    while( std::cin >> s1 >> s2 )
    {
        int ret = cal_lcs( s1, s2 );
        std::cout << ret << std::endl;
    }
    return 0;
}

int cal_lcs( const std::string& X, const std::string Y )
{
    int lengthX = X.size();
    int lengthY = Y.size();

    std::vector< std::vector<int> > dp( lengthX + 1, std::vector<int>( lengthY + 1, int() ) );

    int ans = 0;
    for( int i = 0; i < lengthX + 1; ++i )
        dp[i][0] = 0;
    for( int i = 0; i < lengthY + 1; ++i )
        dp[0][i] = 0;


    for( int i = 1; i < lengthX + 1; ++i )
    {
        for( int j = 1; j < lengthY + 1; ++j )
        {
            if( X[i-1] == Y[j-1] )
                dp[i][j] = dp[i-1][j-1] + 1;
            else
                dp[i][j] = std::max( dp[i-1][j], dp[i][j-1] );

            ans = std::max( ans, dp[i][j] );
        }
    }

    return ans;
}

对于内存开辟使用指针的形式去控制，给出如下代码实现

代码1

int cal_lcs( const std::string& s1, const std::string& s2 )
{
    int sz1 = s1.size();
    int sz2 = s2.size();
    if( !sz1 || !sz2 )
        return 0;

    int** dp = new int*[sz1 + 1];
    if( !dp )
        return -1;
    for( int i = 0; i < sz1 + 1; ++i )
    {
        dp[i] = new int[sz2 + 1];
        if( !dp[i] )
            return -1;
    }

    for( int i = 0; i < sz1 + 1; ++i )
    {
        dp[i][0] = 0;
    }
    for( int i = 0; i < sz2 + 1; ++i )
    {
        dp[0][i] = 0;
    }

    for( int i = 1; i < sz1 + 1; ++i )
    {
        for( int j = 1; j < sz2 + 1; ++j )
        {
            dp[i][j] = (s1[i-1] == s2[j-1])?dp[i-1][j-1] + 1:std::max( dp[i-1][j], dp[i][j-1] );
        }
    }
    int ans = dp[sz1][sz2];
    for( int i = 0; i < sz1 + 1; ++i )
    {
        delete [] dp[i];
    }
    delete [] dp;

    return ans;
}

其他

下面这种二维vector的用法，我是第一次看到。之前我都是用指针开辟二维数组实现的。毕竟在我看ACC的时候，笔者也是建议既然c++有这么强大的STL，解决问题的时候还是首选STL。对于指针一般是用在接近底层的内存管理当中使用，所以以后c里面对于二维数组的使用，放到c++里面都可以用这种形式的vector去实现。

std::vector< std::vector<int> > dp( lengthX + 1, std::vector<int>( lengthY + 1, int() ) );