最长公共子序列和最长递增子序列模板

一：最长公共子序列：（连续）

 

模板：

/*
输入两个字符串序列，求最长公共子序列及其长度。
*/
#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
using namespace std;
const int maxn =105;
char x[maxn],y[maxn];/*两个字符串序列*/
int c[maxn][maxn];/*记录序列Xi和Yj的最长公告子序列的长度。Xi={x1,x2,...,xi};Yj={y1,y2,...,yj}*/
int b[maxn][maxn];/*记录c[i][j]的值时由哪个子问题的解得到的。记录路径*/

void lcsTable(){/*建立一张路径长度表*/
    int m=strlen(x)+1;
    int n=strlen(y)+1;
    for(int i=0;i<m;i++)
        c[i][0]=0;
    for(int j=0;j<n;j++)
        c[0][j]=0;
    for(int i=1;i<m;i++){
        for(int j=1;j<n;j++){
            if(x[i-1]==y[j-1]){
                c[i][j]=c[i-1][j-1]+1;
                b[i][j]=1;
            }else if(c[i-1][j]>c[i][j-1]){
                c[i][j]=c[i-1][j];
                b[i][j]=2;
            }else{
                c[i][j]=c[i][j-1];
                b[i][j]=3;
            }
        }
    }
}

void lcs(int i,int j){/*i,j分别为序列x和序列y的长度*/
    if(i==0 || j==0) return;
    if(b[i][j]==1){
        lcs(i-1, j-1);
        printf("%c ",x[i-1]);
    }else if(b[i][j]==2)
        lcs(i-1,j);
    else
        lcs(i,j-1);
}

int main()
{
   // freopen("input.txt","r",stdin);
    while(~scanf("%s%s",x,y)){
        lcsTable();
        printf("%d\n",c[strlen(x)][strlen(y)]);/*最长公子序列的长度*/
        lcs(strlen(x),strlen(y));/*最长公子序列*/
    }
    return 0;
}

二：最长递增子序列：（非连续）

模板：

/*
第一行输入n，表示序列的长度。
第二行输入序列。
求：最长递增子序列及其长度。
*/
#include<cstdio>
#include<iostream>
using namespace std;
const int maxn=1005;
int a[maxn];/*输入序列*/
int b[maxn];/*与a[]对应，以a[i]结尾的递增序列的长度*/
int c[maxn];/*记录路径*/
int n;/*a[]的长度*/
int len,k;/*最长递增序列的长度及其最后一个元素的位置*/

void LIS(){
    b[0]=1;
    c[0]=-1;
    for(int i=1;i<n;i++){
        int m=0;
        for(int j=0;j<i;j++){
            if(a[i]>a[j] && b[j]>m){
                m=b[j];
                c[i]=j;
            }
        }
        b[i]=m+1;
    }
    len=b[0],k=0;
    for(int i=1;i<n;i++){
        if(b[i]>len){
            len=b[i];
            k=i;
        }
    }
}

void traceBack(int i){
    if(i==-1) return;
    traceBack(c[i]);
    printf("%d ",a[i]);
}

int main()
{
  //  freopen("input.txt","r",stdin);
    while(~scanf("%d",&n)){
        for(int i=0;i<n;i++)
            scanf("%d",&a[i]);
        LIS();
        printf("%d\n",len);/*最长子序列的长度*/
        traceBack(k);/*最长子序列的其中一条*/
    }
    return 0;
}

用二分法优化下：

/*
第一行输入n，k，表示序列的长度，子序列元素在母序列中的下标间隔距离。
第二行输入序列。
求：最长递增子序列的长度。
*/
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=100005;
int n,k;/*序列长度；k为子序列元素在母序列中的下标间隔距离，一般为0*/
int a[maxn],dp[maxn],c[maxn];
/*a数组存放序列
  dp记录在i点时最长的递增子序列长度
  c数组为每次查找时候的标记，记录路径
*/

int bs(int t)
{
    int l = 1,r = n;
    while(l<=r)
    {
        int mid = (l+r)/2;
        if(t>c[mid])
            l = mid+1;
        else
            r = mid-1;
    }
    return l;
}

int LIS()
{
    int ans = 0;
    for(int i =1; i<=n; i++)
    {
        dp[i] = bs(a[i]);
        if(dp[i]>ans)/*更新最长长度*/
            ans = dp[i];
        int j = i-k;/*k为相邻元素的间隔距离，一般为0*/
        if(j>0 && c[dp[j]]>a[j])/*查找标记*/
            c[dp[j]] = a[j];
    }
    return ans;
}

int main()
{
  //  freopen("input.txt","r",stdin);
    while(~scanf("%d%d",&n,&k)){
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            c[i]=10000000;
        }
        printf("%d\n",LIS());
    }
}

例题:小明系列问题——小明序列