注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

ydc的博客

 
 
 

日志

 
 

bzoj 1030  

2013-02-05 13:16:40|  分类: bzoj |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

文本生成器

AC自动机+dp。

首先用补集转化思想,答案就是26^m-不存在一个的方案数

不存在一个的方案数就是AC自动机+DP

任意一个串,对应AC自动机的某一个结点

然后就可以转移了,跟08年湖南省选那题很像,不过这道题不需要矩阵加速

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#define MOD 10007
#define MAXN 80
#define MAXL 110
#define MAXSIZE 26
using namespace std;
int n,m,tot,total,F[MAXL][MAXN*MAXL];
char word[MAXL];
struct AC
{
    int son[MAXSIZE],next,last;
    bool mark;
}tree[MAXN*MAXL];
queue<int>q;
int power(int n,int p)
{
    int ans=1;
    for(int i=n;i;i>>=1,p=p*p%MOD)
        if(i&1)
            ans=ans*p%MOD;
    return ans;
}
void Insert(char a[])
{
    int root=0;
    for(int i=0;a[i];i++)
    {
        if(tree[root].son[a[i]-'A']==0)  tree[root].son[a[i]-'A']=++tot;
        root=tree[root].son[a[i]-'A'];
    }
    tree[root].mark=true;
}
void read()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%s",word);
        Insert(word);
    }
}
void Bfs()
{
    for(int i=0;i<MAXSIZE;i++)
        if(tree[0].son[i])
            q.push(tree[0].son[i]);
    int p;
    while(!q.empty())
    {
        p=q.front(),q.pop();
        for(int i=0;i<MAXSIZE;i++)
            if(tree[p].son[i])
            {
                int u=tree[p].son[i],v=tree[p].next;
                for(;v&&!tree[v].son[i];v=tree[v].next);
                tree[u].next=v=tree[v].son[i],q.push(u),tree[u].last=tree[v].mark?v:tree[v].last;
            }
    }
}
void Dp()
{
    F[0][0]=1;
    for(int i=0;i<m;i++)
        for(int j=0;j<=tot;j++)
            if(F[i][j])
                for(int k=0;k<MAXSIZE;k++)
                {
                    int v=j;
                    for(;v&&!tree[v].son[k];v=tree[v].next);
                    v=tree[v].son[k];
                    if(!tree[v].mark&&!tree[v].last)  F[i+1][v]=(F[i+1][v]+F[i][j])%MOD;
                }
    for(int i=0;i<=tot;i++)
        total=(total+F[m][i])%MOD;
}
int main()
{
    read();
    Bfs();
    Dp();
    printf("%d\n",(power(m,26)-total+MOD)%MOD);
    return 0;
}



  评论这张
 
阅读(714)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017