[HAOI2010]最长公共子序列

一、题目

点此看题

二、解法

设$f[i][j]$为$A$串的前$i$个字符和$B$串的前$j$个字符的最大匹配长度，转移不难。

然后方案数我们就再维护一个$g[i][j]$，从哪里转移来就加上哪里的$g$，但是这里要注意一种情况，当$f[i][j]=f[i-1][j-1]$时，说明$f[i][j]=f[i][j-1]=f[i-1][j]$，而且$f[i-1][j]$和$f[i][j-1]$必定是包含$f[i-1][j-1]$的选取方案的，所以如果出现这种情况我们需要减去$g[i-1][j-1]$。

最后要注意的就是要写滚动数组，时间复杂度$O(nm)$，贴个代码。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int M = 5005;
const int MOD = 100000000;
int read()
{
    int x=0,flag=1;char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int n,m,c,f[2][M],g[2][M];
char a[M],b[M];
int main()
{
    scanf("%s",a+1),n=strlen(a+1)-1;
    scanf("%s",b+1),m=strlen(b+1)-1;
    for(int i=0;i<=m;i++)
        g[0][i]=1;
    g[1][0]=1;
    for(int i=1;i<=n;i++)
    {
        c^=1;
        for(int j=1;j<=m;j++)
        {
            g[c][j]=0;
            f[c][j]=max(f[c^1][j],f[c][j-1]);
            if(a[i]==b[j]) f[c][j]=max(f[c][j],f[c^1][j-1]+1);
            if(a[i]==b[j] && f[c][j]==f[c^1][j-1]+1) g[c][j]+=g[c^1][j-1];
            if(f[c][j]==f[c^1][j]) g[c][j]+=g[c^1][j];
            if(f[c][j]==f[c][j-1]) g[c][j]+=g[c][j-1];
            if(f[c^1][j-1]==f[c][j]) g[c][j]-=g[c^1][j-1];
            g[c][j]%=MOD;
        }
    }
    printf("%d\n%d\n",f[c][m],g[c][m]);
}