最长公共子序列,最长公共子串,最长单调子序列等问题入门

本文详细阐述了子串和子序列的概念区别,探讨了最长公共子序列问题及其动态规划解决方案,介绍了最大子序列和的问题及算法实现,并讨论了最长公共子串的两种求解方法。


首先来说子串和字符串的区别:

子串!=子序列

例如:一个字符串   awbcdewgh

他的子串:  awbc.   awbcd   awbcde   ....很多个子串  但是都是连续在一起  

他的子序列: abc  .  abcd    abcde  ...  很多个子序列  但是子序列中的字符在字符串中不一定是连在一起的, 但是  子序列一定是单调的, (即字符之间ASCII单调递增或单调递减)

最大子序列

最大子序列是要找出由数组成的一维数组中和最大的连续子序列。比如{5,-3,4,2}的最大子序列就是 {5,-3,4,2},它的和是8,达到最大;而 {5,-6,4,2}的最大子序列是{4,2},它的和是6。你已经看出来了,找最大子序列的方法很简单,只要前i项的和还没有小于0那么子序列就一直向后扩展,否则丢弃之前的子序列开始新的子序列,同时我们要记下各个子序列的和,最后找到和最大的子序列。

int maxSubSum(const vector<int> & arr,int &begin,int &end){
    int maxSum=0;
    int currSum=0;
    int newbegin=0;
    for(int i=0;i<arr.size();++i){
        currSum+=arr[i];
        if(currSum>maxSum){
            maxSum=currSum;
            begin=newbegin;
            end=i;
        }
        if(currSum<0){
            currSum=0;
            newbegin=i+1;
        }
    }
    return maxSum;
}



两字符序列的最长公共字符子序列

问题描述:字符序列的子序列是指从给定字符序列中随意地(不一定连续)去掉若干个字符(可能一个也不去掉)后所形成的字符

序列。令给定的字符序列X=“x0,x1,…,xm-1,序列Y=“y0,y1,…,yk-1X的子序列,存在X的一个严格递增下标序列<i

0,i1,…,ik-1>,使得对所有的j=0,1,…,k-1,有xij=yj。例如,X=“ABCBDAB”,Y=“BCDB”是X的一个子序列。

考虑最长公共子序列问题如何分解成子问题,设A=“a0,a1,…,am-1B=“b0,b1,…,bm-1,并Z=“z0,z1,…,zk-

1为它们的最长公共子序列。不难证明有以下性质:

(1) 如果am-1=bn-1,则zk-1=am-1=bn-1,且“z0,z1,…,zk-2是“a0,a1,…,am-2和“b0,b1,…,bn-2的一个最

长公共子序列;

(2) 如果am-1!=bn-1,则若zk-1!=am-1,蕴涵“z0,z1,…,zk-1是“a0,a1,…,am-2和“b0,b1,…,bn-1的一个最

长公共子序列;

(3) 如果am-1!=bn-1,则若zk-1!=bn-1,蕴涵“z0,z1,…,zk-1是“a0,a1,…,am-1和“b0,b1,…,bn-2的一个最

长公共子序列。

这样,在找A和B的公共子序列时,如有am-1=bn-1,则进一步解决一个子问题,找“a0,a1,…,am-2和“b0,b1,…,bm-2

一个最长公共子序列;如果am-1!=bn-1,则要解决两个子问题,找出“a0,a1,…,am-2和“b0,b1,…,bn-1的一个最长公

共子序列和找出“a0,a1,…,am-1和“b0,b1,…,bn-2的一个最长公共子序列,再取两者中较长者作为A和B的最长公共子

序列。

求解:

引进一个二维数组c[][],用c[i][j]记录X[i]与Y[j] 的LCS 的长度,b[i][j]记录c[i][j]是通过哪一个子问题的值求得的,以决定

搜索的方向。

我们是自底向上进行递推计算,那么在计算c[i,j]之前,c[i-1][j-1],c[i-1][j]与c[i][j-1]均已计算出来。此时我们根据X[i] =

 Y[j]还是X[i] != Y[j],就可以计算出c[i][j]。

问题的递归式写成:



分析:
由于每次调用至少向上或向左(或向上向左同时)移动一步,故最多调用(m + n)次就会遇到i = 0或j = 0的情况,此时开始返回。返回时与递归调用时方向相反,步数相同,故算法时间复杂度为Θ(m + n)。

#include <stdio.h>  
#include <string.h>  
#define MAXLEN 100  
  
void LCSLength(char *x, char *y, int m, int n, int c[][MAXLEN], int b[][MAXLEN])  
{  
    int i, j;  
      
    for(i = 0; i <= m; i++)  
        c[i][0] = 0;  
    for(j = 0; j <= n; j++)  
        c[0][j] = 0;  
    for(i = 1; i<= m; i++)  
    {  
        for(j = 1; j <= n; j++)  
        {  
            if(x[i-1] == y[j-1])  //上面第一步(1) 
            {  
                c[i][j] = c[i-1][j-1] + 1;  
                b[i][j] = 0;  
            }  
            else if(c[i-1][j] >= c[i][j-1])  //(2) 
            {  
                c[i][j] = c[i-1][j];  
                b[i][j] = 1;  
            }  
            else  //(3) 
            {  
                c[i][j] = c[i][j-1];  
                b[i][j] = -1;  
            }  
        }  
    }  
}  
  
void PrintLCS(int b[][MAXLEN], char *x, int i, int j)  //可以输出子序列
{  
    if(i == 0 || j == 0)  
        return;  
    if(b[i][j] == 0)  
    {  
        PrintLCS(b, x, i-1, j-1);  
        printf("%c ", x[i-1]);  
    }  
    else if(b[i][j] == 1)  
        PrintLCS(b, x, i-1, j);  
    else  
        PrintLCS(b, x, i, j-1);  
}  
  
int main(int argc, char **argv)  
{  
    char x[MAXLEN] = {"ABCBDAB"};  
    char y[MAXLEN] = {"BDCABA"};  
    int b[MAXLEN][MAXLEN];  
    int c[MAXLEN][MAXLEN];  
    int m, n;  
      
    m = strlen(x);  
    n = strlen(y);  
      
    LCSLength(x, y, m, n, c, b);  
    PrintLCS(b, x, m, n);  
      
    return 0;  
}  

求最大长度的核心代码:

 for(i = 0; i <= m; i++)    
        c[i][0] = 0;    
    for(j = 0; j <= n; j++)    
        c[0][j] = 0;    
    for(i = 1; i<= m; i++)       
        for(j = 1; j <= n; j++)    
        {    
            if(x[i-1] == y[j-1])  //上面第一步(1)      
                c[i][j] = c[i-1][j-1] + 1;      
              
            else (c[i-1][j] >= c[i][j-1])  //(2)      
                c[i][j] = max(c[i-1][j] >= c[i][j-1]) ;    
    }  
	printf("%d\n",c[m][n]);  





#include <cstdio>  
#include <cstring>  
const int maxn=10001;  
char s[maxn];  
int dp[maxn],Max;  
void LICS()  
{  
    int len;  
    memset(dp,0,sizeof(dp));  
    len=strlen(s);  
    for(int i=0;i<len;i++)  
    {  
        dp[i]=1;//给定一个数组求的时候,初始值就是1,一个数组的最大序列肯定会有一个字符;  
        for(int j=0;j<i;j++)  
        {  
            if(s[i]>s[j] && dp[i]<1+dp[j])// 递推公式,如果这个位置比前面的字符都大,就加入到递增序列中来  
                dp[i]=1+dp[j];  
        }  
    }  
    Max=0;  
    for(int i=0;i<len;i++)//求出最大值  
        if(Max<dp[i])  
          Max=dp[i];  
}  
int main()  
{  
    int t;  
    scanf("%d",&t);  
    while(t--)  
    {  
        scanf("%s",s);  
        LICS();  
        printf("%d\n",Max);  
    }  
    return 0;  
}  



最长公共子串:

方法一:

     Longest Common Substring和Longest Common Subsequence是有区别的

     X = <a, b, c, f, b, c>

     Y = <a, b, f, c, a, b>

     X和Y的Longest Common Sequence为<a, b, c, b>,长度为4

     X和Y的Longest Common Substring为 <a, b>长度为2

    其实Substring问题是Subsequence问题的特殊情况,也是要找两个递增的下标序列

    <i1, i2, ...ik> 和 <j1, j2, ..., jk>使

     xi1 == yj1

    xi2 == yj2

    ......

    xik == yjk

    与Subsequence问题不同的是,Substring问题不光要求下标序列是递增的,还要求每次

   递增的增量为1, 即两个下标序列为:

   <i, i+1, i+2, ..., i+k-1> 和 <j, j+1, j+2, ..., j+k-1>

    类比Subquence问题的动态规划解法,Substring也可以用动态规划解决,令

    c[i][j]表示Xi和Yi的最大Substring的长度,比如

   X = <y, e, d, f>

   Y = <y, e, k, f>

   c[1][1] = 1

   c[2][2] = 2

   c[3][3] = 0

   c[4][4] = 1

   动态转移方程为:

   如果xi == yj, 则 c[i][j] = c[i-1][j-1]+1

   如果xi ! = yj,  那么c[i][j] = 0

   最后求Longest Common Substring的长度等于

   max{  c[i][j],  1<=i<=n, 1<=j<=m}

/**  
找出两个字符串的最长公共连续子串的长度 
** author :liuzhiwei   
** data   :2011-08-16 
**/   
#include "stdio.h"  
#include "string.h"  
#include "stdlib.h"  
  
int longest_common_substring(char *str1, char *str2)  
{  
    int i,j,k,len1,len2,max,x,y;  
    len1 = strlen(str1);  
    len2 = strlen(str2);  
    int **c = new int*[len1+1];  
    for(i = 0; i < len1+1; i++)  
        c[i] = new int[len2+1];  
    for(i = 0; i < len1+1; i++)  
        c[i][0]=0;        //第0列都初始化为0  
    for(j = 0; j < len2+1; j++)  
        c[0][j]=0;        //第0行都初始化为0   
    max = -1;  
    for(i = 1 ; i < len1+1 ; i++)  
    {  
        for(j = 1; j < len2+1; j++)  
        {  
            if(str1[i-1]==str2[j-1])     //只需要跟左上方的c[i-1][j-1]比较就可以了  
                c[i][j]=c[i-1][j-1]+1;  
            else                         //不连续的时候还要跟左边的c[i][j-1]、上边的c[i-1][j]值比较,这里不需要  
                c[i][j]=0;  
            if(c[i][j]>max)  
            {  
                max=c[i][j];  
                x=i;  
                y=j;  
            }  
        }  
    }  
  
    //输出公共子串  
    char s[1000];  
    k=max;  
    i=x-1,j=y-1;  
    s[k--]='\0';  
    while(i>=0 && j>=0)  
    {  
        if(str1[i]==str2[j])  
        {  
            s[k--]=str1[i];  
            i--;  
            j--;  
        }  
        else       //只要有一个不相等,就说明相等的公共字符断了,不连续了  
            break;  
    }  
    printf("最长公共子串为:");  
    puts(s);  
    for(i = 0; i < len1+1; i++)         //释放动态申请的二维数组  
        delete[] c[i];  
    delete[] c;  
    return max;  
}  
int main(void)  
{  
    char str1[1000],str2[1000];  
    printf("请输入第一个字符串:");  
    gets(str1);  
    printf("请输入第二个字符串:");  
    gets(str2);  
    int len = longest_common_substring(str1, str2);  
    printf("最长公共连续子串的长度为:%d\n",len);  
    system("pause");  
    return 0;  
}  


方法二:

      将字符串s1和s2分别写在两把直尺上面(我依然用s1,s2来表示这两把直尺),然后将s1固定,s2的头部和s1的尾部对齐,然后逐渐移动直尺s2,比较重叠部分的字符串中的公共子串的长度,直到直尺s2移动到s1的头部。在这个过程中求得的最大长度就是s1、s2最大子串的长度。

     下图是求解过程的图示(下图有点错误,应该是将s2从右往左移动),蓝色部分表示重叠的字符串,红色的部分表示重叠部分相同的子串

      其中s1="shaohui",s2="ahui",最后求得的结果为3


/**  
找出两个字符串的最长公共连续子串的长度 
** author :liuzhiwei   
** data   :2011-08-16 
**/   
#include "stdio.h"  
#include "string.h"  
#include "stdlib.h"  
  
int longest_common_substring(char *str1, char *str2)  
{  
    int i,len1,len2,len,s1_start,s2_start,idx,curmax,max;  
    len1 = strlen(str1);  
    len2 = strlen(str2);  
    len = len1 + len2;  
    max = 0;  
    for(i = 0 ; i < len ; i++)  
    {  
        s1_start = s2_start = 0;  
        if(i < len1)  
            s1_start = len1 - i;    //每次开始匹配的起始位置  
        else  
            s2_start = i - len1;  
        curmax = 0;  
        for(idx = 0 ; ( s1_start + idx < len1 ) && ( s2_start + idx < len2 ); idx++ )  
        {  
            if(str1[s1_start+idx]==str2[s2_start+idx])  
                curmax++;  
            else     //只要有一个不相等,就说明相等的公共字符断了,不连续了,要保存curmax与max中的最大值,并将curmax重置为0  
            {  
                max = curmax > max ? curmax : max;  
                curmax = 0;  
            }  
        }  
        max = curmax > max ? curmax : max;  
    }  
    return max;  
}  
int main(void)  
{  
    char str1[1000],str2[1000];  
    printf("请输入第一个字符串:");  
    gets(str1);  
    printf("请输入第二个字符串:");  
    gets(str2);  
    int len = longest_common_substring(str1, str2);  
    printf("最长公共连续子串的长度为:%d\n",len);  
    system("pause");  
    return 0;  
}  
稍微改动一下,便可以输出公共子串了,就是要保存一下连续公共子串最后一个字符在其中一个字符串中的下标位置:

/**  
找出两个字符串的最长公共连续子串的长度 
** author :liuzhiwei   
** data   :2011-08-16 
**/   
#include "stdio.h"  
#include "string.h"  
#include "stdlib.h"  
  
int longest_common_substring(char *str1, char *str2)  
{  
    int i,k,len1,len2,len,s1_start,s2_start,idx,curmax,max;  
    len1 = strlen(str1);  
    len2 = strlen(str2);  
    len = len1 + len2;  
    max = 0;  
    for(i = 0 ; i < len ; i++)  
    {  
        s1_start = s2_start = 0;  
        if(i < len1)  
            s1_start = len1 - i;    //每次开始匹配的起始位置  
        else  
            s2_start = i - len1;  
        curmax = 0;  
        for(idx = 0 ; ( s1_start + idx < len1 ) && ( s2_start + idx < len2 ); idx++ )  
        {  
            if(str1[s1_start+idx]==str2[s2_start+idx])  
                curmax++;  
            else     //只要有一个不相等,就说明相等的公共字符断了,不连续了,要保存curmax与max中的最大值,并将curmax重置为0  
            {  
                //max = curmax > max ? curmax : max;  
                if(curmax > max)  
                {  
                    max = curmax;  
                    k = s1_start+idx-1;      //保存连续子串长度增加时连续子串最后一个字符在str1字符串中的下标位置,便于输出公共连续子串  
                }  
                curmax = 0;  
            }  
        }  
        //max = curmax > max ? curmax : max;  
        if(curmax > max)  
        {  
            max = curmax;  
            k = s1_start+idx-1;  
        }  
    }  
  
    //输出公共子串  
    char s[1000];  
    for(i=0;i<max;i++)  
        s[i]=str1[k-max+1+i];     //公共字串在str1中的下标起始位置为k-max+1,结束位置为k  
    s[i]='\0';  
    printf("最长公共子串为:");  
    puts(s);  
    return max;  
}  
int main(void)  
{  
    char str1[1000],str2[1000];  
    printf("请输入第一个字符串:");  
    gets(str1);  
    printf("请输入第二个字符串:");  
    gets(str2);  
    int len = longest_common_substring(str1, str2);  
    printf("最长公共连续子串的长度为:%d\n",len);  
    system("pause");  
    return 0;  
}  








评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值