C语言 · 最长公共子序列 · 最长字符序列
如果存在一个严格递增而且长度大于0的下标序列{i1,i2……ik},使得对所有的j=1,2,……k,有x(ij)=z(j),那么我们称Z是X的字符子序列。而且,如果Z既是X的字符子序列又是Y的字符子序列,那么我们称Z为X和Y的公共字符序列。
在我们今天的问题中,我们希望计算两个给定字符序列X和Y的最大长度的公共字符序列,这里我们只要求输出这个最大长度公共子序列对应的长度值。
举例来说,字符序列X=abcd,Y=acde,那么它们的最大长度为3,相应的公共字符序列为acd。
aedfhb

/*
递归思路:
当数组a和b对应位置字符相同时,则直接求解下一个位置;
否则,取两种情况中的较大数值。
*/
#include<stdio.h>
#include<string.h>
char a[],b[];//定义字符串数组
int lena,lenb,lenz=;
int maxlong(int i,int j){
if(i>=lena || j>=lenb) return ;//出口
if(a[i] == b[j]) return +maxlong(i+,j+);//直接求解下一个位置
else //取两种情况中的较大数值
return maxlong(i+,j) > maxlong(i,j+) ? maxlong(i+,j) : maxlong(i,j+);
}
int main(){
scanf("%s%s",&a,&b);//输入字符串
lena=strlen(a);//获取字符串长度
lenb=strlen(b);
printf("%d",maxlong(,));//输出子串长度
return ;
}
动态规划:时间复杂度大大降低。
#include<stdio.h>
#include<string.h>
char a[],b[];
char num[][];//记录中间结果的数组
//动态规划采用二维数组来标识中间计算结果,避免重复的计算来提高效率
void LCS(int lena,int lenb){
for(int i=;i<=lena;i++){
for(int j=;j<=lenb;j++){
if(a[i-]==b[j-]){//注意这里的下标是i-1与j-1
num[i][j]=num[i-][j-]+;
}
else{
num[i][j]=num[i][j-] > num[i-][j] ? num[i][j-] : num[i-][j];
}
}
}
}
int main(){
scanf("%s%s",&a,&b);//输入字符串
int lena = strlen(a);//获取字符串长度
int lenb = strlen(b);
memset(num,,sizeof(num));//数组赋初值
LCS(lena,lenb);
printf("%d",num[lena][lenb]);
return ;
}
求两字符串的最长子串,有两种情况:一种是不要求在原串中连续(即上述情况),另一种要求子串在原字符串中是连续的。
下面转载了网上的方法来做另一种情况的。
用矩阵来记录中间结果:例如原串"bab"和"caba",则数组如下:
b a b
c 0 0 0
a 0 1 0
b 1 0 1
a 0 1 0
则矩阵的斜对角线最长的那个就是我们找的最长公共子串——求最长的由1组成的斜对角线:当要在矩阵是填1时让它等于其左上角元素加1,如下:
b a b
c 0 0 0
a 0 1 0
b 1 0 2
a 0 2 0
这样矩阵中的最大元素就是 最长公共子串的长度。
在构造这个二维矩阵的过程中由于得出矩阵的某一行后其上一行就没用了,所以实际上在程序中可以用一维数组来代替这个矩阵(降低空间复杂度)。
以下代码来自网络:
1 #include<iostream>
2 #include<cstring>
3 #include<vector>
4 using namespace std;
5 //str1为横向,str2这纵向
6 const string LCS(const string& str1,const string& str2){
7 int xlen=str1.size(); //横向长度
8 vector<int> tmp(xlen); //保存矩阵的上一行
9 vector<int> arr(tmp); //当前行
10 int ylen=str2.size(); //纵向长度
11 int maxele=0; //矩阵元素中的最大值
12 int pos=0; //矩阵元素最大值出现在第几列
13 for(int i=0;i<ylen;i++){
14 string s=str2.substr(i,1);
15 arr.assign(xlen,0); //数组清0
16 for(int j=0;j<xlen;j++){
17 if(str1.compare(j,1,s)==0){
18 if(j==0)
19 arr[j]=1;
20 else
21 arr[j]=tmp[j-1]+1;
22 if(arr[j]>maxele){
23 maxele=arr[j];
24 pos=j;
25 }
26 }
27 }
28 tmp.assign(arr.begin(),arr.end());
29 }
30 string res=str1.substr(pos-maxele+1,maxele);
31 return res;
32 }
33 int main(){
34 string str1("21232523311324");
35 string str2("312123223445");
36 string lcs=LCS(str1,str2);
37 cout<<lcs<<endl;
38 return 0;
39 }
最长公共子序列应用
有这样的题目:
回文串是指正反都相同的字符串,比如abba,他正着看反着看都相同,给你一个字符串,长度小于1000,让你删除最少的字符使得该字符串成为回文串,例子:agddgca,删除字符c就可以组成一个回文串。问最少删除几个字符?
想到可以借助最长公共子序列来解决,我们那上面那个例子看一下,我们把这个字符转倒序:acgddga,可以看到他们的最长公共子序列刚好是最长回文串。其实倒过来的字符串和原来的字符串比较,就是那前面的跟后面的字符串比较,如果一样,那么他们就是回文了。
所以求解这道题就比较简单了:
1)先将字符串倒序排列
2)求两个字符串飞串的最长公共子序列
3)用字符串的长度减去最长公共子序列就得到结果
C语言 · 最长公共子序列 · 最长字符序列的更多相关文章
- 算法设计 - LCS 最长公共子序列&&最长公共子串 &&LIS 最长递增子序列
出处 http://segmentfault.com/blog/exploring/ 本章讲解:1. LCS(最长公共子序列)O(n^2)的时间复杂度,O(n^2)的空间复杂度:2. 与之类似但不同的 ...
- 《算法导论》读书笔记之动态规划—最长公共子序列 & 最长公共子串(LCS)
From:http://my.oschina.net/leejun2005/blog/117167 1.先科普下最长公共子序列 & 最长公共子串的区别: 找两个字符串的最长公共子串,这个子串要 ...
- 最长公共子序列&最长公共子串
首先区别最长公共子串和最长公共子序列 LCS(计算机科学算法:最长公共子序列)_百度百科 最长公共子串,这个子串要求在原字符串中是连续的.而最长公共子序列则并不要求连续. 最长公共子序列: http ...
- 动态规划——最长公共子序列&&最长公共子串
最长公共子序列(LCS)是一类典型的动归问题. 问题 给定两个序列(整数序列或者字符串)A和B,序列的子序列定义为从序列中按照索引单调增加的顺序取出若干个元素得到的新的序列,比如从序列A中取出 A ...
- 简单动态规划——最长公共子序列&&最长回文子序列&&最长上升||下降子序列
最长公共子序列,顾名思义当然是求两个字符串的最长公共子序列啦,当然,这只是一道非常菜的动规,所以直接附上代码: #include<iostream> #include<cstdio& ...
- [BZOJ2423][HAOI2010]最长公共子序列
[BZOJ2423][HAOI2010]最长公共子序列 试题描述 字符序列的子序列是指从给定字符序列中随意地(不一定连续)去掉若干个字符(可能一个也不去掉)后所形成的字符序列.令给定的字符序列X=“x ...
- [Python]最长公共子序列 VS 最长公共子串[动态规划]
前言 由于原微软开源的基于古老的perl语言的Rouge依赖环境实在难以搭建,遂跟着Rouge论文的描述自行实现. Rouge存在N.L.S.W.SU等几大子评估指标.在复现Rouge-L的函数时,便 ...
- 动态规划之最长公共子序列(LCS)
转自:http://segmentfault.com/blog/exploring/ LCS 问题描述 定义: 一个数列 S,如果分别是两个或多个已知数列的子序列,且是所有符合此条件序列中最长的,则 ...
- 准备NOIP2017 最长公共子序列(模版)
一些概念: (1)子序列: 一个序列A = a1,a2,--an,中任意删除若干项,剩余的序列叫做A的一个子序列.也可以认为是从序列A按原顺序保留任意若干项得到的序列.例如: 对序列 1,3,5, ...
随机推荐
- Oracle 12C -- Plug in a Non-CDB as a PDB
1.备份non-CDB数据库2.关闭non-CDB数据库 SQL> shutdown immediate; 3.将non-CDB至于只读状态 SQL> startup open read ...
- 关于ViewPager设置属性页setCurrentItem会阻塞主线程ANR总结
关于android开发设置View Pager的直接跳转页set CurrentItem会阻塞主线程ANR. 根据网上解决的说法,分析源码: if (mFirstLayout) { // We don ...
- MFC添加自定义窗口消息
原文链接: http://www.cnblogs.com/smartvessel/archive/2011/07/18/2109472.html 1. 在头文件stdafx.h中增加一个自定义消息宏 ...
- 在vmware fusion上安装linux系统然后再安装vmware tools出现问题
/usr/lib/vmware-tools] In which directory do you want to install the documentation files? [/usr/shar ...
- SQL Server中关于跟踪(Trace)那点事(转载)
前言 一提到跟踪俩字,很多人想到警匪片中的场景,同样在我们的SQL Server数据库中“跟踪”也是无处不在的,如果我们利用好了跟踪技巧,就可以针对某些特定的场景做定向分析,找出充足的证据来破案. 简 ...
- 用python参加Kaggle的经验总结【转】
用python参加Kaggle的经验总结 转载自:http://www.jianshu.com/p/32def2294ae6,作者 JxKing 最近挤出时间,用python在kaggle上试了 ...
- Java实现打印功能
用java实现打印,java.awt中提供了一些打印的API,要实现打印,首先要获得打印对象,然后继承Printable实现接口方法print,以便打印机进行打印,最后用Graphics2D直接输出直 ...
- 微信小程序 confirm(删除提示)提示框,询问框,小程序操作成功提示后跳转
微信小程序删除处理 没有 confrim 那怎么实现这个效果呢 可以使用小程序里的模态框 代码: wx.showModal({ title: '提示', content: '确定要删除吗?', suc ...
- MySQL 5.6学习笔记(数据库基本操作,查看和修改表的存储引擎)
1. 数据库基本操作 1.1 查看数据库 查看数据库列表: mysql> show databases; +--------------------+ | Database | +------ ...
- 配置并使用Android支持的库
原文链接:http://android.eoe.cn/topic/android_sdk Android Support Library(支持库)提供了包含一个API库的JAR文件,当你的应用运行在A ...