• 土地是自然存在的地球的一部分,不是人类劳动成果,没有价值,买房者不应支付土地费,房价之所以虚高,是因为买房者支付了不该支付的土地费。 2019-10-15
  • 月薪过万白领辞职回乡养鸡 亏数十万后回城找工作 2019-10-15
  • 尽管过去计划经济时,由于局限性,工作上有过失误,有过短缺,但没有浪费。不存在笑博士说的做出来的东西无人用的问题。何况那时,主要还是通过人工作的计划。在今天,高科 2019-09-17
  • 即使中美贸易归零,中国也不会屈服美国的大棒政策 2019-09-17
  • 宁夏文联面向全国重奖征文 喜迎自治区成立六十周年 2019-09-10
  • 大兴首开保利·熙悦林语 又一限价房将入市! ——凤凰网房产北京 2019-09-10
  • 大山深处4名孩子的交通安全课 2019-09-06
  • 新疆坚决打好污染防治攻坚战 2019-09-06
  • 分级基金将黯然离场 10只分级B拉响下折“警报” 2019-09-05
  • 刘嘉玲戴礼帽系领带帅过鲜肉 大长腿走路带风霸气在侧漏 2019-09-05
  • 阶级不是“等级差别”,而是私有制基础上剥削和被剥削两大社会集团。 2019-08-28
  • 邮储银行山西省分行零售信贷结余突破200亿元 2019-08-24
  • 智能服装让生活更有趣 2019-08-24
  • 探访广西高考评卷现场:评卷场所警戒线“包围” 2019-08-23
  • 一镜到底绝密视频!10个故事,一个不一样的上合峰会 2019-08-21
  • 意甲桑普多利亚vs:[题解]POJ2778 DNA Sequence

    都灵桑普多利亚关系 www.dqyww.com 3/8/2017来源:ASP.NET技巧人气:10221

    Description

    It’s well known that DNA Sequence is a sequence only contains A, C, T and G, and it’s very useful to analyze a segment of DNA Sequence,For example, if a animal’s DNA sequence contains segment ATC then it may mean that the animal may have a genetic disease. Until now scientists have found several those segments, the PRoblem is how many kinds of DNA sequences of a species don’t contain those segments.

    Suppose that DNA sequences of a species is a sequence that consist of A, C, T and G,and the length of sequences is a given integer n.

    Input

    First line contains two integer m (0 <= m <= 10), n (1 <= n <=2000000000). Here, m is the number of genetic disease segment, and n is the length of sequences.

    Next m lines each line contain a DNA genetic disease segment, and length of these segments is not larger than 10.

    Output

    An integer, the number of DNA sequences, mod 100000.

    Sample Input

    4 3 AT AC AG AA

    Sample Output

    36

    Solution

    题目大意:给定m个字符串和一个长度n,求不包含这m个字符串的长度为n的字符串个数mod 100000的值,字符串中只含’A’、’C’、’T’、’G’。 做法: 首先构造一下AC自动机我们就会发现,题目要求的就是在AC自动机上跑n次同时不经过任何一个表示字符串的点的方案总数。仔细想一想会发现,有两种点不合法:一是直接表示字符串的点,二是fail边连向表示字符串的点的点。所以我们就可以DP,在构造AC自动机的时候给不合法的点打上标记。用fi,j表示跑了i次到达j节点的合法路径总数,dp如下:

    f[0][1]=1; //初始化 for(int i=1;i<=n;i++){ //n为输入的字符串长度 for(int j=1;j<=num;j++){ //num为AC自动机中节点总数,1节点是根 if(!T[j].flag&&f[i-1][j]){ //如果该节点合法并且当前点的答案不为0(如果当前节点答案为0就没必要更新了,更新等于没更新) for(int k=1;k<=4;k++){ //枚举4种字母 int p=j; while(!T[p].next[k])p=T[p].fail; //如果不匹配就沿着失配边走(为了防止死循环,已经提早把T[0].next全部赋为1) (f[i][T[p].next[k]]+=f[i-1][j])%=mod; //更新其他点的答案 } } } }

    但是这样肯定不行,n有2?109,所以我们可以用矩阵乘法优化当前DP。矩阵还是容易推的,实在不懂看看代码就行了。 代码:

    #include<cstdio> #include<cstring> #include<algorithm> using namespace std; typedef long long LL; const int maxn=110,mod=100000; struct node{ int next[5],fail,flag; }T[maxn]; struct matrix{ LL a[maxn][maxn]; int l,r; matrix(){memset(a,0,sizeof a);} }temp,c,I; int n,m,num=1,que[maxn],ans,head,tail,Hash[256]; char ch[12]; void Insert(char *s){ int p=1,len=strlen(s); for(int i=0;i<len;i++){ if(!T[p].next[Hash[s[i]]])T[p].next[Hash[s[i]]]=++num; p=T[p].next[Hash[s[i]]]; } T[p].flag=true; } void Bfs(){ //构造AC自动机的fail边,顺便标记不合法节点 memset(que,head=tail=0,sizeof que); que[++tail]=1; while(head<tail){ int x=que[++head]; for(int i=1;i<=4;i++){ if(T[x].next[i]){ int p=T[x].fail; while(!T[p].next[i])p=T[p].fail; T[T[x].next[i]].fail=T[p].next[i]; if(T[T[p].next[i]].flag)T[T[x].next[i]].flag=true; que[++tail]=T[x].next[i]; } } } } matrix mul(matrix x,matrix y){ matrix re; re.l=x.l;re.r=y.r; for(int i=1;i<=re.l;i++){ for(int j=1;j<=re.r;j++){ for(int k=1;k<=x.r;k++){ (re.a[i][j]+=x.a[i][k]*y.a[k][j])%=mod; } } } return re; } matrix pow(matrix x,int y){ matrix re=I; while(y){ if(y&1)re=mul(re,x); x=mul(x,x);y>>=1; } return re; } int main(){ Hash['A']=1;Hash['C']=2;Hash['T']=3;Hash['G']=4; for(int i=1;i<=4;i++)T[0].next[i]=1; scanf("%d%d",&n,&m); for(int i=1;i<=n;i++){ scanf("\n%s",ch); Insert(ch); } Bfs();temp.l=1;c.l=c.r=I.l=I.r=temp.r=num; for(int i=1;i<=num;i++)I.a[i][i]=1; for(int i=1;i<=num;i++){ //转移矩阵 if(T[i].flag)continue; for(int j=1;j<=4;j++){ int p=i; while(!T[p].next[j])p=T[p].fail; c.a[i][T[p].next[j]]++; } } temp.a[1][1]=1; temp=mul(temp,pow(c,m)); for(int i=1;i<=num;i++){ //统计答案 if(!T[i].flag){ (ans+=temp.a[1][i])%=mod; } } printf("%d\n",ans); return 0; }

  • 土地是自然存在的地球的一部分,不是人类劳动成果,没有价值,买房者不应支付土地费,房价之所以虚高,是因为买房者支付了不该支付的土地费。 2019-10-15
  • 月薪过万白领辞职回乡养鸡 亏数十万后回城找工作 2019-10-15
  • 尽管过去计划经济时,由于局限性,工作上有过失误,有过短缺,但没有浪费。不存在笑博士说的做出来的东西无人用的问题。何况那时,主要还是通过人工作的计划。在今天,高科 2019-09-17
  • 即使中美贸易归零,中国也不会屈服美国的大棒政策 2019-09-17
  • 宁夏文联面向全国重奖征文 喜迎自治区成立六十周年 2019-09-10
  • 大兴首开保利·熙悦林语 又一限价房将入市! ——凤凰网房产北京 2019-09-10
  • 大山深处4名孩子的交通安全课 2019-09-06
  • 新疆坚决打好污染防治攻坚战 2019-09-06
  • 分级基金将黯然离场 10只分级B拉响下折“警报” 2019-09-05
  • 刘嘉玲戴礼帽系领带帅过鲜肉 大长腿走路带风霸气在侧漏 2019-09-05
  • 阶级不是“等级差别”,而是私有制基础上剥削和被剥削两大社会集团。 2019-08-28
  • 邮储银行山西省分行零售信贷结余突破200亿元 2019-08-24
  • 智能服装让生活更有趣 2019-08-24
  • 探访广西高考评卷现场:评卷场所警戒线“包围” 2019-08-23
  • 一镜到底绝密视频!10个故事,一个不一样的上合峰会 2019-08-21
  • 相关推荐:北京pk10高手计划群 psv有什么好玩的游戏 北京pk10骗局全过程 欢乐牛牛 重庆时时龙虎玩法 时时彩一期一计划大小单双 彩票极速赛车冠军定胆 cc天下彩图文资讯 骰子猜大小单双技巧 麻将规则打法 重庆时时彩网址 时时彩压龙虎对着压 下载重庆时时 免费金币二八杠游戏下载 重庆时时彩走势 可靠时时彩计划软件