程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> C++ >> C++入門知識 >> poj 2778 AC自動機+DP+矩陣快速冪

poj 2778 AC自動機+DP+矩陣快速冪

編輯:C++入門知識

[cpp] 
<pre name="code" class="cpp"><pre name="code" class="cpp">#include<cstdio> 
#include<queue> 
#include<cstring> 
#include<cstring> 
using namespace std; 
typedef __int64 type; 
const int kind=4;  //每個節點的子節點的個數上限 
const int mod = 100000; 
const int size=109; //轉移矩陣的行大小 
class AC_auto 

    public: 
    int tot; 
    type Mar[size][size],ans[size][size]; 
 
    struct Node 
    { 
        int key,cnt; 
        struct Node *next[kind],*fail; 
        Node() 
        { 
            key=0; 
            for(int i=0;i<kind;i++) 
            next[i]=0; 
        } 
    }*root,node[size]; 
 
    inline Node * new_Node() 
    { 
        Node *p=&node[tot];  //靜態方式,可以改為動態方式 new Node() 
        p->cnt=tot; 
        p->key=0;  //表示本節點是不是單詞的結尾 
        memset(p->next,0,sizeof(p->next)); 
        tot++; 
        return p; 
    } 
 
    AC_auto(){  tot=0,root=new_Node(); } 
 
    int Index(char c) 
    { 
        switch(c) 
        { 
            case 'A':return 0; 
            case 'C':return 1; 
            case 'G':return 2; 
            case 'T':return 3; 
        } 
    } 
    void insert(char *s)//構造字典樹 
    { 
        int i=0; 
        Node *p=root; 
        while(s[i]) 
        { 
            if(p->next[Index(s[i])]==0) 
            { 
                p->next[Index(s[i])]=new_Node(); 
            } 
            p=p->next[Index(s[i])]; 
            i++; 
        } 
        p->key++; 
    } 
 
    void build_fail() 
    { 
        root->fail=root; 
        queue<Node *> q; 
        q.push(root); 
        while(!q.empty()) 
        { 
            Node *t=q.front();q.pop(); 
            Node *p=0; 
            for(int i=0;i<kind;i++) 
            { 
                if(t->next[i]!=0) 
                { 
                    if(t==root)  t->next[i]->fail = root; 
                    else 
                    { 
                        p=t->fail; //從父節點的失敗指針開始 
                        while(p!=root&&p->next[i]==0) p=p->fail;//循環失敗指針,一直搜索非根節點有相同索引的節點 
                        if(p->next[i])   t->next[i]->fail=p->next[i]; //找到,失敗指針指向它 
                        else  t->next[i]->fail=root;  //無相同索引,失敗指針指向根節點 
                    } 
                    if(t->next[i]->fail->key!=0)  t->next[i]->key++; 
                    q.push(t->next[i]); 
                } 
                else 
                { 
                    if(t==root) t->next[i]=root;//如果根節點無此索引,失敗指針指向根,表示此狀態可以轉換為根所代表的狀態 
                    t->next[i]=t->fail->next[i]; 
                } 
            } 
        } 
    } 
    void get_Mar() 
    { 
        memset(Mar,0,sizeof(Mar)); 
        for(int i=0;i<tot;i++) 
        if(node[i].key==0) 
        { 
            for(int j=0;j<kind;j++) 
            { 
                if(node[i].next[j]->key==0) 
                { 
                    Mar[i][node[i].next[j]->cnt]++; 
                } 
            } 
        } 
    } 
    void MatrixMultiply(type b[][size], type c[][size], int sz)//矩陣乘 
    { 
        int i,j,k; 
        type temp[size][size] = {0}; 
        for (i=0; i<sz; i++) 
        { 
            for (j=0; j<sz; j++) 
            { 
                for (k=0; k<sz; k++) 
                { 
                    temp[i][j] += b[i][k]*c[k][j]; 
                    if (temp[i][j]>=mod) 
                        temp[i][j] %= mod; 
                } 
            } 
        } 
        for (i=0; i<sz; i++) 
        { 
            for (j=0; j<sz; j++) 
            { 
                b[i][j] = temp[i][j]; 
            } 
        } 
    } 
 
    void MatrixPow(type t[][size], type a[][size], int sz, int n) //矩陣的快速冪 
    { 
        while(n>0) 
        { 
            if (n&1) MatrixMultiply(t, a, sz); 
            MatrixMultiply(a, a, sz); 
            n >>= 1; 
        } 
    } 
    void solve(int n) 
    { 
        memset(ans,0,sizeof(ans)); 
        for(int i=0;i<size;i++) 
        ans[i][i]=1; 
        MatrixPow(ans,Mar,tot,n); 
        type res = 0; 
        for(int i=0; i<tot; i++) 
        { 
            if (node[i].key==0) 
            { 
                res += ans[0][i]; 
                if (res>=mod) 
                    res %= mod; 
            } 
        } 
        printf("%I64d\n",res); 
    } 
 
    int query(char *s) 
    { 
        int i=0,sum=0; 
        Node *p=root; 
        while(s[i]) 
        { 
            while(p->next[Index(s[i])]==0&&p!=root)p=p->fail;//某個非根節點匹配失敗,一直循環失敗指針直到有某個字典樹分支可以匹配 
            p=(p->next[Index(s[i])]==0)?root:p->next[Index(s[i])];//一直到根也未找到滿足條件的分支,則從根開始匹配 
            Node *t=p; 
            while(t!=root&&t->key!=-1)//循環失敗指針,記錄匹配次數 
            { 
                sum+=t->key; 
                t->key=-1; 
                t=t->fail; 
            } 
            i++; 
        } 
        return sum; 
    } 
}; 
char keyword[19],str[1000009]; 
int main() 

    int m,n; 
    while(scanf("%d%d",&m,&n)!=EOF) 
    { 
        AC_auto x; 
        while(m--) 
        { 
            scanf("%s",keyword); 
            x.insert(keyword); 
        } 
        x.build_fail(); 
        x.get_Mar(); 
        x.solve(n); 
    } 


作者:wsniyufang

  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved