1. 程式人生 > >最長公共子序列、最長上升子序列、最長公共上升子序列

最長公共子序列、最長上升子序列、最長公共上升子序列

1.最長公共子序列:

核心程式碼:

for(int i=1;i<n;i++)
    {
        for(int j=1;j<m;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]);
            }
        }
    }

hdu5248

2.求最長上升子序列

樹狀陣列優化:nlogn

#include<bits/stdc++.h>
using namespace std;
const int maxn=100005;
char a[maxn];
int len;
int bit[maxn],v[maxn],A[maxn],L[maxn];
int lowbit(int x)
{
    return x&(-x);
}
void add(int x,int d)
{
    for(int i=x;i<=len;i+=lowbit(i))
    {
        if
(bit[i]<d)bit[i]=d; } } int ask(int x) { int res=0; for(int i=x;i>=1;i-=lowbit(i)) { if(res<bit[i])res=bit[i]; } return res; } int main() { int t,n; scanf("%d",&t); while(t--) { scanf("%d",&n); for(int i=1;i<=n;i++) { scanf
("%d",&A[i]); v[i]=A[i]; } len=unique(v+1,v+1+n)-v-1; for(int i=1;i<=n;i++) { int x=lower_bound(v,v+len+1,A[i])-v;// int sum=ask(x-1)+1;//用樹狀陣列查詢比A[i]小的個數再加上1即是以A[i]結尾的最長上升序列 L[i]=sum; //詢問的數在i之前插入。 add(x,sum); } for(int i=1;i<=n;i++) { if(i!=1)printf(" "); printf("%d",L[i]); } printf("\n"); } return 0; }

以上程式碼可以用滾動陣列進行空間優化

3.最長公共上升子序列:(O(N3

解析:f[i][j]表示以b[j]結尾,字串a[i]之前的公共上升子序列最大長度;
顯然:f[i][j]>=f[i1][j];
遞推:若a[i]!=b[j]f[i][j]=f[i1[j];
a[i]==b[j]f[i][j]=max(f[k][j])+1;(1<=k<=j-1&&b[j]>b[k])
核心程式碼

 for(int i=1;i<n;i++)
    {
        for(int j=1;j<m;j++)
        {
            if(a[i]!=b[j])f[i][j]=f[i-1][j];
            else
            {
                int mn=0;
                for(int k=1;k<j;k++)
                {
                    if(b[j]>b[k])
                    mn=max(mn,f[i-1][k]);
                }
                f[i][j]=mn+1;
                ans=max(f[i][j],ans);
            }
        }
    }

4.最長公共上升子序列:(O(N^2))

核心程式碼:
優化:通過記錄對於當前a[i]與b[1-m]的匹配過程中的最大上升子序列mm,邊記錄,邊更新,降低一維

for(int i=1;i<n;i++)
    {
        int mm=0;
        for(int j=1;j<m;j++)
        {
            if(a[i]>b[j])mm=max(mm,f1[j]);//優化
            else if(a[i]==b[j])
                f1[j]=mm+1,ans2=max(ans2,f1[j]);
        }
    }

測試源程式:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int main()
{
    char a[100];
    char b[100];
    a[0]='0';b[0]='0';
    int dp[105][105];
    int f[105][105];
    int f1[105];
    memset(dp,0,sizeof dp);
    scanf("%s%s",a+1,b+1);
    int n=strlen(a);
    int m=strlen(b);
    for(int i=1;i<n;i++)
    {
        for(int j=1;j<m;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]);
            }
        }
    }
    memset(f,0,sizeof f);
    int ans=0;
    for(int i=1;i<n;i++)
    {
        for(int j=1;j<m;j++)
        {
            if(a[i]!=b[j])f[i][j]=f[i-1][j];
            else
            {
                int mn=0;
                for(int k=1;k<j;k++)
                {
                    if(b[j]>b[k])
                    mn=max(mn,f[i-1][k]);
                }
                f[i][j]=mn+1;
                ans=max(f[i][j],ans);
            }
        }
    }
    memset(f1,0,sizeof f1);
    int ans2=0;
    for(int i=1;i<n;i++)
    {
        int mm=0;
        for(int j=1;j<m;j++)
        {
            if(a[i]>b[j])mm=max(mm,f1[j]);
            else if(a[i]==b[j])
                f1[j]=mm+1,ans2=max(ans2,f1[j]);
        }
    }
    cout<<"最長上升子序列:"<<dp[n-1][m-1]<<endl;
    cout<<"最長公共上升子序列n*n*n:"<<ans<<endl;
    cout<<"最長公共上升子序列n*n:"<<ans2<<endl;
    return 0;
}