求x=y=的最長公共子串行.
蠻幹法:
找出x的所有子串行,共2^m個(m是x的元素個數),對每乙個自序列檢視是否在y中,需要n次(n是y的元素個數)
總的時間複雜度: o(n*2^m)
動態規劃演算法:
xi=yj=
設 zk=是xi和yi的最長公共子串行.
則:(1)當xi=yj時,則zk=xi=yj,且zk-1是xi-1與yj-1的最長公共子串行.
(2)若xi != yj,zk != xi, 則zk是xi-1 與yj的最長公共子串行
(3)若xi != yj,zk != yj, 則zk是xi 與yj-1的最長公共子串行
遞推關係:
使用c[i][j]表示xi,yj的最長公共子串行的長度:
c[i][j] =c[i-1][j-1] + 1 如果 i,j>0,xi=yj
=max 如果 i,j>0,xi != yj
c[0][j]=c[i-1][0] 對一切1<=i<=m,1<=j<=n
下面是實現:
package com.動態規劃;
public class lcs
int a_length = a.length;
int b_length = b.length;
if(a_length==0||b_length==0)
this.a = new char[a_length+1];
this.b = new char[b_length+1];
//為了後面計算方便,將兩個序列右移一位
system.arraycopy(a, 0, this.a, 1, a.length);
system.arraycopy(b, 0, this.b, 1, b.length);
//關鍵**
length = new int[a_length+1][b_length+1];
trace = new int[a_length+1][b_length+1];
//初始化
for(int i=0;i=length[i][j-1])else}}
}public void printlcs()
private void printlcs(int m, int n)
//表示a[m]==b[n],輸出
if(trace[m][n]==1)else if(trace[m][n]==2)else
}public static void main(string args) ;
char b = ;
//char a = ;
//char b = ;
lcs lcs = new lcs(a,b);
lcs.printlcs();
}}
結果:
b c b a
最長公共子串行 最長公共子串
1 最長公共子串行 採用動態規劃的思想,用乙個陣列dp i j 記錄a字串中i 1位置到b字串中j 1位置的最長公共子串行,若a i 1 b j 1 那麼dp i j dp i 1 j 1 1,若不相同,那麼dp i j 就是dp i 1 j 和dp i j 1 中的較大者。class lcs el...
最長公共子串行 最長公共子串
1.區別 找兩個字串的最長公共子串,這個子串要求在原字串中是連續的。而最長公共子串行則並不要求連續。2 最長公共子串 其實這是乙個序貫決策問題,可以用動態規劃來求解。我們採用乙個二維矩陣來記錄中間的結果。這個二維矩陣怎麼構造呢?直接舉個例子吧 bab 和 caba 當然我們現在一眼就可以看出來最長公...
最長公共子串 最長公共子串行
子串要求連續 子串行不要求連續 之前的做法是dp求子序列 include include include using namespace std const int inf 0x3f3f3f3f const int mod 1000000007 string s1,s2 int dp 1010 10...