最长公共子串 动态规划_最长公共子串 DNA序列

最长公共子串 动态规划_最长公共子串 DNA序列原题链接题目描述给定两个字符串str1和str2,输出连个字符串的最长公共子序列。如过最长公共子序列为空,则输出-1。输入描述:输出包括两行,第一行代表字符串str1,第二行代表str2。( 1<= length(str1),length(str2)<= 5000)输出描述:输出一行,代表他们最长公共子序列。如果公共子序列的长度为空,则输出-1。示例1输入1A2C3D4B56B1D23CA45B6A输出123456说明”123456″和“12C4B6”都是最长公共


Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

原题链接

题目描述
给定两个字符串str1和str2,输出连个字符串的最长公共子序列。如过最长公共子序列为空,则输出-1。
输入描述:
输出包括两行,第一行代表字符串str1,第二行代表str2。( 1<= length(str1),length(str2)<= 5000)
输出描述:
输出一行,代表他们最长公共子序列。如果公共子序列的长度为空,则输出-1。
示例1
输入

1A2C3D4B56
B1D23CA45B6A

输出

123456

说明
“123456”和“12C4B6”都是最长公共子序列,任意输出一个。
备注:

时间复杂度O(nm)O(n∗m),空间复杂度O(nm)O(n∗m)。(n,m分别表示两个字符串长度)

#include<bits/stdc++.h>
#define x first
#define y second
#define send string::nops
using namespace std;
typedef long long ll;
const int N = 1e4 + 10;
const int M = 3 * N;
const int INF = 0x3f3f3f3f;
typedef pair<int,int> PII;
int f[N][N],path[N][N];//path用来记录dp路径
int main(){ 
   
    string str1,str2;
    cin>>str1>>str2;
    int n = str1.size(),m = str2.size();
    str1.insert(0,1,' ');
    str2.insert(0,1,' ');
// cout<<str1<<endl;
    for(int i = 1;i < str1.size();i ++){ 
   
        for(int j = 1;j < str2.size();j ++){ 
   
            if(str1[i] == str2[j])f[i][j] = f[i - 1][j - 1] + 1,path[i][j] = m * (i - 1) + j - 1;
            else { 
   
                f[i][j] = max(f[i - 1][j],f[i][j - 1]);
                if(f[i - 1][j] > f[i][j - 1])path[i][j] = m * (i - 1) + j;
                else path[i][j] = m * i + j - 1;
            }
        }
    }
    int a = str1.size() - 1,b = str2.size() - 1;
    vector<char>res;
    while(a > 0 && b > 0){ 
   
// cout<<a <<"--"<<b<<endl;
// cout<<f[a][b]<<endl;
// cout<<f[a - 1][b - 1]<<endl;
// cout<<f[a - 1][b]<<endl;
// cout<<f[a][b - 1]<<endl;
        int tt = path[a][b];
        int x = tt / m,y = tt % m;
// cout<<x<<"--"<<y<<endl;
        if(x == a - 1 && y == b - 1 && !(x == 0 && y == 0))res.push_back(str1[a]);
        a = x,b = y;
    }
    reverse(res.begin(),res.end());
    if(res.size() == 0){ 
   
        cout<<-1<<endl;
        return 0;
    }
    for(int i = 0;i < res.size();i ++)cout<<res[i];
    return 0;
}

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/168791.html原文链接:https://javaforall.net

赞 (0)
全栈程序员-站长的头像全栈程序员-站长

相关推荐

  • 因特网是一种什么网_广域网是因特网吗

    因特网是一种什么网_广域网是因特网吗以太网详解(一)

    2022年10月18日
  • 一篇文带你从0到1了解建站及完成CMS系统编写

    一篇文带你从0到1了解建站及完成CMS系统编写学习目标了解搭建一般网站的简便方式了解最原始一般站点搭建了解内容管理站点搭建了解权限设计及完成了解使用设计模式减少代码冗余了解前端拖拽页面生成及生成了解自定义数据的创建了解动态生成的前端页如何绑定自定义数据开发环境Windows7*64SP1php5.6apache/nginxthinkphp5.1mysqlphpstudy2018sqlyoglayoutit声明文章为从0到1了解内容管理系统搭建与编写,由于一篇文章内容篇幅过长,文章内容经过压缩,该项目中相

    2022年6月5日
  • Oracle新建索引、删除索引和查看索引

    Oracle新建索引、删除索引和查看索引1、创建索引createindex索引名on表名(列名);2、删除索引dropindex索引名;3、创建组合索引createindex索引名on表名(列名1,列名2,列名3,…);查看目标表中已添加的索引–在数据库中查找表名select*fromuser_tableswheretable_namelike‘tablename%’;–查看该表的所有索引select*fromall_indexeswheretable_name=‘tablen

    2022年9月4日
  • webstrom最新激活码[在线序列号]

    webstrom最新激活码[在线序列号],https://javaforall.net/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

    2022年3月20日
  • 树莓派连接WIFI显示无效参数「建议收藏」

    树莓派连接WIFI显示无效参数「建议收藏」把你的wifi名字从中文改成英文或者数字就ok了

    2022年5月1日
  • RGBD融合原理及实践[通俗易懂]

    RGBD融合原理及实践[通俗易懂]RGBD融合原理及实践前言原理部分实践前言好久没更新博客了,主要是因为懒,最近有些得闲,决定纪录下之前的工作。RT,RGBD数据融合其实就是将3D摄像机的RGB与Depth数据做融合显示的过程,做法也不难理解,就是将depthcamera与rgbcamera的像素对应起来即可。原理部分原理部分主要借鉴这篇博文,详细的公式在这就不作重复了,贴张图吧。从上面的博客或图片可以看出,…

    2022年9月17日

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号