CF587F-Duff is Mad【AC自动机,根号分治】

正题

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


题目大意

给出nnn个字符串sssqqq次询问给出l,r,kl,r,kl,r,k要求输出sl..rs_{l..r}sl..rsks_ksk中出现了多少次。

1≤n,q,∑∣si∣≤1051\leq n,q,\sum |s_i|\leq 10^51n,q,si105


解题思路

考虑一个比较暴力的做法,先把所有的构出一棵ACACAC自动机,一个串SSS有后缀TTT当且仅当在failfailfail树上SSS的节点在TTT的子树内。

所以暴力的做法就是统计sks_ksk的节点有多少在sl..rs_{l..r}sl..r的终止节点子树内的节点。

需要优化,因为有L=∑∣si∣≤105L=\sum |s_i|\leq 10^5L=si105,所以可以考虑根号分治,设定值T=LT=\sqrt{L}T=L

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

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

对于长度小于等于TTTkkk,我们差分后按照询问端点排序,从左到右扫过nnn个串每次子树内的权值加一,然后每次暴力扫描sks_ksk的所有节点统计权值就好了。
使用树状数组时间复杂度O(nnlog⁡n)O(n\sqrt n\log n)O(nnlogn),已经可以通过本题了。但是如果把树状数组改成分块就可以平衡到O(nn)O(n\sqrt n)O(nn),但是我懒


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;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318609.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

牛客题霸 [最长重复子串] C++题解/答案

牛客题霸 [最长重复子串] C题解/答案 题目描述 一个重复字符串是由两个相同的字符串首尾拼接而成&#xff0c;例如abcabc便是长度为6的一个重复字符串&#xff0c;而abcba则不存在重复字符串。 给定一个字符串&#xff0c;请编写一个函数&#xff0c;返回其最长的重复字符子串…

codeforces1552 D. Array Differentiation(思维+暴力)

D. Array Differentiation 因为相减的顺序可以变化&#xff0c;所以这个环中数的正负性以及相减顺序其实是没有影响的&#xff0c;那么我们可以规定一个方向&#xff0c;然后再枚举所有可能的正负性。 设这环中点分别是 v1,v2,⋯vkv_1,v_2,⋯v_kv1​,v2​,⋯vk​&#xff0c;那…

粉刷木板(ybtoj-单调队列)

题目描述 解析 头疼 定义dp[i]:只用前i块板的最大价值 对于新加入的一个木匠&#xff1a; 不难写出dp转移式&#xff1a; dp[i]max(dp[k](i-k)*p)k表示开始刷的前一个 其中i>s i-k<l; 要是这么转移会是n^2m 所以我就不废了。。。 其实离正解很接近了 把上面移一下项&am…

P6563-[SBCOI2020]一直在你身旁【dp,单调队列】

正题 题目链接:https://www.luogu.com.cn/problem/P6563 题目大意 长度为nnn的序列aia_iai​&#xff0c;现在有一个随机[1,n][1,n][1,n]的整数&#xff0c;每次你可以花费aia_iai​询问这个数字是否大于iii&#xff0c;求猜出所有数至少要多少花费。 T≤500,∑n≤7000T\leq …

牛客题霸 [字符串的排列] C++题解/答案

牛客题霸 [字符串的排列] C题解/答案 题解&#xff1a; stl真好用&#xff0c;emmm。。。 代码&#xff1a; class Solution { public:vector<string> Permutation(string str) {if (str.empty()) return {};sort(str.begin(), str.end());vector<string>ans;i…

POJ - 3415 Common Substrings(长度不小于K的公共子串个数)

Common Substrings 后缀数组单调栈 题解1 题解2 题解3 #include<cstdio> #include<cstring> #include<iostream> using namespace std; typedef long long ll; // sa[i]: 排名是i位的是第几个后缀 // rk[i]: 第i个后缀的排名是多少 // height[i]: sa[i]与s…

跳房子(ybtoj-单调队列)

文章目录题目描述解析代码thanks for reading&#xff01;题目描述 洛谷传送门 跳房子&#xff0c;也叫跳飞机&#xff0c;是一种世界性的儿童游戏&#xff0c;也是中国民间传统的体育游戏之一。 跳房子的游戏规则如下&#xff1a; 在地面上确定一个起点&#xff0c;然后在起…

使用Dapper持久化IdentityServer4

最近研究dotnet core,微软将IdentityServer4作为推荐的服务授权和验证的组件,其独立性特别适合微服务或者分布式的服务扩展验证,所以非常受广大dotnet开发人员的青睐.默认的IdentityServer4默认使用内存对象的验证和授权,而在IdentityServer的官方推荐只有Entity Framework cor…

P5470-[NOI2019]序列【模拟费用流】

正题 题目链接:https://www.luogu.com.cn/problem/P5470 题目大意 两个长度为nnn的序列a,ba,ba,b&#xff0c;求出它们两个长度为KKK的子序列&#xff0c;且这两个子序列至少有LLL个位置下标相等。 求最大化两个子序列的和。 T≤10,1≤n≤2105,∑n≤106T\leq 10,1\leq n\leq…

牛客题霸 [数组中未出现的最小正整数] C++题解/答案

牛客题霸 [数组中未出现的最小正整数] C题解/答案 题目描述 给定一个无序数组arr&#xff0c;找到数组中未出现的最小正整数 例如arr [-1, 2, 3, 4]。返回1 arr [1, 2, 3, 4]。返回5 [要求] 时间复杂度为O(n)O(n)&#xff0c;空间复杂度为O(1)O(1) 题解&#xff1a; 如果…

再不学习我们就out了

前不久我们组来了个Graduate Developer&#xff0c;刚毕业&#xff0c;经验不多&#xff0c;有一次闹了个乌龙&#xff0c;把Stage数据库直接删掉了……好在Azure有备份&#xff0c;不然就算Stage没有重要数据&#xff0c;也得花点时间重建&#xff0c;其他的开发、测试工作都得…

矩阵快速幂一篇通

文章目录概述快速幂解析代码矩阵运算定义加法乘法单位矩阵一、斐波拉契&#xff08;基础模板&#xff09;题目描述解析代码二、行为方案&#xff08;实际应用&#xff09;题目描述解析代码三、矩阵求和&#xff08;子矩阵作为矩阵元素&#xff09;题目描述解析代码四、最短路径…

SPOJ687 Repeats(重复次数最多的连续子串)

Repeats hihoCoder 1419 后缀数组四重复旋律4(重复次数最多的连续子串) #include<bits/stdc.h> using namespace std;// sa[i]: 排名是i位的是第几个后缀 // rk[i]: 第i个后缀的排名是多少 // height[i]: sa[i]与sa[i-1] const int N50010; char s[N]; int rk[N],sa[N…

玉米田(加加强版)【插头dp】

前言 水解警告&#xff0c;数据水勉强卡过的 正题 题目大意 n∗mn*mn∗m的网格里面有些格子被禁止&#xff0c;现在求选取若干个不相邻的格子的方案数。 1≤n≤120,1≤m≤211\leq n\leq 120,1\leq m\leq 211≤n≤120,1≤m≤21 解题思路 听说是插头dpdpdp然后想了一下觉得比插…

牛客题霸 [将字符串转化为整数] C++题解/答案

牛客题霸 [将字符串转化为整数] C题解/答案 题目描述 实现函数 atoi 。函数的功能为将字符串转化为整数 提示&#xff1a;仔细思考所有可能的输入情况。这个问题没有给出输入的限制&#xff0c;你需要自己考虑所有可能的情况。 题解&#xff1a; 题目很简单&#xff0c;但是…

用ABP入门DDD

前言ABP框架一直以来都是用DDD&#xff08;领域驱动设计&#xff09;作为宣传点之一。但是用过ABP的人都知道&#xff0c;ABP并不是一个严格遵循DDD的开发框架&#xff0c;又或者说&#xff0c;它并没有完整实现DDD的所有概念。但是反过来说&#xff0c;认真学过DDD的人会发现&…

多重背包的二进制优化(ybtoj-宝物筛选)

文章目录题目描述解析朴素算法代码二进制优化代码thanks for reading!题目描述 解析 朴素算法 首先考虑朴素算法 把数量为num的物体拆成num个子物体 其价值与重量是原物体的1&#xff0c;2&#xff0c;3…num倍 然后当成独立的物体求就行了 注意应该先枚举重量&#xff0c;再…

P6800-[模板]Chirp Z-Transform【NTT】

正题 题目链接:https://www.luogu.com.cn/problem/P6800 题目大意 给出一个nnn此多项式PPP&#xff0c;对于k∈[0,m−1]k\in[0,m-1]k∈[0,m−1]所有的求P(ck)P(c^k)P(ck) 输出答案对998244353998244353998244353取模 1≤n,m≤1061\leq n,m\leq 10^61≤n,m≤106 解题思路 g(n)…

AtCoder Grand Contest 021 D - Reversed LCS(区间dp)

D - Reversed LCS 繁凡さん 设 f[l,r,k]f [ l , r , k ]f[l,r,k] 表示区间 [l,r][ l , r ][l,r] 中修改 kkk 次能得到的最长回文子序列的长度&#xff0c;直接区间DP转移即可。 字符串的 最长回文子序列&#xff08;lps&#xff09; 长度等于其自身与反转的 最长公共子序列&am…

牛客题霸 [滑动窗口的最大值] C++题解/答案

牛客题霸 [滑动窗口的最大值] C题解/答案 题目描述 给定一个数组和滑动窗口的大小&#xff0c;找出所有滑动窗口里数值的最大值。例如&#xff0c;如果输入数组{2,3,4,2,6,2,5,1}及滑动窗口的大小3&#xff0c;那么一共存在6个滑动窗口&#xff0c;他们的最大值分别为{4,4,6,…