[日常摸鱼]poj2778 DNA Sequence
2024-09-08 02:06:42
这题太神啦
题意:求长度为$n$的不包含给定DNA序列的DNA序列个数,给定的不超过10个
构建出Trie图,用$danger[i]$来表示不能走到$i$,对于DNA序列结尾的结点$danger$设为1,构建$fail$指针的时候对于一个结点$i$的某个后缀如果$danger$为1那么$danger[i]$也应该为1.
然后根据Trie图再构造出对应的转移矩阵,自乘$n$次统计答案.
为什么自乘$n$次就行了?[感性理解.jpg]
人傻自带大常数?
hhh
#include<cstdio>
#include<cstring>
typedef long long lint;
const lint MOD=100000;
const int N=105;
struct matrix
{
lint w[N][N];
matrix(){memset(w,0,sizeof(w));}
};
int n,m,cnt,head,tail,ans;
int tr[N][5],fail[N],q[N],map[300];
bool danger[N];char s[N];
inline void insert(char *c)
{
int len=strlen(c+1),k=0;
for(register int i=1;i<=len;i++)
{
int t=map[(int)c[i]];
if(!tr[k][t])tr[k][t]=++cnt;
k=tr[k][t];
}
danger[k]=1;
}
inline void build()
{
for(register int i=0;i<4;i++)if(tr[0][i])q[tail++]=tr[0][i],fail[tr[0][i]]=0;
while(head<tail)
{
int k=q[head++];
for(register int i=0;i<4;i++)
{
if(!tr[k][i])tr[k][i]=tr[fail[k]][i];
else
{
fail[tr[k][i]]=tr[fail[k]][i];
danger[tr[k][i]]|=danger[fail[tr[k][i]]];
q[tail++]=tr[k][i];
}
}
}
}
inline matrix mul(matrix a,matrix b)
{
matrix res;
for(register int i=0;i<=cnt;i++)
for(register int k=0;k<=cnt;k++)if(a.w[i][k])
for(register int j=0;j<=cnt;j++)
{
lint temp=(lint)(a.w[i][k]*b.w[k][j]);
res.w[i][j]=(res.w[i][j]+temp);
if(res.w[i][j]>MOD)res.w[i][j]%=MOD;
}
return res;
}
inline matrix powmod(matrix a,int b)
{
matrix res;for(register int i=0;i<=cnt;i++)res.w[i][i]=1;
for(;b;b>>=1,a=mul(a,a))if(b&1)res=mul(res,a);
return res;
}
int main()
{
//freopen("input.in","r",stdin);
map['A']=0;map['C']=1;map['G']=2;map['T']=3;
scanf("%d%d",&m,&n);
for(register int i=1;i<=m;i++)
{
scanf("%s",s+1);insert(s);
}
build();matrix res;
for(register int i=0;i<=cnt;i++)if(!danger[i])
{
for(register int j=0;j<4;j++)
if(!danger[tr[i][j]])res.w[i][tr[i][j]]++;
}
res=powmod(res,n);
for(register int i=0;i<=cnt;i++)
{
ans=ans+res.w[0][i];
while(ans>MOD)ans-=MOD;
}
printf("%d",ans);
return 0;
}
最新文章
- AlloyRenderingEngine之Shape
- WCF学习之旅——第一个WCF示例(一)
- Python学习路程day20
- linux进程间通信-管道
- Codeforces Round #344 (Div. 2)(按位或运算)
- bigDecimal 使用小结
- QTY N.W G.W
- 完整的站内搜索Demo(Lucene.Net+盘古分词)
- visual studio 2010 无法连接到ASP.NET Development Server
- Java缓存框架
- php+mysql 除了设置主键防止表单提交内容重复外的另一种方法
- 【春华秋实】深入源码理解.NET Core中Startup的注册及运行
- vuex state使用
- Maven学习第3期---m2eclipse使用
- 牛客练习赛38 D 题 出题人的手环 (离散化+树状数组求逆序对+前缀和)
- Matconvnet环境配置一些坑
- BZOJ1433 [ZJOI2009]假期的宿舍 二分图匹配 匈牙利算法
- vscode 调试 TypeScript
- ios黑科技
- 【Python】unittest-5
热门文章
- NO.A.0010——Windows常用快捷键使用教程
- 「CF645E」 Intellectual Inquiry
- Jinja2语法自动补全配置
- 【问题记录】— web页面调用本地程序
- IPSec传输模式/隧道模式下ESP报文的装包与拆包过程
- 流量控制--6.Classful Queuing Disciplines (qdiscs)
- 磁盘冗余阵列之RAID5、RAID10
- 笔记本无法连接校园网,windows诊断显示校园网之未响应
- JZOJ 【NOIP2017提高A组模拟9.14】捕老鼠
- JZOJ2020年8月12日提高组反思