BZOJ1014 JSOI2008 火星人prefix


Description

  火星人最近研究了一种操作:求一个字串两个后缀的公共前缀。比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字符予以标号:序号: 1 2 3 4 5 6 7 8 9 10 11 字符 m a d a m i m a d a m 现在,火星人定义了一个函数LCQ(x, y),表示:该字符串中第x个字符开始的字串,与该字符串中第y个字符开始的字串,两个字串的公共前缀的长度。比方说,LCQ(1, 7) = 5, LCQ(2, 10) = 1, LCQ(4, 7) = 0 在研究LCQ函数的过程中,火星人发现了这样的一个关联:如果把该字符串的所有后缀排好序,就可以很快地求出LCQ函数的值;同样,如果求出了LCQ函数的值,也可以很快地将该字符串的后缀排好序。 尽管火星人聪明地找到了求取LCQ函数的快速算法,但不甘心认输的地球人又给火星人出了个难题:在求取LCQ函数的同时,还可以改变字符串本身。具体地说,可以更改字符串中某一个字符的值,也可以在字符串中的某一个位置插入一个字符。地球人想考验一下,在如此复杂的问题中,火星人是否还能够做到很快地求取LCQ函数的值。

Input

  第一行给出初始的字符串。第二行是一个非负整数M,表示操作的个数。接下来的M行,每行描述一个操作。操作有3种,如下所示

1、询问。语法:Qxy,x,y均为正整数。功能:计算LCQ(x,y)限制:1<=x,y<=当前字符串长度。

2、修改。语法:Rxd,x是正整数,d是字符。功能:将字符串中第x个数修改为字符d。限制:x不超过当前字符串长度。

3、插入:语法:Ixd,x是非负整数,d是字符。功能:在字符串第x个字符之后插入字符d,如果x=0,则在字符串开头插入。限制:x不超过当前字符串长度

Output

  对于输入文件中每一个询问操作,你都应该输出对应的答案。一个答案一行。

Sample Input

madamimadam

7

Q 1 7

Q 4 8

Q 10 11

R 3 a

Q 1 7

I 10 a

Q 2 11

Sample Output

5

1

0

2

1

HINT

1、所有字符串自始至终都只有小写字母构成。

2、M<=150,000

3、字符串长度L自始至终都满足L<=100,000

4、询问操作的个数不超过10,000个。

对于第1,2个数据,字符串长度自始至终都不超过1,000

对于第3,4,5个数据,没有插入操作。


我们发现有修改和添加字符的操作,就想到了平衡树

又因为我比较喜欢非旋转Treap

然后就写了一波

注意这里的Treap不能按照字符的大小来建树,我们应该根据字符串中字符的先后顺序来建树,添加字符就直接在第k个位置后面添加就好了,修改比较常规,直接先分裂再合并就行了

现在我们考虑怎么求字符串的两个后缀的最长公共前缀,首先如果没有修改和添加操作我们肯定会选择后缀数组,但现在我们有了修改和添加操作,所以我们需要换一种方式思考。在一般情况下,字符串的hash值是否相同可以判断字符串是否相同,所以我们可以用hash,但是如果这样暴力判断,时间效率依旧不优秀,达到O(n2logn)" role="presentation">O(n2logn)O(n2logn),那么我们可以发现这个hash段的长度是具有二分性的,所以我们就可以用二分来查询,就可以将时间优化成O(nlog2n)" role="presentation">O(nlog2n)O(nlog2n)

那么我们怎么维护hash呢,首先,在Treap上进行向上更新操作的时候可以通过合并左右区间hash值来维护,那么我们就需要与处理base偏移量(pow数组),然后通过

hash[t]=hash[ls[t]]+pow[siz[ls[t]]]∗val[t]+pow[siz[ls[t]]+1]∗hash[rs[t]]" role="presentation">hash[t]=hash[ls[t]]+pow[siz[ls[t]]]∗val[t]+pow[siz[ls[t]]+1]∗hash[rs[t]]hash[t]=hash[ls[t]]+pow[siz[ls[t]]]∗val[t]+pow[siz[ls[t]]+1]∗hash[rs[t]]

来计算当前节点为根的子树hash值,其中val是当前节点所代表的字符的编号(字符-‘a’+1)

然后我们在二分的时候直接剥离出从x、y开始长度为len的子树就好了


#include<bits/stdc++.h>
using namespace std;
#define N 1000010
#define Base 37
typedef pair<int,int> pi;
struct Treap{
int tot,root,n;
int siz[N],key[N],ls[N],rs[N],val[N];
unsigned int hash[N],pow[N];
void init(){
tot=root=0;
pow[0]=1;
for(int i=1;i<=N;i++)
pow[i]=pow[i-1]*Base;
}
int new_treap_point(int vl){
int newt=++tot;
siz[newt]=1;
hash[newt]=val[newt]=vl;
key[newt]=rand();
ls[newt]=rs[newt]=0;
return newt;
}
void update(int t){
siz[t]=siz[ls[t]]+siz[rs[t]]+1;
hash[t]=hash[ls[t]]+pow[siz[ls[t]]]*val[t]+pow[siz[ls[t]]+1]*hash[rs[t]];
}
int merge(int a,int b){
if(!a)return b;
if(!b)return a;
if(key[a]<key[b]){
rs[a]=merge(rs[a],b);
update(a);return a;
}else{
ls[b]=merge(a,ls[b]);
update(b);return b;
}
}
pi split(int t,int k){
if(!t)return pi(0,0);
pi ans;
if(siz[ls[t]]>=k){
pi tmp=split(ls[t],k);
ls[t]=tmp.second;
update(t);
ans.first=tmp.first;
ans.second=t;
}else{
pi tmp=split(rs[t],k-siz[ls[t]]-1);
rs[t]=tmp.first;
update(t);
ans.first=t;
ans.second=tmp.second;
}
return ans;
}
void insert(int k,int vl){
pi x=split(root,k);
int t=new_treap_point(vl);
root=merge(merge(x.first,t),x.second);
}
void erase(int k){
pi x=split(root,k);
root=merge(split(x.first,k-1).first,x.second);
}
void modify(int k,int vl){
erase(k);
insert(k-1,vl);
}
int queryhash(int t,int len){
if(t+len-1>siz[root])return -1;
pi x=split(root,t-1);
pi y=split(x.second,len);
int ans=hash[y.first];
root=merge(x.first,merge(y.first,y.second));
return ans;
}
int query(int x,int y){
int l=0,r=N,ans=0;
while(l<=r){
int mid=(l+r)>>1;
int px=queryhash(x,mid),py=queryhash(y,mid);
if(px==-1||py==-1){r=mid-1;continue;}
if(px==py)ans=mid,l=mid+1;
else r=mid-1;
}
return ans;
}
}treap;
int n,m;
char s[N],c[5];
int main(){
srand(time(0));
treap.init();
scanf("%s",s+1);
n=strlen(s+1);
treap.n=n;
for(int i=1;i<=n;i++)
treap.root=treap.merge(treap.root,treap.new_treap_point(s[i]-'a'+1));
scanf("%d",&m);
while(m--){
scanf("%s",c);
if(c[0]=='Q'){
int x,y;scanf("%d%d",&x,&y);
printf("%d\n",treap.query(x,y));
}else if(c[0]=='I'){
int x;scanf("%d%s",&x,c);
treap.insert(x,c[0]-'a'+1);
}else{
int x;scanf("%d%s",&x,c);
treap.modify(x,c[0]-'a'+1);
}
}
return 0;
}

BZOJ1014 JSOI2008 火星人prefix 【非旋转Treap】*的更多相关文章

  1. [bzoj1014][JSOI2008]火星人prefix_非旋转Treap_hash_二分

    火星人prefix bzoj-1014 JSOI-2004 题目大意:给定一个字符串,支持三种操作:1.查询:两个后缀之间的$LCP$:2.单点修改:3.插入一个字符. 注释:$1\le n\le 1 ...

  2. [Bzoj1014][JSOI2008]火星人prefix(无旋Treap&hash)

    题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1014 因为涉及到增加和修改,所以后缀数组就被pass掉了,想到的就是平衡树维护hash值 ...

  3. [BZOJ1014][JSOI2008]火星人prefix

    [BZOJ1014][JSOI2008]火星人prefix 试题描述 火星人最近研究了一种操作:求一个字串两个后缀的公共前缀.比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字 ...

  4. 2018.06.28 BZOJ1014 [JSOI2008]火星人prefix(非旋treap+hash)

    [JSOI2008]火星人prefix Time Limit: 10 Sec Memory Limit: 162 MB Submit: 8951 Solved: 2860 Description 火星 ...

  5. bzoj千题计划106:bzoj1014 [JSOI2008]火星人prefix

    http://www.lydsy.com/JudgeOnline/problem.php?id=1014 两个后缀的最长公共前缀:二分+hash 带修改带插入:splay维护 #include< ...

  6. Jewel Magic UVA - 11996 || bzoj1014: [JSOI2008]火星人prefix

    Jewel Magic UVA - 11996 这是一道用splay/非旋treap做的题(这里用的是非旋treap) 1/2/3是splay/非旋treap的常规操作.对于操作4,可以用哈希法求LC ...

  7. [BZOJ1014] [JSOI2008] 火星人prefix (splay & 二分答案)

    Description 火星人最近研究了一种操作:求一个字串两个后缀的公共前缀.比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字符予以标号:序号: 1 2 3 4 5 6 7 ...

  8. bzoj1014: [JSOI2008]火星人prefix splay+hash+二分

    Description 火星人最近研究了一种操作:求一个字串两个后缀的公共前缀.比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字符予以标号:序号: 1 2 3 4 5 6 7 ...

  9. bzoj1014: [JSOI2008]火星人prefix(splay+hash+二分)

    题目大意:一个字符串三个操作:①求两个后缀的LCP②插入一个字符③修改一个字符. 前几天刚学了hash+二分求lcp,就看到这题. 原来splay还能这么用?!原来splay模板这么好写?我以前写的s ...

随机推荐

  1. 教你如何使用理解懒Redis是更好的Redis

    前言 大家都知道 Redis 是单线程的.真正的内行会告诉你,实际上 Redis 并不是完全单线程,因为在执行磁盘上的特定慢操作时会有多线程.目前为止多线程操作绝大部分集中在 I/O 上以至于在不同线 ...

  2. [java]Stream API——collect、reduce、orElse(x)

    一.collect 1.R collect(Supplier supplier, BiConsumer accumulator, BiConsumer combiner) supplier:一个能创造 ...

  3. javascript模块化编写

    目录: 1. 开篇语 2. 对象形式写法 3. 立即执行函数写法 4. prototype写法 5. 总结 1. 开篇语 现在我们写代码不再是一个人包办所有的活儿,都是在多人合作的情况下完成的.我们只 ...

  4. 树 & 二叉树

    2018-01-04 19:13:46 一.树 在计算机科学中,树(英语:tree)是一种数据结构,用来模拟具有树状结构性质的数据集合.它是由n(n>0)个有限节点组成一个具有层次关系的集合.把 ...

  5. 实例化后的map的默认值

    public class map默认值 { public static void main(String[] args) { Map<String, String> resultMap = ...

  6. Project facet Java 1.8 is not supported by target runtime Apache Tomcat v7.0.

    找到项目下的“.settings”

  7. SPOJ 5152 Brute-force Algorithm EXTREME && HDU 3221 Brute-force Algorithm 快速幂,快速求斐波那契数列,欧拉函数,同余 难度:1

    5152. Brute-force Algorithm EXTREME Problem code: BFALG Please click here to download a PDF version ...

  8. 九、dbms_ddl(提供了在PL/SQL块中执行DDL语句的方法)

    1.概述 作用:提供了在PL/SQL块中执行DDL语句的方法,并且也提供了一些DDL的特殊管理方法. 2.包的组成 1).alter_compile说明:用于重新编译过程.函数和包语法:dbms_dd ...

  9. python获取文件路径, 文件名, 后缀名

    def get_filePath_fileName_fileExt(fileUrl): """ 获取文件路径, 文件名, 后缀名 :param fileUrl: :ret ...

  10. LeetCode OJ:Add and Search Word - Data structure design(增加以及搜索单词)

    Design a data structure that supports the following two operations: void addWord(word) bool search(w ...