BZOJ3881[Coci2015]Divljak——AC自动机+树状数组+LCA+dfs序+树链的并
题目描述
输入
输出
样例输入
a
bc
abc
5
1 abca
2 1
1 bca
2 2
2 3
样例输出
2
1
提示
这道题和bzoj2434都是很好的AC自动机练习题,都利用了fail树的性质来解决字符串问题。首先要知道一点:如果x串是y串的子串,那x串一定是y串一个前缀的后缀。我们知道AC自动机上每个点表示这个点到根节点的字符串。这道题的大体思路是把S集合建成AC自动机,每添加一个T集合中的字符串P就把他在AC自动机上跑一遍,把所有遍历的点及它们fail指针能达到的所有点的答案数都加1(遍历的点自然是这个串的子串,fail指针能达到的点都是遍历的点的后缀,当然也是加入的串的子串),然后查询。利用fail树的性质可以发现,每次插入一个串P,除了遍历的点,其他答案要加1的点都是这些点在fail树上的祖先,因此只要用树上差分将遍历的点+1就把所有答案要加的点都加了。但这样有的点就会被多加好几次,每次插入却最多只能把每个点加一次,因此,要把遍历的点按dfs序排序,然后将相邻两个点的lca到根节点答案都-1(也就是用差分将lca-1)就去重了(因为按dfs序相邻两个点最近,lca也最深)。查询时用树状数组在dfs序上区间求和就行了。
最后附上代码。
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m;
int sum;
int num;
int cnt;
int tot;
int opt;
int g[100010];
int k[2000010];
int v[2000010];
int l[2000010];
int r[2000010];
int d[2000010];
char s[2000010];
int to[2000010];
int head[2000010];
int next[2000010];
int fail[2000010];
int f[23][2000010];
int a[2000010][26];
void add(int x,int y)
{
tot++;
next[tot]=head[x];
head[x]=tot;
to[tot]=y;
}
void change(int x,int val)
{
for(int i=x;i<=num;i+=i&-i)
{
v[i]+=val;
}
}
int ask(int x)
{
int res=0;
for(int i=x;i;i-=i&-i)
{
res+=v[i];
}
return res;
}
void build(char *s,int c)
{
int now=0;
int len=strlen(s);
for(int i=0;i<len;i++)
{
int x=s[i]-'a';
if(!a[now][x])
{
a[now][x]=++cnt;
}
now=a[now][x];
}
g[c]=now;
}
void getfail()
{
queue<int>q;
for(int i=0;i<26;i++)
{
if(a[0][i])
{
q.push(a[0][i]);
fail[a[0][i]]=0;
}
}
while(!q.empty())
{
int now=q.front();
q.pop();
add(fail[now],now);
f[0][now]=fail[now];
for(int i=0;i<26;i++)
{
if(a[now][i])
{
fail[a[now][i]]=a[fail[now]][i];
q.push(a[now][i]);
}
else
{
a[now][i]=a[fail[now]][i];
}
}
}
}
bool cmp(int x,int y)
{
return l[x]<l[y];
}
void dfs(int x)
{
l[x]=++num;
d[x]=d[f[0][x]]+1;
for(int i=1;i<=21;i++)
{
f[i][x]=f[i-1][f[i-1][x]];
}
for(int i=head[x];i;i=next[i])
{
dfs(to[i]);
}
r[x]=num;
}
int lca(int x,int y)
{
if(d[x]<d[y])
{
swap(x,y);
}
int dep=d[x]-d[y];
for(int i=0;i<=21;i++)
{
if(((1<<i)&dep)!=0)
{
x=f[i][x];
}
}
if(x==y)
{
return x;
}
for(int i=21;i>=0;i--)
{
if(f[i][x]!=f[i][y])
{
x=f[i][x];
y=f[i][y];
}
}
return f[0][x];
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%s",s);
build(s,i);
}
getfail();
dfs(0);
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
scanf("%d",&opt);
if(opt==1)
{
scanf("%s",s);
int len=strlen(s);
int now=0;
int x=0;
for(int j=0;j<len;j++)
{
now=a[now][s[j]-'a'];
k[++x]=now;
}
sort(k+1,k+1+x,cmp);
for(int j=1;j<=x;j++)
{
change(l[k[j]],1);
} for(int j=1;j<x;j++)
{
change(l[lca(k[j],k[j+1])],-1);
}
}
else
{
scanf("%d",&sum);
printf("%d\n",ask(r[g[sum]])-ask(l[g[sum]]-1));
}
}
}
BZOJ3881[Coci2015]Divljak——AC自动机+树状数组+LCA+dfs序+树链的并的更多相关文章
- Codeforces 1111E DP + 树状数组 + LCA + dfs序
题意:给你一颗树,有q次询问,每次询问给你若干个点,这些点可以最多分出m组,每组要满足两个条件:1:每组至少一个点,2:组内的点不能是组内其它点的祖先,问这样的分组能有多少个? 思路:https:// ...
- BZOJ 3881 [COCI2015]Divljak (Trie图+Fail树+树链的并+树状数组维护dfs序)
题目大意: Alice有n个字符串S_1,S_2...S_n,Bob有一个字符串集合T,一开始集合是空的. 接下来会发生q个操作,操作有两种形式: “1 P”,Bob往自己的集合里添加了一个字符串P. ...
- (好题)树状数组+离散化+DFS序+离线/莫队 HDOJ 4358 Boring counting
题目传送门 题意:给你一棵树,树上的每个节点都有树值,给m个查询,问以每个点u为根的子树下有多少种权值恰好出现k次. 分析:首先要对权值离散化,然后要将树形转换为线形,配上图:.然后按照右端点从小到大 ...
- 【BZOJ-3881】Divljak AC自动机fail树 + 树链剖分+ 树状数组 + DFS序
3881: [Coci2015]Divljak Time Limit: 20 Sec Memory Limit: 768 MBSubmit: 508 Solved: 158[Submit][Sta ...
- 【BZOJ4999】This Problem Is Too Simple! 离线+树状数组+LCA
[BZOJ4999]This Problem Is Too Simple! Description 给您一颗树,每个节点有个初始值. 现在支持以下两种操作: 1. C i x(0<=x<2 ...
- 【BZOJ】2819: Nim(树链剖分 / lca+dfs序+树状数组)
题目 传送门:QWQ 分析 先敲了个树链剖分,发现无法AC(其实是自己弱,懒得debug.手写栈) 然后去学了学正解 核心挺好理解的,$ query(a) $是$ a $到根的异或和. 答案就是$ l ...
- Ultra-QuickSort(树状数组求逆序对数)
Ultra-QuickSort 题目链接:http://poj.org/problem?id=2299 Time Limit: 7000MS Memory Limit: 65536K Total ...
- hdu4605 树状数组+离散化+dfs
Magic Ball Game Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- hdu2838树状数组解逆序
离散化和排序后的序号问题搞得我实在是头痛 不过树状数组解逆序和偏序一类问题真的好用 更新:hdu的数据弱的真实,我交上去错的代价也对了.. 下面的代码是错的 /* 每个点的贡献度=权值*在这个点之前的 ...
随机推荐
- 3分钟学会做智能插座(DIY)
转载请注明:@小五义http://www.cnblogs.com/xiaowuyiQQ群:64770604 感谢博达科技提供的技术支持,博达科技新出了turnip智能插座,通过微信控制,实现了语音控制 ...
- SkylineGlobe 邻近度(Proximity)分析JavaScript源代码
邻近度(Proximity)描述了地理空间中两个地物距离相近的程度,是空间分析的一个重要手段. <html xmlns="http://www.w3.org/1999/xhtml&qu ...
- BZOJ1758 WC2010 重建计划 二分答案、点分治、单调队列
传送门 看到平均数最大,自然地想到二分答案.那么我们的$check$函数就是要求:是否存在一条长度在$[L,U]$的路径,满足其权值和$\geq 0$. 看到长度在$[L,U]$,自然地想到点分治求解 ...
- Mac安装使用MongoDB
Mac 下安装 MongoDB 一般有两种方法,一种是通过源码安装,一种是直接使用 homebrew ,个人推荐使用 homebrew ,简单粗暴. 一.安装 homebrew : /usr/bin/ ...
- Luogu P2827 蚯蚓
看到题目就可以想到直接开的堆模拟的过程了吧,这个还是很naive的 注意在用堆做的时候也是要明智一点的,对于蚯蚓长度的相加肯定不能直接遍历并加上,还是可以差分一下的 其实说白了就是把集体加->单 ...
- springboot @Value 获取计算机中绝对路径文件的内容
默认情况下使用 @Value("aaa.txt") private Resource txtResource; 这样获取到的是项目classpath 下的 aaa.txt 如果想获 ...
- SpringBoot日记——Spring的安全配置-登录认证与授权
安全是每个项目开发中都需要考虑的,比如权限控制,安全认证,防止漏洞攻击等. 比较常见的安全框架有:Apache的shiro.Spring Security等等,相信用shiro的用户群体更多,而sec ...
- Jmeter(三十三)_JsonPath表达式提取响应
我们在用jmeter做接口测试的时候,有的时候会遇到一些复杂的json响应.比如多层list嵌套时的取值 一个简单的例子: $..Name:列出所有省份 $..Province[0].Name 提取P ...
- centos7.2部署vnc服务记录
不做过多介绍了,下面直接记录下centos7系统下安装配置vncserver的操作记录 0)更改为启动桌面或命令行模式 获取当前系统启动模式 [root@localhost ~]# systemctl ...
- div+css实现圆形div以及带箭头提示框效果
.img{ width:90px; height:90px; border-radius:45px; margin:0 40%; border:solid rgb(100,100,100) 1px;& ...