最長公共子序列、最長上升子序列、最長公共上升子序列
阿新 • • 發佈:2019-01-31
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.求最長上升子序列
樹狀陣列優化:
#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()
解析:f[i][j]表示以b[j]結尾,字串a[i]之前的公共上升子序列最大長度;
顯然:
遞推:若
若(f[k][j])+1;(<=j-1&&)
核心程式碼
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;
}