BZOJ 3230 后缀数组+ST
思路:
首先我们已经会了后缀数组求本质不同的子串个数
这道题跟那个差不多
首先我们可以知道按字典序排好的每个后缀之前包含多少本质不同的字串
就是sigma(n-sa[i]+1-ht[i]+bi[i-1])
在这上面二分就可以求出来原串是什么了
就把两个后缀的LCP和原串的长度取个min即可
再把串倒过来也差不多这么搞一遍就好啦
//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define int long long
const int N=100050;
int n,q,A[N],B[N],cntA[N],cntB[N],rk[N],ht[N],sa[N],tsa[N],f[N][20],bi[N],base[N];
char s[N];
struct Ask{int x,y,bx,by,backx,backy,lim,ans;}ask[N];
void SA(){
memset(cntA,0,sizeof(cntA));
for(int i=1;i<=n;i++)cntA[s[i]]++;
for(int i='b';i<='z';i++)cntA[i]+=cntA[i-1];
for(int i=n;i;i--)sa[cntA[s[i]]--]=i;
rk[sa[1]]=1;
for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
for(int l=1;rk[sa[n]]<n;l<<=1){
memset(cntA,0,sizeof(cntA));
memset(cntB,0,sizeof(cntB));
for(int i=1;i<=n;i++)
cntA[A[i]=rk[i]]++,
cntB[B[i]=(i+l<=n?rk[i+l]:0)]++;
for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1];
for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
rk[sa[1]]=1;
for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]);
}
for(int i=1,j=0;i<=n;i++){
j=j?j-1:0;
while(s[i+j]==s[sa[rk[i]-1]+j])j++;
ht[rk[i]]=j;
}
}
void init_rmq(){
for(int i=1;i<=n;i++)f[i][0]=ht[i];
for(int j=1;j<=18;j++)
for(int i=1;i+(1<<(j-1))<=n;i++)
f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
int LCP(int x,int y){
int t=base[y-x+1];
return min(f[x][t],f[y-(1<<t)+1][t]);
}
signed main(){
scanf("%lld%lld%s",&n,&q,s+1);
base[0]=-1;
for(int i=1;i<=n;i++)base[i]=base[i>>1]+1;
SA(),init_rmq();
for(int i=1;i<=n;i++)bi[i]=n-sa[i]+1-ht[i]+bi[i-1];
for(int i=1;i<=q;i++){
scanf("%lld%lld",&ask[i].x,&ask[i].y);
if(ask[i].x>ask[i].y)swap(ask[i].x,ask[i].y);
ask[i].bx=lower_bound(bi+1,bi+1+n,ask[i].x)-bi,
ask[i].by=lower_bound(bi+1,bi+1+n,ask[i].y)-bi;
int lena=n-(bi[ask[i].bx]-ask[i].x)-sa[ask[i].bx]+1,
lenb=n-(bi[ask[i].by]-ask[i].y)-sa[ask[i].by]+1;
ask[i].backx=n-(lena+sa[ask[i].bx]-1)+1;
ask[i].backy=n-(lenb+sa[ask[i].by]-1)+1;
if(ask[i].by>n)ask[i].lim=-1;
else ask[i].lim=min(lena,lenb);
if(~ask[i].lim){
int t=min(ask[i].lim,ask[i].bx!=ask[i].by?LCP(ask[i].bx+1,ask[i].by):N);
ask[i].ans+=t*t;
}
}
for(int i=1;i<=n/2;i++)swap(s[i],s[n-i+1]);
SA(),init_rmq();
for(int i=1;i<=n;i++)bi[i]=n-sa[i]+1-ht[i]+bi[i-1];
for(int i=1;i<=q;i++){
ask[i].bx=lower_bound(bi+1,bi+1+n,ask[i].x)-bi,
ask[i].by=lower_bound(bi+1,bi+1+n,ask[i].y)-bi;
if(~ask[i].lim){
int tempx=rk[ask[i].backx],tempy=rk[ask[i].backy];
if(tempx>tempy)swap(tempx,tempy);
int t=min(ask[i].lim,tempx<tempy?LCP(tempx+1,tempy):N);
printf("%lld\n",ask[i].ans+1LL*t*t);
}
else puts("-1");
}
}
BZOJ 3230 后缀数组+ST的更多相关文章
- BZOJ 4516 后缀数组+ST+set
写了一半 没了啊啊啊 重新写的 思路: 先不考虑后缀自动机 (我不会啊) 那这道题只能用后缀数组了 先把原串倒一下 后缀->前缀 相当于每回在前面加了一个字母 求不同的子串个数 首先 正常的求子 ...
- SPOJ 687 Repeats(后缀数组+ST表)
[题目链接] http://www.spoj.com/problems/REPEATS/en/ [题目大意] 求重复次数最多的连续重复子串的长度. [题解] 考虑错位匹配,设重复部分长度为l,记s[i ...
- POJ 3693 Maximum repetition substring(后缀数组+ST表)
[题目链接] poj.org/problem?id=3693 [题目大意] 求一个串重复次数最多的连续重复子串并输出,要求字典序最小. [题解] 考虑错位匹配,设重复部分长度为l,记s[i]和s[i+ ...
- BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay
BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay Description 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1.2 拼凑起来形成一个魔 ...
- UVA10829 L-Gap Substrings(后缀数组+ST表)
后缀数组+ST表. 代填的坑. \(Code\ Below:\) #include <bits/stdc++.h> #define ll long long using namespace ...
- bzoj千题计划312:bzoj2119: 股市的预测(后缀数组+st表)
https://www.lydsy.com/JudgeOnline/problem.php?id=2119 题意:将给定数组差分后,求ABA形式的字串个数,要求|B|=m,|A|>0 1.后缀数 ...
- BZOJ 4556 [Tjoi2016&Heoi2016]字符串 ——后缀数组 ST表 主席树 二分答案
Solution 1: 后缀数组暴力大法好 #include <map> #include <cmath> #include <queue> #include &l ...
- POJ3693 Maximum repetition substring [后缀数组 ST表]
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9458 Acc ...
- 【BZOJ-4310】跳蚤 后缀数组 + ST表 + 二分
4310: 跳蚤 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 180 Solved: 83[Submit][Status][Discuss] De ...
随机推荐
- Codeforces Round #447
QAQ #include<stdio.h> #include<string.h> #include<stdlib.h> #include<vector> ...
- JS网站图集相册特效
JS网站图集相册特效是一款可以直接使用鼠标进行前后导航,也可以通过缩略图来切换图片. 在线演示本地下载
- 使用纯 CSS 实现 Google Photos 照片列表布局
文章太长,因为介绍了如何一步一步进化到最后接近完美的效果的,不想读的同学可以直接跳到最后一个大标题之后看代码.demo及原理就好,或者也可以直接看下面这个链接的源代码. 不过还是建议顺序读下去,因为后 ...
- js-字符串方法
字符串 遍历字符串 方法:(类似数组) 使用for 或 for… in 结果:得到字符串中的每个字符 查找字符 ² charAt(索引值) 注: 超出索引值范围时,则返回空字符 ² ch ...
- Spring DATA MongoDB @DBref查询,or和and联合查询
@DBref文档关联,在按该类型查询的时候,在字段名后加上关联表的字段名即可,如: Criteria.where("bloggroup.$id"), $id代表关联表的oid字段. ...
- Caffe: Caffe的Python接口
官方参考:http://caffe.berkeleyvision.org/installation.html 官方介绍是这样的: Python The main requirements are nu ...
- vue项目布局
1.底部有分类布局 类似这种底部有分类的,点击四个tap分别道不同的页面这样的,每个页面都是一个路由,把底部作为一个组件在每一个页面中引入就行.组件就是公用的,能公用的就写成组件.如下 { path: ...
- 【转】ORACLE SQL基础—DDL语言 礼记八目 2017-12-23 21:26:21
原文地址:https://www.toutiao.com/i6502733303550837261/ SQL语言分为:DDL数据定义语言,DML数据操纵语言,DCL是数据库控制语言,TC事务控制语言 ...
- Smallest Common Multiple FreeCodeCamp
题目:找出能被两个给定参数和它们之间的连续数字整除的最小公倍数. 范围是两个数字构成的数组,两个数字不一定按数字顺序排序. 分析:首先题目的意思求一个连续数列的所有数字的最小公倍数,这连续的数字序列 ...
- nyoj11-奇偶数分离
奇偶数分离 时间限制:3000 ms | 内存限制:65535 KB 难度:1 描述 有一个整型偶数n(2<= n <=10000),你要做的是:先把1到n中的所有奇数从小到大输出,再 ...