[BZOJ2555]SubString LCT+后缀自动机
2555: SubString
Time Limit: 30 Sec Memory Limit: 512 MB
Submit: 3253 Solved: 975
[Submit][Status][Discuss]
Description
懒得写背景了,给你一个字符串init,要求你支持两个操作
(1):在当前字符串的后面插入一个字符串
(2):询问字符串s在当前字符串中出现了几次?(作为连续子串)
你必须在线支持这些操作。
Input
第一行一个数Q表示操作个数
第二行一个字符串表示初始字符串init
接下来Q行,每行2个字符串Type,Str
Type是ADD的话表示在后面插入字符串。
Type是QUERY的话表示询问某字符串在当前字符串中出现了几次。
为了体现在线操作,你需要维护一个变量mask,初始值为0
读入串Str之后,使用这个过程将之解码成真正询问的串TrueStr。
询问的时候,对TrueStr询问后输出一行答案Result
然后mask = mask xor Result
插入的时候,将TrueStr插到当前字符串后面即可。
HINT:ADD和QUERY操作的字符串都需要解压
Output
Sample Input
A
QUERY B
ADD BBABBBBAAB
Sample Output
HINT
40 % 的数据字符串最终长度 <= 20000,询问次数<= 1000,询问总长度<= 10000
100 % 的数据字符串最终长度 <= 600000,询问次数<= 10000,询问总长度<= 3000000
新加数据一组--2015.05.20
Source
显然我们要求当前子串子树的后缀结束结点个数。
要在线实现,用lct维护sam
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define maxn 1200000
using namespace std;
int mask;
char st[];
string chars;
void gets(int mask) {
scanf("%s",st);
chars=st;
for(int j=;j<chars.length();j++) {
mask=(mask*+j)%chars.length();
char t=chars[j];
chars[j]=chars[mask];
chars[mask]=t;
}
}
struct data {
int last,cnt;
int son[maxn][],link[maxn],step[maxn];
data() {last=cnt=;}
struct lct {
int fa[maxn],s[maxn][],v[maxn],tag[maxn];
bool isroot(int x) {return s[fa[x]][]!=x&&s[fa[x]][]!=x;}
void add(int x,int y){if(x) v[x]+=y,tag[x]+=y;}
void pushdown(int x) {
if(tag[x]) {
add(s[x][],tag[x]);add(s[x][],tag[x]);
tag[x]=;
}
}
void push(int x) {
if(!isroot(x)) push(fa[x]);
pushdown(x);
}
void rorate(int x) {
int y=fa[x],z=fa[y];
bool l=(s[y][]==x),r=l^;
if(!isroot(y)) s[z][s[z][]==y]=x;
fa[x]=z;fa[y]=x;s[y][l]=s[x][r];
fa[s[x][r]]=y;s[x][r]=y;
}
void splay(int x) {
push(x);
while(!isroot(x)) {
int y=fa[x],z=fa[y];
if(!isroot(y)) {
if(s[y][]==x^s[z][]==y) rorate(x);
else rorate(y);
}
rorate(x);
}
}
void access(int x) {for(int t=;x;t=x,x=fa[x]) {splay(x);s[x][]=t;}}
void link(int x,int y) {fa[x]=y;access(y);splay(y);add(y,v[x]);}
void cut(int x) {access(x);splay(x);add(s[x][],-v[x]);s[x][]=fa[s[x][]]=;}
}t;
void extend(int c) {
int p=last,np=last=++cnt;step[np]=step[p]+;t.v[np]=;
while(p&&!son[p][c]) son[p][c]=np,p=link[p];
if(!p) link[np]=,t.link(np,);
else {
int q=son[p][c];
if(step[q]==step[p]+) link[np]=q,t.link(np,q);
else {
int nq=++cnt;
memcpy(son[nq],son[q],sizeof(son[q]));
link[nq]=link[q];t.link(nq,link[q]);
link[q]=link[np]=nq;
t.cut(q);t.link(np,nq);t.link(q,nq);
step[nq]=step[p]+;
while(son[p][c]==q&&p) son[p][c]=nq,p=link[p];
}
}
}
void add() {
gets(mask);
int l=chars.length();
for(int i=;i<l;i++) extend(chars[i]-'A');
}
int query() {
gets(mask);
int p=,l=chars.length();
for(int i=;i<l;i++) {
p=son[p][chars[i]-'A'];if(!p) return ;
}
t.splay(p);
return t.v[p];
}
}sam;
int main() {
int Q;scanf("%d",&Q);
scanf("%s",st+);
int len=strlen(st+);
for(int i=;i<=len;i++) sam.extend(st[i]-'A');
while(Q--) {
char s[];
scanf("%s",s);
if(s[]=='A') sam.add();
else {
int ans=sam.query();
printf("%d\n",ans);
mask^=ans;
}
}
}
[BZOJ2555]SubString LCT+后缀自动机的更多相关文章
- bzoj2555 substring(LCT 后缀自动机)
/* 动态求right集合的大小 LCT维护parent树即可 注意 由于树是有向的不会换根并且每次操作单一, 于是不需要维护子树和(写起来很麻烦) 直接打标记修改即可 */ #include< ...
- BZOJ2555 SubString 【后缀自动机 + LCT】
题目 懒得写背景了,给你一个字符串init,要求你支持两个操作 (1):在当前字符串的后面插入一个字符串 (2):询问字符串s在当前字符串中出现了几次?(作为连续子串) 你必须在线支持这些操作. 输入 ...
- BZOJ2555: SubString(后缀自动机,LCT维护Parent树)
Description 懒得写背景了,给你一个字符串init,要求你支持两个操作 (1):在当前字符串的后面插入一个字符串 (2):询问字符串s在当前字符串中出现了几次?(作为连续子串) 你必须在线支 ...
- BZOJ2555 SubString【后缀自动机+LCT】
Description 懒得写背景了,给你一个字符串init,要求你支持两个操作 (1):在当前字符串的后面插入一个字符串 (2):询问字符串s在当前字符串中出现了几次?(作为连续子串) 你必须在线支 ...
- BZOJ2555 SubString(后缀自动机+LCT)
询问串放在SAM上不跳fail跑到的节点的|right|即为答案.用LCT维护parent树即可.可以直接维护子树信息,也可以转化为路径加.注意强制在线所使用的mask是作为参数传进去的. #incl ...
- 【BZOJ2555】SubString(后缀自动机,Link-Cut Tree)
[BZOJ2555]SubString(后缀自动机,Link-Cut Tree) 题面 BZOJ 题解 这题看起来不难 每次要求的就是\(right/endpos\)集合的大小 所以搞一个\(LCT\ ...
- spoj 7258 Lexicographical Substring Search (后缀自动机)
spoj 7258 Lexicographical Substring Search (后缀自动机) 题意:给出一个字符串,长度为90000.询问q次,每次回答一个k,求字典序第k小的子串. 解题思路 ...
- 【SPOJ】Longest Common Substring(后缀自动机)
[SPOJ]Longest Common Substring(后缀自动机) 题面 Vjudge 题意:求两个串的最长公共子串 题解 \(SA\)的做法很简单 不再赘述 对于一个串构建\(SAM\) 另 ...
- bzoj 2555: SubString【后缀自动机+LCT】
一直WA--找了半天错的发现居然是解密那里的mask其实是不能动的--传进去的会变,但是真实的那个不会变-- 然后就是后缀自动机,用LCT维护parent树了--注意不能makeroot,因为自动机的 ...
随机推荐
- [转] UnityVS(Visual Studio Tools For Unity)的安装与使用
一些废话 Unity 的开发者们,尤其是微软系的Unity开发者们,用Mono是不是烦死了?你是不是跟我一样,用vs来写代码,用Mono来跟踪调试?好麻烦啊好麻烦. 也许你会说,傻逼你不会用Unity ...
- HDU 1798 Tell me the area
http://acm.hdu.edu.cn/showproblem.php?pid=1798 Problem Description There are two circles in the ...
- x86/x64的stack*****************************TBD
1.push parameter, %rdi,%rsi,%rdx,%rcx,%r8,%r9 用作函数参数,依次对应第1参数,第2参数... 2. push return address 3. push ...
- BZOJ4753 JSOI2016最佳团体(分数规划+树形dp)
看到比值先二分答案.于是转化成一个非常裸的树形背包.直接暴力背包的话复杂度就是O(n2),因为相当于在lca处枚举每个点对.这里使用一种更通用的dfs序优化树形背包写法.https://www.cnb ...
- [Leetcode] Best time to buy and sell stock ii 买卖股票的最佳时机
Say you have an array for which the i th element is the price of a given stock on day i. Design an a ...
- 洛谷P4593 [TJOI2018]教科书般的亵渎 【数学】
题目链接 洛谷P4593 题解 orz dalao upd:经典的自然数幂和,伯努利数裸题 由题我们只需模拟出代价,只需使用\(S(n,k) = \sum\limits_{i = 1}^{n} i^{ ...
- yum命令Header V3 RSA/SHA1 Signature, key ID c105b9de: NOKEY
yum命令Header V3 RSA/SHA1 Signature, key ID c105b9de: NOKEY 博客分类: linux 三种解决方案 我采取第三种方案解决的 第一种: linu ...
- POJ - 1017 贪心训练
Packets Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 59725 Accepted: 20273 Descrip ...
- 在xml文件中引入带有@Configuration的Java类
在xml引入带有@Configuration的Java类,就是把这个带有@Configuration的Java类,当做一个普通的的类用<bean>标签引入: 核心代码如下: @Config ...
- matlab求矩阵的鞍点
function count = andian(a) v = max(a,[],2); count = 0; for i=1:length(v) [r2,c2] = find(a==v(i)); mi ...