HDU 4899 Hero meet devil (状压DP, DP预处理）

题意：给你一个基因序列s（只有A,T,C,G四个字符，假设长度为n），问长度为m的基因序列s1中与给定的基因序列LCS是0，1......n的有多少个？

思路：最直接的方法是暴力枚举长度为m的串，然后再用求LCS的dp。当然我们可以在枚举的时候同时进行dp,但是复杂的仍然为O(4 ^ m)。我们可以观察求LCS 的状态转移方程：dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]) 若s[i] == s1[j] dp[i][j] = max(dp[i - 1][j - 1] + 1)。可以发现，每一行的相邻的状态最多只会差1，那么我们可以用差分的方法转化为状压dp。剩下的部分这两篇博客讲的很清楚了：https://www.cnblogs.com/RabbitHu/p/BZOJ3864.html， https://www.cnblogs.com/owenyu/p/6724616.html

代码：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const LL mod = 1000000007;
const int maxn = 1010;
char mp[4] = {'A', 'T', 'C', 'G'};
int ans[20];
int dp[2][1 << 15];
int trans[1 << 15][4], cnt[1 << 15];
char s[maxn];
void init(int n) {int pre[20], cur[20];memset(pre, 0, sizeof(pre));memset(cur, 0, sizeof(cur));for (int i = 0; i < (1 << n); i++) {if(i)cnt[i] = cnt[i >> 1] + (i & 1);pre[0] = i & 1;for (int j = 1; j < n; j++)pre[j] = pre[j - 1] + (i >> j & 1); for (int k = 0; k < 4; k++) {int now = 0;cur[0] = pre[0];if(mp[k] == s[0]) cur[0] = 1;now |= cur[0];for (int j = 1; j < n; j++) {cur[j] = max(cur[j - 1], pre[j]);if(mp[k] == s[j]) {cur[j] = max(cur[j], pre[j - 1] + 1);}now |= ((cur[j] - cur[j - 1]) << j);}trans[i][k] = now;}}
}
int main() {int T, m, n;scanf("%d", &T);while(T--) {scanf("%s", s);scanf("%d", &m);int n = strlen(s);init(n);memset(dp, 0, sizeof(dp));memset(ans, 0, sizeof(ans));dp[0][0] = 1;    for (int i = 1; i <= m; i++) {memset(dp[i & 1], 0, sizeof(dp[i & 1]));int pre = (i & 1) ^ 1;int now = i & 1;for (int j = 0; j < (1 << n); j++) {for (int k = 0; k < 4; k++) {dp[now][trans[j][k]] = (dp[pre][j] + dp[now][trans[j][k]]) % mod;}}}  for (int i = 0; i < (1 << n); i++)ans[cnt[i]] = (ans[cnt[i]] + dp[m & 1][i]) % mod;for(int i = 0; i <= n; i++)printf("%d\n", ans[i]); }
}

转载于:https://www.cnblogs.com/pkgunboat/p/10440763.html

HDU 4899 Hero meet devil (状压DP, DP预处理）相关推荐

[BZOJ 3864][HDU 4899]Hero meet devil(DP套DP)
题意给你一个只由AGCT组成的字符串S(|S|≤15),对于每个1≤i≤|S| 询问有多少个只由AGCT组成的长度为m(1≤m≤1000)的字符串T,使得LCS(S,T)=i. |S|<=15 ...
HDU 4899 Hero meet devil
题意:一段DNA序列s,只包含字符ATGC,长度不超过15,求有多少种长度为n的DNA序列与s的最长公共子序列长度为0~len. 解法:状压dp. LSC: if(a[i] == b[j]) dp[i ...
HDU 4336 Card Collector（状压 + 概率DP 期望）题解
题意:每包干脆面可能开出卡或者什么都没有,一共n种卡,每种卡每包爆率pi,问收齐n种卡的期望思路:期望求解公式为:$E(x) = \sum_{i=1}^{k}pi * xi + (1 - \sum_ ...
bzoj 3864: Hero meet devil [dp套dp]
3864: Hero meet devil 题意: 给你一个只由AGCT组成的字符串S (|S| ≤ 15),对于每个0 ≤ .. ≤ |S|,问有多少个只由AGCT组成的长度为m(1 ≤ m ≤ ...
#3864. Hero meet devil dp套dp + 状压 + 状态机
传送门文章目录题意: 思路: 题意: 给你一个只包含ACGTACGTACGT的串sss,再给你一个mmm,第iii行输出有多少个长度为mmm且只包含ACGTACGTACGT的串与sss的lcslc ...
HDU 5691 Sitting in Line 状压dp
Sitting in Line 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5691 Description 度度熊是他同时代中最伟大的数学家,一切 ...
BZOJ 3864: Hero meet devil (从dp性质实现dp套dp)
题意:求长度为m的,字符集大小为4的,字符串,中,与字符串S(|S|<=15)的最长公共子序列长度=i的字符串数量.i∈0→∣S∣i \ \in 0 \to |S|i ∈0→∣S∣ 发现这个状态 ...
HDU 3001 Travelling (三进制状压dp)
题意 n(n<=10)n(n个城市,知道每个城市间的旅行费用,但每个城市最多走两遍.问最小花费是多少 . 也就是每个城市可以走两次的tsp问题. 分析最多走两次,三进制0 1 2可满足,即用三 ...
bzoj4565 [HAOI2016]字符合并结论+状压+区间dp
如果k==2的话就是记搜, 但这个题用记搜的思路的话是需要枚举k个断点的, 所以对于枚举断点,就很可能有优化,比如到一个断点,一个决策的最优值于是考虑区间dp模型,相当于是插入一个数,然后看影响. ...

HDU 4899 Hero meet devil (状压DP, DP预处理）

HDU 4899 Hero meet devil (状压DP, DP预处理）相关推荐

最新文章

热门文章