[JSOI2007]文本生成器

题目大意：
　　给定$n(nle60)$个字符串$s_i(|s_i|le100)$，问有多少个长度为$m(mle100)$的字符串$t$满足至少包含一个$s_i$？保证所有字符串均由大写英文字母构成。

思路：
　　建立AC自动机。用$f[i][j]$表示构造到$t$的第$i$位，对应自动机上编号为$j$的结点，不包含任何$s_i$的方案数。显然遇到危险结点需要跳过，否则直接转移即可。答案即为$26^m-sum f[m][i]$。

#include<queue>
#include<cstdio>
#include<cctype>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=60,L=101,S=26,mod=10007;
char s[L];
int f[L][N*L];
class AhoCorasick {
    private:
        bool val[N*L];
        std::queue<int> q;
        int ch[N*L][S],fail[N*L];
        int sz,new_node() {
            return ++sz;
        }
        int idx(const int &c) const {
            return c-'A';
        }
    public:
        void insert(const char s[]) {
            int p=0;
            for(register int i=0;s[i];i++) {
                const int c=idx(s[i]);
                p=ch[p][c]?:ch[p][c]=new_node();
            }
            val[p]=true;
        }
        void get_fail() {
            for(register int c=0;c<S;c++) {
                if(ch[0][c]) q.push(ch[0][c]);
            }
            while(!q.empty()) {
                const int &x=q.front();
                for(register int c=0;c<S;c++) {
                    int &y=ch[x][c];
                    if(!y) {
                        y=ch[fail[x]][c];
                        continue;
                    }
                    fail[y]=ch[fail[x]][c];
                    val[y]|=val[fail[y]];
                    q.push(y);
                }
                q.pop();
            }
        }
        int solve(const int &n,const int &m) const {
            int ret=f[0][0]=1;
            for(register int i=1;i<=m;i++) {
                (ret*=S)%=mod;
                for(register int j=0;j<=sz;j++) {
                    if(!f[i-1][j]) continue;
                    for(register int c=0;c<S;c++) {
                        if(!val[ch[j][c]]) (f[i][ch[j][c]]+=f[i-1][j])%=mod;
                    }
                }
            }
            for(register int i=0;i<=sz;i++) {
                ret=(ret-f[m][i]+mod)%mod;
            }
            return ret;
        }
};
AhoCorasick ac;
int main() {
    const int n=getint(),m=getint();
    for(register int i=0;i<n;i++) {
        scanf("%s",s);
        ac.insert(s);
    }
    ac.get_fail();
    printf("%d
",ac.solve(n,m));
    return 0;
}