hdu2243

背单词，始终是复习英语的重要环节。在荒废了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<N<6,0<L<2^31) 第二行有N个词根，每个词根仅由小写字母组成，长度不超过5。两个词根中间用一个空格分隔开。

 

Output

对于每组数据，请在一行里输出一共可能的单词数目。 由于结果可能非常巨大，你只需要输出单词总数模2^64的值。

 

Sample Input

2 3
aa ab
1 2
a

 

Sample Output

104
52

 

 

和poj那个基本一样 ，，不过求A^1+A^2----A^n 有了花头 //之前一直以为白书上的是最简单的 没想到还有更简单的 。。

 

 

http://www.cnblogs.com/kuangbin/p/3159520.html   //代码看这里吧

 

 

 

还是附上我的代码吧=_=

 

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#include<algorithm>
using namespace std;
typedef unsigned long long LL;
struct Mat
{
    LL mat[36][36];
    int n;
    Mat(int _n)
    {
        n=_n;
        memset(mat,0,sizeof(mat));
    }
    Mat operator *(const Mat &B)const
    {
        Mat C(n);
        for(int k=0; k<n; ++k)
            for(int i=0; i<n; ++i)
            {
                if(mat[i][k]==0)  continue;
                for(int j=0; j<n; ++j)
                {
                    if(B.mat[k][j]==0) continue;
                    C.mat[i][j]=C.mat[i][j]+mat[i][k]*B.mat[k][j];
                }
            }
        return C;
    }
    Mat operator ^(int k)
    {
        Mat C(n);
        for(int i=0; i<n; ++i) C.mat[i][i]=1;
        while(k)
        {
            if(k&1) C=C*(*this),--k;
            k>>=1;
            *this=(*this)*(*this);
        }
        return C;
    }
};
struct AC
{
    int ch[38][26],fail[38],val[38],sz,rt;
    void init()
    {
        sz=rt=0;
        memset(ch[rt],-1,sizeof(ch[rt]));
    }
    void insert(char *str)
    {
        int len=strlen(str),u=rt;
        for(int i=0; i<len; ++i)
        {
            if(ch[u][str[i]-'a']==-1)
            {
                ++sz;
                memset(ch[sz],-1,sizeof(ch[sz]));
                val[sz]=0;
                ch[u][str[i]-'a']=sz;
            }
            u=ch[u][str[i]-'a'];
        }
        val[u]=1;
    }
    void build()
    {
        queue<int>Q;
        int u=rt;
        for(int i=0; i<26; ++i)
        {
            if(ch[u][i]==-1) ch[u][i]=rt;
            else
            {
                fail[ch[u][i]]=rt;
                Q.push(ch[u][i]);
            }
        }
        while(!Q.empty())
        {
            u=Q.front();
            Q.pop();
            val[u]|=val[fail[u]];
            for(int i=0; i<26; ++i)
            {
                if(ch[u][i]==-1) ch[u][i]=ch[fail[u]][i];
                else
                {
                    fail[ch[u][i]]=ch[fail[u]][i];
                    Q.push(ch[u][i]);
                }
            }
        }
    }
    Mat get()
    {
        Mat ret(sz+2);
        for(int i=0; i<=sz; ++i)
            for(int j=0; j<26; ++j)
                if(!val[ch[i][j]]) ++ret.mat[i][ch[i][j]];
        for(int i=0; i<=sz+1; ++i)
            ret.mat[i][sz+1]=1;
        return ret;
    }
} ac;
char s[55];
int main()
{
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        ac.init();
        for(int i=0; i<n; ++i)
        {
            scanf("%s",s);
            ac.insert(s);
        }
        ac.build();
        Mat A=ac.get();
        A=A^m;
        LL res=0;
        for(int i=0; i<A.n; ++i)
            res+=A.mat[0][i];
        --res;
        A=Mat(2);
        A.mat[0][0]=26;
        A.mat[1][1]=A.mat[0][1]=1;
        A=A^m;
        Mat B(2);
        B.mat[0][0]=B.mat[1][0]=1;
        A=A*B;
        LL ans=A.mat[0][0];
        ans-=res+1;
        cout<<ans<<endl;
    }
}