算法导论(LCS最长公共子序列)

真动态规划版(自底而下):

#include<iostream>
using namespace std;

const int maxn = 10001;
int m,n;

char b[maxn][maxn]; //回溯时保存路径
int c[maxn][maxn]; //最长子序列
void LCS(char *x,char *y)
{
    int i,j;
    for(i=0;i<=m;i++)
    {
        c[i][0] = 0;
    }
    for(j=0;j<=n;j++)
    {
        c[0][j] = 0;
    }
    for(i=1;i<=m;i++)
    {
        for(j=1;j<=n;j++)
        {
            if(x[i] == y[j])
            {
                c[i][j] = c[i-1][j-1] + 1;
            }
            else if(c[i-1][j] >= c[i][j-1])
            {
                c[i][j] = c[i-1][j];
            }
            else
            {
                c[i][j] = c[i][j-1];
            }
        }
    }
}

int main()
{
    int i;
    char x[maxn];
    char y[maxn];
    while(cin>>m>>n)
    {
        for(i=1;i<=m;i++)
        {
            cin>>x[i];
        }
        for(i=1;i<=n;i++)
        {
            cin>>y[i];
        }
        LCS(x,y);
        cout<<c[m][n]<<endl;
    }
    return 0;
}

递归版带备忘(自顶而下):

#include<iostream>
#include<cstring>
using namespace std;

const int maxn = 10001;
char b[maxn][maxn]; //回溯时保存路径
int c[maxn][maxn]; //最长子序列
int LCS(char *x,char *y,int m,int n)
{
    int i,j;
    if(m == 0 || n == 0)  //两个串有一个序列为0时那么最长公共序列一定为0
    {
        c[m][n] = 0;
    }
    else if(c[m][n] > 0)
    {
        return c[m][n]; //这里是重点(采用了备忘的方法)，用空间换取时间，因为树里面重复的部分可能很多
    }
    else
    {
        if(x[m] == y[n])  //第一种情况
        {
            c[m][n] = LCS(x,y,m-1,n-1) + 1;
        }
        else //第二种情况
        {
            c[m][n] = max(LCS(x,y,m-1,n),LCS(x,y,m,n-1));
        }
    }
    return c[m][n]; //递归完成后返回值
}

int main()
{
    int i,j;
    char x[maxn];
    char y[maxn];
    int m,n;
    while(cin>>m>>n)
    {
        for(i=1; i<=m; i++)
        {
            cin>>x[i];
        }
        for(i=1; i<=n; i++)
        {
            cin>>y[i];
        }
        for(i=0; i<=m; i++)
        {
            c[i][0] = 0;
        }
        for(j=0; j<=n; j++)
        {
            c[0][j] = 0;
        }
        cout<<LCS(x,y,m,n)<<endl;
    }
    return 0;
}