[algorithm]求最长公共子序列问题
最直白方法:时间复杂度是O(n3), 空间复杂度是常数
reference:http://blog.csdn.net/monkeyandy/article/details/7957263
/**
** copyright@andy
** http://blog.csdn.net/MonkeyAndy
**/
首先介绍动态规划方法的相关知识
动态规划方法的基本思想:
分成若干个子问题,先求解子问题,然后根据子问题的解求得原问题的解。经分解得到的子问题往往不是互相独立的。可重复利用!
其核心思想就是分治,分治方法的特点是分解后的子问题相对独立,可以通过简单的合并算法得到原问题的解!
动态规划方法的应用对象:
优化问题:
- 一个优化问题可能有很多可行解,每个解都有一个求解代价,
- 我们希望选择一个具有最少代价的解。
- 一个优化问题可能有多个优化解。
动态规划方法:
当一个优化问题可分为多个子问题,子问题的解在构造上一级问题的求解过程中被重复使用。 这样可以节省计算时间与空间。
动态规划算法的步骤:
- 分析优化解的结构
- 递归地定义最优解的代价
- 自底向上地计算优化解的代价保存之,并获取构造最优解的信息。
- 根据计算优化解的代价,构造优化解
2. 具有重叠子问题(Overlapping sub-problems)
动态规划算法利用该性质,以求得的子问题的解保存后,自底向上地,重复利用子问题的解获得上级问题的解,节省空间开销和时间开销。
如果不具备重叠子问题的性质动态规划方法的空间开销和时间开销很大。利用动态规划算法没有意义!
设x=(x1, x2, ... xm)是一个序列, z=(z1, z2, ... zm) 是x的一个子序列,如果存在1i1<i2<...<ik<m, 使xij=zj
z=(z1, z2, ... zm)即是x=(x1, x2, ... xm)中删除一些元素以后得到的序列。
2.公共子序列定义
Z是序列X与Y的公共子序列,如果Z是X的子序列且Z是Y的子序列。
3.最长公共子序列问题(LCS)
输入:X=(x1, x2, ... xn),Y=(y1, y2, ...ym)
输出:X与Y的最长公共子序列Z
定义. 设X=(x1 , x2 , ... xn) 是一个序列, Xi表示X的第i个前缀,即 Xi=(x1 , ... xi )。
例. X=(A,B,D,C,A), X1=(A), X2=(A,B), X3=(A,B,D)
定理1(LCS的优化结构)设X=(x1 , ... xm), Y=(y1 , y2 , ... yn)是两个序列,Z=(z1 , z2 , ... zk) 是X与Y的LCS。下列结论成立:
- 如xm= yn, 则zk=xm=yn, Zk-1是Xm-1和Yn-1的LCS,即,LCS(X,Y)=LCS(Xm-1,Yn-1)+xm
- 若xm != yn,且zk!=xm,则Z是Xm-1和Y的LCS,即 LCS(X,Y)=LCS(Xm-1,Y)
- 若xm != yn , 且zk !=yn,则Z是X与Yn-1的LCS,即 LCS(X,Y)=LCS(X,Yn-1)
If xm !=yn , 必须求解Xm-1和Y的LCS以及X和Yn-1 的LCS,长者是X和Y的LCS。C[i,j] 表示Xi与Yj的LCS的长度,则建立求解LCS长度的递归方程
1.基本思想
计算C[i,j] 需先计算C[i-1,j-1]、C[i,j-1]、C[i-1,j]
C[i-1,j-1] | C[ i, j-1] |
C[i-1, j] | C[ i, j] |
先计算出第0行与第0列,然后逐行计算。
2. 算法
- /**
- *copyright@andy
- *date@ 2012/09/08
- * m : x的长度
- * n : y的长度
- * c : xi与yi的LCS的长度
- * b : 用于标记xi和yi的关系
- **/
- void LCSLength( int m, int n, char *x, char *y int **c, Type **b )
- {
- int i, j;
- /*c 初始化*/
- for ( i=1; i<=m; i++) c[ 0 ][ i ]=0;
- for( j=1; j<=n, j++) c[ j ][0 ]=0;
- for ( i=1; i<=m; i++)
- for( j=1; j<=n, j++)
- {
- if(x[i]==y[j]) //xi = yj 时LCS 增1
- {
- c[ i ][ j ]=c[ i-1 ][ j-1]+1;
- b[ i ][ j ]='\'';
- }
- else if ( c[i-1][j] >= c[i][j-1]) //当xi ≠ yj时保存长子序列
- {
- c[ i ][ j ]=c[ i-1 ][ j];
- b[ i ][ j ]='|';
- }
- else
- {
- c[ i ][ j ]=c[ i ][ j-1];
- b[ i ][ j ]='--';
- }
- }
- }
- /**
- *copyright@andy
- *date@ 2012/09/08
- * i : x的长度
- * j : y的长度
- * x : x字符串
- * b : 用于标记xi和yi的关系
- **/
- void LCS( int i, int j, char *x, Type **b )
- {
- if ( i==0 || j==0) return;
- if (b[ i ][ j ] =='\'')//当xi = yj 时,输出xi ,并删除xi 和 yj ,之后在子序列中继续求解LCS
- {
- LCS(i-1, j-1, x, b);
- out<<x[ i ];
- }
- else if ( b[i][j] == '|')//当Xi-1 和 Yj 的LCS不小于Xi 和 Yj-1 的LCS 时
- LCS(i-1, j, x, b); //在Xi-1 和 Yj中继续求解LCS;否则,在Xi 和 Yj-1 中继续求解LCS;
- else
- LCS( i, j-1, x, b);
- }
[algorithm]求最长公共子序列问题的更多相关文章
- HDU 4681 string 求最长公共子序列的简单DP+暴力枚举
先预处理,用求最长公共子序列的DP顺着处理一遍,再逆着处理一遍. 再预处理串a和b中包含串c的子序列,当然,为了使这子序列尽可能短,会以c 串的第一个字符开始 ,c 串的最后一个字符结束 将这些起始位 ...
- HDU 1243 反恐训练营 (动态规划求最长公共子序列)
反恐训练营 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total Subm ...
- Java实现 LeetCode 583 两个字符串的删除操作(求最长公共子序列问题)
583. 两个字符串的删除操作 给定两个单词 word1 和 word2,找到使得 word1 和 word2 相同所需的最小步数,每步可以删除任意一个字符串中的一个字符. 示例: 输入: " ...
- poj1458 求最长公共子序列 经典DP
Common Subsequence Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 45763 Accepted: 18 ...
- codevs 1862 最长公共子序列(求最长公共子序列长度并统计最长公共子序列的个数)
题目描述 Description 字符序列的子序列是指从给定字符序列中随意地(不一定连续)去掉若干个字符(可能一个也不去掉)后所形成的字符序列.令给定的字符序列X=“x0,x1,…,xm-1”,序列Y ...
- 动态规划求最长公共子序列(Longest Common Subsequence, LCS)
1. 问题描述 子串应该比较好理解,至于什么是子序列,这里给出一个例子:有两个母串 cnblogs belong 比如序列bo, bg, lg在母串cnblogs与belong中都出现过并且出现顺序与 ...
- 【dp】求最长公共子序列
[题目描述] 一个给定序列的子序列是在该序列中删去若干元素后得到的序列.确切地说,若给定序列X=<x1,x2,…,xm>X=<x1,x2,…,xm>,则另一序列Z=<z1 ...
- Coincidence (动态规划求最长公共子序列)(王道)
题目描述: Find a longest common subsequence of two strings. 输入: First and second line of each input case ...
- hdu 1159求最长公共子序列
题目描述:给出两个字符串,求两个字符串的公共子序列(不是公共子串,不要求连续,但要符合在原字符串中的顺序) in: abcfbc abfcab programming contest abcd mnp ...
随机推荐
- ContentProvider使用
package activity.cyq.contentprovidertenderer; import android.content.ContentProvider; import android ...
- part 5 Two way databinding in AngularJS
- Hive通过查询语句向表中插入数据注意事项
最近在学习使用Hive(版本0.13.1)的过程中,发现了一些坑,它们或许是Hive提倡的比关系数据库更加自由的体现(同时引来一些问题),或许是一些bug.总而言之,这些都需要使用Hive的开发人员额 ...
- ASP.NET 跨域获取JSON天气数据
前几天做一个门户网站,在首页需要加载气象数据,采用了中央气象局的接口. 刚开始采用JSONP在前台跨域请求数据,没成功~ 后换成在c#后台请求数据返回... 前端代码: $(function () { ...
- iOS定位 - 普通定位(没有地图) - 反地理编码(得到具体位置)
#import <CoreLocation/CoreLocation.h> 使用到的头文件 要引入CoreLocation这个包 <CLLocationManagerDelegate ...
- UI1_UIScrollView
// // AppDelegate.m // UI1_UIScrollView // // Created by zhangxueming on 15/7/10. // Copyright (c) 2 ...
- 工作单元(Unit of Work)
维护受业务事务影响的对象列表,并协调变化的写入和并发问题的解决. 从DB中存取Data时,必须记录增删改动作,以将对DB有影响的数据写会到DB中去. 如果在每次修改对象模型时就对DB进行相应的修改,会 ...
- JAVA多线程通信
JAVA多线程通信 package com.frank.thread; /** * author:pengyan * date:Jun 16, 2011 * file:ProducerAndCusto ...
- 信号处理程序(signal handler)会被重置的信号
首先说明我的系统,CentOS 6.6,内核为2.6.32-504.12.2.el6.i686. 当用signal对某个信号设定信号处理函数的时候,有些信号的处理函数会被重置,有些则不会,这种情况的具 ...
- php var_export与var_dump 输出的不同
var_export必须返回合法的php代码,var_export返回的代码,可以直接当作php代码赋值个一个变量. 而这个变量就会取得和被var_export一样的类型的值. 问题描述: 在跟踪 ...