AC自动机——Aho-Corasick Automaton

这是一个英文版的讲的比较好的AC自动机资料。

http://www.cs.uku.fi/~kilpelai/BSA05/lectures/slides04.pdf

如果不爱看英文，可以看我整理的大致的翻译，再加上点解释说明啥的，建议中英两个版本结合着看，毕竟我翻译的里面可能有些错误。

http://download.csdn.net/download/morgan_xww/4476863

可以预览一下文档中的截图：

我的模板：

class ACAutomaton{public:    static const int MAX_N = 10000 * 50 + 5;        //最大结点数：模式串个数 X 模式串最大长度    static const int CLD_NUM = 26;        //从每个结点出发的最多边数，字符集Σ的大小，一般是26个字母    int n;                          //trie树当前结点总数    int id['z'+1];                  //字母x对应的结点编号为id[x]    int fail[MAX_N];                //fail指针    int tag[MAX_N];                 //根据题目而不同    int trie[MAX_N][CLD_NUM];       //trie树，也就是goto函数    void init()    {        for (int i = 0; i < CLD_NUM; i++)            id['a'+i] = i;    }    void reset()    {        memset(trie[0], -1, sizeof(trie[0]));        tag[0] = 0;        n = 1;    }    //插入模式串s，构造单词树(keyword tree)    void add(char *s)    {        int p = 0;        while (*s)        {            int i = id[*s];            if ( -1 == trie[p][i] )            {                memset(trie[n], -1, sizeof(trie[n]));                tag[n] = 0;                trie[p][i] = n++;            }            p = trie[p][i];            s++;        }        tag[p]++;         //因题而异    }    //构造AC自动机，用BFS来计算每个结点的fail指针，就是构造trie图    void construct()    {        queue<int> Q;        fail[0] = 0;        for (int i = 0; i < CLD_NUM; i++)        {            if (-1 != trie[0][i])            {                fail[trie[0][i]] = 0; //root下的第一层结点的fail指针都指向root                Q.push(trie[0][i]);            }            else            {                trie[0][i] = 0;    //这是阶段一中的第2步            }        }        while ( !Q.empty() )        {            int u = Q.front();            Q.pop();            for (int i = 0; i < CLD_NUM; i++)            {                int &v = trie[u][i];                if ( -1 != v )                {                    Q.push(v);                    fail[v] = trie[fail[u]][i];                    tag[u] += tag[fail[u]];     //因题而异，某些题目中不需要这句话                }                else                {            //当trie[u][i]==-1时，设置其为trie[fail[u]][i]，就构造了trie图                    v = trie[fail[u]][i];                }            }        }    }    //因题而异    //在目标串t中匹配模式串    int solve(char *t)    {        int q = 0, ret = 0;        while ( *t )        {            q = trie[q][id[*t]];            int u = q;            while ( u != 0 )            {                ret += tag[u];                tag[u] = 0;                u = fail[u];            }            t++;        }        return ret;    }} ac;