博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
POJ 2778 DNA Sequence (AC自己主动机 + dp)
阅读量:5825 次
发布时间:2019-06-18

本文共 2951 字,大约阅读时间需要 9 分钟。

DNA Sequence

题意:DNA的序列由ACTG四个字母组成,如今给定m个不可行的序列。问随机构成的长度为n的序列中。有多少种序列是可行的(仅仅要包括一个不可行序列便不可行)。个数非常大。对100000取模。
思路:推荐一个 ,讲的很清楚。

这样的题目。n非常大,首先想到的就是用矩阵来优化。那么怎样构造转移方程呢:首先建立一棵Trie,然后依照AC自己主动机的方式构造fail指针,然后会发现。当一个状态分别加入ACTG之后,会得到还有一个状态。

(详细解释见代码)

代码:
/*ID: wuqi9395@126.comPROG:LANG: C++*/#include#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define INF (1<<30)#define PI acos(-1.0)#define mem(a, b) memset(a, b, sizeof(a))#define rep(i, n) for (int i = 0; i < n; i++)#define debug puts("===============")typedef long long ll;using namespace std;const int maxn = 110;const int maxm = 110;ll mod = 100000;struct Matrix { int n, m; ll a[maxn][maxm]; void clear() { n = m = 0; memset(a, 0, sizeof(a)); } Matrix operator * (const Matrix &b) const { //实现矩阵乘法 Matrix tmp; tmp.clear(); tmp.n = n; tmp.m = b.m; for (int i = 0; i < n; i++) for (int j = 0; j < m; j++) { if (!a[i][j]) continue; for (int k = 0; k < b.m; k++) tmp.a[i][k] += a[i][j] * b.a[j][k], tmp.a[i][k] %= mod; } return tmp; }}A, res;const int maxnode = 11 * 11;const int charset = 4;struct ACAutomaton { int ch[maxnode][charset]; int fail[maxnode]; int Q[maxnode]; int val[maxnode]; int sz; int id(char ch) { if (ch == 'A') return 0; else if (ch == 'C') return 1; else if (ch == 'T') return 2; return 3; } void init() { fail[0] = 0; //for (int i = 0; i < charset; i++) ID[i] = i; } void reset() { sz = 1; memset(ch[0], 0, sizeof(ch[0])); } void Insert(char* s, int key) { int u = 0; for (; *s; s++) { int c = id(*s); if (!ch[u][c]) { memset(ch[sz], 0, sizeof(ch[sz])); val[sz] = 0; ch[u][c] = sz++; } u = ch[u][c]; } val[u] = key; } void Construct () { int *s = Q, *e = Q; for (int i = 0; i < charset; i++) { if (ch[0][i]) { *e++ = ch[0][i]; fail[ch[0][i]] = 0; } } while(s != e) { int u = *s++; if (val[fail[u]]) val[u] = 1; for (int i = 0; i < charset; i++) { int &v = ch[u][i]; if (v) { *e++ = v; fail[v] = ch[fail[u]][i]; } else { v = ch[fail[u]][i]; } } } } /* dp[i][j]表示长度为i。后缀为j的状态 最多就仅仅有10*10个后缀 所以可以通过dp[n][j] = a0 * dp[n-1][0] + ... + ak * dp[n - 1][k]得到状态转移的矩阵 */ void work() { for (int i = 0; i < sz; i++) { for (int j = 0; j < charset; j++) { //对于i状态,通过加入ACTG可以得到新的状态(且之前已经构造过AC自己主动机,ch[i][j]便表示新状态) if (!val[i] && !val[ch[i][j]]) { //两个状态都必须是可行的,转化才有意义 A.a[i][ch[i][j]]++; } } } }} AC;Matrix Matrix_pow(Matrix A, ll k, ll mod) { res.clear(); res.n = res.m = AC.sz; for (int i = 0; i < AC.sz; i++) res.a[i][i] = 1; while(k) { if (k & 1) res = res * A; A = A * A; k >>= 1; } return res;}int main () { int m, n; A.clear(); AC.init(); AC.reset(); char str[15]; scanf("%d%d", &m, &n); for (int i = 0; i < m; i++) { scanf("%s", str); AC.Insert(str, 1); } A.n = A.m = AC.sz; AC.Construct(); //之前的都是AC自己主动机构造部分 AC.work(); //得到状态转移的矩阵 res = Matrix_pow(A, n, mod); int ans = 0; rep(i, AC.sz) ans += res.a[0][i]; printf("%d\n", ans % mod); return 0;}

转载地址:http://aqidx.baihongyu.com/

你可能感兴趣的文章
iptables 配置需要保存
查看>>
.NET各种小问题
查看>>
ApkTool反编译和重新打包
查看>>
OpenState: Programming Platform-independent Stateful OpenFlow Applications Inside the Switch
查看>>
java中的Volatile关键字
查看>>
前端自定义图标
查看>>
sqlserver 取取月初月末和月份间隔
查看>>
Vagrant的一个BUG - 不支持'change_host_name'
查看>>
实验二
查看>>
MongoDB数据库迁移
查看>>
独立开发一个云(PaaS)的核心要素, Go, Go, Go!!!
查看>>
java的继承性
查看>>
tomcat 实例
查看>>
MyBatis使用DEMO及cache的使用心得
查看>>
网站文章如何能自动判定是抄袭?一种算法和实践架构剖析
查看>>
【OpenCV学习】滚动条
查看>>
ofo用科技引领行业进入4.0时代 用户粘性连续8个月远甩摩拜
查看>>
无法拒绝|华为618最高优惠1000元 更有梅西签名球衣奉送
查看>>
乐信Q2季报图解:调整后净利过5亿 同比增长776%
查看>>
兰州青年志愿者“中西合璧”玩快闪 温暖旅客回家路
查看>>