NOIP5424. 【NOIP2017提高A组集训10.25】凤凰院凶真 LCIS 方案

题意：求lcis以及其方案。
原题：千帆渡
n<=5000.
很明显的DP。
设f[i][j]表示做到a的第i位，结尾是b的第j位，这样子便于记录方案。
那么明显有
f[i][j]=max(f[i][k-1])+1(b[j]>b[k]&&j>k)
由于规定了a[i]=b[j]，所以上面的可以变成
f[i][j]=max(f[i][k-1])+1(a[i]>b[k]&&j>k)
至于j>k这个显然。
所以我们可以在dp过程中不断更新k，然后扫完b以后更新。
轻松切，然后至于方案的话记录一下每个状态由什么转移过来，然后由于我们的dp的状态定义是一定以j结尾的，所以这就很简单了。

#include<cstdio>#include<algorithm>#include<cstring>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fd(i,a,b) for(int i=a;i>=b;i--)using namespace std;const int N=5e3+5;int a[N],b[N];int f[N][N],n,m,lis[N];struct node{    int x,y,z;}pre[N][N];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)    {        int k=0,x=0,y=0;        fo(j,1,m)        {            f[i][j]=f[i-1][j];            pre[i][j].z=f[i-1][j];            pre[i][j].x=i-1,pre[i][j].y=j;            if (a[i]>b[j]&&f[i-1][j]>k)            {                k=f[i-1][j];                x=i-1,y=j;            }            if (a[i]==b[j])            {                f[i][j]=k+1;                pre[i][j].z=k;                pre[i][j].x=x;                pre[i][j].y=y;            }        }    }    int ans=0,p=0;    fo(j,0,m)    if (ans<f[n][j])ans=f[n][j],p=j;    printf("%d\n",ans);    int minx=1e9,tot=0;    while(f[n][p])    {        if (f[n][p]<minx&&f[n][p]>0)        {            minx=f[n][p];            lis[++tot]=b[p];        }        int xx=pre[n][p].x;        int yy=pre[n][p].y;        n=xx,p=yy;    }    fd(i,tot,1)printf("%d ",lis[i]);    printf("\n");}