[TJOI2013]单词(AC自动机,树形dp)

题目

给定N个单词,这N个单词形成一篇文章(单词间隔断).统计每个单词出现次数.

数据范围:N\leq 200, \sum len\leq 10^{6}

题解

可以看出每个单词的贡献是独立的,也是对整体有影响的,便想到用Trie树来统计贡献.

如果插入一个单词(如aaa),则每层计数的cnt+1.这样就统计了aaa所有前缀的贡献.如aab没有计算ab.

考虑AC自动机的fail指针的转移:fail(x)=maxlen\{P_k|\ P_k \sqsupset x \}(即P_kx的后缀)

如果在插入所有单词后建立fail指针,那么可以知道一个单词如aab的所有后缀(如果后缀有贡献)会通过fail指针链接起来

形成了一棵树,进行树形dp即可求出所有的贡献

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <queue>
using namespace std;

struct Trie {
  int next[26];
  int fail, cnt, sum;
} trie[500005];
int cnt_trie;
void insert(string );
void getfail();

struct edge {
  int next, to;
} e[500005];
int head[500005], cnt_adj;
void add_edge(int, int);
bool vis[500005];
void getsum(int);
int find(string);
int N;
string s[205];

int main(){
  cin >> N;
  int i;
  for(i = 1; i <= N; i++){
    cin >> s[i];
    insert(s[i]);
  }
  getfail();
  getsum(0);
  for(i = 1; i <= N; i++)
    cout << find(s[i]) << endl;
  return 0;
}

void insert(string s){
  int i, len = s.size();
  int id, u = 0;
  for(i = 0; i < len; i++){
    id = s[i] - 'a';
    if(!trie[u].next[id])
      trie[u].next[id] = ++cnt_trie;
    u = trie[u].next[id];
    trie[u].cnt++;
  }
}

void getfail(){
  int i;
  queue<int> Q;
  for(i = 0; i < 26; i++)
    if(trie[0].next[i])
      Q.push(trie[0].next[i]);

  while(!Q.empty()){
    auto u = Q.front(); Q.pop();

    for(i = 0; i < 26; i++){
      if(trie[u].next[i]){
        Q.push(trie[u].next[i]);
        trie[trie[u].next[i]].fail = trie[trie[u].fail].next[i];
      }
      else
        trie[u].next[i] = trie[trie[u].fail].next[i];
    }
    add_edge(trie[u].fail, u);
  }
}

void getsum(int u){
  int i;
  if(vis[u]) return;
  vis[u] = true;
  trie[u].sum = trie[u].cnt;
  for(i = head[u]; i; i = e[i].next){
    int v = e[i].to;
    getsum(v), trie[u].sum += trie[v].sum;
  }
}

int find(string s){
  int i, len = s.size();
  int id, u = 0;
  for(i = 0; i < len; i++){
    id = s[i] - 'a';
    u = trie[u].next[id];
  }
  return trie[u].sum;
}

void add_edge(int x, int y){
  e[++cnt_adj].to = y;
  e[cnt_adj].next = head[x];
  head[x] = cnt_adj;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值