动态规划(四)--最长公共子序列
来源:互联网 发布:coc防空火箭数据2016 编辑:程序博客网 时间:2024/05/11 00:12
最长公共子序列问题
一种相似度的概念;一个给定的序列的子序列是将序列中零个或多个元素去掉之后得到的结果。
定义:给定一个序列X=<x1,x2,。。。,xm>,另一个序列Z=<z1,z2,...,zk>满足如下条件时称为X
的子序列。即存在一个严格递增的X的下标序列<i1,i2,...,ik>,对所有j=1,2,...,k,满足xij=zj,
问题描述:给定两个序列X=<x1,x2,。。。,xm>和Y=<y1,y2,...,yn>,求X和Y的长度最长的公共子序列。
根据第一篇的理论,
Step 1:刻画最长公共子序列的特征
LCS的最优子结构
定理:
令X=<x1,x2,。。。,xm>和Y=<y1,y2,...,yn>为两个序列,Z=<z1,z2,...,zk>为X和Y的任意LCS。
1.若Xm = Yn,则zk = xm=yn且Zk-1是Xm-1和Yn-1的一个LCS。
2.若Xm不等于Yn,那么Zk不等于Xm,意味着Z是Xm-1和Y的一个LCS。
3.若Xm不等于Yn,那么Zk不等于Xm,意味着Z是Xm和Y的一个LCS。
上述定力告诉我们,两个序列的LCS包含两个序列的前缀的LCS,因此LCS问题具有最优子结构。
Step 2 : 一个递归解
在求解X=<x1,x2,。。。,xm>和Y=<y1,y2,...,yn>的一个LCS时,我们需要求解一个或两个子问题,如果Xm=Yn,我们应该求解Xm-1和Yn-1的一个LCS。将Xm=Yn追加到这个LCS的末尾,就得到X和Y的一个LCS。若Xm不等于Yn,我们必须求解两个子问题:求Xm-1和Y的一个LCS与X和Yn-1的一个LCS;两个LCS较长者即为X和Y的一个LCS。
定义c[i][j]表示Xi和Yj的LCS的长度可得以下公式
{ 0 若i=0或j=0
c[i][j] = {c[i-1][j-1] + 1 若i,j>0且Xi=Yj
{max(c[i][j-1],c[i-1][j]) 若i,j>0且xi≠Yj
Step 3:计算LCS的长度及构造LCS
自底向上方法,先生成一张表,再根据表中元素所指的方向迭代出最长子序列
/** * @author BiangHoo * * 2013年9月11日 */public class LCS {public static void main(String[] args) {String X[] = {"A","B","C","B","D","A","B"};String Y[] = {"B","D","C","A","B","A"};LCS_Length(X,Y);}static void display(String array[][]){for(int i=0;i<array.length;i++){for(int j=0;j<array[0].length;j++){System.out.print(array[i][j]+" ");}System.out.print("\n");}System.out.print("\n"+"the LCS is: ");}static void LCS_Length(String X[],String Y[]){int xlen = X.length;int ylen = Y.length;String b[][] = new String[xlen+1][ylen+1];int c[][] = new int[xlen+1][ylen+1];for(int i=0;i<xlen+1;i++){c[i][0] = 0;}for(int i=0;i<ylen+1;i++){c[0][i] = 0;}for(int i=1;i<xlen+1;i++){//bottom to topfor(int j=1;j<ylen+1;j++){if(X[i-1]==Y[j-1]){c[i][j] = c[i-1][j-1]+1;b[i][j] ="arrow";}else if(c[i-1][j] >= c[i][j-1]){c[i][j] = c[i-1][j];b[i][j] = "up";}else{c[i][j] = c[i][j-1];b[i][j] = "left";}}}display(b);Print_LCS(b,X,xlen,ylen);}static void Print_LCS(String [][] b,String[]X,int i,int j){if(b[i][j] == null){return ;}if(b[i][j] == "arrow"){Print_LCS(b,X,i-1,j-1);System.out.print(X[i-1]+" ");}else if(b[i][j] == "up"){Print_LCS(b,X,i-1,j);}else{Print_LCS(b,X,i,j-1);}}}
- 动态规划(四)--最长公共子序列
- 动态规划(四)最长公共子序列
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列(动态规划)
- 最长公共子序列&&最长公共子串---[动态规划]
- 动态规划-最长公共子序列、最长公共子串
- 动态规划之最长公共子序列
- 动态规划 ------- 最长公共子序列
- 动态规划实现最长公共子序列
- 【动态规划】最长公共子序列LCS
- 动态规划--最长公共子序列
- 通过Terracotta实现Tomcat集群(2)
- 关于cannot find -ltcl -ltk的问题
- JAVA学习笔记---类对象排序
- 构造函数与析构函数中抛出异常
- Java学习总结之聊天室项目
- 动态规划(四)--最长公共子序列
- 第93天的交易(2013-9-12)(+50)(71.08)(数据单总盈利:737.35)(贵金属总盈利:-537.92)
- Flash导致弹出的div被隐藏
- 字节流(day20)
- 外部添加的类和SDI程序View类之间的相互访问
- Attribute在.net编程中的应用
- Ubuntu 软件包管理详解
- Linux rm命令
- p115 E2