【AC自动机+矩阵快速幂】poj2778 DNA Sequence

原创 2017年08月13日 13:38:00

渐渐开始变难orz
先老老实实建树;
然后建ac自动机,注意两个:
1. 若ch[I][j] ==0 ,某个节点没有某一字符的儿子,考虑将其指向其失配节点的儿子:ch[I][j] = ch[f[I]][j];Tire数就变成Tire图了。

2.若val[f[I]]==1,某节点失配节点是危险节点,则将此节点标为危险节点。

然后构造矩阵mp,如果a,b节点相连且b不是危险节点,mp[a][b]++;

然后把矩阵mp扔去快速幂,求mp[0][i](I>=0&&I<=sz )的和,就是答案;

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
const int maxn = 101;
const int mod = 100000;

int n, l;

struct Matrix{
    ll mat[101][101];
    int r, c;
    Matrix(int x){ memset(mat, 0, sizeof(mat)); r = x; c = x; }

    Matrix operator*(const Matrix &b) const{
        Matrix res(b.r);
        for (int i = 0; i <= b.r; i++){
            for (int j = 0; j <= b.r; j++){
                for (int k = 0; k <= b.r; k++){
                    res.mat[i][j] = (res.mat[i][j] + mat[i][k] * b.mat[k][j]) % mod;
                }
            }
        }
        return res;
    }
};

Matrix Qpow(Matrix x, ll n){
    Matrix res(x.c), tmp = x;
    for (int i = 0; i <= x.c; i++) res.mat[i][i] = 1;

    while (n){
        if (n & 1) res = res *tmp;
        n >>= 1;
        tmp = tmp*tmp;
    }

    return res;
}

struct Trie{

    int ch[maxn][26];
    int val[maxn];
    int last[maxn];
    int f[maxn];
    int sz;
    int idx(char x){
        if (x == 'A') return 0;
        if (x == 'T') return 1;
        if (x == 'C') return 2;
        if (x == 'G') return 3;
    }

    void insert(char *s){
        int u = 0, n = strlen(s);
        for (int i = 0; i < n; i++){
            int c = idx(s[i]);
            if (!ch[u][c]){
                ch[u][c] = sz++;
            }
            u = ch[u][c];
        }
        val[u] = 1;
    }

    void init(){
        memset(ch, 0, sizeof(ch));
        memset(val, 0, sizeof(val));
        sz = 1;
        for (int i = 0; i < n; i++){
            char str[55]; scanf("%s", str);
            insert(str);
        }
    }

    void build(){
        queue<int> q;
        memset(f, 0, sizeof(f));
        memset(last, 0, sizeof(last));
        for (int i = 0; i < 26; i++){
            int u = ch[0][i];
            if (u){ q.push(u); }
        }
        while (!q.empty()){
            int r = q.front(); q.pop();
            for (int i = 0; i < 26; i++){
                int u = ch[r][i];
                if (!u) {
                    ch[r][i] = ch[f[r]][i];
                    continue;
                };
                q.push(u);

                int v = f[r];
                while (v && !ch[v][i]) v = f[v];
                f[u] = ch[v][i];
                last[u] = val[f[u]] ? f[u] : last[f[u]];
                if (val[f[u]] == 1) val[u] = 1;
            }
        }
    }

    Matrix getM(){
        Matrix mp(sz);
        for (int i = 0; i < sz; i++){
            for (int j = 0; j < 4; j++){
                if (!val[i] && !val[ch[i][j]])
                    mp.mat[i][ch[i][j]]++;
            }
        }
        return mp;
    }
};

Trie ac;
int main(){
    while (~scanf("%d%d", &n, &l)){
        ac.init();
        ac.build();
        Matrix m = ac.getM();
        m = Qpow(m, l);
        ll ans = 0;
        for (int i = 0; i < ac.sz; i++){
            ans = (ans + m.mat[0][i]) % mod;
        }
        cout << ans << endl;
    }
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

poj 2778 DNA Sequence 【ac自动机 + dp + 矩阵快速幂】

这题基本上是纯模版,开始wa了,然后和模版对照了一下,居然是建fail指针的时候搞错了。。。少写了一个q.push(to) #include #include #include #include...
  • zz_1215
  • zz_1215
  • 2012年05月17日 11:01
  • 704

POJ2778 DNA Sequence AC自动机+矩阵快速幂

题目描述:给出m(m        这道题目深层次地利用了AC自动机,但是并没有用到AC自动机中的find(query)函数,而是用到了BFS序建立失败指针走向,以及为字典树中的空节点加上...
  • jijijix
  • jijijix
  • 2017年02月17日 15:32
  • 310

poj2778 DNA Sequence(AC自动机+矩阵快速幂 )

大神附图的题解:http://blog.csdn.net/morgan_xww/article/details/7834801 【题解】 将所有病毒串建立成字典树,并标记词尾结点,以下称"非法结点" ...
  • cjk_cjk
  • cjk_cjk
  • 2015年06月24日 12:44
  • 301

POJ2778 DNA Sequence AC自动机+DP+矩阵快速幂

Problem Address:http://poj.org/problem?id=2778【前言】这道题算是出自Trie主题的。算是一道比较综合的题目。【思路】(1)首先用所给的字符串构造出一个AC...

POJ 2778 DNA Sequence (AC自动机 + 矩阵快速幂)

题目链接:DNA Sequence 解析:AC自动机 + 矩阵加速(快速幂)。 这个时候AC自动机 的一种状态转移图的思路就很透彻了,AC自动机就是可以确定状态的转移。 ...

poj2778 DNA Sequence AC自动机+矩阵快速幂

DNA Sequence Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 11721   Accepted: 44...
  • corncsd
  • corncsd
  • 2014年09月11日 20:47
  • 272

POJ 2778 DNA Sequence (AC自动机+矩阵快速幂)

题目大意有m个字符串,求所有长度为n的字符串中,不包含这m个串的个数。这些串全部以A,C,G,T组成。分析 把m个串构造成AC自动机 然后用M[i][j]表示结点i到结点j只走1步,可行的方法数。那么...

[POJ 2778] DNA Sequence (AC自动机+trie图+矩阵快速幂)

链接POJ 2778题意给出一些关键词,问长度为n的不含有这些关键词中任何一个的字符串有多少个。 关键词和字符串均由’A’、’C’、’T’、’G’四个字符构成。题解一道关于trie图的题目。首先补充...

【POJ2778】【AC自动机】【矩阵快速幂】17.2.9 T2 DNA Sequence 题解

DNA Sequence 基因序列可以看做是只含大写ACGT 的字符串。 现有n 个癌症序列,每个癌症序列有至多10 个字符。 现请问不含癌症序列的长度为len 的基因有多少个?答案对1000...

E - DNA Sequence POJ - 2778 (AC自动机+图论小知识+矩阵快速幂)

这道题目很有趣,我觉得很好玩。虽然不会做,看了题解才懂得。 刚刚开始看这道题目,以为是dp,但是所求的字符长度n的大小有20亿, 就直接放弃了。然后看了题解, 看了半天才看懂。 首先构造AC自动机...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【AC自动机+矩阵快速幂】poj2778 DNA Sequence
举报原因:
原因补充:

(最多只允许输入30个字)