[DP] bzoj2423: [HAOI2010]最长公共子序列

bzoj2423: [HAOI2010]最长公共子序列:https://www.lydsy.com/JudgeOnline/problem.php?id=2423

DP
神啊（DP都是神奇的东西）
第一问直接N^2求
第二问怎么办呢？
设f2[i][j]表示a[ ]匹配到第i个 b[ ]匹配到第j个 的最长子序列个数
对于每一位f2[i][j]
如果a[i]==b[i] 加上f1[i-1][j-1]的个数
如果a[i]==b[i]并且f1[i-1][j-1]+1==f1[i][j] 说明当前这个位置是最长的子序列的结尾 加上f[i-1][j-1]的个数
如果当前这一位是从f1[i-1][j]/f1[i][j-1]继承来的 那就加上
但是如果它和f1[i-1][j-1]相等的话 那么f1[i-1][j-1]==f1[i-1][j]==f1[i][j-1]==f1[i][j] 那一定会加重复 然后减掉就好了
开滚动数组然后O(N^2)跑就ok

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
long long f1[2][5100],f2[2][5100];
int m=100000000;
char s1[5100],s2[5100];
int main()
{
    int len1,len2;
    scanf("%s",s1+1); scanf("%s",s2+1);
    len1=strlen(s1+1);len2=strlen(s2+1);
    len1--;len2--;
    int x=0,y=1;
    for (int i=0;i<=len2;i++) f2[0][i]=1;
    f2[1][0]=1;
    for (int i=1;i<=len1;i++)
    {
        x^=1;y^=1;
        for (int j=1;j<=len2;j++)
        {
            f1[x][j]=max(f1[x][j-1],f1[y][j]);//继承 
            f2[x][j]=0;
            if (s1[i]==s2[j]) f1[x][j]=max(f1[x][j],f1[y][j-1]+1);
            if (s1[i]==s2[j]&&f1[x][j]==f1[y][j-1]+1) f2[x][j]+=f2[y][j-1];
            if (f1[x][j]==f1[y][j]) f2[x][j]+=f2[y][j];
            if (f1[x][j]==f1[x][j-1]) f2[x][j]+=f2[x][j-1];
            if (f1[x][j]==f1[y][j-1]) f2[x][j]-=f2[y][j-1];
            f2[x][j]=(f2[x][j]+m)%m;
        }
    }
    printf("%lld\n%lld\n",f1[x][len2],f2[x][len2]);
    return 0;
}