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];
else
trie[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]++;
}
else
a.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和矩阵
---------------------------------

上一篇:nios II--实验1——hello_world硬件部分


下一篇:淘宝天猫关键词SEO优化