题意
给出一个长度为n的基因序列S(仅包含字母ACGT的串)和一个数m。
对于每个x=0~n,求有多少个长度为m的基因序列T满足LCS(S,T)=x。
n,m<=15
题解
外层:
f(i,g)
表示T考虑到前i位,状态为s的方案数。
这个内层的状态g是一个长度为n的序列,表示LCS(i,1~n)这个的第二维。
可以想到这个序列是不降的且每次最多加1,所以只需要差分后用一个二进制数保存即可。
内层状态的转移需要预处理出来,
nxt[g][k]=gnow
表示在g这个状态时,在T后多加一个字母k得到的新状态
gnow
。具体转移的话就是和经典的二维DP求LCS一样的道理。
然后外层每次枚举T后面放那个字母,然后简单转移即可。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const char ltr[4]={'A','C','G','T'};
const int MOD=1000000007;
int _test,n,m,nxt[(1<<15)+5][4],f[2][(1<<15)+5],ans[20],res[(1<<15)+5];
int g1[20],g2[20];
char st[20];
void Pre(){
for(int s=0;s<=(1<<n)-1;s++){
g1[0]=0; for(int i=1;i<=n;i++) g1[i]=g1[i-1]+((s>>i-1)&1); res[s]=g1[n];
for(int k=0;k<=3;k++){
g2[0]=0;
for(int i=1;i<=n;i++){
g2[i]=max(g2[i-1],g1[i]);
if(ltr[k]==st[i]) g2[i]=max(g2[i],g1[i-1]+1);
}
nxt[s][k]=0; for(int i=1;i<=n;i++) nxt[s][k]|=((g2[i]-g2[i-1])<<i-1);
}
}
}
int main(){
freopen("bzoj3864.in","r",stdin);
freopen("bzoj3864.out","w",stdout);
scanf("%d",&_test);
while(_test--){
scanf("%s%d",st+1,&m); n=strlen(st+1);
Pre();
memset(f,0,sizeof(f)); f[0][0]=1;
for(int i=0;i<=m-1;i++){
memset(f[(i&1)^1],0,sizeof f[(i&1)^1]);
for(int s=0;s<=(1<<n)-1;s++)
for(int k=0;k<=3;k++) f[(i&1)^1][nxt[s][k]]=(f[(i&1)^1][nxt[s][k]]+f[i&1][s])%MOD;
}
memset(ans,0,sizeof ans);
for(int s=0;s<=(1<<n)-1;s++) ans[res[s]]=(ans[res[s]]+f[m&1][s])%MOD;
for(int i=0;i<=n;i++) printf("%d\n",ans[i]);
}
return 0;
}