hdu 1159:Common Subsequence（动态规划）

Common Subsequence

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 18765    Accepted Submission(s): 7946

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

 

Source

Southeastern Europe 2003

 

Recommend

Ignatius   |   We have carefully selected several similar problems for you:  1087 1176 1058 1069 1421

---

 

之前做过好多次，一直不解其意，最近重温一遍。现在写下解题心得。

这道题的目的是求出a字符串和b字符串的最长公共子序列，用到动态规划。

动态规划的解法：

　　先定义两个字符数组存储两个字符串

—— char a[1000]、b[1000];

　　然后再定义一个二维数组，存储求解最终问题过程中产生的所有子问题的解

—— int dp[1001][1001];

最长公共子序列的状态转移方程为：

if(a[i]==b[j])  

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

else 

    dp[i][j]=dp[i-1][j]>dp[i][j-1]?dp[i-1][j]:dp[i][j-1];

根据以上写出程序即可。

另外摘取别人的一段对动态规划的解释：

【动态规划法】

　　经常会遇到复杂的问题不能简单的分解成几个子问题，而会分解出一系列的子问题。简单的采用把大问题分解成子问题，并综合所有子问题的解求出大问题的解的方法，问题求解耗时会按问题规模呈幂级数增加。

　　为了节约重复求相同子问题的时间，引入一个数组，不管他们是否对最终解有用，把所有子问题的解存于数组中，这就是动态规划法所采用的基本做法。

 

网易公开课的《算法导论》也有详细的讲解：

[第15集] 动态规划，最长公共子序列 

 

下面给出代码：

【C++】

 #include <iostream>

 using namespace std;
 int dp[][];
 int main()
 {
     //dp[i][j]代表着a取前i个字符和b取前j个字符时的最长公共子序列的大小
     char a[],b[];
     while(cin>>a>>b){
         int i,j;
         int al,bl;
         for(i=;a[i]!='\0';i++);    //计算a、b字符串长度
         for(j=;b[j]!='\0';j++);
         al=i;bl=j;

         for(i=;i<=al;i++)   //dp[][]初始化
             dp[i][]=;
         for(i=;i<=bl;i++)
             dp[][i]=;

         for(i=;i<=al;i++)   //计算dp[][]
             for(j=;j<=bl;j++){
                 if(a[i-]==b[j-])
                     dp[i][j]=dp[i-][j-]+;
                 else
                     dp[i][j] = dp[i-][j] > dp[i][j-] ? dp[i-][j] : dp[i][j-];
             }

         cout<<dp[al][bl]<<endl;
     }
     return ;
 }

【C】

 #include <stdio.h>
 #include <stdlib.h>
 int dp[][];
 int main()
 {
     char a[],b[];
     while(scanf("%s%s",a,b)!=EOF){
         int i,j;
         int al,bl;
         for(i=;a[i]!='\0';i++);
         for(j=;b[j]!='\0';j++);
         al=i;bl=j;
         for(i=;i<=al;i++)
             dp[i][]=;
         for(j=;j<=bl;j++)
             dp[][j]=;
         for(i=;i<=al;i++)
             for(j=;j<=bl;j++){
                 if(a[i-]==b[j-])
                     dp[i][j] = dp[i-][j-]+;
                 else
                     dp[i][j] = dp[i-][j] > dp[i][j-] ? dp[i-][j] : dp[i][j-];
             }
         printf("%d\n",dp[al][bl]);
     }
     return ;
 }

Freecode : www.cnblogs.com/yym2013