传送门

(bzoj上不去我也很无奈啊)

题意:维护一个字符串,支持后面加字符串,给定串询问出现次数。强制在线。

数据范围:暴力跑不过

前置知识:重量平衡树

众所周知,平衡树都用了一些策略保证平衡。

平衡有三种:严格平衡,均摊平衡,期望平衡

其中严格平衡和期望平衡的(当然不是很严格,规模上是平衡的就可以了)称为重量平衡树。(该定义存疑)

这种平衡树每次操作对树的结构影响很小(最大的子树为log级别)。这意味着你可以用它维护子树规模的集合。

也就是说可以在里面套一个动态开点线段树维护子树的所有权值。

常见的有替罪羊树和treap。

(fhq好像也可以,但好像没人用……)

以维护子树集合为例,替罪羊树每次重构的时候暴力更新集合。

treap就比较玄学,旋转完后暴力更新,期望是log的……


回到这道题,我们需要一个叫后缀平衡树的黑科技。

其实就是用重量平衡树维护后缀数组。

首先把所有串倒序,注意是假装一来就是倒的,显然不影响答案,只是改成了在前面加。

把字符一个一个添加,这样每添加一个只会多一个后缀,只需要把这个后缀插入平衡树。

注意:把添加的字符按顺序记下来,这样得到的串是输入顺序,但实际上是倒着存的。

现在的问题是如何快速比较两个后缀。

由于新来的后缀不在树上,但之前的树上的后缀是可以快速比较的。

所以我们比较第一个字符,如果一样就丢掉然后在树上比较。

显然可以O(logn)O(logn)O(logn)做。一个技巧,对于每个节点维护(l,r)(l,r)(l,r),其权值为val=l+r2val=\frac{l+r}2val=2l+r​,左儿子(l,val)(l,val)(l,val),右儿子(val,r)(val,r)(val,r),可以直接用权值O(1)O(1)O(1)比较,因为深度是logloglog的,精度可以保证。 重构时暴力更新。

询问的时候查比SSS小的后缀数,再在后面加一个很大的字符再查一次,两个减一下。由于深度是logloglog的,所以直接暴力比较复杂度是O(∣S∣logn)O(|S|logn)O(∣S∣logn)

莫名觉得比后缀数组好写是怎么回事

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
//#include <windows.h>
#include <cassert>
#define MAXN 3600005
using namespace std;
const double alpha=0.7;
char c[MAXN];
int ch[MAXN][2],siz[MAXN],rt,tot;
double L[MAXN],R[MAXN];
inline double val(const int& x){return x?L[x]+R[x]:-1e6;}
inline bool cmp(const int& x,const int& y){return c[x]<c[y]||(c[x]==c[y]&&val(x-1)<val(y-1));}
int p[MAXN];
void dfs(int u)
{if (ch[u][0]) dfs(ch[u][0]);p[++p[0]]=u;if (ch[u][1]) dfs(ch[u][1]);
}
void build(int& x,int l,int r,double vl,double vr)
{if (l>r) return (void)(x=0);int mid=(l+r)>>1;x=p[mid];siz[x]=r-l+1;L[x]=vl;R[x]=vr;double vmid=(vl+vr)/2;build(ch[x][0],l,mid-1,vl,vmid);build(ch[x][1],mid+1,r,vmid,vr);
}
void insert(int& x,int i,double l,double r)
{if (!x) return (void)(siz[x=i]=1,L[x]=l,R[x]=r);if (siz[ch[x][0]]>siz[x]*alpha||siz[ch[x][1]]>siz[x]*alpha){p[0]=0;dfs(x);build(x,1,p[0],l,r);}++siz[x];L[x]=l;R[x]=r;double mid=(l+r)/2;if (cmp(i,x)) insert(ch[x][0],i,l,mid);else insert(ch[x][1],i,mid,r);
}
inline void insert(char i){c[++tot]=i;insert(rt,tot,0,1e6);}
int de;
int query(int x,char *s,int len)
{if (!x) return 0;int i;for (i=0;s[i]==c[x-i]&&i<len;i++);if (i==len||s[i]<c[x-i]) return query(ch[x][0],s,len);else{de=x;int res=siz[ch[x][0]]+1+query(ch[x][1],s,len);return res;}
}
char s[MAXN];
int len,mask;
inline void read(int type)
{#ifndef ONLINE_JUDGE type=0;#endifscanf("%s",s);len=strlen(s);if (type){int t=mask;for (int i=0;i<len;i++) swap(s[i],s[t=(t*131+i)%len]); }
}
int main()
{int q;scanf("%d",&q);read(0);for (int i=0;i<len;i++)   insert(s[i]);while (q--){read(0);if (s[0]=='A'){read(1);for (int i=0;i<len;i++) insert(s[i]);}else{read(1);reverse(s,s+len);int ans=-query(rt,s,len);int t=de;s[len]='Z'+1;ans+=query(rt,s,len+1);printf("%d\n",ans);mask^=ans;}}return 0;
}

【bzoj2555】Substring【后缀平衡树入门】相关推荐

  1. bzoj4768: wxh loves substring //后缀平衡树

    bzoj4768: 2555加强版之wxh loves substring 题意 给出一个字符串,要求资瓷: 在末尾添加/删除字符: 询问一个串的出现次数. 原串长与变化长度之和<=800000 ...

  2. 牛客-139 I. Substring(后缀数组 or 后缀自动机)

    牛客-139 I. Substring(后缀数组 or 后缀自动机) 题目链接 题意 一个由{a,b,c}\{a, b, c\}{a,b,c}组成的字符串SSS,求S子串的最大的集合,使得集合里的字符 ...

  3. 3682: Phorni 后缀平衡树 线段树

    国际惯例的题面: 考虑如果没有强制在线我们能怎么水掉这个题,先构造出字符串,各种方法求一下后缀数组,然后线段树维护区间rank最小的位置即可. 然而他要求强制在线,支持插入后缀,并比较后缀大小(求ra ...

  4. [学习笔记]后缀平衡树

    后缀数组+平衡树=后缀平衡树 支持动态插入字符(只能往前插入),即插入一个后缀,维护所有后缀的排名 插入后缀找到位置?平衡树上二分 法一: 哈希+二分,太慢 法二: 第一个字符不同,已经可以比较,否则 ...

  5. HDU 5769 Substring(后缀数组)

    Substring Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total ...

  6. Bzoj2555 SubString

    Time Limit: 30 Sec  Memory Limit: 512 MB Submit: 2430  Solved: 719 Description 懒得写背景了,给你一个字符串init,要求 ...

  7. HDU多校4 - 6988 Display Substring(后缀自动机+二分)

    题目链接:点击查看 题目大意:给出一个长度为 nnn 的字符串 sss,每个字母都有一个权值,现在要求所有本质不同子串中权值和第 kkk 大的权值 题目分析:如果没有本质不同,那有一个很简单的二分套二 ...

  8. HDU - 5769 Substring(后缀数组)

    题目链接:点击查看 题目大意:给出一个字符串 s 和一个字符 ch,问字符串 s 中有多少个本质不同的子串包含字符 ch 题目分析:其实就是求本质不同的子串加了一点条件,对于每个sa[ i ],其本质 ...

  9. HihoCoder - 1441 后缀自动机一·基本概念(模拟,后缀自动机入门好题)

    题目链接:点击查看 题目大意:给出后缀自动机的定义以及各部分的参数,要求模拟后缀自动机的各个部分 题目分析:直接按照后缀自动机的定义模拟就好了,因为字符串的长度只有50,所以可以用n^3的算法暴力,还 ...

最新文章

  1. 管理和配置客户端访问服务器
  2. 乔布斯亲笔签名的Apple II手册 拍卖了511万元
  3. 2021高考一诊成绩查询,双流棠湖中学高2021届成都市“一诊”成绩已揭晓
  4. tortoisegitpull 并合_tortoiseGIT 本地分支创建合并
  5. python调用java的jar包_python调用java的jar包报错127
  6. 【记】jQuery中的选择器:visible对visibility:hidden的处理
  7. 第一弹!安利10个让你爽到爆的IDEA必备插件!
  8. java集合框架(set)
  9. Java/JSP中调用SQL Server存储过程完整示例收藏
  10. WPF制作简易串口调试助手(上位机部分)
  11. 如何快速深入理解监控知识? | 技术干货
  12. java计算机毕业设计易医就医购药交互平台源码+系统+mysql数据库+lw文档
  13. 格志AK890打印驱动安装
  14. 将鼠标放在A标签上时浏览器左下角会显示跳转后界面地址,如何隐藏
  15. word文档里怎么撤回上一步的操作
  16. Win10 VS2019+QT/OpenCV/灰点相机/函数信号发生器 配置及其使用
  17. 五种线粒体呼吸链复合体研究工具
  18. 什么是结构化数据、半结构化数据、非结构化数据
  19. windows 安装vagrant reload 失败; No Virtualbox Guest Additions installation found.
  20. getUserProfile:fail can only be invoked by user TAP gesture

热门文章

  1. sql年月日24小时制_24小时制的「无码」真人秀,令人叫绝
  2. 如何可以让球的时速达到100公里?
  3. 中国古代数学有多牛?
  4. 当你死后尸体如何处理?两种新玩法了解一下!
  5. 国外的程序猿可以工作到退休而国内的为什么这么短命(思维认知)
  6. SVM支持向量机绘图
  7. 大一计算机课实训总结1000字,大一计算机实训报告.doc
  8. json 查询 java_怎样在java中查询mysql得到如下的json格式的结果?
  9. pythonjam进不去怎么办_教程看了一遍又一遍,自学Python还是连门都进不去?
  10. 前目的地罗伯森是谁_距离目的地只剩10公里,开车师傅却在高速公路上睡着了...