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和矩阵
---------------------------------
data矩阵poj 2778 DNA Sequence的更多相关文章
- POJ 2778 DNA Sequence(AC自动机+矩阵加速)
DNA Sequence Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9899 Accepted: 3717 Desc ...
- 线性代数(矩阵乘法):POJ 2778 DNA Sequence
DNA Sequence Description It's well known that DNA Sequence is a sequence only contains A, C, T and ...
- POJ 2778 DNA Sequence (ac自动机+矩阵快速幂)
DNA Sequence Description It's well known that DNA Sequence is a sequence only contains A, C, T and G ...
- poj 2778 DNA Sequence AC自动机DP 矩阵优化
DNA Sequence Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 11860 Accepted: 4527 Des ...
- POJ 2778 DNA Sequence (AC自己主动机 + dp)
DNA Sequence 题意:DNA的序列由ACTG四个字母组成,如今给定m个不可行的序列.问随机构成的长度为n的序列中.有多少种序列是可行的(仅仅要包括一个不可行序列便不可行).个数非常大.对10 ...
- poj 2778 DNA Sequence AC自动机
DNA Sequence Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 11860 Accepted: 4527 Des ...
- poj 2778 DNA Sequence ac自动机+矩阵快速幂
链接:http://poj.org/problem?id=2778 题意:给定不超过10串,每串长度不超过10的灾难基因:问在之后给定的长度不超过2e9的基因长度中不包含灾难基因的基因有多少中? DN ...
- POJ 2778 DNA Sequence(AC自动机+矩阵快速幂)
题目链接:http://poj.org/problem?id=2778 题意:有m种DNA序列是有疾病的,问有多少种长度为n的DNA序列不包含任何一种有疾病的DNA序列.(仅含A,T,C,G四个字符) ...
- POJ 2778 DNA Sequence(AC自动机+矩阵)
[题目链接] http://poj.org/problem?id=2778 [题目大意] 给出一些字符串,求不包含这些字符串的长度为n的字符串的数量 [题解] 我们将所有串插入自动机计算match,对 ...
随机推荐
- Android DES AES MD5加密
AES加密: <span style="font-size:18px;">package com.example.encrypdate.util; import jav ...
- installshield 32位打包和64位打包的注意事项
原文:installshield 32位打包和64位打包的注意事项 32/64位问题要把握几点:1. 明确你的产品是否需要区分32/64位2. 明确你的产品中是否有32/64位的服务注册3. 了解In ...
- linux下C语言中的flock函数使用方法 .
表头文件 #include<sys/file.h> 定义函数 int flock(int fd,int operation); 函数说明 flock()会依參数operation所指 ...
- lua迭代器和仿制药for
不管是什么样的结构,你只需要同意遍历集合可以称为迭代器的所有元素.lua常用来形容叙事功能迭代器.个元素.每个迭代器都须要保存一些状态来知道当前处于什么位置和怎样进行下一次迭代. 对于这种任务.闭包提 ...
- HTTP2协议之HPACK--之头部压缩规范介绍
接下来打算把HTTP2协议的头部压缩算法给翻译下,敬请等候. HPACK - Header Compression for HTTP/2 HPACK:HTTP/2头部压缩 概要说明 这个规范定义了HP ...
- [翻译]如何编写GIMP插件(一)
近期想尝试编写gimp插件,在gimp官网看到了三篇简明教程,顺便翻译了下,由于本人英文,计算机知识有限,文中难免有warning,error出现,欢迎指正. <How to write a G ...
- Katana介绍以及使用
Katana介绍以及使用 接上篇OWIN产生的背景以及简单介绍,在了解了OWIN规范的来龙去脉后,接下来看一下Katana这个OWIN规范的实现,并看看如何使用在我们的Web开发中. 阅读目录: 一. ...
- Tabs - 选项卡插件
接上篇Tabs - 选项卡插件 其中12)Yet (E)Another Tab Interface没有依赖任何javascript框架,以作补充 9)Flipping C ...
- html5学习(一)--canvas画时钟
利用空余时间学习一下html5. <!doctype html> <html> <head></head> <body> <canva ...
- java入门学习(十二)运算语句 if switch
这两天在网上做兼职,耽误了些博客见谅哈 欢迎来我的博客:www.taomaipin.com java中的运算语句而且频繁用到的无法就是条件语句和循环语句,包括if,for,while,switch,b ...