[hdu2243][ac自动机]考研路茫茫——单词情结

考研路茫茫——单词情结

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 6436 Accepted Submission(s): 2209

Problem Description
背单词，始终是复习英语的重要环节。在荒废了3年大学生涯后，Lele也终于要开始背单词了。
一天，Lele在某本单词书上看到了一个根据词根来背单词的方法。比如”ab”,放在单词前一般表示”相反，变坏，离去”等。

于是Lele想，如果背了N个词根，那这些词根到底会不会在单词里出现呢。更确切的描述是：长度不超过L，只由小写字母组成的，至少包含一个词根的单词，一共可能有多少个呢？这里就不考虑单词是否有实际意义。

比如一共有2个词根 aa 和 ab ，则可能存在104个长度不超过3的单词，分别为
(2个) aa,ab,
(26个)aaa,aab,aac…aaz,
(26个)aba,abb,abc…abz,
(25个)baa,caa,daa…zaa,
(25个)bab,cab,dab…zab。

这个只是很小的情况。而对于其他复杂点的情况，Lele实在是数不出来了，现在就请你帮帮他。

Input
本题目包含多组数据，请处理到文件结束。
每组数据占两行。
第一行有两个正整数N和L。(0

sol::

在我前面做了那些ac自动机的基础上，这题考的就是矩乘了，这题坑点是模数，而且m<2^31，如果你+1的话就ll了

#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;
typedef unsigned long long ll;
int n,m;
inline int read()
{
    char c;
    int res,flag=0;
    while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
    res=c-'0';
    while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
    return flag?-res:res;
}
const int N=40;

int tot,l,top[N][26],fail[N];
int go[N][26];
char sr[11];
bool mark[N];
inline void insert(int &x,int y)
{
    if(!x)
    {
        x=++tot;
        fail[x]=0;
        mark[x]=0;
        for(int i=0;i<=25;++i) top[x][i]=go[x][i]=0;
    }
    if(y>l)
    {
        mark[x]=1;
        return;
    }
    insert(go[x][sr[y]-'a'],y+1);
}
int q[N];
inline void build()
{
    int t=0,w,u,v;
    q[w=1]=1;
    for(int i=0;i<=25;++i) top[0][i]=go[0][i]=1;
    while(t<w)
    {
        u=q[++t];
        mark[u]=mark[u]|mark[fail[u]];

        for(int ch=0;ch<=25;++ch)
        if(go[u][ch])
        {
            top[u][ch]=go[u][ch];
            q[++w]=go[u][ch];
            v=fail[u];
            while(!go[v][ch]) v=fail[v];
            v=go[v][ch];
            fail[go[u][ch]]=v;
        }
        else top[u][ch]=top[fail[u]][ch];
    }
}
int rt;
ll tmp[N][N],g[N][N],f[N][N];
inline void mul(ll a[N][N],ll b[N][N])
{
    memset(tmp,0,sizeof(tmp));
    for(int i=1;i<=tot;++i)
    for(int j=1;j<=tot;++j)
    for(int k=1;k<=tot;++k)
    tmp[i][j]+=a[i][k]*b[k][j];
    for(int i=1;i<=tot;++i)
    for(int j=1;j<=tot;++j)
    a[i][j]=tmp[i][j];
}
inline void ksm(ll t)
{
    while(t)
    {
        if(t&1) mul(f,g);
        mul(g,g);
        t>>=1;
    }
}
int main()
{
//    freopen("2243.in","r",stdin);
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        tot=rt=0;
        for(int i=1;i<=n;++i)
        {
            scanf("%s",sr+1);
            l=strlen(sr+1);
            insert(rt,1);
        }
        build();
        for(int i=1;i<=tot+1;++i)
        for(int j=1;j<=tot+1;++j)
        f[i][j]=g[i][j]=0;
        f[1][1]=1;
        for(int j=1;j<=tot;++j)
        if(!mark[j])
        for(int ch=0;ch<=25;++ch)
        g[j][top[j][ch]]++;
        else g[j][j]+=26,g[j][tot+1]++;
        ++tot;
        g[tot][tot]=1;
        ksm((ll)m+1);
        printf("%I64u\n",f[1][tot]);
    }
}