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;
}