题意:给一个长度<=15的串S(仅包含ACGT),对于所有的i(0<=i<=|S|),求有多少长度为M的串T满足串T和串S的最长公共子序列长度为i。
这个题让我想起了之前做过的一个数位dp,问有多少个数的数字组成的最长上升子序列长度为x。那个题中显然最长上升子序列不超过10,我们用状态压缩来模拟那个做lis时的栈即可。
这个题|S|<=15,一样用状态压缩来表示那些字母可以包含在lcs中。注意求lcs的经典n^2的dp实际是在线的,即添加一个T中的字符就可以求出与S形成lcs。我们对所有状态拓展出它加上一个字符后新形成的状态即可,注意需要将二进制状态还原并暴力跑一次lcs再转为新状态。
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define erp(i,a,b) for(int i=a;i>=b;--i)
#define LL long long
using namespace std;
const int mo = 1000000007;
const char str[] = {'A','C','G','T'};
inline void up(int&x,