题目链接：https://darkbzoj.cf/problem/3864

这个是DP套DP的典型例题。

常规LCS的求法：
设dp[i][j]表示到a[i],b[j]的LCS。
那么if \(a[i]!=b[j]\),then \(dp[i][j]=max(dp[i][j-1],dp[i-1][j])\)
if \(a[i]==b[j]\),then \(dp[i][j]=dp[i-1][j-1]+1\)

现在我们知道字符串a，但是不知道b。但是很显然的，对于b字符串，每增加一位（枚举），必然会有对相应LCS的贡献。
我们设f[i][s]表示当前处理到字符串b的第i位，且当前它与a的LCS的状态为s。g[s][0/1/2/3]表示在当前LCS状态为s的时候，在后面添加'A','T','C','G’之后的LCS情况。
那么可以列出递推式\(f[i][g[s][0/1/2/3]]+=f[i-1][s]\)
\(g\)数组的处理可以写到函数里，进行预处理。

然后就是要注意。。。。空间的问题，这题空间卡得太丧心病狂了。。。。。qwqwqwq

代码如下：

#include<cstring>
#include<cstdio>
#include<algorithm>
#define mod 1000000007
using namespace std;
int n,m,t;
int g[(1<<15)+2][4],f[1001][(1<<15)+2],ans[111],pre[16],now[16],a[16];
char ss[16];
inline int calc(int x)
{int cur_ans=0;for(int i=20;i>=0;i--)if(x&(1<<i))cur_ans++;return cur_ans;
}
inline int trans(char x)
{if(x=='A') return 0;else if(x=='T') return 1;else if(x=='C') return 2;else return 3;
}
inline int solve(int s,int k)
{memset(pre,0,sizeof(pre));memset(now,0,sizeof(now));for(int i=0;i<n;i++) pre[i+1]=pre[i]+((s>>i)&1);for(int i=1;i<=n;i++){if(k==a[i]) now[i]=max(now[i],pre[i-1]+1);now[i]=max(max(now[i],now[i-1]),pre[i]);}int cur_ans=0;for(int i=0;i<n;i++) cur_ans+=(now[i+1]-now[i])*(1<<i);return cur_ans;
}
int main()
{scanf("%d",&t);while(t--){memset(f,0,sizeof(f));memset(g,0,sizeof(g));memset(ans,0,sizeof(ans));scanf("%s",ss+1);n=strlen(ss+1);for(int i=1;i<=n;i++) a[i]=trans(ss[i])+1;int maxx=1<<n;scanf("%d",&m);for(int i=0;i<=maxx;i++)for(int j=1;j<=4;j++)g[i][j]=solve(i,j);f[0][0]=1;for(int i=1;i<=m;i++)for(int j=0;j<=maxx;j++)for(int k=1;k<=4;k++)f[i][g[j][k]]=(f[i][g[j][k]]+f[i-1][j])%mod;for(int i=0;i<=maxx;i++)ans[calc(i)]=(ans[calc(i)]+f[m][i])%mod;for(int i=0;i<=n;i++)printf("%d\n",ans[i]);}return 0;
}