[TJOI2013]单词

题目大意：
　　给定由$n$个单词组成的文章$(sum|s_i|leq10^6)$，统计每个单词出现的次数。

思路：
　　这一题乍一眼看上去像是用所有单词构造一个AC自动机，然后再用每个字符串去匹配，每次跳转失配指针统计答案。这样确实能在BZOJ上A掉，但是在洛谷上却被卡成90分（考虑一些很长的字符串仅由同一种字母构成，且长度依次递增，这样匹配每个串是$O(n^2)$的）。正确的算法应该是用所有的失配指针构造Fail树，每次匹配不跳转指针，只在当前结点上记录匹配的次数，最后再在Fail树上推一遍。

#include<list>
#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
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=201,L=1e6+1,S=26;
int cnt[N],p[N];
char s[L];
class AhoCorasick {
    private:
        std::queue<int> q;
        int ch[L][S],fail[L],sum[L];
        std::list<int> list[L],e[L];
        int sz,new_node() {
            return ++sz;
        }
        int idx(const char &c) const {
            return c-'a';
        }
        void dfs(const int &x) {
            for(std::list<int>::const_iterator i=e[x].begin();i!=e[x].end();i++) {
                const int &y=*i;
                dfs(y);
                sum[x]+=sum[y];
            }
        }
    public:
        void insert(const char s[],const int &id) {
            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();
            }
            list[p].push_front(id);
        }
        void get_fail() {
            for(register int c=0;c<S;c++) {
                if(ch[0][c]) {
                    q.push(ch[0][c]);
                    e[0].push_front(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];
                    e[ch[fail[x]][c]].push_front(y);
                    q.push(y);
                }
                q.pop();
            }
        }
        void find(const char s[],const int &n) {
            for(register int i=0,p=0;i<n;i++) {
                sum[p=ch[p][idx(s[i])]]++;
            }
        }
        void stat() {
            dfs(0);
            for(register int i=1;i<=sz;i++) {
                for(register std::list<int>::const_iterator j=list[i].begin();j!=list[i].end();j++) {
                    cnt[*j]+=sum[i];
                }
            }
        }
};
AhoCorasick ac;
int main() {
    const int n=getint();
    for(register int i=0;i<n;i++) {
        p[i]=strlen(s);
        scanf("%s",&s[p[i]]);
        ac.insert(&s[p[i]],i);
    }
    p[n]=strlen(s);
    ac.get_fail();
    for(register int i=0;i<n;i++) {
        ac.find(&s[p[i]],p[i+1]-p[i]);
    }
    ac.stat();
    for(register int i=0;i<n;i++) {
        printf("%d
",cnt[i]);
    }
    return 0;
}