标签：
主要利用状态之间的转换吧，可以把各个状态之间的转换成矩阵，利用矩阵乘法来找出有多少条路径。题解转自：
http://blog.csdn.net/morgan_xww/article/details/7834801
2 1 0 0 1
2 1 1 0 0
1 1 0 1 1
2 1 0 0 1
2 1 0 0 1
M[i,j]表示从结点i到j只走一步有几种走法。
那么M的n次幂就表示从结点i到j走n步有几种走法。
注意：危险结点要去掉，也就是去掉危险结点的行和列。结点3和4是单词结尾所以危险，结点2的fail指针指向4，当匹配”AC”时也就匹配了”C”，所以2也是危险的。
矩阵变成M：
2 1
2 1
计算M[][]的n次幂，然后 Σ(M[0,i]) mod 100000 就是答案。
由于n很大，可以使用二分来计算矩阵的幂
自己写的代码：
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <string.h>
#include <queue>
#include <cmath>
#include <vector>
#define LL __int64
using namespace std;
const int dictsize=4;
const int QN=1000;
const int root=0;
const LL MOD=100000;
int ID[130];
int head,tail;
int que[QN];
struct Node {
int fail,next[dictsize];
bool tag;
void initial(){
fail=-1,tag=false;
for(int i=0;i<dictsize;i++) next[i]=-1;
}
}trie[150];
int tot,n,l;
char str[15];
struct Matrix{
LL mat[105][105];
};
Matrix a,per;
Matrix operator *(Matrix a,Matrix b){
Matrix c;
for(int i=0;i<=tot;i++){
for(int j=0;j<=tot;j++){
c.mat[i][j]=0;
for(int k=0;k<=tot;k++)
c.mat[i][j]=(c.mat[i][j]+a.mat[i][k]*b.mat[k][j])%MOD;
}
}
return c;
}
void Insert_trie(){
int p=0,i=0;
while(str[i]){
if(trie[p].next[ID[str[i]]]==-1) trie[p].next[ID[str[i]]]=++tot;
p=trie[p].next[ID[str[i]]];
i++;
}
trie[p].tag=true;
}
void build_ac(){
que[tail++]=root;
int i,tmp,p;
while(head!=tail){
tmp=que[head++];
p=-1;
for(int i=0;i<dictsize;i++){
if(trie[tmp].next[i]!=-1){
if(tmp==root) trie[trie[tmp].next[i]].fail=root;
else{
p=trie[tmp].fail;
while(p!=-1){
if(trie[p].next[i]!=-1){
trie[trie[tmp].next[i]].fail=trie[p].next[i];
break;
}
p=trie[p].fail;
}
if(p==-1){
trie[trie[tmp].next[i]].fail=root;
}
}
if(trie[trie[trie[tmp].next[i]].fail].tag){
trie[trie[tmp].next[i]].tag=true;
}
que[tail++]=trie[tmp].next[i];
}
else{ //trie[tmp].next[i]==-1
if(tmp==root) trie[tmp].next[i]=root;
else{
p=trie[tmp].fail;
while(p!=-1){
if(trie[p].next[i]!=-1){
trie[tmp].next[i]=trie[p].next[i];
break;
}
p=trie[p].fail;
}
if(p==-1){
trie[tmp].next[i]=root;
}
}
}
}
}
}
void build_map(){
memset(a.mat,0,sizeof(a.mat));
for(int i=0;i<=tot;i++){
for(int j=0;j<dictsize;j++){
if(!trie[i].tag&&!trie[trie[i].next[j]].tag)
a.mat[i][trie[i].next[j]]++;
}
}
}
Matrix multi(int k){
Matrix ans,p=a;
memset(ans.mat,0,sizeof(ans.mat));
for(int i=0;i<=tot;i++)
ans.mat[i][i]=1;
while(k){
if(k&1) ans=ans*p;
k>>=1;
p=p*p;
}
return ans;
}
int main(){
ID[‘A‘]=0,ID[‘G‘]=1,ID[‘C‘]=2,ID[‘T‘]=3;
while(scanf("%d%d",&n,&l)!=EOF){
tot=head=tail=0;
for(int i=0;i<110;i++)
trie[i].initial();
for(int i=0;i<n;i++){
scanf("%s",str);
Insert_trie();
}
build_ac();
build_map();
Matrix ans=multi(l);
LL print=0;
for(int i=0;i<=tot;i++)
print=(print+ans.mat[0][i])%MOD;
printf("%I64d\n",print);
}
return 0;
}
标签：
原文地址：http://www.cnblogs.com/jie-dcai/p/4308813.html
