CF 163E //AC自动机fail树上建线段树 || 建树状数组+差分

最新推荐文章于 2022-05-06 15:31:35 发布
最新推荐文章于 2022-05-06 15:31:35 发布
阅读量1k 收藏 2
点赞数 1
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41730604/article/details/102330495
版权

题目

题意

n ( 1 e 6 ) n(1e6) n(1e6)个询问， k k k个模式串(总长度 ( 1 e 6 ) (1e6) (1e6))。
询问：
1：给定串 s s s求， s s s串里模式串出现几次。
2：将第 i i i个串删除(若已经删除无视此操作。
2：将第 i i i个串恢复(若已经存在无视此操作。

思路

建立AC自动机，在fail树上，用当前结点去更新后代。
可以在fail树上先序 d f s 序 dfs序 dfs序上维护线段树(区间更新，单点操作，此功能可以用树状数组差分维护）。
将所有模式串尾的结点 v i v_i vi​， [ d f n ( v i ) , d f s ( v i ) + s o n ( v i ) − 1 ] [dfn(v_i),dfs(v_i)+son(v_i)-1] [dfn(vi​),dfs(vi​)+son(vi​)−1]这一段+1,表示这个模式串作为后缀对其他结点的贡献。
询问时，逐位询问，每次询问到一个结点，得到的就是所有模式串可以作为当前后缀个数，每次都加上当前后缀模式串个数。
(一顿口胡
线段树：

/*   Author : Rshs
 *   Data : 2019-10-07-20.20
 */
#include<bits/stdc++.h>
using namespace std;
#define FI first
#define SE second
#define LL long long
#define MP make_pair
#define PII pair<int,int>
#define SZ(a) (int)a.size()
const double pai = acos(-1);
const double eps = 1e-10;
const LL mod = 1e9+7;
const int MXN = 1e6+5;
/*******segTree******************************************************************/
int lz[MXN<<2],T[MXN<<2];
inline void pushUp(int rt){
    T[rt]=T[rt<<1]+T[rt<<1|1];
}
void pushDown(int rt,int le){
    if(lz[rt]){
        lz[rt<<1]+=lz[rt];
        lz[rt<<1|1]+=lz[rt];
        T[rt<<1]+=lz[rt]*(le-(le>>1));
        T[rt<<1|1]+=lz[rt]*(le>>1);
        lz[rt]=0;
    }
}
void build(int l,int r,int rt){
    lz[rt]=0;
    if(l==r){
        T[rt]=0;
        return;
    }
    int m=(l+r)>>1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
    pushUp(rt);
}
void update(int L,int R,LL c,int l,int r,int rt){
    if(l>=L&&r<=R){
        lz[rt]+=c;
        T[rt]+=c*(LL)(r-l+1);
        return;
    }
    pushDown(rt,r-l+1);
    int m=(l+r)>>1;
    if(m>=L) update(L,R,c,l,m,rt<<1);
    if(m<R)  update(L,R,c,m+1,r,rt<<1|1);
    pushUp(rt);
}
LL treeQuery(int L,int R,int l,int r,int rt){
    if(l>=L&&r<=R)
        return T[rt];
    pushDown(rt,r-l+1);
    int m=(l+r)>>1;
    LL re=0;
    if(m>=L) re+=treeQuery(L,R,l,m,rt<<1);
    if(m<R)  re+=treeQuery(L,R,m+1,r,rt<<1|1);
    return re;
}
/*******segTree*********************************************************************/
/*******A C*****************************************************************/
int tr[MXN][26],col[MXN],fail[MXN],pos[MXN],nonum;
int dfn[MXN],cson[MXN],dfs_clock;
vector<int>g[MXN];
void ins(char *x,int ii){ //����
    int le=strlen(x),p=0;
    for(int i=0;i<le;i++){
        int c=x[i]-'a';
        if(!tr[p][c]) tr[p][c]=++nonum;
        p=tr[p][c];
    }
    pos[ii]=p;
    col[p]++; //�ַ���x����++
}
void ACGet(){
    queue<int>q;
    memset(fail,0,sizeof(fail));
    for(int i=0;i<26;i++)if(tr[0][i]) q.push(tr[0][i]);//���ַ����,ֱ�ӽ���(0)�����Ϊ�˱���ָ���Լ�
    while(!q.empty()){
        int now=q.front();q.pop();
        for(int i=0;i<26;i++){
            if(tr[now][i]){
                fail[tr[now][i]]=tr[fail[now]][i];//ָ��ǰһ��fail(��ͬ���׺)
                q.push(tr[now][i]);
            }
            else tr[now][i]=tr[fail[now]][i]; //ģʽ��β����һ���ߣ�����ƥ��
        }
    }
}
void ACDfs(int u,int pr){
    dfn[u]=++dfs_clock;
    cson[u]=1;
    for(auto v:g[u]){
        if(v==pr)continue;
        ACDfs(v,u);
        cson[u]+=cson[v];
    }
}
int ACQuery(char *t){
    int p=0;
    LL re=0;
    for(int i=0;t[i];i++){
        p=tr[p][t[i]-'a'];
        re+=treeQuery(dfn[p],dfn[p],1,dfs_clock,1);
    }
    return re;
}
/*******A C******************************************************************/
char sss[MXN];
int us[MXN];
int main(){
    int k,n;cin>>n>>k;
    for(int i=1;i<=k;i++){
        scanf("%s",sss);ins(sss,i);
    }
    ACGet();
    for(int i=1;i<=nonum;i++){
        g[fail[i]].push_back(i);g[i].push_back(fail[i]);
    }
    ACDfs(0,-1);
    for(int i=1;i<=nonum;i++){
        if(col[i]) update(dfn[i],dfn[i]+cson[i]-1,1,1,dfs_clock,1);
    }
    for(int i=1;i<=k;i++) us[i]=1;
    while(n--){
        char op;cin>>op;
        if(op=='+'){
            int x;scanf("%d",&x);
            if(us[x]==0){
                update(dfn[pos[x]],dfn[pos[x]]+cson[pos[x]]-1,1,1,dfs_clock,1);
                us[x]=1;
            }

        }
        if(op=='-'){
            int x;scanf("%d",&x);
            if(us[x]==1){
                update(dfn[pos[x]],dfn[pos[x]]+cson[pos[x]]-1,-1,1,dfs_clock,1);
                us[x]=0;
            }
        }
        if(op=='?'){
            scanf("%s",sss);
            LL ans=ACQuery(sss);
            cout<<ans<<'\n';
        }
    }
    return 0;
}

树状数组+差分：

/*   Author : Rshs
 *   Data : 2019-10-07-20.20
 */
#include<bits/stdc++.h>
using namespace std;
#define FI first
#define SE second
#define LL long long
#define MP make_pair
#define PII pair<int,int>
#define SZ(a) (int)a.size()
const double pai = acos(-1);
const double eps = 1e-10;
const LL mod = 1e9+7;
const int MXN = 1e6+5;
/*******segTree******************************************************************/
int bit[MXN];
void add(int i,LL c,int nn){ //��ai+c
    while(i<=nn){
        bit[i]+=c;
        i=i+(i&(-i));
    }
}
LL query(int i){ //ѯ��[l,r]�ĺͣ�query(r)-query(l-1)
    LL re=0;
    while(i>0){
        re+=bit[i];
        i=i-(i&(-i));
    }
    return re;
}
/*******segTree*********************************************************************/
/*******A C*****************************************************************/
int tr[MXN][26],col[MXN],fail[MXN],pos[MXN],nonum;
int dfn[MXN],cson[MXN],dfs_clock;
vector<int>g[MXN];
void ins(char *x,int ii){ //����
    int le=strlen(x),p=0;
    for(int i=0;i<le;i++){
        int c=x[i]-'a';
        if(!tr[p][c]) tr[p][c]=++nonum;
        p=tr[p][c];
    }
    pos[ii]=p;
    col[p]++; //�ַ���x����++
}
void ACGet(){
    queue<int>q;
    memset(fail,0,sizeof(fail));
    for(int i=0;i<26;i++)if(tr[0][i]) q.push(tr[0][i]);//���ַ����,ֱ�ӽ���(0)�����Ϊ�˱���ָ���Լ�
    while(!q.empty()){
        int now=q.front();q.pop();
        for(int i=0;i<26;i++){
            if(tr[now][i]){
                fail[tr[now][i]]=tr[fail[now]][i];//ָ��ǰһ��fail(��ͬ���׺)
                q.push(tr[now][i]);
            }
            else tr[now][i]=tr[fail[now]][i]; //ģʽ��β����һ���ߣ�����ƥ��
        }
    }
}
void ACDfs(int u,int pr){
    dfn[u]=++dfs_clock;
    cson[u]=1;
    for(auto v:g[u]){
        if(v==pr)continue;
        ACDfs(v,u);
        cson[u]+=cson[v];
    }
}
int ACQuery(char *t){
    int p=0;
    LL re=0;
    for(int i=0;t[i];i++){
        p=tr[p][t[i]-'a'];
        re+=query(dfn[p]);
    }
    return re;
}
/*******A C******************************************************************/
char sss[MXN];
int us[MXN];
int main(){
    int k,n;cin>>n>>k;
    for(int i=1;i<=k;i++){
        scanf("%s",sss);ins(sss,i);
    }
    ACGet();
    for(int i=1;i<=nonum;i++){
        g[fail[i]].push_back(i);g[i].push_back(fail[i]);
    }
    ACDfs(0,-1);
    for(int i=1;i<=nonum;i++){
        if(col[i]) add(dfn[i],1,dfs_clock),add(dfn[i]+cson[i],-1,dfs_clock);
    }
    for(int i=1;i<=k;i++) us[i]=1;
    while(n--){
        char op;cin>>op;
        if(op=='+'){
            int x;scanf("%d",&x);
            if(us[x]==0){
                add(dfn[pos[x]],1,dfs_clock),add(dfn[pos[x]]+cson[pos[x]],-1,dfs_clock);
                us[x]=1;
            }

        }
        if(op=='-'){
            int x;scanf("%d",&x);
            if(us[x]==1){
                add(dfn[pos[x]],-1,dfs_clock),add(dfn[pos[x]]+cson[pos[x]],1,dfs_clock);
                us[x]=0;
            }
        }
        if(op=='?'){
            scanf("%s",sss);
            LL ans=ACQuery(sss);
            cout<<ans<<'\n';
        }
    }
    return 0;
}
RSHS
关注
  • 1
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
DogLCD:树莓派 DogLCD 驱动程序库
06-04
介绍该库通过 4 线 SPI 串行接口（SI：串行输入，CLK：时钟，RS：寄存器选择和 CSB：芯片选择）支持带有 ST7036 驱动器 IC 的 DOGM081E、DOGM163E 和 DOGM183E 显示器。 您必须至少连接 SI、CLK 和 RS，但可以使用每...
delphi7系统五子棋游戏系统
05-06
5D5994295C5D945C5D945D5D5D5C945C9494949605CF05CFCFCF051062F1F15C 94295D5D94949495945D945D945D105D9494945D94945D94949495CDFECECBCD 5896942A942A2B620AFDFD059EF1FDF1CFCFFDF105FDCF05052BCD95CB965D94 ...
AC自动机算法模板
nEo.大魔王
07-08 820
PS：还不是很会，上午刚学，把仅学的一点先记录下来，肯定有错，学会了再来纠正（逃。 测试题目：http://acm.hdu.edu.cn/showproblem.php?pid=2222 想学自动AC机，忽然发现当初最开始学习时听到的一句话如今真的学到这里了：AC自动机fail树上dfs序建可持久化线段树。嘤 AC自动机 = KMP + Trie时间复杂度O（（n+m）*L） 具体思想，...
uva11481
lyc
01-27 820
此题纠结了我一段时间了。今天在acseed的提示下我ac了。。 首先题目很容易就可以想到首先算出前面C（m，k）种方案。那么就要考虑其他的n-k个位置都要不在其位置上。 那么我首先想到了错排。但是不知道为什么如果这么算的话是错误的。以5，3，2为例，答案是12，但是用错排算却是6错的。我就不知道为什么了。 求大神指教我上面的想法为什么错了啊。 然后被提醒说是容斥原理，首先第一步还是C（m，
Mike and Friends（CodeForces - 547E，AC自动机fail树dfs序上建可持久化线段树）
The___Flash的博客
05-06 1012
一、题目链接 Mike and Friends 二、题目大意 有 nnn 个字符串 s1,s2,⋯ ,sns_1,s_2,\cdots,s_ns1​,s2​,⋯,sn​ 和 qqq 次询问. 定义函数 f(S,T)f(S,T)f(S,T) 为字符串 SSS 在字符串 TTT 中的出现次数. 每次询问给出三个整数 l,r,kl,r,kl,r,k 表示询问 ∑i=lrf(si,sk)\sum_{i=l}^{r}{f(s_i,s_k)}∑i=lr​f(si​,sk​). 1≤n≤2×105,1≤q≤5×105,1
bzoj3881 [Coci2015]Divljak（AC自动机+fail树+dfs序+树状数组+树链剖分）
Bfk的博客
02-05 970
bzoj3881 [Coci2015]Divljak 原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=3881 题意： Alice有n个字符串S1,S2...Sn" role="presentation" style="position: relative;">S1,S2...SnS1,S2...SnS_1,S_2...S_n，Bob
【模板】AC自动机
Fighting_Peter的博客
01-24 291
ACM模板 目录插入以及构建AC自动机【模板】AC自动机（二次加强版）ac自动机fail树上dfs序建可持久化线段树 插入以及构建AC自动机 #include<queue> #include<string> const int N=200010; struct node { int chd[26],fail,cnt; }tree[N]; void insert(string s) { int p=0; for(int i=0;i<s.size();i+
洛谷p5726答案C语言,洛谷2414（构建ac自动机fail树dfs序后遍历Trie树维护bit及询问答案）...
weixin_35810223的博客
05-19 315
要点这是一道蔡队题，看我标题行事任意询问y串上有多少个x串，暴力找每个节点是不是结尾肯定是炸的，考虑本质：如果某节点是x的结尾，根据ac自动机的性质，x一定是此(子)串后缀。又有每个Trie节点的fail只指向另一个节点，故有fail树的概念。问题就变成了“对于串x的尾节点，在fail树中它的子树中有多少个点是在y串上”。解决方法是巧妙的。离线记录查询的信息。然后搜索原Trie树，遇到尾节点就扫描...
[bzoj2434][ac自动机][线段树合并][Noi2011]阿狸的打字机
Gzb1128的栖息地
02-07 249
2434: [Noi2011]阿狸的打字机 Time Limit: 10 Sec Memory Limit: 256 MB Submit: 3803 Solved: 2081 [Submit][Status][Discuss] Description 阿狸喜欢收藏各种稀奇古怪的东西，最近他淘到一台老式的打字机。打字机上只有28个按键，分别印有26个小写英文字母和’B’、’P’两个字母...
UVALive5092Permutation Counting(DP)
Rain
03-31 397
题目链接：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=3093 题目大意： 问1-n的数组成的排列，其中ai > i的个数恰好等于k的排列有多少种。 分析： dp[i][j]表示n个数组成的排列有多少满足条件k的个数
可持久化线段树&主席树、dfs序 Task Highways (aut)
SUNSET
08-12 378
Highways Memory limit: 128 MB Byteland is a small country, with  cities connected by  bidirectional roads. From each city there is just one way of reaching every other city using the road
HDU - 4117 GRE Words （AC自动机+Fail树+线段树+DFS序+DP）
LP_Cong
11-20 3245
GRE Words Problem Description Recently George is preparing for the Graduate Record Examinations (GRE for short). Obviously the most important thing is reciting the words. Now George is workin
CodeForces - 1437G Death DBMS(AC自动机fail树上树链剖分建线段树/暴跳fail)
Falcon的博客
10-28 539
题目链接：点击查看 题目大意：给出 n 个模式串，每个模式串初始时的权值为 0，然后有 m 次操作： 1 i x：将第 i 个模式串的权值修改为 x 2 s：给出一个字符串 s，询问字符串 s 作为主串时，所有可以与其匹配的模式串中，权值的最大值 题目分析：记录一下第一次单切 div2 的 G 题吧，不知道是 cf 题目质量变捞了还是我变强了，比赛结束前 20 分钟看到 G 出的比 F 还多，反正在挂机也是无聊，就开了 G 题，发现就是 AC 自动机的裸题，写了一发暴跳 fail 的代码，交上去果不其
BZOJ 2905: 背单词 AC自动机+fail树+dfs序+线段树
EM LGH
08-30 305
Description 给定一张包含N个单词的表，每个单词有个价值W。要求从中选出一个子序列使得其中的每个单词是后一个单词的子串，最大化子序列中W的和。 Input 第一行一个整数TEST，表示数据组数。接下来TEST组数据，每组数据第一行为一个整数N。接下来N行，每行为一个字符串和一个整数W。 Output TEST行，每行一个整数，表示...
AC自动机模板 （模板题 bzoj 4327）
执子之手,与子偕老
08-07 246
题目：这里写链接内容 可持久化版本待续 //AC自动机 //功能，一堆串和一个串进行匹配，看哪些位置可以被匹配。 //fail指针。当前串最长的后缀等于某个节点所表示的前缀。 //fail树。fail指针所形成的树。fail树上的所有祖先都是该节点的子串 //常见套路： //记录节点状态（当前匹配到哪个节点），DP，快速幂优化 //利用所有后缀的前缀是所有子串的性质 //直接利用fail树和t...
CF 163E. e-Government ac自动机+fail树+树状数组
weixin_30341745的博客
09-16 90
E. e-Government 题目： 　　给出n个字符串，表示n个人名，有两种操作： 　　?string ，统计字符串string中出现的属于城市居民的次数。 　　+id，把编号为id的人变为城市居民，如果已经是忽略。 　　-id，把编号为id的人变为不是城市居民，如果已经不是的话忽略。 　　现有m个操作，对于?输出结果。 分析： 　　很容易想到建立ac自动机，+-操作倒简...
P5357 【模板】AC自动机（二次加强版）（树上差分+AC自动机）
满月照亮的路
04-24 167
#include<bits/stdc++.h> using namespace std; #define debug puts("YES"); #define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++) #define ll long long #define lrt int l,int r,int rt #define lson l,m...
bzoj 3881: [Coci2015]Divljak AC自动机+树链的并+树状数组
beginend
11-08 368
题意Alice有n个字符串S_1,S_2…S_n，Bob有一个字符串集合T，一开始集合是空的。 接下来会发生q个操作，操作有两种形式： “1 P”，Bob往自己的集合里添加了一个字符串P。 “2 x”，Alice询问Bob，集合T中有多少个字符串包含串S_x。（我们称串A包含串B，当且仅当B是A的子串） Bob遇到了困难，需要你的帮助。 1 <= n,q <= 100000； Alice
cf 163e e-Government
lyc
10-02 1034
此题是那次做了GRE的题目后，相同做法的题。 先接着上次说的说说fail树吧：考虑把fail的指针反向建边，由于对于fail节点每个节点所指向的fail唯一，根据fail的原理，我们会从某个节点返回到root，所以fail的反向建边后性质就成了个fail树，这个树表示的是对于如果访问到某个子节点，那么他的祖先节点必然可以作为子串出现，那么鉴于此，我们要统计子串访问多少次，就是该节点返回到root
Fail树--学习笔记
niiick
02-28 2378
先看下面这个问题 BZOJ3172 [Tjoi2013]单词 Time Limit: 10 Sec Memory Limit: 512 MB Description 某人读论文，一篇论文是由许多单词组成。但他发现一个单词会在论文中出现很多次，现在想知道每个单词分别在论文中出现多少次。 Input 第一个一个整数N,表示有多少个单词，接下来N行每行一个单词。每个单词由小写字母组成，N&amp;amp;amp;amp;lt;=20...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值