分析:
YY了一个做法,没想到真的能过。。。(太乱来了)
观察常规的最长公共子序列问题的DP方式:DP[i][j]表示第1个串匹到i位置,第2个串匹到j位置的最长公共子序列长度。
然后有个很显然的性质:
D
P
[
i
]
[
j
−
1
]
≤
D
P
[
i
]
[
j
]
≤
D
P
[
i
]
[
j
−
1
]
+
1
DP[i][j-1]\leq DP[i][j]\leq DP[i][j-1]+1
DP[i][j−1]≤DP[i][j]≤DP[i][j−1]+1
所以可以把 D P [ i ] [ 0 ] 、 D P [ i ] [ 1 ] 、 D P [ i ] [ 2 ] … … D P [ i ] [ ∣ S ∣ ] DP[i][0]、DP[i][1]、DP[i][2]……DP[i][|S|] DP[i][0]、DP[i][1]、DP[i][2]……DP[i][∣S∣]都表示出来,根据上面的性质,差分一下就可以了。
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define SF scanf
#define PF printf
#define MAXN 1010
#define MAXK (1<<15)
#define MOD 1000000007
using namespace std;
int tran[MAXK][4];
char s[MAXN];
int len,m;
int pres[MAXN],newp[MAXN];
char letter[6]={'A','G','C','T'};
void prepare(){
for(int i=0;i<(1<<len);i++){
for(int j=0;j<len;j++)
pres[j]=(i>>j)&1;
for(int j=1;j<len;j++)
pres[j]+=pres[j-1];
for(int nxt=0;nxt<4;nxt++){
for(int j=0;j<len;j++){
newp[j]=pres[j];
if(j!=0)
newp[j]=max(newp[j],newp[j-1]);
if(s[j]==letter[nxt]){
if(j!=0)
newp[j]=max(newp[j],pres[j-1]+1);
else
newp[j]|=1;
}
}
int mask=newp[0];
for(int j=1;j<len;j++)
mask|=(1<<j)*(newp[j]-newp[j-1]);
tran[i][nxt]=mask;
}
}
}
int count_one(int x){
int sum=0;
while(x){
if(x&1)
sum++;
x>>=1;
}
return sum;
}
int dp[2][MAXK],ans[MAXN];
int main(){
int t;
SF("%d",&t);
while(t--){
SF("%s",s);
len=strlen(s);
prepare();
SF("%d",&m);
int now=0;
memset(dp,0,sizeof dp);
dp[now][0]=1;
for(int i=1;i<=m;i++){
memset(dp[now^1],0,sizeof dp[now^1]);
for(int mask=0;mask<(1<<len);mask++)
for(int nxt=0;nxt<4;nxt++){
int mask1=tran[mask][nxt];
dp[now^1][mask1]=(dp[now^1][mask1]+dp[now][mask])%MOD;
}
now^=1;
}
memset(ans,0,sizeof ans);
for(int i=0;i<(1<<len);i++)
(ans[count_one(i)]+=dp[now][i])%=MOD;
for(int i=0;i<=len;i++)
PF("%d\n",ans[i]);
}
}