data矩阵poj 2778 DNA Sequence

最近研究data矩阵，稍微总结一下，以后继续补充：

    ac自动机处理字符串，dp计算谜底，用矩阵来减速

每日一道理
巴尔扎克说过“不幸，是天才的进升阶梯，信徒的洗礼之水，弱者的无底深渊”。风雨过后，眼前会是鸥翔鱼游的天水一色；走出荆棘，前面就是铺满鲜花的康庄大道；登上山顶，脚下便是积翠如云的空蒙山色。 在这个世界上，一星陨落，黯淡不了星空灿烂，一花凋零，荒芜不了整个春天。人生要尽全力度过每一关，不管遇到什么困难不可轻言放弃。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
const long long mod=100000;
struct
{int next[30],tmp,fail,lon;
}trie[111];
int lon;
int m,n;
struct matrix
{long long data[111][111];int size;matrix(int tmp,int lon){size=lon;for(int i=0;i<=size;i++)for(int j=0;j<=size;j++)data[i][j]=tmp;}matrix operator * (const matrix &xx) const{matrix ans(0,size);for(int p=0;p<=size;p++)for(int q=0;q<=size;q++)for(int i=0;i<=size;i++){ans.data[p][q]+=data[p][i]*xx.data[i][q];ans.data[p][q]%=mod;}return(ans);}
};
void trieini()
{memset(trie,0,sizeof(trie));lon=0;
}
void insert(char s[])
{int t=0;int n=strlen(s+1);for(int i=1;i<=n;i++){if(trie[t].next[s[i]-'A']==0){trie[t].next[s[i]-'A']=++lon;trie[lon].lon=i;}t=trie[t].next[s[i]-'A'];if(i==n)trie[t].tmp++;}
}void getfail()
{int root=0;queue <int> q;q.push(root);trie[root].fail=root;while(!q.empty()){int t=q.front();q.pop();for(int i=0;i<26;i++){if(trie[t].next[i]){if(t==root){int u=trie[t].next[i];trie[u].fail=root;q.push(u);continue;}int u=trie[t].next[i];int tmp=trie[t].fail;while(tmp!=root&&trie[tmp].next[i]==0)tmp=trie[tmp].fail;if(trie[tmp].next[i])trie[u].fail=trie[tmp].next[i];elsetrie[u].fail=root;q.push(u);}}}for(int i=1;i<=lon;i++){int t=i;while(t!=root){int u=trie[t].fail;trie[t].tmp+=trie[u].tmp;t=u;}}
}void find(matrix &a)
{int root=0;for(int k=0;k<=lon;k++){if(trie[k].tmp) continue;for(int i=0;i<26;i++){if(i+'A'!='A')if(i+'A'!='C')if(i+'A'!='T')if(i+'A'!='G')continue;int t=k;while(t!=root&&trie[t].next[i]==0)t=trie[t].fail;if(trie[t].next[i]){int u=trie[t].next[i];if(!trie[u].tmp)a.data[k][u]++;}elsea.data[k][0]++;}}
}int main()
{while(scanf("%d %d",&m,&n)!=EOF){trieini();char s[15];for(int i=1;i<=m;i++){scanf("%s",s+1);insert(s);}getfail();matrix a(0,lon);find(a);matrix ans(a);n--;while(n){if(n&1)ans=ans*a;a=a*a;n/=2;}int answer=0;for(int i=0;i<=lon;i++){answer+=ans.data[0][i];answer%=mod;}cout<<answer<<endl;}return 0;
}

文章结束给大家分享下程序员的一些笑话语录： IBM和波音777
　　波音777是有史以来第一架完全在电脑虚拟现实中设计制造的飞机，所用的设备完全由IBM公司所提供。试飞前，波音公司的总裁非常热情的邀请IBM的技术主管去参加试飞，可那位主管却说道:“啊，非常荣幸，可惜那天是我妻子的生日，So..”..
　　波音公司的总载一听就生气了:“胆小鬼，我还没告诉你试飞的日期呢!”

--------------------------------- 原创文章 By
data和矩阵
---------------------------------