tjut 4899

#include<stdio.h>  
#include<iostream>  
#include<algorithm>  
#include<string>  
#include<string.h>  
#include<math.h>  
#include<limits.h>  
#include<time.h>  
#include<stdlib.h>  
#include<map>  
#include<queue>  
#include<set>  
#include<stack>  
#include<vector>  
#define LL long long  
using namespace std;  
const int mod = 1e9 + 7;  
char dna[5] = "ATGC";  
LL dp[2][70000], prt[20];  
int res[20], res1[20], ans[70000][5];  
int main()  
{  
    int T;  
    while(~scanf("%d", &T))  
        while(T--)  
        {  
            memset(res, 0, sizeof res);  
            memset(ans, 0, sizeof ans);  
            memset(dp, 0, sizeof dp);  
            string s;  
            cin >> s;  
            int n;  
            cin >> n;  
            int len = s.size();  
            int maxx = 1 << len;  
            for(int i = 0; i < maxx; i++)  
            {  
                res[0] = 0;  
                int tmp = i;  
                int cnt = 1;  
                for(int j = 0; j < len; j++)//转化当前状态为正常dp  
                {  
                    if(tmp & 1)  
                        res[cnt] = res[cnt - 1] + 1;  
                    else  
                        res[cnt] = res[cnt - 1];  
                    cnt++;  
                    tmp >>= 1;  
                }  
                for(int j = 0; j < 4; j++)//计算分别加4个字符后状态变化  
                {  
                    res1[0] = 0;  
                    for(int k = 1; k <= len; k++)  
                    {  
                        if(s[k - 1] == dna[j])  
                            res1[k] = res[k - 1] + 1;  
                        else  
                            res1[k] = max(res1[k - 1], res[k]);  
                    }  
                    for(int k = len; k > 0; k--)//记录i状态后添加字符j后的状态  
                    {  
                        ans[i][j] <<= 1;  
                        if(res1[k] > res1[k - 1])  
                            ans[i][j] += 1;  
                    }  
                }  
            }  
            dp[0][0] = 1;  
            for(int i = 0; i < n; i++)  
            {  
                memset(dp[(i + 1) & 1], 0, sizeof dp[(i + 1) & 1]);  
                for(int j = 0; j < maxx; j++)  
                {  
                    for(int k = 0; k < 4; k++)  
                    {  
                        dp[(i + 1) & 1][ans[j][k]] += dp[i & 1][j];  
                        if(dp[(i + 1) & 1][ans[j][k]] > mod)  
                            dp[(i + 1) & 1][ans[j][k]] %= mod;  
                    }  
                }  
            }  
            memset(prt, 0, sizeof prt);  
            for(int i = 0; i < maxx; i++)  
            {  
                int cnt = 0;  
                int tmp = i;  
                while(tmp)  
                {  
                    if(tmp & 1)  
                        cnt++;  
                    tmp >>= 1;  
                }  
                prt[cnt] += dp[n & 1][i];  
                if(prt[cnt] > mod)  
                    prt[cnt] %= mod;  
            }  
            for(int i = 0; i <= len; i++)  
                printf("%lld\n", prt[i]);  
        }  
    return 0;  
}  

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值