正题

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


题目大意

给出一个长度为\(n\)的字符串\(S\)。\(q\)次询问给出一个串\(T\)和一个区间\([L,R]\),求\(T\)有多少个本质不同的子串不是\(S_{L\sim R}\)的子串。

\(1\leq n\leq 5\times 10^5,1\leq Q\le 10^5,\sum|T|\leq 10^6\)


解题思路

因为给了很多\(L=1,R=n\)的部分分所以应该是提示我们先从这个方面考虑。

这个部分比较简单,考虑改为求有多少个本质不同的子串在\(S_{L\sim R}\)中出现过,因为是本质不同的子串,我们可以先建一个\(S\)的\(SAM\)和一个\(T\)的\(SAM\)。

然后把\(T\)串拿到\(S\)的\(SAM\)上面跑,然后每次跑出来的一个匹配长度记为\(len\)。对于\(T\)的\(SAM\)上的每一个节点我们记录一个\(pos\)表示这个节点属于的长度位置,然后跑到这个位置的\(len\)就是能够匹配的长度了,记为\(ans\)。然后答案就是\(max\{len_i-max\{ans_{pos_i},len_{fa_i}\},0\}\)(防匹配长度超出\([len_{fa_i},len_i]\)的范围)

这样一次的时间复杂度就是\(O(|T|)\)的了。

然后考虑带区间的怎么做,我们需要保证我们在\(S\)的\(SAM\)上面跳的时候需要保证这些节点都是属于\(S_{L\sim R}\)的自动机上的,而且要保证我们提取出来的\(len\)也是在那个上面的。

其实如果这个节点的\(endpos\)类里面有\(L\sim R\)的信息就好了,这个用线段树合并维护一下\(endpos\)类的信息就可以了。

时间复杂度\(O((n+\sum |T|)\log n)\)


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,线段树合并】的更多相关文章

  1. [NOI2018]你的名字(SAM+线段树合并)

    考虑l=1,r=n的68分,对S和T建SAM,对T的SAM上的每个节点,计算它能给答案带来多少贡献. T上节点x代表的本质不同的子串数为mx[x]-mx[fa[x]],然后需要去掉所代表子串与S的最长 ...

  2. NOI2018 你的名字——SAM+线段树合并

    题目链接在这里洛谷/LOJ 题目大意 有一个串\(S\),每次询问给你一个串\(T\),两个数\(L\)和\(R\),问你\(T\)有多少个本质不同的子串不是\(S[L,R]\)的子串 SOLUTIO ...

  3. 【BZOJ5417】[NOI2018]你的名字(线段树,后缀自动机)

    [BZOJ5417][NOI2018]你的名字(线段树,后缀自动机) 题面 BZOJ 洛谷 题解 首先考虑\(l=1,r=|S|\)的做法,对于每次询问的\(T\)串,暴力在\(S\)串的\(SAM\ ...

  4. 【NOI2018】你的名字(SAM & 线段树合并)

    Description Hint Solution 不妨先讨论一下无区间限制的做法. 首先"子串"可以理解为"前缀的后缀",因此我们定义一个 \(\lim(i) ...

  5. UOJ#395. 【NOI2018】你的名字 字符串,SAM,线段树合并

    原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ395.html 题解 记得同步赛的时候这题我爆0了,最暴力的暴力都没调出来. 首先我们看看 68 分怎么做 ...

  6. luogu4770 [NOI2018]你的名字 (SAM+主席树)

    对S建SAM,拿着T在上面跑 跑的时候不仅无法转移要跳parent,转移过去不在范围内也要跳parent(注意因为范围和长度有关,跳的时候应该把长度一点一点地缩) 这样就能得到对于T的每个前缀,它最长 ...

  7. CF1037H Security——SAM+线段树合并

    又是一道\(SAM\)维护\(endpos\)集合的题,我直接把CF700E的板子粘过来了QwQ 思路 如果我们有\([l,r]\)对应的\(SAM\),只需要在上面贪心就可以了.因为要求的是字典序比 ...

  8. 洛谷P4482 [BJWC2018]Border 的四种求法 字符串,SAM,线段树合并,线段树,树链剖分,DSU on Tree

    原文链接https://www.cnblogs.com/zhouzhendong/p/LuoguP4482.html 题意 给定一个字符串 S,有 q 次询问,每次给定两个数 L,R ,求 S[L.. ...

  9. Codeforces 700E. Cool Slogans 字符串,SAM,线段树合并,动态规划

    原文链接https://www.cnblogs.com/zhouzhendong/p/CF700E.html 题解 首先建个SAM. 一个结论:对于parent树上任意一个点x,以及它所代表的子树内任 ...

随机推荐

  1. 异步编程async体会

    namespace 异步编程{ class Program { static void Main(string[] args) { Console.WriteLine("mian this ...

  2. QT5学习:分割窗口类的使用

    分割窗口在应用程序中经常用到,它可以灵活分布窗口布局,经常用于类似文件资源管理器的窗口设计中,然后抱着这样的想法简单的实现了下 [cpp]  view plain copy   //main.cpp ...

  3. Semaphore 类 的使用理解C#

    示例 下面的代码示例创建一个信号量,其最大计数为3,初始计数为零. 该示例启动五个线程,这会阻止等待信号量. 主线程使用 Release(Int32) 方法重载将信号量计数增加到其最大值,从而允许三个 ...

  4. uwp 的work project 的 取消闹钟

    private void initalAlarmHanle() { string cancelAlarm = "CancelAlarmEvent"; ConnectionManag ...

  5. 关于Ubuntu18.04 linux系统下使用Tim QQ 微信

    先配上张图 步骤: 1.1 :需要安装环境deepin-wine 1.1:(你把他理解为jdk就好,没有jdk无法运行java程序,同理没有deepin-wine环境无法运行腾讯产品) 1.2 :去哪 ...

  6. c++制表符

    1 #include<iostream> 2 using namespace std; 3 int main() 4 { 5 //制表符是8个字符,当在一个制表符开始的位置(明白什么位置是 ...

  7. 使用spring向service里面注入dao不成功。

    因为原来的程序没有使用spring.后来加spring的时候action有个地方的new没有改!!! new了个新的实现层 不是spring管理的对象.

  8. Ubuntu16.04 Linux 下无痛安装、配置Gogs

    本文在Win7+VMware的ubuntu 16.04中测试,安装Gogs,Install from binary. 准备工作: sudo apt-get install git sudo addus ...

  9. 【SpringMVC】RESTFul简介以及案例实现

    RESTful 概念 REST:Representational State Transfer,表现层资源状态转移. 资源 资源是一种看待服务器的方式,即,将服务器看作是由很多离散的资源组成.每个资源 ...

  10. centos7 netstat

    netstat 是控制台命令,它可以显示路由表.实际的网络连接以及每一个网络接口设备的状态信息.Netstat 用于显示与 IP . TCP . UDP 和 ICMP 协议相关的统计数据,一般用于检验 ...