这题太神啦

题意:求长度为$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;
}

最新文章

  1. AlloyRenderingEngine之Shape
  2. WCF学习之旅——第一个WCF示例(一)
  3. Python学习路程day20
  4. linux进程间通信-管道
  5. Codeforces Round #344 (Div. 2)(按位或运算)
  6. bigDecimal 使用小结
  7. QTY N.W G.W
  8. 完整的站内搜索Demo(Lucene.Net+盘古分词)
  9. visual studio 2010 无法连接到ASP.NET Development Server
  10. Java缓存框架
  11. php+mysql 除了设置主键防止表单提交内容重复外的另一种方法
  12. 【春华秋实】深入源码理解.NET Core中Startup的注册及运行
  13. vuex state使用
  14. Maven学习第3期---m2eclipse使用
  15. 牛客练习赛38 D 题 出题人的手环 (离散化+树状数组求逆序对+前缀和)
  16. Matconvnet环境配置一些坑
  17. BZOJ1433 [ZJOI2009]假期的宿舍 二分图匹配 匈牙利算法
  18. vscode 调试 TypeScript
  19. ios黑科技
  20. 【Python】unittest-5

热门文章

  1. NO.A.0010——Windows常用快捷键使用教程
  2. 「CF645E」 Intellectual Inquiry
  3. Jinja2语法自动补全配置
  4. 【问题记录】— web页面调用本地程序
  5. IPSec传输模式/隧道模式下ESP报文的装包与拆包过程
  6. 流量控制--6.Classful Queuing Disciplines (qdiscs)
  7. 磁盘冗余阵列之RAID5、RAID10
  8. 笔记本无法连接校园网,windows诊断显示校园网之未响应
  9. JZOJ 【NOIP2017提高A组模拟9.14】捕老鼠
  10. JZOJ2020年8月12日提高组反思