求S串与T串的 最长公共子序列 的 长度 及其 个数.
动态规划递推式：
[f(i,j)=maxleft{ f(i-1,j), f(i,j-1)
ight} quad (S_i
eq T_j)
]
[f(i,j)=maxleft{ ext{ $^{(varphi)}$ } f(i-1,j), f(i,j-1), ext{ $^{(lambda)}$ } f(i-1,j-1)+1
ight} quad (S_i= T_j)
]
[g(i,j)= g(i-1,j) + g(i,j-1) - g(i-1, j-1) quad left(S_i
eq T_j ext{ and }f(i,j)=f(i-1,j-1)
ight)
]
[g(i,j)= g(i-1,j) + g(i,j-1) + g(i-1, j-1) quad left(S_i= T_j ext{ and }(lambda)
ight)
]
滚动数组压缩空间.
#include <cstdio>
#include <cstring>
#include <algorithm>
const int mod=100000000;
char S[5003], T[5003];
int sl, tl;
int f[2][5003], g[2][5003];
int main() {
scanf("%s", S+1);
scanf("%s", T+1);
sl=strlen(S+1)-1;
tl=strlen(T+1)-1;
for (register int i=0; i<=tl; ++i) g[0][i]=1;
for (register int i=1; i<=sl; ++i) {
g[i&1][0]=1;
for (register int j=1; j<=tl; ++j) {
g[i&1][j]=0;
f[i&1][j]=std::max(f[i&1^1][j], f[i&1][j-1]);
if (S[i]==T[j]) f[i&1][j]=std::max(f[i&1][j], f[i&1^1][j-1]+1);
if (f[i&1][j]==f[i&1^1][j]) g[i&1][j]+=g[i&1^1][j];
if (f[i&1][j]==f[i&1][j-1]) g[i&1][j]+=g[i&1][j-1];
if (S[i]==T[j] && f[i&1][j]==f[i&1^1][j-1]+1) g[i&1][j]+=g[i&1^1][j-1];
if (S[i]!=T[j] && f[i&1][j]==f[i&1^1][j-1]) g[i&1][j]-=g[i&1^1][j-1];
g[i&1][j]=(g[i&1][j] + mod) % mod;
}
}
printf("%d
%d
", f[sl&1][tl], g[sl&1][tl]);
return 0;
}