luogu P2516 [HAOI2010]最长公共子序列

传送门

首先那个(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
%d
",f[la][m],g[la][m]);
  return 0;
}