#include#include#include#include#include#include#include#include#include#include#defineLLlongl" />

亚洲免费在线-亚洲免费在线播放-亚洲免费在线观看-亚洲免费在线观看视频-亚洲免费在线看-亚洲免费在线视频

poj 2778 DNA Sequence

系統 1915 0

http://poj.org/problem?id=2778

AC自動機 + 矩陣連乘

代碼:

      #include<iostream>

#include<cmath>

#include<cstdio>

#include<string>

#include<cstring>

#include<vector>

#include<stack>

#include<queue>

#include<set>

#include<map>

#include<algorithm>



#define LL long long



using namespace std;



const int INF=0x3f3f3f3f;

const int MOD=100000;

const LL LMOD=100000;

const int N=1005;

const int M=1005;

const int K=4;

const int MATRIXN=5;

const int MATRIXM=115;

const int MATRIXH=115;

struct nodeTrie

{

    int v;

    int fail;

    int next[K];

    void initialize()

    {

        v=0;

        fail=-1;

        memset(next,-1,sizeof(next));

    }

}trie[M];

int cnt,root;

char s[N];

LL ma[MATRIXN][MATRIXM],mb[MATRIXM][MATRIXH],mc[MATRIXM][MATRIXH];

int getNewNode()

{

    ++cnt;

    trie[cnt].initialize();

    return cnt;

}

void addWord(int p,char *s,int k)

{

    if(s[0]=='\0') return ;

    for(int i=0;s[i]!='\0';++i)

    {

        if(trie[p].next[s[i]-'A']==-1)

        trie[p].next[s[i]-'A']=getNewNode();

        p=trie[p].next[s[i]-'A'];

    }

    (trie[p].v)=k;

}

void init(int n)

{

    cnt=0;

    root=getNewNode();

    for(int i=0;i<4;++i)

    {

        s[0]='A'+i;s[1]='\0';

        addWord(root,s,0);

    }

    while(n--)

    {

        gets(s);

        for(int i=0;s[i]!='\0';++i)

        {

            if(s[i]=='C')

            s[i]='B';

            else if(s[i]=='G')

            s[i]='C';

            else if(s[i]=='T')

            s[i]='D';

        }

        addWord(root,s,1);

    }

}

void bfs(int p)

{

    trie[p].fail=root;

    queue<int>qt;

    qt.push(p);

    while(!qt.empty())

    {

        int y;

        int x=qt.front();qt.pop();

        for(int i=0;i<K;++i)

        if(trie[x].next[i]!=-1)

        {

            qt.push(trie[x].next[i]);

            if(x==root)

            {trie[trie[x].next[i]].fail=root;continue;}

            y=trie[x].fail;

            while(y!=root&&trie[y].next[i]==-1)

            y=trie[y].fail;

            if(trie[y].next[i]!=-1)

            trie[trie[x].next[i]].fail=trie[y].next[i];

            else

            trie[trie[x].next[i]].fail=root;

        }

    }

}

void makeMatrix(LL a[][MATRIXM],int n)

{

    for(int i=1;i<=n;++i)

    for(int j=1;j<=n;++j)

    a[i][j]=0;

    for(int i=1;i<=n;++i)

    for(int j=0;j<4;++j)

    {//cout<<i<<endl;

        int p=i;

        int q=0;

        while(p!=root&&trie[p].next[j]==-1)

        p=trie[p].fail;

        if(trie[p].next[j]!=-1)

        q=trie[p].next[j];

        int w=q;

        bool flag=true;

        while(w!=root)

        {

            if(trie[w].v>0)

            {flag=false;break;}

            w=trie[w].fail;

        }

        if(trie[w].v>0)

        flag=false;

        if(flag==true)

        a[i][q]=1;

    }

}

void matrixMul(LL a[][MATRIXM],LL b[][MATRIXM],int n,int m,int h)

{

    for(int i=1;i<=n;++i)

    for(int j=1;j<=h;++j)

    {

        LL tmp=0;

        for(int l=1;l<=m;++l)

        tmp=(tmp+a[i][l]*b[l][j])%LMOD;

        mc[i][j]=tmp;

    }

    for(int i=1;i<=n;++i)

    for(int j=1;j<=h;++j)

    a[i][j]=mc[i][j];

}

int main()

{

    //freopen("data.in","r",stdin);

    int n,m;

    while(scanf("%d %d ",&n,&m)!=EOF)

    {

        init(n);

        bfs(root);

        int MN=1,MM=cnt,MH=cnt;

        makeMatrix(mb,MM);

        memset(ma,0,sizeof(ma));

        ma[1][1]=1;

        while(m>0)

        {

            if((m&1)==1)

            matrixMul(ma,mb,MN,MM,MH);

            matrixMul(mb,mb,MM,MM,MH);

            m=(m>>1);

        }

        LL ans=0;

        for(int i=1;i<=MN;++i)

        for(int j=1;j<=MM;++j)

        ans=(ans+ma[i][j])%LMOD;

        cout<<ans<<endl;

    }

    return 0;

}


    

poj 2778 DNA Sequence


更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 久久只有精品视频 | 久青草视频在线观看 | 久久99热久久国产精品 | 曰批免费视频播放在线看片二 | 亚洲人成在线精品不卡网 | 91论坛在线 | 奇米在线播放 | 亚洲综合中文网 | 精品国产你懂的在线观看 | 精品99久久 | 在线观看不卡视频 | 亚洲综合在线另类色区奇米 | 日韩在线a视频免费播放 | 男女一级做片a性视频 | 日韩中文字幕在线有码视频网 | 国产激情对白一区二区三区四 | 国产情精品嫩草影院88av | 黄网站www | 亚洲日本aⅴ片在线观看香蕉 | 91大学生视频 | 日韩精品一区二区三区乱码 | 国产一二视频 | 狠狠色噜噜狠狠狠米奇9999 | 99成人国产精品视频 | 亚洲综合网在线观看 | 国产成年女一区二区三区 | 起碰成人免费公开网视频 | 日本精品久久久中文字幕 | 中文字幕一区视频 | 欧美日韩激情在线一区 | 性感美女香蕉视频 | 一级毛片一级毛片a毛片欧美 | 热久久99精品这里有精品 | 奇米第| 亚洲综合成人网在线观看 | 奇米影视第四色7777 | 特级一级毛片 | 综合爱 | 四虎永久在线 | 国产在线拍揄自揄视频不卡99 | 国产精品原创 |