【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...

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

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

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

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

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

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

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

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

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

DNA Sequence Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 6232   A...

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

先构建AC自动机,然后通过ACzidongji

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

题目链接:点击打开链接 题意描述:给m字符串,这m个字符串是带有病毒的DNA。然后问一个长度为n的字符串不带有任何病毒有多少中可能?所有字符串之后ACGT这几个字符串组成 解题思路:AC自动机+矩...

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

DNA Sequence Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 16494   Accepted: 63...

100道动态规划——28 POJ 2778 DNA Sequence AC自动机+DP+矩阵快速幂 很不错嘛

一眼看出是AC自动机上的DP,因为数据量的问题显然不能记忆化搜索或者是递推,然后就不知道怎么做了         天知道还要联系到矩阵快速幂,图论的东西好久都没写过了=_=         定义一...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

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