jzoj5424【NOIP2017提高A组集训10.25】凤凰院凶真

题目

有一个有n个数的集合a，m个数的集合b
求它们的最长公共上升子序列
n,m<=5000

题解

设f[i][j]表示a考虑到第i个数，b考虑到第j个数，并且b[j]一定要取的方案数
那么有f[i][j]=f[i-1][j]
当a[i]=b[j]时，可以从f[i][j]还可以由f[i-1][k](b[k]< b[j])转移而来
方案的话顺便记录一下就好了

贴代码

#include<iostream>#include<algorithm>#include<cstdio>#include<cmath>#include<cstring>#define fo(i,a,b) for(i=a;i<=b;i++)#define fo1(i,b,a) for(i=b;i>=a;i--)using namespace std;const int maxn=5e3+5;int f[maxn][maxn],g[maxn][maxn][2];bool bz[maxn][maxn];int a[maxn],b[maxn],an[maxn];int i,j,k,l,n,m,x,y,z,ans,t1,t2,t3,t4,now;int main(){    freopen("okarin.in","r",stdin);    freopen("okarin.out","w",stdout);    scanf("%d",&n);    fo(i,1,n) scanf("%d",&a[i]);    scanf("%d",&m);    fo(i,1,m) scanf("%d",&b[i]);    fo(i,1,n){        x=0; y=0; z=0;        fo(j,1,m){            f[i][j]=f[i-1][j];            g[i][j][0]=g[i-1][j][0];            g[i][j][1]=g[i-1][j][1];            if (bz[i-1][j]==true){                g[i][j][0]=i-1;                g[i][j][1]=j;            }            if (b[j-1]<a[i] && f[i-1][j-1]>x){                x=f[i-1][j-1]; y=i-1; z=j-1;            }            if (a[i]==b[j] && x+1>f[i][j]){                f[i][j]=x+1;                g[i][j][0]=g[y][z][0];                g[i][j][1]=g[y][z][1];                if (bz[y][z]==true){                    g[i][j][0]=y;                    g[i][j][1]=z;                }                bz[i][j]=true;            }            if (f[i][j]>ans){                ans=f[i][j];                t1=i;                t2=j;            }        }    }    printf("%d\n",ans);    while (t1>0 && t2>0){        an[++now]=a[t1];        t3=t1; t4=t2;        t1=g[t3][t4][0];        t2=g[t3][t4][1];    }    fo1(i,now,1) printf("%d ",an[i]);    return 0;}