思路:

首先我们已经会了后缀数组求本质不同的子串个数

这道题跟那个差不多

首先我们可以知道按字典序排好的每个后缀之前包含多少本质不同的字串

就是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的更多相关文章

  1. BZOJ 4516 后缀数组+ST+set

    写了一半 没了啊啊啊 重新写的 思路: 先不考虑后缀自动机 (我不会啊) 那这道题只能用后缀数组了 先把原串倒一下 后缀->前缀 相当于每回在前面加了一个字母 求不同的子串个数 首先 正常的求子 ...

  2. SPOJ 687 Repeats(后缀数组+ST表)

    [题目链接] http://www.spoj.com/problems/REPEATS/en/ [题目大意] 求重复次数最多的连续重复子串的长度. [题解] 考虑错位匹配,设重复部分长度为l,记s[i ...

  3. POJ 3693 Maximum repetition substring(后缀数组+ST表)

    [题目链接] poj.org/problem?id=3693 [题目大意] 求一个串重复次数最多的连续重复子串并输出,要求字典序最小. [题解] 考虑错位匹配,设重复部分长度为l,记s[i]和s[i+ ...

  4. BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay

    BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay Description 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1.2 拼凑起来形成一个魔 ...

  5. UVA10829 L-Gap Substrings(后缀数组+ST表)

    后缀数组+ST表. 代填的坑. \(Code\ Below:\) #include <bits/stdc++.h> #define ll long long using namespace ...

  6. bzoj千题计划312:bzoj2119: 股市的预测(后缀数组+st表)

    https://www.lydsy.com/JudgeOnline/problem.php?id=2119 题意:将给定数组差分后,求ABA形式的字串个数,要求|B|=m,|A|>0 1.后缀数 ...

  7. BZOJ 4556 [Tjoi2016&Heoi2016]字符串 ——后缀数组 ST表 主席树 二分答案

    Solution 1: 后缀数组暴力大法好 #include <map> #include <cmath> #include <queue> #include &l ...

  8. POJ3693 Maximum repetition substring [后缀数组 ST表]

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 9458   Acc ...

  9. 【BZOJ-4310】跳蚤 后缀数组 + ST表 + 二分

    4310: 跳蚤 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 180  Solved: 83[Submit][Status][Discuss] De ...

随机推荐

  1. java 抽象工厂模式简单实例

    抽象工厂模式:提供一个创建一系列的相关的或者依赖的对象的接口,无需指定它们的具体实现类,具体的时间分别在子类工厂中产生. 类似于工厂模式:隔离了具体类的生产实现,使得替换具体的工厂实现类很容易.包含有 ...

  2. mysql5.5和5.6版本更新内容

    mysql 5.5,5.6 比5.1改进地方: 1,5.5默认存储引擎为innodb2,5.5增加cpu多核处理能力:innodb_read_io_threads innodb_write_io_th ...

  3. <script runat=server>与<%%>,<%=%>与<%response.write%>

    我想问一下:在语句<script runat="server"> </script>中编写后台代码和在后台.cs文件中编写后台代码有什么不同,执行效率会不会 ...

  4. python版本及ML库

    一:关于Python版本的选择问题 关于Python的选择问题:要看学术界能不能把科学库迁移到Python3. 1:多个版本共用: 最近发现SciPy的最高版本是3.2,只能是退而求其次,不使用最新版 ...

  5. 【前端分享】jQuery.lazyload详解(转)

    jQuery实现图片延迟加载,不知道是否可以节省带宽呢?有人知道吗?这究竟只是一个视觉特效还是真的能延迟加载减少服务器的请求呢? <script type="text/javascri ...

  6. 一个单元格占两行三列的HTML代码为

    主要是这两个属性: colspan 单元格占多少列 rowspan 单元格占多少行 <table width="200" border="1">&l ...

  7. JS 20180416课时训练

    <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...

  8. 路飞学城Python-Day96

    51-数据库表关系之一对多 为什么有多表的概念? 一个Book有单表的时候只能查询Book表中的内容,要想查询更多的出版社的信息,需要再将单表的内容不断更新 导致单表会存储大量的重复的信息,浪费了大量 ...

  9. K3 新单到老单关联字段的添加

    新单到老单字段的添加分为两种: 一种为文本字段信息的关联,新单与老单字段的信息均为文本字段: 另一种为基础资料信息的关联,新单与老单均为基础资料字段信息.       K3 WISE 11.0中存储老 ...

  10. python数据分析------文本挖掘(jieba)

    1.import jieba jieba的cut函数有三个模式:全模式.精准模式.搜索引擎模式 1 精确模式,试图将句子最精确地切开,适合文本分析: 2 全模式,把句子中所有的可以成词的词语都扫描出来 ...