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 ...
随机推荐
- hiho一下 第174周
题目1 : Dice Possibility 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 What is possibility of rolling N dice ...
- linux下使用Android studio启动模拟器时提示 waiting for target device to come online 的问题
方法来自:http://stackoverflow.com/questions/42612468/how-can-i-get-more-information-about-waiting-for-ta ...
- 开源作品-ThinkPHP在线分析工具(单文件绿色版)-TPLogAnalysis_PHP_1_0
TPLogAnalysis_PHP_1_0 前言:项目开发基于ThinkPHP框架,但是在调试程序的时候,没有一款日志可视化分析工具.在网络也找不到任何相关的TP日志分析工具.求人不如求己,于是决定抽 ...
- HDL之Bitwise operation
1 Verilog 1.1 Bitwise operator Bitwise operators perform a bit wise operation on two operands. They ...
- ubuntu 安装 OpenCV-CUDA
参考链接:http://www.cnblogs.com/platero/p/3993877.html 官方指导:https://help.ubuntu.com/community/OpenCV 0.这 ...
- Socket 数据包顺序的问题
今天遇到一个问题,到现在还未查明原因,记录一下,留后续跟踪. 基于Netty的Socket通讯问题,Server在向Client发送数据时,假设数据原顺序为1,2,3,4... 但到了客户端顺序可能 ...
- 高手的C++学习忠告,虚心学习下~~[转载]
1.把C++当成一门新的语言学习(和C没啥关系!真的.): 2.看<Thinking In C++>,不要看<C++变成死相>: 3.看<The C++ Programm ...
- Java导入Excel文件页面实现JS
Excel导入: 页面创建导入按钮,如: 代码: <button class="layui-btn layui-btn-small layui-btn-primary ajax-all ...
- 关于配置websocket,nginx转发https至wss问题
在本地测试通过的socket,再放到现在的有nginx代理之后发现会报:failed: Error in connection establishment: net::ERR_NAME_NOT_RES ...
- RobotFrameWork+APPIUM实现对安卓APK的自动化测试----第一篇【安装】
文章来源http://blog.csdn.net/deadgrape/article/details/50563119 前言:关于RobotFrameWork+APPIUM实现对安卓APK的自动化测试 ...