模板统计LA 4670 Dominating Patterns

本篇文章是一篇关于模板统计的帖子

    题目链接：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&category=29&page=show_problem&problem=2671

    统计文本中模板单词的涌现频率。

    AC自动机模板题：

    每日一道理
人生是洁白的画纸，我们每个人就是手握各色笔的画师；人生也是一条看不到尽头的长路，我们每个人则是人生道路的远足者；人生还像是一块神奇的土地，我们每个人则是手握农具的耕耘者；但人生更像一本难懂的书，我们每个人则是孜孜不倦的读书郎。

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <map>
#include <queue>
#include <algorithm>
using namespace std;

#define SIGMA_SIZE 26
#define MAXNODE 20000
#define TEXT_SIZE 1000005
#define P_SIZE 80
#define P_NUM 160

map<string,int> ms;
struct AcCorsickAutomata
{
    int cnt[P_NUM];
    int sz;
    int ch[MAXNODE][SIGMA_SIZE];
    int f[MAXNODE];
    int val[MAXNODE];
    int last[MAXNODE];
    void init()
    {
        sz = 1;
        memset(ch[0],0,sizeof(ch[0]));
        memset(cnt,0,sizeof(cnt));
        ms.clear();
    }
    int idx(char c)
    {
        return c - 'a';
    }
    void insert(char *s,int v)
    {
        int u = 0,n = strlen(s);
        for(int i=0; i<n; i++)
        {
            int c = idx(s[i]);
            if(!ch[u][c])
            {
                memset(ch[sz],0,sizeof(ch[sz]));
                val[sz] = 0;
                ch[u][c] = sz++;
            }
            u = ch[u][c];
        }
        val[u] = v;
        ms[string(s)] = v;
    }

    void print(int j)
    {
        if(j)
        {
            cnt[val[j]]++;
            print(last[j]);
        }
    }
    void find(char *T)
    {
        int n = strlen(T);
        int j = 0;
        for(int i=0; i<n; i++)
        {
            int c = idx(T[i]);
            while(j && !ch[j][c]) j = f[j];
            j = ch[j][c];
            if(val[j]) print(j);
            else if(last[j]) print(last[j]);
        }
    }
    void getFail()
    {
        queue<int> q;
        f[0] = 0;
        for(int c = 0; c<SIGMA_SIZE; c++)
        {
            int u = ch[0][c];
            if(u)
            {
                f[u] = 0;
                q.push(u);
                last[u] = 0;
            }
        }
        while(!q.empty())
        {
            int r = q.front();
            q.pop();
            for(int c = 0; c<SIGMA_SIZE; c++)
            {
                int u = ch[r][c];
                if(!u) continue;
                q.push(u);
                int v = f[r];
                while(v && !ch[v][c]) v = f[v];
                f[u] = ch[v][c];
                last[u] = val[f[u]] ? f[u] : last[f[u]];
            }
        }
    }
};
AcCorsickAutomata ac;
char text[TEXT_SIZE],p[P_NUM][P_SIZE];
int main()
{
/*#ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
#endif*/
    int n;
    while(scanf(" %d",&n)!=EOF && n!=0)
    {
        ac.init();
        for(int i=1; i<=n; i++)
        {
            scanf(" %s",p[i]);
            ac.insert(p[i],i);
        }
        scanf(" %s",text);
        ac.getFail();
        ac.find(text);
        int best = -1;
        for(int i=1; i<=n; i++)
        {
            if(ac.cnt[i] > best) best = ac.cnt[i];
        }
        printf("%d\n",best);
        for(int i=1; i<=n; i++)
        {
            if(ac.cnt[ms[string(p[i])]] == best) printf("%s\n",p[i]);
        }
    }
    return 0;
}

    

文章结束给大家分享下程序员的一些笑话语录： 开发时间
　　项目经理: 如果我再给你一个人,那可以什么时候可以完工?程序员: 3个月吧!项目经理: 那给两个呢?程序员: 1个月吧!
项目经理: 那100呢?程序员: 1年吧!
项目经理: 那10000呢?程序员: 那我将永远无法完成任务.

--------------------------------- 原创文章 By
模板和统计
---------------------------------