【XSY1551】往事 广义后缀数组 线段树合并
题目大意
给你一颗trie树,令\(s_i\)为点\(i\)到根的路径上的字符组成的字符串。求\(max_{u\neq v}(LCP(s_u,s_v)+LCS(s_u,s_v))\)
\(LCP=\)最长公共前缀,\(LCS=\)最长公共后缀
\(1\leq n\leq 200000\),字符集为\(\{0\ldots 300\}\)
题解
我们先看看这个\(LCP(s_u,s_v)\)怎么求
广义后缀自动机不行。广义后缀树可能可以,但我不会。广义后缀数组可以。然后我就开始手推广义后缀数组
广义后缀数组:和后缀数组类似,求出\(s_i\)的排名以及\(LCP(s_{sa_{}i-1},s_{sa_i})\)
实现也和后缀数组类似,倍增,把两段\(2^{i-1}\)的信息合并成\(2^i\)的信息。另外还要保存\(s_i\)的长度为\(2^j\)的前缀在所有字符串的长度为\(2^j\)的前缀中的排名。
求完\(sa_i\)和\(rk_i\)后,我们用二分+哈希求\(LCP(s_{sa_{i-1}},s_{sa_i})\)。上面保存下来的所有长度为\(2^j\)的前缀的排名可以当哈希值来用(考场上我写了哈希)。然后用st表来维护区间最小值。现在我们可以\(O(1)\)求出\(LCP(s_u,s_v)\)了
考虑以\(x\)为根的子树,若\(u,v\)在\(x\)的不同子树内,则\(LCS(s_u,s_v)=d_x-1\)。这里\(d_x\)为点\(x\)的深度,根的深度为1。
因为\(LCP(s_{sa_{i-1}},s_{sa_i})\geq LCP(s_{sa_x},s_{sa_y})~(x\leq i-1<i\leq y)\),所以我们只用求出以\(rk_i\)为关键字排序后相邻两个点的\(LCP\)。这个可以用线段树维护。然后线段树合并即可。
时间复杂度:\(O(n\log n)\)
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
struct list
{
int v[200010];
int t[200010];
int w[200010];
int h[200010];
int n;
list()
{
n=0;
memset(h,0,sizeof h);
}
void clear()
{
n=0;
memset(h,0,sizeof h);
}
void add(int x,int y,int z)
{
n++;
v[n]=y;
t[n]=h[x];
w[n]=z;
h[x]=n;
}
};
list l,l2;
int f[200010][20];
ll hs[200010][20];
int d[200010];
int e[200010];
int sa[200010];
int rk[200010];
int sx[200010];
int sy[200010];
int b[200010];
ll mod=1000000007;
ll pw[200010];
int ht[200010];
int st[200010][20];
int lo[200010];
int rt[200010];
int n;
int getmin(int x,int y)
{
int t=lo[y-x+1];
return min(st[x][t],st[y-(1<<t)+1][t]);
}
int query(int x,int y)
{
// x=rk[x];
// y=rk[y];
if(x==y)
return 0x3fffffff;
if(x>y)
swap(x,y);
return getmin(x+1,y);
}
namespace seg
{
struct p
{
int s,first,last,sz;
p()
{
s=first=last=sz=0;
}
};
int cnt;
int ls[4000010];
int rs[4000010];
p s[4000010];
void init()
{
memset(ls,0,sizeof ls);
memset(rs,0,sizeof rs);
cnt=0;
}
p mt(p a,p b)
{
if(!a.sz)
return b;
if(!b.sz)
return a;
p c;
c.sz=a.sz+b.sz;
c.first=a.first;
c.last=b.last;
c.s=max(max(a.s,b.s),query(a.last,b.first));
return c;
}
int insert(int p,int x,int l,int r)
{
if(!p)
p=++cnt;
if(l==r)
{
s[p].sz=1;
s[p].first=s[p].last=x;
return p;
}
int mid=(l+r)>>1;
if(x<=mid)
ls[p]=insert(ls[p],x,l,mid);
else
rs[p]=insert(rs[p],x,mid+1,r);
s[p]=mt(s[ls[p]],s[rs[p]]);
return p;
}
int merge(int x,int y)
{
if(!x||!y)
return x+y;
ls[x]=merge(ls[x],ls[y]);
rs[x]=merge(rs[x],rs[y]);
s[x]=mt(s[ls[x]],s[rs[x]]);
return x;
}
};
int ans=0;
void solve(int x)
{
rt[x]=seg::insert(rt[x],rk[x],1,n);
int i;
for(i=l.h[x];i;i=l.t[i])
{
solve(l.v[i]);
rt[x]=seg::merge(rt[x],rt[l.v[i]]);
}
if(seg::s[rt[x]].sz>1)
ans=max(ans,d[x]-1+seg::s[rt[x]].s);
}
int main()
{
memset(rt,0,sizeof rt);
seg::init();
int i,j,x,y;
pw[0]=311;
for(i=1;i<=50;i++)
pw[i]=pw[i-1]*pw[i-1]%mod;
freopen("recollection.in","r",stdin);
freopen("recollection.out","w",stdout);
memset(f,0,sizeof f);
scanf("%d",&n);
d[1]=1;
e[1]=310;
hs[1][0]=310;
for(i=2;i<=n;i++)
{
scanf("%d%d",&x,&y);
l.add(x,i,y);
f[i][0]=x;
d[i]=d[x]+1;
e[i]=y+1;
hs[i][0]=e[i];
for(j=1;j<=18;j++)
{
f[i][j]=f[f[i][j-1]][j-1];
hs[i][j]=(hs[f[i][j-1]][j-1]*pw[j-1]+hs[i][j-1])%mod;
}
}
int sz=310,k,o;
for(i=1;i<=sz;i++)
b[i]=0;
for(i=1;i<=n;i++)
b[sx[i]=e[i]]++;
for(i=2;i<=sz;i++)
b[i]+=b[i-1];
for(i=n;i>=1;i--)
sa[b[sx[i]]--]=i;
for(j=1;(1<<j)<=n;j++)
{
l2.clear();
k=0;
for(i=1;i<=n;i++)
if(d[sa[i]]>(1<<(j-1)))
l2.add(f[sa[i]][j-1],sa[i],0);
for(i=1;i<=n;i++)
for(o=l2.h[sa[i]];o;o=l2.t[o])
sy[++k]=l2.v[o];
for(i=1;i<=n;i++)
if(d[i]<=(1<<(j-1)))
sy[++k]=i;
for(i=1;i<=sz;i++)
b[i]=0;
for(i=1;i<=n;i++)
b[sx[sy[i]]]++;
for(i=2;i<=sz;i++)
b[i]+=b[i-1];
for(i=n;i>=1;i--)
sa[b[sx[sy[i]]]--]=sy[i];
k=0;
swap(sx,sy);
for(i=1;i<=n;i++)
if(i!=1&&sy[sa[i]]==sy[sa[i-1]]&&((d[sa[i]]<=(1<<(j-1))&&d[sa[i-1]]<=(1<<(j-1)))||(d[sa[i]]>(1<<(j-1))&&d[sa[i-1]]>(1<<(j-1))&&sy[f[sa[i]][j-1]]==sy[f[sa[i-1]][j-1]])))
sx[sa[i]]=k;
else
sx[sa[i]]=++k;
if(k>=n)
break;
sz=k;
}
for(i=1;i<=n;i++)
rk[sa[i]]=i;
ht[0]=0;
for(i=2;i<=n;i++)
{
x=sa[i-1];
y=sa[i];
int now=0;
for(j=18;j>=0;j--)
if(d[x]-1>=(1<<j)&&d[y]-1>=(1<<j)&&hs[x][j]==hs[y][j])
{
now+=(1<<j);
x=f[x][j];
y=f[y][j];
}
ht[i]=now;
}
for(i=1;i<=n;i++)
st[i][0]=ht[i];
for(j=1;j<=18;j++)
for(i=1;i+(1<<j)-1<=n;i++)
st[i][j]=min(st[i][j-1],st[i+(1<<(j-1))][j-1]);
lo[1]=0;
for(i=2;i<=n;i++)
lo[i]=lo[i/2]+1;
solve(1);
printf("%d\n",ans);
return 0;
}
【XSY1551】往事 广义后缀数组 线段树合并的更多相关文章
- cf666E. Forensic Examination(广义后缀自动机 线段树合并)
题意 题目链接 Sol 神仙题Orz 后缀自动机 + 线段树合并 首先对所有的\(t_i\)建个广义后缀自动机,这样可以得到所有子串信息. 考虑把询问离线,然后把\(S\)拿到自动机上跑,同时维护一下 ...
- CF 666E Forensic Examination——广义后缀自动机+线段树合并
题目:http://codeforces.com/contest/666/problem/E 对模式串建广义后缀自动机,询问的时候把询问子串对应到广义后缀自动机的节点上,就处理了“区间”询问. 还要处 ...
- 【CF666E】Forensic Examination - 广义后缀自动机+线段树合并
广义SAM专题的最后一题了……呼 题意: 给出一个长度为$n$的串$S$和$m$个串$T_{1\cdots m}$,给出$q$个询问$l,r,pl,pr$,询问$S[pl\cdots pr]$在$T_ ...
- Codeforces 666E Forensic Examination(广义后缀自动机+线段树合并)
将所有串(包括S)放一块建SAM.对于询问,倍增定位出该子串所在节点,然后要查询的就是该子串在区间内的哪个字符串出现最多.可以线段树合并求出该节点在每个字符串中的出现次数. #include<b ...
- Codeforces.666E.Forensic Examination(广义后缀自动机 线段树合并)
题目链接 \(Description\) 给定串\(S\)和\(m\)个串\(T_i\).\(Q\)次询问,每次询问\(l,r,p_l,p_r\),求\(S[p_l\sim p_r]\)在\(T_l\ ...
- [Luogu5161]WD与数列(后缀数组/后缀自动机+线段树合并)
https://blog.csdn.net/WAautomaton/article/details/85057257 解法一:后缀数组 显然将原数组差分后答案就是所有不相交不相邻重复子串个数+n*(n ...
- BZOJ3413: 匹配(后缀自动机 线段树合并)
题意 题目链接 Sol 神仙题Orz 后缀自动机 + 线段树合并... 首先可以转化一下模型(想不到qwq):问题可以转化为统计\(B\)中每个前缀在\(A\)中出现的次数.(画一画就出来了) 然后直 ...
- 模板—字符串—后缀自动机(后缀自动机+线段树合并求right集合)
模板—字符串—后缀自动机(后缀自动机+线段树合并求right集合) Code: #include <bits/stdc++.h> using namespace std; #define ...
- 【BZOJ4556】[TJOI2016&HEOI2016] 字符串(后缀自动机+线段树合并+二分)
点此看题面 大致题意: 给你一个字符串\(s\),每次问你一个子串\(s[a..b]\)的所有子串和\(s[c..d]\)的最长公共前缀. 二分 首先我们可以发现一个简单性质,即要求最长公共前缀,则我 ...
随机推荐
- 广播频道-BroadcastChannel
BroadcastChannel,就字面意思来言,叫做“广播频道”,官方文档说,该API是用于同源不同页面之间完成通信的功能. 1. 概况 它与window.postMessage的区别就是,Broa ...
- H5 25-CSS三大特性之层叠性
25-CSS三大特性之层叠性 我是段落 <!DOCTYPE html> <html lang="en"> <head> <meta cha ...
- javaweb之Cookie学习
Cookie简介 HTTP是无状态协议,服务器不能记录浏览器的访问状态,也就是说服务器不能区分中两次请求是否由一个客户端发出.这样的设计严重阻碍的Web程序的设计.如:在我们进行网购时,买了一条裤子, ...
- php 对象赋值后改变成员变量影响赋值对象
话不多说看代码 打印结果 对obj1的操作 直接影响了obj2 , 对obj2的操作 直接影响了obj1
- MyEclipse 配置 Tomcat
安装好Tomcat,MyEclipse 之后,利用这两个工具可以开发部署Web 应用,步骤相对手动部署要简洁的多,这里有一个特别要注意的地方:系统里安装JDK.Tomcat.MyEclipse 的版本 ...
- spark、standalone集群 (2)集群zookeeper 热备
测试 cmd spark-examples-1.6.0-hadoop2.6.0.jar spark 2.0以后 就没有这个 jar.需要下载 ./bin/spark-submit -- ...
- Oracle NVL空值处理函数
--NVL空值处理函数 --需求:显示价格表中业主类型ID为1的价格记录 如果上限值为null,则显示9999999 ) from dual; select * from t_pricetable ) ...
- mysql [assword expired
mysql 5.6 在使用Navicat在其他机器上进行远程登录数据库时 会出现 password expired ,需要重新设置一下密码. SET PASSWORD FOR 'root'@'%' = ...
- php变量详解
变量是用于存储信息的"容器". 定义一个变量的语法: $变量名 = 值; 使用变量的例子: <?php $x=5; $y=6; $z=$x+$y; echo $z; ?> ...
- WPF实现滚动显示的TextBlock
在我们使用TextBlock进行数据显示时,经常会遇到这样一种情况就是TextBlock的文字内容太多,如果全部显示的话会占据大量的界面,这是我们就会只让其显示一部分,另外的一部分就让其随着时间的推移 ...