P4770-[NOI2018]你的名字【SAM,线段树合并】
正题
题目链接:https://www.luogu.com.cn/problem/P4770
题目大意
给出一个长度为nnn的字符串SSS。qqq次询问给出一个串TTT和一个区间[L,R][L,R][L,R],求TTT有多少个本质不同的子串不是SL∼RS_{L\sim R}SL∼R的子串。
1≤n≤5×105,1≤Q≤105,∑∣T∣≤1061\leq n\leq 5\times 10^5,1\leq Q\le 10^5,\sum|T|\leq 10^61≤n≤5×105,1≤Q≤105,∑∣T∣≤106
解题思路
因为给了很多L=1,R=nL=1,R=nL=1,R=n的部分分所以应该是提示我们先从这个方面考虑。
这个部分比较简单,考虑改为求有多少个本质不同的子串在SL∼RS_{L\sim R}SL∼R中出现过,因为是本质不同的子串,我们可以先建一个SSS的SAMSAMSAM和一个TTT的SAMSAMSAM。
然后把TTT串拿到SSS的SAMSAMSAM上面跑,然后每次跑出来的一个匹配长度记为lenlenlen。对于TTT的SAMSAMSAM上的每一个节点我们记录一个pospospos表示这个节点属于的长度位置,然后跑到这个位置的lenlenlen就是能够匹配的长度了,记为ansansans。然后答案就是max{leni−max{ansposi,lenfai},0}max\{len_i-max\{ans_{pos_i},len_{fa_i}\},0\}max{leni−max{ansposi,lenfai},0}(防匹配长度超出[lenfai,leni][len_{fa_i},len_i][lenfai,leni]的范围)
这样一次的时间复杂度就是O(∣T∣)O(|T|)O(∣T∣)的了。
然后考虑带区间的怎么做,我们需要保证我们在SSS的SAMSAMSAM上面跳的时候需要保证这些节点都是属于SL∼RS_{L\sim R}SL∼R的自动机上的,而且要保证我们提取出来的lenlenlen也是在那个上面的。
其实如果这个节点的endposendposendpos类里面有L∼RL\sim RL∼R的信息就好了,这个用线段树合并维护一下endposendposendpos类的信息就可以了。
时间复杂度O((n+∑∣T∣)logn)O((n+\sum |T|)\log n)O((n+∑∣T∣)logn)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=1e6+10;
int n,q,ql,qr,p[N],c[N],rt[N],ans[N],pos[N];
char s[N];
struct SegTree{int cnt,w[N<<5],ls[N<<5],rs[N<<5];int Change(int x,int L,int R,int pos){int p=++cnt;w[p]=max(w[x],pos);if(L==R)return p;int mid=(L+R)>>1;if(pos<=mid)ls[p]=Change(ls[x],L,mid,pos),rs[p]=rs[x];else rs[p]=Change(rs[x],mid+1,R,pos),ls[p]=ls[x];return p;}int Ask(int x,int L,int R,int l,int r){if(!x)return 0;if(L==l&&R==r)return w[x];int mid=(L+R)>>1;if(r<=mid)return Ask(ls[x],L,mid,l,r);if(l>mid)return Ask(rs[x],mid+1,R,l,r);return max(Ask(ls[x],L,mid,l,mid),Ask(rs[x],mid+1,R,mid+1,r));}int Merge(int x,int y){if(!x||!y)return x+y;int p=++cnt;w[p]=max(w[x],w[y]);ls[p]=Merge(ls[x],ls[y]);rs[p]=Merge(rs[x],rs[y]);return p;}
}R;
struct SAM{int cnt,last,ch[N][26],len[N],fa[N];void init(){memset(ch[1],0,sizeof(ch[1]));last=cnt=1;}int Insert(int c){int p=last,np=last=++cnt;len[np]=len[p]+1;memset(ch[np],0,sizeof(ch[np]));for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;if(!p)fa[np]=1;else{int q=ch[p][c];if(len[p]+1==len[q])fa[np]=q;else{int nq=++cnt;len[nq]=len[p]+1;memcpy(ch[nq],ch[q],sizeof(ch[nq]));fa[nq]=fa[q];fa[q]=fa[np]=nq;pos[nq]=pos[q];for(;p&&ch[p][c]==q;p=fa[p])ch[p][c]=nq;}}last=np;return np;}void Build(){for(int i=1;i<=cnt;i++)c[len[i]]++;for(int i=1;i<=n;i++)c[i]+=c[i-1];for(int i=1;i<=cnt;i++)p[c[len[i]]--]=i;for(int i=cnt;i>=1;i--){int x=p[i];rt[fa[x]]=R.Merge(rt[fa[x]],rt[x]);}return;}void next(int &x,int &l,int c){while(x){if(ch[x][c]){int maxl=R.Ask(rt[ch[x][c]],1,n,1,qr)-ql+1;if(len[fa[x]]<maxl){l=min(l+1,maxl);x=ch[x][c];return;}}x=fa[x];l=len[x];}l=0;x=1;return;}ll calc(){ll prt=0;for(int i=2;i<=cnt;i++)prt+=max(len[i]-max(ans[pos[i]],len[fa[i]]),0);return prt;}
}S,T;
void work(char *s){int m=strlen(s+1);T.init();for(int i=1;i<=m;i++){int x=T.Insert(s[i]-'a');pos[x]=i;}int x=1,l=0;for(int i=1;i<=m;i++){int c=s[i]-'a';S.next(x,l,c);ans[i]=l;}printf("%lld\n",T.calc());
}
signed main()
{scanf("%s",s+1);n=strlen(s+1);S.init();for(int i=1;i<=n;i++){int x=S.Insert(s[i]-'a');rt[x]=R.Change(rt[x],1,n,i);}S.Build();scanf("%d",&q);while(q--){scanf("%s",s+1);scanf("%d%d",&ql,&qr);work(s);}return 0;
}
P4770-[NOI2018]你的名字【SAM,线段树合并】相关推荐
- CF1037H Security——SAM+线段树合并
又是一道\(SAM\)维护\(endpos\)集合的题,我直接把CF700E的板子粘过来了QwQ 思路 如果我们有\([l,r]\)对应的\(SAM\),只需要在上面贪心就可以了.因为要求的是字典序比 ...
- CF700E-Cool Slogans【SAM,线段树合并,dp】
正题 题目链接:https://www.luogu.com.cn/problem/CF700E 题目大意 给出一个字符串SSS,求一个最大的kkk使得存在kkk个字符串其中s1s_1s1是SSS的子 ...
- 洛谷P4482 [BJWC2018]Border 的四种求法 字符串,SAM,线段树合并,线段树,树链剖分,DSU on Tree...
原文链接https://www.cnblogs.com/zhouzhendong/p/LuoguP4482.html 题意 给定一个字符串 S,有 q 次询问,每次给定两个数 L,R ,求 S[L.. ...
- 2021CCPC华为云挑战赛:HDU 7091 重叠的子串(SAM + 线段树合并)
重叠的子串 给定一个长度为n(1≤∣s∣≤105)n(1 \le \mid s \mid \le 10 ^ 5)n(1≤∣s∣≤105)的只由小写字母构成的字符串sss,有m,(1≤m≤106)m, ...
- CF666E-Forensic Examination【广义SAM,线段树合并】
正题 题目链接:https://www.luogu.com.cn/problem/CF666E 解题思路 给出一个串SSS和nnn个串TiT_iTi.mmm次询问Sa∼bS_{a\sim b}Sa∼ ...
- YbtOJ#532-往事之树【广义SAM,线段树合并】
正题 题目链接:https://www.ybtoj.com.cn/problem/532 题目大意 给出nnn个点的一个TrieTrieTrie树,定义SxS_xSx表示节点xxx代表的字符串 求m ...
- CF204E-Little Elephant and Strings【广义SAM,线段树合并】
正题 题目链接:https://www.luogu.com.cn/problem/CF204E 题目大意 nnn个字符串的一个字符串集合,对于每个字符串求有多少个子串是这个字符串集合中至少kkk个字符 ...
- UOJ #395 BZOJ 5417 Luogu P4770 [NOI2018]你的名字 (后缀自动机、线段树合并)
NOI2019考前做NOI2018题.. 题目链接: (bzoj) https://www.lydsy.com/JudgeOnline/problem.php?id=5417 (luogu) http ...
- P4770:你的名字(SAM、线段树合并)
文章目录 前言 解析 前言 1000A快乐!!!awa 没有想象中的那么恶心. 解析 先考虑每次都询问 [1,n][1,n][1,n] 如何做. 正难则反,用T所有本质不同串数量减去是S串子串又是T的 ...
最新文章
- vue-i18n国际化实例
- python之装饰器篇
- Excel 下来公式 内容却一样
- 【随笔】工程师都是性情中人
- 三篇论文之Google MapReduce中文版
- filter在ios不兼容 js_jq的animate属性不支持transform 和 ios手机上box-shadow兼容性问题 和...
- Java 算法 学做菜
- 【HTML+CSS网页设计与布局 从入门到精通】第15章-表格格式:边距/居中/边框/折叠
- ndnsim r语言 ubuntu_Ubuntu14.04下配置ndnSIM-nom-rapid-car2car
- BZOJ3828 : [Poi2014]Criminals
- Ros代码迁移,编译头文件报错解决办法
- 刚装完系统的简单优化
- [swift 进阶]读书笔记-C2P3_Set、 C2P4_Range
- 和数传媒:区块链博弈论机制设计是未来关键
- AT24CXX、DS1302、红外HS0038 20170610 周六
- ThoughtWorks面试经历——武汉java开发
- 【IDE】IAR for ARM官网下载链接(包含一些历史版本)
- 浏览器工作原理(超级详解!!!!)
- chapter5. Java数据结构与java算法
- Android端恶意锁屏勒索应用分析
热门文章
- composer 查看php 版本_最常用的PHP版本:PHP 7.3取代7.2
- a*算法的时间复杂度_数据结构(1)——算法和时间复杂度
- java 线程访问控件_C#多线程与跨线程访问界面控件的方法
- python新闻评论分析_从新闻文章中提取评论
- java虚拟机工作原理图_Java虚拟机工作原理
- 技能高考本科计算机类,技能高考多少分上本科
- MBR的Linux分区机制启动过程,linux系统启动流程(MBR)
- 算法题目——读书知识点统计问题(POJ-3320)(尺取法)
- Millenium Leapcow POJ - 2111 (千禧年跳牛)(贪心找最长路径,记忆化)
- Zookeeper实践与应用--分布式锁实现