题目大意

给定一个由AGCT组成的串\(t\), 求对于所有的\(L \in [1, |t|]\), 有多少个由AGCT组成的串\(s\)满足\(LCS(s, t) = L\).

Solution

传说中的DP套DP.

我们用\(f_{i, j}\)表示\(s\)的前\(i\)位与\(t\)的前\(j\)位的最长公共子序列, 则我们有

\[f_{i, j} = \max \begin{cases} f_{i - 1, j} \\ f_{i, j - 1} \\ f_{i - 1, j - 1} + [s_i = t_j] \end{cases}
\]

则\(LCS(s, t) = f[|S|][|T|]\).

逐位考虑\(s\), 假设当前到\(s\)的第\(i\)位, 我们用\(F\)表示状态: \(F = \{ f[i][1], f[i][2], \cdots, f[i][|t|] \}\). 考虑当\(i\)变成\(i + 1\)时, \(F\)会怎么变化: 对于确定的\(t\), \(F\)的变化只与\(s[i + 1]\)有关, 因此我们令\(c =s[i + 1]\), 用\(T(F, c)\)表示当\(s[i + 1] = c\)时\(F\)会变成怎么样.

我们用\(g[i][F]\)表示\(s\)的前\(i\)位与\(t\)的最长公共子序列状态为\(F\)的串\(s\)的数量, 则对于每一个\(c\), 有\(g[i + 1][T(F, c)] += g[i][F]\).

我们发现\(F\)不容易被记录, 又因为注意到\(0 \le F[i] - F[i - 1] \le 1\)且\(F[0] = 0\), 因此我们用\(S\)来表示\(\{ F[i] - F[i - 1] \}\)的状态集合. 这样我们就可以轻易地进行状态压缩了.

考虑怎么统计答案:

\[ans[k] = \sum_{F} g[|T|][F] \times [F[|T|] = k]
\]

#include <cstdio>
#include <cstring>
#include <algorithm> using namespace std;
const int N = 15, MOD = (int)1e9 + 7;
int t[N + 7];
int n, m;
int main()
{ #ifndef ONLINE_JUDGE freopen("gene.in", "r", stdin);
freopen("gene.out", "w", stdout); #endif int T; scanf("%d\n", &T);
for (int cs = 0; cs < T; ++ cs)
{
static char str[N + 7]; scanf("%s", str + 1);
n = strlen(str + 1); scanf("%d\n", &m);
for (int i = 1; i <= n; ++ i)
if (str[i] == 'A') t[i] = 1;
else if (str[i] == 'T') t[i] = 2;
else if (str[i] == 'G') t[i] = 3;
else if (str[i] == 'C') t[i] = 4;
static int trans[(1 << N) + 7][7];
for (int j = 0; j < 1 << n; ++ j) for (int c = 1; c <= 4; ++ c)
{
static int a[N + 7], b[N + 7];
memset(a, 0, sizeof a);
for (int k = 0; k < n; ++ k) a[k + 1] = a[k] + (j >> k & 1);
b[0] = 0;
for (int k = 1; k <= n; ++ k) b[k] = max(max(a[k], b[k - 1]), a[k - 1] + (c == t[k]));
int stt = 0;
for (int k = 0; k < n; ++ k) if (b[k + 1] - b[k]) stt |= 1 << k;
trans[j][c] = stt;
}
static int f[(1 << N) + 7], g[(1 << N) + 7];
memset(f, 0, sizeof f); f[0] = 1;
for (int i = 0; i < m; ++ i)
{
memset(g, 0, sizeof g);
for (int j = 0; j < 1 << n; ++ j) if (f[j]) for (int c = 1; c <= 4; ++ c)
g[trans[j][c]] = (g[trans[j][c]] + f[j]) % MOD;
swap(f, g);
}
static int ans[N + 7];
memset(ans, 0, sizeof ans);
for (int i = 0; i < 1 << n; ++ i)
{
int cnt = 0;
for (int tmp = i; tmp; tmp >>= 1) if (tmp & 1) ++ cnt;
ans[cnt] = (ans[cnt] + f[i]) % MOD;
}
for (int i = 0; i <= n; ++ i) printf("%d\n", ans[i]);
}
}

BZOJ 3864 Hero Meets Devil的更多相关文章

  1. bzoj 3864: Hero meet devil [dp套dp]

    3864: Hero meet devil 题意: 给你一个只由AGCT组成的字符串S (|S| ≤ 15),对于每个0 ≤ .. ≤ |S|,问 有多少个只由AGCT组成的长度为m(1 ≤ m ≤ ...

  2. BZOJ 3864 Hero meet devil 超详细超好懂题解

    题目链接 BZOJ 3864 题意简述 设字符集为ATCG,给出一个长为\(n(n \le 15)\)的字符串\(A\),问有多少长度为\(m(m \le 1000)\)的字符串\(B\)与\(A\) ...

  3. bzoj 3864: Hero meet devil

    bzoj3864次元联通们 第一次写dp of dp (:з」∠) 不能再颓废啦 考虑最长匹配序列匹配书转移 由于dp[i][j]的转移可由上一行dp[i-1][j-1],dp[i-1][j],dp[ ...

  4. bzoj 3864: Hero meet devil(dp套dp)

    题面 给你一个只由\(AGCT\)组成的字符串\(S (|S| ≤ 15)\),对于每个\(0 ≤ .. ≤ |S|\),问 有多少个只由\(AGCT\)组成的长度为\(m(1 ≤ m ≤ 1000) ...

  5. BZOJ 3864 Hero meet devil (状压DP)

    最近写状压写的有点多,什么LIS,LCSLIS,LCSLIS,LCS全都用状压写了-这道题就是一道状压LCSLCSLCS 题意 给出一个长度为n(n<=15)n(n<=15)n(n< ...

  6. 【BZOJ3864】Hero meet devil DP套DP

    [BZOJ3864]Hero meet devil Description There is an old country and the king fell in love with a devil ...

  7. bzoj千题计划241:bzoj3864: Hero meet devil

    http://www.lydsy.com/JudgeOnline/problem.php?id=3864 题意: 给你一个DNA序列,求有多少个长度为m的DNA序列和给定序列的LCS为0,1,2... ...

  8. HDU 4899 Hero meet devil(状压DP)(2014 Multi-University Training Contest 4)

    Problem Description There is an old country and the king fell in love with a devil. The devil always ...

  9. bzoj3864: Hero meet devil

    Description There is an old country and the king fell in love with a devil. The devil always asks th ...

随机推荐

  1. (原、整)BSP的江湖传说

    @author:黑袍小道 查看随缘,当苦无妨,良人可归.     引言 为什么叫江湖传说,因为实现了第一人是卡马克,就这么简单.(不接受那啥) Quake3:http://www.mralligato ...

  2. pin

    sjhh@123456 Michael zhang zhangxiaocong_2011@yeah.net

  3. HDU 4714 Tree2cycle 找规律

    假设最少删除的边的个数为cost,显然,最终答案即为cost+cost+1 (因为删除一条边,就会增加一个链,所以删除cost条边后,就会有cost+1条链,将这cost+1条链连接起来的代价为cos ...

  4. Android记事本08

    昨天: Anr问题异常的原因和解决方案. 今天: Activity数据传递之通用方式. Activity数据传递之静态变量. Activity数据传递之全局变量. 遇到的问题: 无.

  5. 软工实践 - 第十六次作业 Alpha 冲刺 (7/10)

    队名:起床一起肝活队 组长博客:https://www.cnblogs.com/dawnduck/p/10013959.html 作业博客:班级博客本次作业的链接 组员情况 组员1(队长):白晨曦 过 ...

  6. user namespace 2

    昨天又试了一下, echo 0 1000 1 > /proc/<PID>/uid_map 想着直接把一个进程的uid设置成root,这样他就可以在自己的user namespace ...

  7. Scala 基础(2)—— 基本数据结构

    1. Scala 的面向对象 在学习 Java 的时候,我们说 Java 是一门面向对象的语言,然而 Java 其实并没有完全遵守“一切皆对象”这一准则. 例如:Java 的8种基本数据类型 & ...

  8. Codeforces 585D Lizard Era: Beginning | 折半搜索

    参考这个博客 #include<cstdio> #include<algorithm> #include<cstring> #include<map> ...

  9. Delivering Goods UVALive - 7986(最短路+最小路径覆盖)

    Delivering Goods UVALive - 7986(最短路+最小路径覆盖) 题意: 给一张n个点m条边的有向带权图,给出C个关键点,问沿着最短路径走,从0最少需要出发多少次才能能覆盖这些关 ...

  10. CF911F Tree Destruction 解题报告

    CF911F Tree Destruction 题意翻译 给你一棵树,每次挑选这棵树的两个叶子,加上他们之间的边数(距离),然后将其中一个点去掉,问你边数(距离)之和最大可以是多少. 输入输出格式 输 ...