传送门

首先那个\(O(n^2)\)的dp都会吧,不会自己找博客或者问别人,或是去做模板题(误)

对以下内容不理解的,强势推荐flash的博客

我们除了原来记录最长上升子序列的\(f_{i,j}\),再记\(g_{i,j}\)表示到\(i,j\)时的最长上升子序列个数,同时设两个字符串为\(A,B\)

若\(A_i=B_j\) ,则有\(f_{i,j}=f_{i-1,j-1}+1,g_{i,j}=g_{i-1,j-1}\)

否则\(f_{i,j}=max(f_{i-1,j},f_{i,j-1}),g_{i,j}\)的话看能否从\(f_{i-1,j}\)或\(f_{i,j-1}\)转移,如果可以就加上对应的\(g\)

注意,如果\(f_{i,j}=f_{i-1,j-1}\)(等价于\(f_{i,j}=f_{i-1,j}=f_{i,j-1}\)),那么\(g_{i-1,j-1}=g_{i-1,j}=g_{i,j-1}=\frac{1}{2}g_{i,j}\),所以要减去\(g_{i-1,j-1}\)

至于为什么,yyb都没写,我也不好说.,反正要么感性理解(雾),要么打表手玩.还不理解就戳flash吧qwq

// luogu-judger-enable-o2

#include<bits/stdc++.h>

#define LL long long

#define il inline

#define re register

#define db double

#define max(a,b) ((a)>(b)?(a):(b))

using namespace std;

const int N=5000+10,mod=100000000;

il LL rd()

{

    re LL x=0,w=1;re char ch=0;

    while(ch<'0'||ch>'9') {if(ch=='-') w=-1;ch=getchar();}

    while(ch>='0'&&ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}

    return x*w;

}

int f[2][N],g[2][N],n,m;

char cc[N],ss[N];

int main()

{

  scanf("%s%s",cc+1,ss+1);

  n=strlen(cc+1)-1,m=strlen(ss+1)-1;

  g[0][0]=g[1][0]=1;

  for(int j=1;j<=m;j++) g[0][j]=1;

  int nw=1,la=0;

  for(int i=1;i<=n;i++)

    {

      for(int j=1;j<=m;j++)

        {

    	  f[nw][j]=g[nw][j]=0;

          if(cc[i]==ss[j]) f[nw][j]=f[la][j-1]+1,g[nw][j]=g[la][j-1];

          else f[nw][j]=max(f[nw][j-1],f[la][j]);

          if(f[nw][j]==f[la][j]) g[nw][j]+=g[la][j];

          if(f[nw][j]==f[nw][j-1]) g[nw][j]+=g[nw][j-1];

          if(f[nw][j]==f[la][j-1]) g[nw][j]-=g[la][j-1];

          g[nw][j]%=mod;

        }

      swap(nw,la);

    }

  printf("%d\n%d\n",f[la][m],g[la][m]);

  return 0;

}