LCS(最长公共子序列)及其O(n)空间优化，O(nlogn)时间复杂度优化

n^2 的版本

int LCS_n_2(vector<char> &a,vector<char> &b)
{vector<vector<int>> dp(a.size(),vector<int>(b.size(), 0));for(int i = 0;i<a.size();++i) if (a[i] == b[0]) dp[i][0] = 1;for(int j = 0;j<b.size();++j) if (a[0] == b[j]) dp[0][j] = 1;for(int i = 1;i<a.size();++i)for(int j = 1;j<b.size();++j){if(a[i] == b[j]) dp[i][j] = dp[i-1][j-1]+1;else dp[i][j] = max(dp[i-1][j],dp[i][j-1]);}return dp[a.size()-1][b.size()-1];
}

此时注意到其实每次算dp只需要三个数，dp[i][j],dp[i-1][j],dp[i][j-1]，那么只要保存一行的dp，后续直接更新即可，于是空间复杂度可以优化为 O ( N ) O(N) O(N)

template<typename vec>
int LCS_n(vec &a,vec &b)
{vector<int> dp(b.size()+1,0);int last = 0;for(int i = 1;i<=a.size();++i){last = 0;for(int j = 1;j<=b.size();++j){if(a[i-1] == b[j-1]){int tmp = last;last = dp[j-1] + 1;if(j-1>0)dp[j-1] = tmp;}else{int tmp = last;last = max(last,dp[j]);if(j-1>0)dp[j-1] = tmp;}//show(dp);}dp[b.size()] = last;}return last;
}

那么其实还有一种可以把时间复杂度优化到 O ( n l o g n ) O(nlogn) O(nlogn)的方法。

总的来说，就是把LCS转化成LIS，然后用LIS的NlogN算法来求解。
实现如下：（引用）
假设有两个序列 s1[ 1~6 ] = { a, b, c , a, d, c }, s2[ 1~7 ] = { c, a, b, e, d, a, b }。
记录s1中每个元素在s2中出现的位置, 再将位置按降序排列, 则上面的例子可表示为：loc( a）= { 6, 2 }, loc( b ) = { 7, 3 }, loc( c ) = { 1 }, loc( d ) = { 5 }。
将s1中每个元素的位置按s1中元素的顺序排列成一个序列s3 = { 6, 2, 7, 3, 1, 6, 2, 5, 1 }。
在对s3求LIS得到的值即为求LCS的答案。

可以参考下列实现

int LCS_n_logn(vector<char> &a,vector<char> &b)
{unordered_map<char,vector<int> >dp;for(int i = 0;i<a.size();++i){if(dp.count(a[i]) == 0)dp[a[i]].push_back(-1);}for(int i = 0;i<b.size();++i){if(dp.count(b[i]) == 1)dp[b[i]].push_back(i+1);}vector<int> dp2;for(int i = 0;i<a.size();++i){if(dp[a[i]].size() > 1){for(int j = dp[a[i]].size() - 1; j > 0; --j)dp2.push_back(dp[a[i]][j]);}}return LIS(dp2);
}

LCS(最长公共子序列)及其O(n)空间优化，O(nlogn)时间复杂度优化相关推荐

LCS最长公共子序列和LIS最长上升子序列——例题剖析
一.LCS最长公共子序列最长公共子序列(LCS)问题算法详解+例题(转换成LIS,优化为O(nlogn),看不懂你来打我) longest comment subsequence 模板题 longe ...
算法设计 - LCS 最长公共子序列最长公共子串 LIS 最长递增子序列
出处 http://segmentfault.com/blog/exploring/ 本章讲解: 1. LCS(最长公共子序列)O(n^2)的时间复杂度,O(n^2)的空间复杂度: 2. 与之类似但不 ...
动态规划求解LCS最长公共子序列问题c++
文章目录求最长公共子序列最长公共子序列问题刻画最长公共子序列的特征递归公式 C++代码求最长公共子序列长度打印任意一组最长公共子序列模板题求最长公共子序列最长公共子序列问题 (Lon ...
LCS(最长公共子序列)
题意描述求两个字符串的最长公共子序列的长度动态规划用二维数组 C[i][j] 记录串x1x2⋯xi与y1y2⋯yj 的 LCS长度,则可得到状态转移方程: 代码实现: #include<i ...
LCS/最长公共子序列/最长公共子串实现 Python/Java
参考 http://blog.csdn.net/u012102306/article/details/53184446 http://blog.csdn.net/hrn1216/article/det ...
大白话讲解LCS(最长公共子序列)
今天看了七月在线算法课.再一次认识了LCS,现在整理记录: LCS(Longest Common Subsequence)最长公共子序列. 一个序列S任意删除若干个字符得到新序列T,那么T叫做S的子序 ...
LCS最长公共子序列（最优线性时间O(n)）
这篇日志主要为了记录这几天的学习成果. 最长公共子序列根据要不要求子序列连续分两种情况. 只考虑两个串的情况,假设两个串长度均为n. 一,子序列不要求连续. (1)动态规划(O(n*n)) (转自:h ...
POJ 1458 Common Subsequence DP LCS 最长公共子序列
最长公共子序列,照抄<算法设计与分析导论>P138-140 设输入的两个字符串分别为a1,a2,```,am(串a) b1,b2,````,bn(串b) 设d(i,j)为字符串a1,a2, ...
LCS最长公共子序列
例如 b c d d e和 a c e e d e的公共子串为c d e. 如果使用暴力,复杂度太高会直接超时.就需要使用动态规划 dp[i][j]表示a串第i个结尾,b串第j个结尾的最长公共子串的数 ...

LCS(最长公共子序列)及其O(n)空间优化，O(nlogn)时间复杂度优化

LCS(最长公共子序列)及其O(n)空间优化，O(nlogn)时间复杂度优化

LCS(最长公共子序列)及其O(n)空间优化，O(nlogn)时间复杂度优化相关推荐

最新文章

热门文章