DHU：1159 Common Subsequence

Common Subsequence

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)

Total Submission(s): 33697    Accepted Submission(s): 15333

Problem Description

A subsequence of a given sequence is the given sequence with some elements (possible none) left out. Given a sequence X = <x1, x2, ..., xm> another sequence Z = <z1, z2, ..., zk> is a subsequence of X if there exists a strictly increasing sequence <i1, i2, ..., ik> of indices of X such that for all j = 1,2,...,k, xij = zj. For example, Z = <a, b, f, c> is a subsequence of X = <a, b, c, f, b, c> with index sequence <1, 2, 4, 6>. Given two sequences X and Y the problem is to find the length of the maximum-length common subsequence of X and Y. 

The program input is from a text file. Each data set in the file contains two strings representing the given sequences. The sequences are separated by any number of white spaces. The input data are correct. For each set of data the program prints on the standard output the length of the maximum-length common subsequence from the beginning of a separate line. 

Sample Input

abcfbc abfcab

programming contest 

abcd mnp

Sample Output

4

2

0

解题思路：用动态规划求最长公共子序列。

3、1 

问题描述•

若给定序列X={x1,x2,…,xm}，则另一序列Z={z1,z2,…,zk}，是X的子序列是指存在一个严格递增下标序列{i1,i2,…,ik}使得对于所有j=1,2,…,k有：zj=xij。例如，序列Z={B，C，D，B}是序列X={A，B，C，B，D，A，B}的子序列，相应的递增下标序列为

{2，3，5，7}。

•公共子序列的定义：

给定2个序列X和Y，当另一序列Z既是X的子序列又是Y的子序列时，称Z是序列X和Y的公共子序列。

3、2 最长公共子序列的结构

设序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}的最长公共子序列为Z={z1,z2,…,zk}，

(1)若xm=yn，则zk=xm=yn，且zk-1是X   m -1和Y n - 1 的最长公共子序列。 

(2)若xm≠yn且zk≠xm,则z是xm-1和Y的最长公共子序列。

(3)若xm≠yn且zk≠yn，                                                   则 z是X和yn-1的最长公共子序列。

由此可见，2个序列的最长公共子序列包含了这2个序列的前缀的最长公共子序列。因

此，最长公共子序列问题具有最优子结构性质。

3、3 子问题的递归结构  由最长公共子序列问题的最优子结构性质建立子问题最优值的递归关系。用c[i][j]记录序列和的最长公共子序列的长度。其中， Xi={x1,x2,…,xi}；Yj={y1,y2,…,yj}。当i=0或j=0时，空序列是Xi和Yj的最长公共子序列。故此时C[i][j]=0。其他情况下，由最优子结构性质可建立递归关系如下：

当i=0,j=0;则c[i][j]=0;

当I,j>0,xi==yj,c[i][j]=c[i-1][j-1]+1;

当I,j>0,xi!=yj,c[i][j]=max(c[i-1][j],c[i][j-1]);

题目代码：

#include<iostream>

#include <cstdio>

#include <cstring>

//#include <queue>

//#include <stack>

//#include <algorithm>

#define max(a,b)a>b?a:b

using namespace std;

int main()

{

    char s1[1500];

    char s2[1500];

    int dp[1500][1500];

    int i,j;

    while(scanf("%s%s",s1,s2)!=EOF)

    {

        memset(dp,0,sizeof(dp));

        int len1=strlen(s1);

        int len2=strlen(s2);

        for(i=1;i<=len1;i++)

            for(j=1;j<=len2;j++)

            {

                if(s1[i-1]==s2[j-1])

                    dp[i][j]=dp[i-1][j-1]+1;

                else

                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);

            }

        printf("%d\n",dp[len1][len2]);

    }

    return 0;

}