UOJ#395. 【NOI2018】你的名字 字符串,SAM,线段树合并
原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ395.html
题解
记得同步赛的时候这题我爆0了,最暴力的暴力都没调出来。
首先我们看看 68 分怎么做
——求两个串的本质不同的公共子串个数。
它是一个模板题,然而我当时并不会,甚至连SAM都忘了怎么写QAQ。
再简化一下:如何求一个串的本质不同的子串个数。
给串建一个SAM,把所有节点代表的字符串个数(也就是 Max(x) - Max(fa(x)) 加起来就好了。
回到上一个问题。
假设这两个串分别是 S,T 。对 T 建个SAM。
对于T的SAM,考虑对于它的任何一个节点 x ,算出 x 的 Right 集合代表的所有前缀与 S 的所有前缀的 LCS 的最大值(也就是这个节点代表的状态能在 S 上匹配的最长长度),设为 val(x)。然后对于所有 x 把 $(1,val(x)] \cap (Max(fa(x)),Max(x)]$ 的长度加起来就好了。
那么如何求那个最长的匹配长度?对 S 建一个 SAM,然后用 T 在 S 的 SAM 上走一遍,找到 T 的每一个前缀的 最长的是 S 的子串的后缀 然后 T 的 SAM 上的一个节点的 val 就是他在 parent 树上的所有后代节点的 Max 。
由于 S 的 SAM 可以预先建好,所以询问一个 T 串的复杂度是 $O(|T|)$ 的。
那么 S 有 [L,R] 的限制呢?
线段树合并预处理一下 S 的 SAM 的每一个节点的 Right 集合。
修改一下求最长的匹配长度的过程,保证走转移边的时候在 [L,R] 中有匹配。
注意这里有一个易错点:我们匹配失败跳 father 的时候,不能直接 len' = Max(father) ,只能不断减一。原因是在 len 不断减一的过程中可能会找到匹配,而直接跳 father 会漏过这个匹配。然而出题人数据出的很水,没注意到这个东西还是有96分!
至此,我们得到了一个 $O((|S|+\sum |T|)\log |S|)$ 的做法。
但是,由于在 SAM 上遍历节点暴力跳祖先的复杂度是 $O(n\sqrt n)$ 的,然后加个线段树合并多个 $\log$ ,总复杂度 $O(n\sqrt n \log n)$ 的可以通过原题数据……wft??(UOJ Hack数据过不去的)
代码
#include <bits/stdc++.h>
#define clr(x) memset(x,0,sizeof (x))
using namespace std;
typedef long long LL;
LL read(){
LL x=0,f=0;
char ch=getchar();
while (!isdigit(ch))
f|=ch=='-',ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return f?-x:x;
}
const int N=500005*4;
int n,m,q;
char s[N];
struct Node{
int Next[26],fa,Max,pos;
};
namespace Seg{
const int S=N*35;
int ls[S],rs[S],cnt=0;
void Ins(int &rt,int L,int R,int x){
if (!rt)
rt=++cnt;
if (L==R)
return;
int mid=(L+R)>>1;
if (x<=mid)
Ins(ls[rt],L,mid,x);
else
Ins(rs[rt],mid+1,R,x);
}
int Merge(int a,int b,int L,int R){
if (!a||!b)
return a+b;
int rt=++cnt;
if (L<R){
int mid=(L+R)>>1;
ls[rt]=Merge(ls[a],ls[b],L,mid);
rs[rt]=Merge(rs[a],rs[b],mid+1,R);
}
return rt;
}
int Query(int rt,int L,int R,int xL,int xR){
if (!rt||R<xL||L>xR||xL>xR)
return 0;
if (xL<=L&&R<=xR)
return 1;
int mid=(L+R)>>1;
return Query(ls[rt],L,mid,xL,xR)
|Query(rs[rt],mid+1,R,xL,xR);
}
}
namespace SAM{
Node t[N];
int root,last,size;
int rt[N],id[N];
void Init(){
while (size){
clr(t[size].Next);
t[size].fa=t[size].Max=t[size].pos=rt[size]=0;
size--;
}
root=last=size=1;
}
void extend(int c,int ps){
int p=last,np=++size,q,nq;
t[np].Max=t[p].Max+1,t[np].pos=ps;
Seg::Ins(rt[np],1,n,ps);
for (;p&&!t[p].Next[c];p=t[p].fa)
t[p].Next[c]=np;
if (!p)
t[np].fa=1;
else {
q=t[p].Next[c];
if (t[p].Max+1==t[q].Max)
t[np].fa=q;
else {
nq=++size;
t[nq]=t[q],t[nq].Max=t[p].Max+1,t[nq].pos=ps;
t[np].fa=t[q].fa=nq;
for (;p&&t[p].Next[c]==q;p=t[p].fa)
t[p].Next[c]=nq;
}
}
last=np;
}
void Sort(){
static int tax[N];
for (int i=0;i<=size;i++)
tax[i]=0;
for (int i=1;i<=size;i++)
tax[t[i].Max]++;
for (int i=1;i<=size;i++)
tax[i]+=tax[i-1];
for (int i=1;i<=size;i++)
id[tax[t[i].Max]--]=i;
}
void build(){
Sort();
for (int i=size;i>1;i--){
int x=id[i],f=t[x].fa;
rt[f]=Seg::Merge(rt[f],rt[x],1,n);
}
}
}
namespace sam{
Node t[N];
int root,last,size;
int id[N],val[N];
void Init(){
while (size){
clr(t[size].Next);
t[size].fa=t[size].Max=t[size].pos=val[size]=0;
size--;
}
root=last=size=1;
}
void extend(int c,int ps){
int p=last,np=++size,q,nq;
t[np].Max=t[p].Max+1,t[np].pos=ps;
for (;p&&!t[p].Next[c];p=t[p].fa)
t[p].Next[c]=np;
if (!p)
t[np].fa=1;
else {
q=t[p].Next[c];
if (t[p].Max+1==t[q].Max)
t[np].fa=q;
else {
nq=++size;
t[nq]=t[q],t[nq].Max=t[p].Max+1,t[nq].pos=ps;
t[np].fa=t[q].fa=nq;
for (;p&&t[p].Next[c]==q;p=t[p].fa)
t[p].Next[c]=nq;
}
}
last=np;
}
void Sort(){
static int tax[N];
for (int i=0;i<=size;i++)
tax[i]=0;
for (int i=1;i<=size;i++)
tax[t[i].Max]++;
for (int i=1;i<=size;i++)
tax[i]+=tax[i-1];
for (int i=1;i<=size;i++)
id[tax[t[i].Max]--]=i;
}
LL solve(){
Sort();
LL ans=0;
for (int i=size;i>1;i--){
int x=id[i],f=t[x].fa;
val[f]=max(val[x],val[f]);
ans+=max(0,t[x].Max-max(t[f].Max,val[x]));
}
return ans;
}
}
int main(){
scanf("%s",s+1);
n=strlen(s+1),q=read();
SAM::Init();
for (int i=1;i<=n;i++)
SAM::extend(s[i]-'a',i);
SAM::build();
SAM::t[0].Max=-1;
while (q--){
scanf("%s",s+1);
m=strlen(s+1);
sam::Init();
int L=read(),R=read();
int x=1,len=0;
for (int i=1;i<=m;i++){
int c=s[i]-'a',nowx=sam::size+1;
sam::extend(c,i);
while (x){
int nx=SAM::t[x].Next[c];
if (nx&&Seg::Query(SAM::rt[nx],1,n,L+len,R))
break;
if ((--len)==SAM::t[SAM::t[x].fa].Max)
x=SAM::t[x].fa;
}
if (!x)
x=1,len=0;
else {
x=SAM::t[x].Next[c];
sam::val[nowx]=++len;
}
}
printf("%lld\n",sam::solve());
}
return 0;
}
UOJ#395. 【NOI2018】你的名字 字符串,SAM,线段树合并的更多相关文章
- bzoj5417/luoguP4770 [NOI2018]你的名字(后缀自动机+线段树合并)
bzoj5417/luoguP4770 [NOI2018]你的名字(后缀自动机+线段树合并) bzoj Luogu 给出一个字符串 $ S $ 及 $ q $ 次询问,每次询问一个字符串 $ T $ ...
- 【NOI2018】你的名字(SAM & 线段树合并)
Description Hint Solution 不妨先讨论一下无区间限制的做法. 首先"子串"可以理解为"前缀的后缀",因此我们定义一个 \(\lim(i) ...
- BZOJ5417[Noi2018]你的名字——后缀自动机+线段树合并
题目链接: [Noi2018]你的名字 题目大意:给出一个字符串$S$及$q$次询问,每次询问一个字符串$T$有多少本质不同的子串不是$S[l,r]$的子串($S[l,r]$表示$S$串的第$l$个字 ...
- 洛谷P4482 [BJWC2018]Border 的四种求法 字符串,SAM,线段树合并,线段树,树链剖分,DSU on Tree
原文链接https://www.cnblogs.com/zhouzhendong/p/LuoguP4482.html 题意 给定一个字符串 S,有 q 次询问,每次给定两个数 L,R ,求 S[L.. ...
- P4770-[NOI2018]你的名字【SAM,线段树合并】
正题 题目链接:https://www.luogu.com.cn/problem/P4770 题目大意 给出一个长度为\(n\)的字符串\(S\).\(q\)次询问给出一个串\(T\)和一个区间\([ ...
- Codeforces 700E. Cool Slogans 字符串,SAM,线段树合并,动态规划
原文链接https://www.cnblogs.com/zhouzhendong/p/CF700E.html 题解 首先建个SAM. 一个结论:对于parent树上任意一个点x,以及它所代表的子树内任 ...
- loj#2059. 「TJOI / HEOI2016」字符串 sam+线段树合并+倍增
题意:给你一个子串,m次询问,每次给你abcd,问你子串sa-b的所有子串和子串sc-d的最长公共前缀是多长 题解:首先要求两个子串的最长公共前缀就是把反过来插入变成最长公共后缀,两个节点在paren ...
- BZOJ.5417.[NOI2018]你的名字(后缀自动机 线段树合并)
LOJ 洛谷 BZOJ 考虑\(l=1,r=|S|\)的情况: 对\(S\)串建SAM,\(T\)在上面匹配,可以得到每个位置\(i\)的后缀的最长匹配长度\(mx[i]\). 因为要去重,对\(T\ ...
- luogu4770 [NOI2018]你的名字 后缀自动机 + 线段树合并
其实很水的一道题吧.... 题意是:每次给定一个串\(T\)以及\(l, r\),询问有多少个字符串\(s\)满足,\(s\)是\(T\)的子串,但不是\(S[l .. r]\)的子串 统计\(T\) ...
随机推荐
- python学习day11 函数Ⅲ (内置函数与lambda表达式)
函数Ⅲ(内置函数&lambda表达式) 1.函数小高级 函数可以当做变量来使用: def func(): print(123) func_list = [func, func, func] # ...
- java替换ascii表字符
如下: //处理特殊字符 public String dealSpecialXml(String xml){ String result = ""; //result = xml. ...
- (二叉树 BFS) leetcode102. Binary Tree Level Order Traversal
Given a binary tree, return the level order traversal of its nodes' values. (ie, from left to right, ...
- A Graph Partitioning Game Theoretical Approach for the VNF Service Chaining Problem
文章名称:A Graph Partitioning Game Theoretical Approach for the VNF Service Chaining Problem 发表时间:2017 期 ...
- 网页三剑客之CSS
1.CSS概述 CSS中文简称层叠样式表(英文全称:Cascading Style Sheets),用来控制页面的表现,即使页面更好看的语言. 2.CSS基本语法和页面引用 2.1 css的定义方法 ...
- Client-Side Template Injection with AngularJS
<html> <head> <meta charset="utf-8"> <script src="https://cdn.bo ...
- django - 总结 - cnblog
1.头像预览 -------方法1-------- 点击头像------>点击input img和input重合; img在label,input-->display:none $(&qu ...
- HTML(五)HTML表格
标准表格 <table border="1"> <caption>Monthly savings</caption> <tr> &l ...
- mysql常见的问题
1.为什么选择某一个版本 各个版本之间的区别及优缺点 首先,服务器特性 mysql percona mysql mariaDB 开源 开源 开源 支持分区表 支持分区表 支持分区表 innodb Xt ...
- 《我是一只IT小小鸟读后感》
在我步入大学前,并未了解何为IT,真是毫无知晓.由于种种原因最终还是选择了软件工程专业,是 对是错,是福是祸,不知该不该去考虑,但即已然 选择了这条路,便得付出一些努力,这个世界总 是有许多在默默努力 ...