如标题所言,我已经沉迷于AC自动机无法自拔了。。。
这又是一道AC自动的题,红红火火恍恍惚惚



这题目做起来真舒服
简单概括一下:\(AC\)自动机\(fail\)树上树链剖分\(+\)树状数组
这种类型的题其实还蛮多的,比如这道:[HDU 4117] GRE Words
$
$
首先把\(AC\)自动机建出来,然后在所有子节点连一条由\(fail\)指向该点的边,这样一棵\(fail\)树就出来了。
题目问的是:求\(x\)在\(y\)中出现多少次,把问题放到\(fail\)树上来,就变成了:求从根到\(y\)的的节点中(这里指的是\(dfs\)序从根到\(y\)),有多少个在\(x\)的子树内
那么这个东西就很好求了,像普通的树链剖分题那样,用线段树就能维护,但是这道题只要查\(root\)到\(y\),\(bit\)也适用且常数要更小些
具体做法:我们离线来做这道题,把\(y\)相同的询问放到一起来处理
考虑这样几个做法:

  1. 每次遇到\('P'\)则统计答案当前点\(y\)的所有询问的答案;
  2. 遇到\('B'\)则将当前点的\(dfn\)从\(bit\)(或线段树)中删除;
  3. 否则往下跳,并将该节点插入\(bit\)(或线段树);
//made by Hero_of_Someone
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<vector>
#define N (100010)
using namespace std;

char s[N];

struct Trie{
  int size,root;
  int son[N][26],fail[N];
  int val[N],n,fa[N],ans[N];

  void init(){ size=1,root=0; }

  void insert(){
    int cur=root;
    for(int i=0;s[i];i++){
      if(s[i]=='P') val[++n]=cur;
      else if(s[i]=='B') cur=fa[cur];
      else{
    int id=s[i]-'a';
    if(!son[cur][id]) son[cur][id]=size++;
    fa[son[cur][id]]=cur,cur=son[cur][id];
      }
    }
  }

  int num,head[N],nxt[N],to[N];
  void add(int u,int v){
    nxt[++num]=head[u];to[num]=v;head[u]=num;
  }

  void build(){
    int que[N];
    int hd=0,tl=0;
    for(int i=0;i<26;i++)
      if(son[root][i]){
    que[tl++]=son[root][i];
    fail[son[root][i]]=root;
      }
      else son[root][i]=root;

    while(hd<tl){
      int cur=que[hd++];
      for(int i=0;i<26;i++){
    int Son=son[cur][i];
    if(Son){
      int f=fail[cur];
      while(f&&!son[f][i]) f=fail[f];
      fail[Son]=son[f][i];
      que[tl++]=Son;
    }
    else son[cur][i]=son[fail[cur]][i];
      }
    }
    for(int i=1;i<size;i++) add(fail[i],i);
  }

  int st[N],ed[N],tim,t[N];
  void dfs(int x){
    st[x]=++tim;
    for(int i=head[x];i;i=nxt[i]) dfs(to[i]);
    ed[x]=tim;
  }

  int lowbit(int x){return x&-x;}
  void up(int x,int val){for(;x<=tim;x+=lowbit(x))t[x]+=val;}
  int Q(int x){int ret=0;for(;x;x-=lowbit(x))ret+=t[x];return ret;}

  struct node{ int x,id; }; vector<node>p[N];

  void Ans(){
    int m; scanf("%d",&m);
    for(int i=1;i<=m;i++){
      int x,y; scanf("%d%d",&x,&y);
      p[val[y]].push_back((node){x,i});
    }
    int x=0;
    for(int i=0;s[i];i++){
      if(s[i]=='P')
    for(int j=0,l=p[x].size();j<l;j++){
      int X=p[x][j].x,ID=p[x][j].id;
      ans[ID]=Q(ed[val[X]])-Q(st[val[X]]-1);
    }
      else if(s[i]=='B') up(st[x],-1),x=fa[x];
      else x=son[x][s[i]-'a'],up(st[x],1);
    }
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
  }

}AC;

void init(){ AC.init(); scanf("%s",s); AC.insert(); AC.build(); }

void work(){ AC.dfs(0); AC.Ans(); }

int main(){ init(); work(); return 0; }

沉迷AC自动机无法自拔之:[BZOJ2434] [Noi2011] 阿狸的打字机的更多相关文章

  1. 沉迷AC自动机无法自拔之:穿越广场 square

    如标题所言,我已经沉迷于AC自动机无法自拔了... 这又是一道AC自动的题,红红火火恍恍惚惚 穿越广场 [问题描述] L 国的仪仗队要穿越首都广场了.首都广场可以看做是一块 N*M 的矩形网格,仪仗队 ...

  2. BZOJ2434 [Noi2011]阿狸的打字机 【AC自动机 + fail树 + 树状数组】

    2434: [Noi2011]阿狸的打字机 Time Limit: 10 Sec  Memory Limit: 256 MB Submit: 3610  Solved: 1960 [Submit][S ...

  3. BZOJ2434: [NOI2011]阿狸的打字机(AC自动机+dfs序+树状数组)

    [NOI2011]阿狸的打字机 题目链接:https://www.luogu.org/problemnew/show/P2414 题目背景 阿狸喜欢收藏各种稀奇古怪的东西,最近他淘到一台老式的打字机. ...

  4. BZOJ2434 [Noi2011]阿狸的打字机(AC自动机 + fail树 + DFS序 + 线段树)

    题目这么说的: 阿狸喜欢收藏各种稀奇古怪的东西,最近他淘到一台老式的打字机.打字机上只有28个按键,分别印有26个小写英文字母和'B'.'P'两个字母.经阿狸研究发现,这个打字机是这样工作的: 输入小 ...

  5. BZOJ2434: [Noi2011]阿狸的打字机(AC自动机 树状数组)

    Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 4140  Solved: 2276[Submit][Status][Discuss] Descript ...

  6. BZOJ2434[Noi2011]阿狸的打字机——AC自动机+dfs序+树状数组

    题目描述 阿狸喜欢收藏各种稀奇古怪的东西,最近他淘到一台老式的打字机.打字机上只有28个按键,分别印有26个小写英文字母和'B'.'P'两个字母. 经阿狸研究发现,这个打字机是这样工作的: l 输入小 ...

  7. 沉迷AC自动机无法自拔之:[UVALive 4126] Password Suspects

    图片加载可能有点慢,请跳过题面先看题解,谢谢 一看到这么多模式串就非常兴奋,又是\(AC\)自动机 题目就是要求:经过 \(n\) 个节点,把所有单词都遍历一遍的方案数,和那道题差不多嘛 所以这样设: ...

  8. 沉迷AC自动机无法自拔之:[UVA 11468] Substring

    图片加载可能有点慢,请跳过题面先看题解,谢谢 这个鬼题目,上一波套路好了 先用题目给的模板串建\(AC\)自动机,把单词结尾标记为 \(val=1\),然后在建好的\(AC\)自动机上跑 \(dp\) ...

  9. BZOJ2434:[NOI2011]阿狸的打字机(AC自动机,线段树)

    Description 阿狸喜欢收藏各种稀奇古怪的东西,最近他淘到一台老式的打字机.打字机上只有28个按键,分别印有26个小写英文字母和'B'.'P'两个字母. 经阿狸研究发现,这个打字机是这样工作的 ...

随机推荐

  1. dpkg打包与解包

    1.打包 dpkg -b 2.解包 2.1 dpkg -X 解出包内容 2.2 dpkg -e 输出包控制信息

  2. 【强化学习】python 实现 q-learning 例三(例一改写)

    本文作者:hhh5460 本文地址:https://www.cnblogs.com/hhh5460/p/10139738.html 例一的代码是函数式编写的,这里用面向对象的方式重新撸了一遍.好处是, ...

  3. 使用while和read命令读取文件内容

    转:使用while和read命令读取文件内容 1.准备数据文件 $cat a.txt 200:2 300:3 400:4 500:5 2.用while循环从文件中读取数据 #!/bin/ksh whi ...

  4. JS 数据处理技巧及小算法汇总

    前言: 金秋九月的最后一天,突然发现这个月博客啥也没更新,不写点什么总觉得这个月没啥长进,逆水行舟,不进则退,前进的路上贵在坚持,说好的每个月至少一到两篇,不能半途而废!好多知识写下来也能加深一下自身 ...

  5. item 4: 知道怎么去看推导的类型

    本文翻译自modern effective C++,由于水平有限,故无法保证翻译完全正确,欢迎指出错误.谢谢! 博客已经迁移到这里啦 对于推导类型结果的查看,根据不同的软件开发阶段,你想知道的信息的不 ...

  6. Webpack 2 视频教程

    这是我免费发布的高质量超清「Webpack 2 视频教程」. Webpack 作为目前前端开发必备的框架,Webpack 发布了 2.0 版本,此视频就是基于 2.0 的版本讲解的. 这个基本就是目前 ...

  7. MongoDB日常运维操作命令小结

    总所周知,MongoDB是一个NoSQL非数据库系统,即一个数据库可以包含多个集合(Collection),每个集合对应于关系数据库中的表:而每个集合中可以存储一组由列标识的记录,列是可以自由定义的, ...

  8. 太白教你学python---博客分类目录

    太白非技术类随笔(持续更新中...猛击这里!!!) python基础 python基础一 pytcharm安装详细教程 python基础二 python基础数据类型 Python最详细,最深入的代码块 ...

  9. gitblit 配置图文详解

    Windows平台下Git服务器搭建 前提是确保存在JDK环境. 第一步:下载Gitblit.下载地址:http://www.gitblit.com/ 第二步:解压缩下载的压缩包即可,无需安装. 第三 ...

  10. Individual P1: Summary

    经过5个小时成功把simple mode写差不多了..orz 也是蛮拼的. 开始毫无头绪,本能地开始从度娘搜索‘c# 单词统计’= =看了两段代码也算是见过c#的人了.差不多花了我1小时的时间. 然后 ...