【算法实验四】--【动态规划】--最长公共子序列

1041.最长公共子序列

时限：1000ms 内存限制：200000K  总时限：3000ms

描述

一个给定序列的子序列是在该序列中删去若干元素后得到的序列。确切地说，若给定序列X=<x1, x2,…, xm>，则另一序列Z=<z1, z2,…, zk>是X的子序列是指存在一个严格递增的下标序列 <i1, i2,…, ik>，使得对于所有j=1,2,…,k有：

Xij ＝ Zj

如果一个序列S即是A的子序列又是B的子序列，则称S是A、B的公共子序列。
求A、B所有公共子序列中最长的序列的长度。

 

输入

输入共两行，每行一个由字母和数字组成的字符串，代表序列A、B。A、B的长度不超过200个字符。

 

输出

一个整数，表示最长各个子序列的长度。
格式：printf("%d\n");

 

输入样例

programming
contest

 

输出样例

2

解析：这是一个比较典型的动态规划，但思路和做法却非常简单，先来说一下大概逻辑。对于长度分别为m,n的两个数组x,y，如果x[m]==y[n]的话，目前相同度flag=1，那么比较剩余的x[m-1]和y[n-1]即可，如果x[m]!=y[n]的话，就比较剩余x[m-1]，y[n]或者x[m]，y[n-1]的相似度即可。这就是基本思想。代码实现的时候是反着来的。申请一个记录相同的数组res，如果x[i]和y[j]相同的话就把res[i][j]记为1，然后在比较下面，由于序列是严格递增，所以下一个相同的必须在当前已经找到的行和列后面，具体实现见代码。需要说明的是，下面的代码还有输出相同元素的功能，如果不需要去掉即可。代码如下：

代码一：

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
int main()
{
    char a[101],b[101];
    gets(a);
    gets(b);
    int l=strlen(a),n=strlen(b);
    int m[101][101];
    memset(m,0,sizeof(m));
    for(int i=0;i<=l;i++)
        m[i][0]=0;
    for(int i=0;i<=n;i++)
        m[0][i]=0;
    for(int i=l;i>0;i--)
        a[i]=a[i-1];
    for(int i=n;i>0;i--)
        b[i]=b[i-1];
    for(int i=1;i<=l;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(a[i]==b[j])
                m[i][j]=m[i-1][j-1]+1;
            else
            {
                if(m[i-1][j]>=m[i][j-1])
                    m[i][j]=m[i-1][j];
                else
                    m[i][j]=m[i][j-1];
            }
        }
    }
    printf("%d\n",m[l][n]);
}
注：以上代码是针对这个题的，下面的代码是拓展之后可以记录路径的代码。

代码二：

#include <stdio.h>
#include <stdlib.h>
#include<string.h>
int res[100][100],dir[100][100];
void LCS(char x[],char y[],int m,int n)
{
    int i,j;
    for(i=0;i<=m;i++)
        res[i][0]=0;
    for(j=0;j<=n;j++)
        res[0][j]=0;
    for(i=1;i<=m;i++)
    {
        for(j=1;j<=n;j++)
        {
            if(x[i]==y[j])
            {
                res[i][j]=res[i-1][j-1]+1;
                dir[i][j]=1;
            }
            else if(res[i-1][j]>=res[i][j-1])
            {
                res[i][j]=res[i-1][j];
                dir[i][j]=2;
            }
            else
            {
                res[i][j]=res[i][j-1];
                dir[i][j]=3;
            }
        }
    }
}
void Print(char x[],char y[],int m,int n)
{
    int i=m,j=n;
    while(i>=0&&j>=0&&dir[i][j]!=0)
    {
        if(dir[i][j]==1)
        {
            printf("%c\n",x[i]);
            i--;j--;
        }

        else if(dir[i][j]==2)
            i--;
        else
            j--;
    }
}
int main()
{
    char x[100],y[100];
    int m,n;
    gets(x);
    gets(y);
    m=strlen(x);
    n=strlen(y);
    for(int i=m;i>=0;i--)
        x[i]=x[i-1];
    for(int i=n;i>=0;i--)
        y[i]=y[i-1];
    LCS(x,y,m,n);
    printf("%d\n",res[m][n]);
    Print(x,y,m,n);
    return 0;
}