正题

题目链接:https://www.luogu.com.cn/problem/CF587F


题目大意

给出\(n\)个字符串\(s\)。\(q\)次询问给出\(l,r,k\)要求输出\(s_{l..r}\)在\(s_k\)中出现了多少次。

\(1\leq n,q,\sum |s_i|\leq 10^5\)


解题思路

考虑一个比较暴力的做法,先把所有的构出一棵\(AC\)自动机,一个串\(S\)有后缀\(T\)当且仅当在\(fail\)树上\(S\)的节点在\(T\)的子树内。

所以暴力的做法就是统计\(s_k\)的节点有多少在\(s_{l..r}\)的终止节点子树内的节点。

需要优化,因为有\(L=\sum |s_i|\leq 10^5\),所以可以考虑根号分治,设定值\(T=\sqrt{L}\)。

对于长度大于\(T\)的\(k\),我们所有的询问统一差分变为前缀问题处理。因为这样的串的个数不会超过\(\sqrt L\)个,所以可以每次暴力\(O(n)\)处理。

我们把所有\(s_k\)路径上的节点权值加一,然后暴力扫描\(n\)个串,每次终止节点统计子树内的权值和。

对于长度小于等于\(T\)的\(k\),我们差分后按照询问端点排序,从左到右扫过\(n\)个串每次子树内的权值加一,然后每次暴力扫描\(s_k\)的所有节点统计权值就好了。

使用树状数组时间复杂度\(O(n\sqrt n\log n)\),已经可以通过本题了。但是如果把树状数组改成分块就可以平衡到\(O(n\sqrt n)\),但是我懒


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<cmath>
#define ll long long
#define mp(x,y) make_pair(x,y)
#define lowbit(x) (x&-x)
using namespace std;
const ll N=1e5+10;
struct node{
ll to,next;
}a[N];
ll n,m,T,tot,dfr,L,top[N],ls[N],l[N],t[N],ans[N];
ll cnt,ch[N][26],fail[N],pos[N],rfn[N],ed[N];
vector<pair<ll,ll> > v[N],V[N];
char *s[N],st[N];queue<ll> q;
void Change(ll x,ll val){
while(x<=dfr){
t[x]+=val;
x+=lowbit(x);
}
return;
}
ll Ask(ll x){
ll ans=0;
while(x){
ans+=t[x];
x-=lowbit(x);
}
return ans;
}
void addl(ll x,ll y){
a[++tot].to=y;
a[tot].next=ls[x];
ls[x]=tot;return;
}
ll Insert(char *s,ll l){
ll x=0;
for(ll i=0;i<l;i++){
ll c=s[i]-'a';
if(!ch[x][c])ch[x][c]=++cnt;
x=ch[x][c];
}
return x;
}
void GetFail(){
for(ll i=0;i<26;i++)
if(ch[0][i])q.push(ch[0][i]);
while(!q.empty()){
ll x=q.front();q.pop();
for(ll i=0;i<26;i++){
if(!ch[x][i])ch[x][i]=ch[fail[x]][i];
else{
fail[ch[x][i]]=ch[fail[x]][i];
q.push(ch[x][i]);
}
}
}
for(ll i=1;i<=cnt;i++)addl(fail[i],i);
return;
}
void dfs(ll x,ll fa){
rfn[x]=++dfr;
for(ll i=ls[x];i;i=a[i].next){
ll y=a[i].to;
if(y==fa)continue;
dfs(y,x);
}
ed[x]=dfr;
return;
}
signed main()
{
scanf("%lld%lld",&n,&m);
for(ll i=1;i<=n;i++){
s[i]=st+top[i];
scanf("%s",s[i]);l[i]=strlen(s[i]);
top[i+1]=top[i]+l[i];pos[i]=Insert(s[i],l[i]);
}
L=top[n+1];T=sqrt(L);
GetFail();
dfs(0,0);
for(ll i=1;i<=m;i++){
ll L,R,K;scanf("%lld%lld%lld",&L,&R,&K);
if(l[K]<=T)v[R].push_back(mp(K,i)),v[L-1].push_back(mp(K,-i));
else V[K].push_back(mp(R,i)),V[K].push_back(mp(L-1,-i));
}
for(ll p=1;p<=n;p++)
if(l[p]>T){
memset(t,0,sizeof(t));
sort(V[p].begin(),V[p].end());
ll x=0;
for(ll i=0;i<l[p];i++){
x=ch[x][s[p][i]-'a'];
Change(rfn[x],1);
}
ll z=0,sum=0;
while(z<V[p].size()&&!V[p][z].first)z++;
for(ll i=1;i<=n;i++){
sum+=Ask(ed[pos[i]])-Ask(rfn[pos[i]]-1);
while(z<V[p].size()&&V[p][z].first<=i){
ll id=V[p][z].second,op=1;
if(id<0)id=-id,op=-op;
ans[id]+=sum*op;z++;
}
}
}
memset(t,0,sizeof(t));
for(ll p=1;p<=n;p++){
Change(rfn[pos[p]],1);
Change(ed[pos[p]]+1,-1);
for(ll i=0;i<v[p].size();i++){
ll k=v[p][i].first,id=v[p][i].second,op=1;
if(id<0)id=-id,op=-op;ll x=0;
for(ll j=0;j<l[k];j++){
x=ch[x][s[k][j]-'a'];
ans[id]+=op*Ask(rfn[x]);
}
}
}
for(ll i=1;i<=m;i++)
printf("%lld\n",ans[i]);
return 0;
}

CF587F-Duff is Mad【AC自动机,根号分治】的更多相关文章

  1. [CF587F]Duff is Mad[AC自动机+根号分治+分块]

    题意 给你 \(n\) 个串 \(s_{1\cdots n}\) ,每次询问给出 \(l,r,k\) ,问在 \(s_{l\cdots r}\) 中出现了多少次 \(s_k\) . \(n,q,\su ...

  2. 【CF587F】Duff is Mad AC自动机+分块

    [CF587F]Duff is Mad 题意:给出n个串$s_1,s_2..s_n$,有q组询问,每次给出l,r,k,问你编号在[l,r]中的所有串在$s_k$中出现了多少次. $\sum|s_i|, ...

  3. CF587F Duff is Mad(AC自动机+树状数组+分块)

    考虑两一个暴力 1 因为询问\([a,b]\)可以拆成\([1,b]\)-\([1,a-1]\)所以把询问离线,然后就是求\([1,x]\)中被\(S_i\)包含的串的数量.考虑当\([1,x-1]- ...

  4. CF587F Duff is Mad

    题目 有趣的思想 首先暴力的话,自然是对每一个询问在\(AC\)自动机上跑一遍\(k\),看看跑出来的节点在\(fail\)树到根的路径上有多少个\(l\)到\(r\)之间的结束标记就好了 我们发现无 ...

  5. Codeforces 587F - Duff is Mad(根号分治+AC 自动机+树状数组)

    题面传送门 第一眼看成了 CF547E-- 话说 CF587F 和 CF547E 出题人一样欸--还有另一道 AC 自动机的题 CF696D 也是这位名叫 PrinceOfPersia 的出题人出的- ...

  6. NOI.AC#2266-Bacteria【根号分治,倍增】

    正题 题目链接:http://noi.ac/problem/2266 题目大意 给出\(n\)个点的一棵树,有一些边上有中转站(边长度为\(2\),中间有一个中转站),否则就是边长为\(1\). \( ...

  7. 【XSY3320】string AC自动机 哈希 点分治

    题目大意 给一棵树,每条边上有一个字符,求有多少对 \((x,y)(x<y)\),满足 \(x\) 到 \(y\) 路径上的边上的字符按顺序组成的字符串为回文串. \(1\leq n\leq 5 ...

  8. Codeforces 710F - String Set Queries(AC 自动机)

    题面传送门 题意:强制在线的 AC 自动机. \(n,\sum|s|\leq 3\times 10^5\) 如果不是强制在线那此题就是道 sb 题,加了强制在线就不那么 sb 了. 这里介绍两种做法: ...

  9. 【题解】Tree-String Problem Codeforces 291E AC自动机

    Prelude 传送到Codeforces:(/ω\)--- (/ω•\) Solution 很水的一道题. 对查询的串建出来AC自动机,然后树上随便跑跑就行了. 为什么要写这篇题解呢? 我第一眼看到 ...

随机推荐

  1. commandBinding 的命令

    <Window x:Class="WpfApplication1.Window29" xmlns="http://schemas.microsoft.com/win ...

  2. Python中的reduce()函数

    reduce()函数也是Python内置的一个高阶函数.reduce()函数接收的参数和 map()类似,一个函数 f,一个list,但行为和 map()不同,reduce()传入的函数 f 必须接收 ...

  3. PipedInputStream and PipedOutputStream example

    必须要有PipedInputStream ,PipedOutputStream 在不同线程,不然死锁 Java.io.PipedOutputStream and java.io.PipedInputS ...

  4. Jemeter压力测试

    Jmeter教程 简单的压力测试 Jmeter是一个非常好用的压力测试工具.  Jmeter用来做轻量级的压力测试,非常合适,只需要十几分钟,就能把压力测试需要的脚本写好. 转载自小坦克:https: ...

  5. vue@cli3 public目录下的静态图片,如何使用在css类文件中(sass可行,纯css不行)

    之前写了一篇vue文件怎么使用的文章,有人问我怎么在css文件中使用public下的文件,这是个好问题,因为我之前都没有研究过 需要解决的2个问题 一开始按照vue文件的使用方式(https://ww ...

  6. LeetCode《买卖股票的最佳时机》系列题目,最详解

    目录 说在前面 引例:只能交易一次 一.动态数组定义 二.状态转移方程 三.初始化 四.优化 无限制买卖 一.动态数组定义 二.状态转移方程 三.初始化 四.优化 交易 2 次,最大利润? 一.动态数 ...

  7. J2EE之DAO设计模式及简单实现

    JAVAEE(Java Enterprise Edition ) 模式 : DAO模式 因此在了解DAO模式之前,我们先来学习一下Java EE的体系结构: (一)JavaEE体系结构 客户端: 客户 ...

  8. 基于Linux系统的网络服务——高速缓存DNS及企业级域名解析服务

    1.DNS域名系统 DNS(Domain Name System,域名系统),因特网上作为域名和IP地址相互映射的一个分布式数据库,能够使用户更方便的访问互联网,而不用去记住能够被机器直接读取的IP数 ...

  9. Java 学习:对象和类

    对象和类 从认识的角度考虑是先有对象后有类.对象,是具体的事物.类,是抽象的,是对对象的抽象. 从代码运行角度考虑是先有类后又对象.类是对象的模板. 对象:对象是类的一个实例,有状态和行为. 类:类是 ...

  10. gimp 缩放图片 python script

    滤镜 -> Python Fu -> 控制台,复制粘贴下面代码. 然后 scale_image(800, 800) 图片自动缩放为 800*800 了 这个代码是傻傻的缩放,你可以把它改成 ...