题链:

http://www.lydsy.com/JudgeOnline/problem.php?id=4453

题解:

后缀数组,离线询问,栈
看了一堆题解才看懂,太弱啦 ~

如果对于一个区间[l,r]有两个后缀i,j,(令 i<j),那么哪个更优呢?
1).如果rank[i]>rank[j],那么i肯定比j优。
2).如果rank[i]<rank[j]&&lcp(i,j)<r?j+1,那么j肯定比i优。
3).如果rank[i]<rank[j]&&lcp(i,j)>=r?j+1,那么i比j更优。
所以,从后向前遍历一遍字符串,每次把遍历到的当前位置固定为上述区间的左端点l,
那么后面的区间就会被分为若干块,每一块里面的位置作为右端点的答案相同。也就是说答案是呈区间分部的。


所以对于每次遍历到的左端点 i,首先求出它所影响的区间范围[i,r']
然后在把左端点为 i的询问回答了。
但是怎样做到上述的两个操作呢:求出影响范围,回答询问
学习了各大博主的解法:
维护一个栈,里面按顺序存储着一个个的区间(之前不是说答案连呈区间分部分,维护的就是这些区间,存储当前区间的左右端点和答案)
然后就可以在栈里面搞事情了。
1).求出影响范围:
判断栈顶区间是否可以被当前的 i点的影响范围覆盖:
可以全部覆盖的话就移除栈顶,然后继续判断。
无法覆盖的话就把答案为 i 的那个区间放在栈顶。
部分覆盖的话就在栈顶区间里二分出准确位置,然后把栈顶区间的缩小,并把答案为 i 的那个区间放在栈顶。
2).回答询问:
直接二分出询问的右端点在栈里面的哪个答案区间就好了。

建议直接看看代码的实现。

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 105000
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
struct Stack{int l,r,p;}stk[MAXN];
struct Query{
int l,r,id;
bool operator <(const Query &rtm)const{
return l>rtm.l;
}
}q[MAXN];
char S[MAXN];
int sa[MAXN],rak[MAXN],hei[MAXN];
int ANS[MAXN],stm[MAXN][20],log2[MAXN];
void build(int N,int M){
static int cc[MAXN],ta[MAXN],tb[MAXN],*x,*y,h,p;
x=ta; y=tb; h=0;
for(int i=0;i<M;i++) cc[i]=0;
for(int i=0;i<N;i++) cc[x[i]=S[i]]++;
for(int i=1;i<M;i++) cc[i]+=cc[i-1];
for(int i=N-1;i>=0;i--) sa[--cc[x[i]]]=i;
for(int k=1;p=0,k<N;k<<=1){
for(int i=N-k;i<N;i++) y[p++]=i;
for(int i=0;i<N;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
for(int i=0;i<M;i++) cc[i]=0;
for(int i=0;i<N;i++) cc[x[y[i]]]++;
for(int i=1;i<M;i++) cc[i]+=cc[i-1];
for(int i=N-1;i>=0;i--) sa[--cc[x[y[i]]]]=y[i];
swap(x,y); y[N]=-1; x[sa[0]]=0; M=1;
for(int i=1;i<N;i++)
x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?M-1:M++;
if(M>=N) break;
}
for(int i=0;i<N;i++) rak[sa[i]]=i;
for(int i=0,j;i<N;i++){
if(h) h--;
if(rak[i]){
j=sa[rak[i]-1];
while(S[i+h]==S[j+h]) h++;
}
stm[rak[i]][0]=hei[rak[i]]=h;
}
for(int k=1;k<=log2[N];k++)
for(int i=(1<<k)-1;i<N;i++)
stm[i][k]=min(stm[i-(1<<(k-1))][k-1],stm[i][k-1]);
}
int LCP(int l,int r){
static int k;
l=rak[l]; r=rak[r];
if(l>r) swap(l,r); l++;
k=log2[r-l+1];
return min(stm[l+(1<<k)-1][k],stm[r][k]);
}
bool check(int i,int j,int r){
static int lcp;
if(rak[i]>rak[j]) return 1;
lcp=LCP(i,j);
return lcp>=r-j+1;
}
void solve(int N,int M){
int top=1,j=1,l,r,mid,pos; bool L,R,fg;
stk[0].l=N;stk[top]=(Stack){N-1,N-1,N-1};
for(;q[j].l==N-1&&j<=M;j++)ANS[q[j].id]=stk[top].p+1;
for(int i=N-2;fg=0,i>=0&&j<=M;i--){
for(;top;top--){
L=check(i,stk[top].p,stk[top].l);
R=check(i,stk[top].p,stk[top].r);
if(L^R) fg=1;
if(!R) break;
}
if(fg){
pos=l=stk[top].l; r=stk[top].r;
while(l<=r){
mid=(l+r)>>1;
if(check(i,stk[top].p,mid)) pos=mid,l=mid+1;
else r=mid-1;
}
stk[top].l=pos+1;
}
else pos=stk[top].l-1;
stk[++top].r=pos;
stk[top].l=stk[top].p=i;
for(;q[j].l==i&&j<=M;j++){
l=1; r=top;
while(l<=r){
mid=(l+r)>>1;
if(stk[mid].l<=q[j].r&&q[j].r<=stk[mid].r) break;
else if(q[j].r<stk[mid].l) l=mid+1;
else r=mid-1;
}
ANS[q[j].id]=stk[mid].p+1;
}
}
}
int main()
{
log2[1]=0; for(int i=2;i<=100000;i++) log2[i]=log2[i>>1]+1;
scanf("%s",S); int N=strlen(S),M;
build(N,300); scanf("%d",&M);
for(int i=1;i<=M;i++) scanf("%d%d",&q[i].l,&q[i].r),q[i].l--,q[i].r--,q[i].id=i;
sort(q+1,q+M+1);
solve(N,M);
for(int i=1;i<=M;i++) printf("%d\n",ANS[i]);
return 0;
}

●BZOJ 4453 cys就是要拿英魂!的更多相关文章

  1. BZOJ 4453: cys就是要拿英魂![后缀数组 ST表 单调栈类似物]

    4453: cys就是要拿英魂! Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 90  Solved: 46[Submit][Status][Discu ...

  2. BZOJ.4453.cys就是要拿英魂!(后缀数组 单调栈)

    BZOJ 求字典序最大,容易想到对原串建后缀数组求\(rk\). 假设当前区间是\([l,r]\),对于在\([l,r]\)中的两个后缀\(i,j\)(\(i<j\)),显然我们不能直接比较\( ...

  3. bzoj 4453 cys就是要拿英魂!——后缀数组+单调栈+set

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4453 询问离线,按R排序. 发现直接用 rk[ ] 的错误情况就是前面的某个位置 j 和自己 ...

  4. bzoj 4453 cys就是要拿英魂! —— 后缀数组+单调栈+set

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4453 这种问题...一般先把询问离线,排序: 区间对后缀排名的影响在于一些排名大而位置靠后的 ...

  5. BZOJ 4453 cys就是要拿英魂!(后缀数组+单调栈+平衡树)

    一开始的时候感觉就是一个主席树裸题. 然后发现自己错了. 首先建出后缀数组. 设\(i<j\) 如果\(rk[i]>rk[j]\)显然i更优. 如果\(rk[i]<rk[j]\)不一 ...

  6. 【BZOJ4453】cys就是要拿英魂! 后缀数组+单调栈+set

    [BZOJ4453]cys就是要拿英魂! Description pps又开始dota视频直播了!一群每天被pps虐的蒟蒻决定学习pps的操作技术,他们把pps在这局放的技能记录了下来,每个技能用一个 ...

  7. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  8. BZOJ 2127: happiness [最小割]

    2127: happiness Time Limit: 51 Sec  Memory Limit: 259 MBSubmit: 1815  Solved: 878[Submit][Status][Di ...

  9. BZOJ 3275: Number

    3275: Number Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 874  Solved: 371[Submit][Status][Discus ...

随机推荐

  1. 乐动力APP案例

    第一部分 调研, 评测 下载软件并使用起来,描述最简单直观的个人第一次上手体验. 这款软件的主界面功能还是比较完善,里面有多个关于运动相关的数据,还有一些推荐健身教程,记录功能也十分不错,其中最难理解 ...

  2. 第四十六条:for-each循环优先于传统的for循环

    for(Elements e : list) { //doSomeThing-- }

  3. AWK读书笔记

    1.awk 'parttern {action}' filename 从文件中逐行读取并匹配parttern,若匹配成功执行action否则读取下一行. parttern和action都可选,若省略p ...

  4. vue 在methods中调用mounted中的方法?

    首先可以在data中先声明一个变量 比如 isShow=' ' mounted 中 ---> methods 中 --->  this.sureDelBox(item) 直接this调用 ...

  5. 10-TypeScript中的接口

    接口是一种规约的约定,从接口继承的类必须实现接口的约定.在高级开发中,通常接口是用于实现各种设计模式的基础,没有接口,设计模式无从谈起. 定义接口: interface ILog{ recordlog ...

  6. Oracle闪回技术

    (一)闪回技术概要 闪回技术是数据库备份与恢复的重要补充手段,主要包括以下7种特性: 特性 原理 数据库支持 闪回查询(Flashback Query) 利用undo表空间中的回退信息,查询过去某个时 ...

  7. Mego(08) - 高级建模

    对于模型建立Mego还提供了一些高级主题 数据库函数映射 我们可以将现有的CLR方法映射到指定数据库的标题函数上,如下所示 public class OrderManageEntities : DbC ...

  8. 改变this不要怕!!!

    在之前的学习和工作中,会不必要的涉及到 改变this的指向问题: 脑子里会想到的是 apply, call, bind 这三个 可逆知道他的区别吗? 1. bind区别于 apply和call 是因为 ...

  9. Let's Encrypt,站点加密之旅

    HTTPS(全称:Hyper Text Transfer Protocol over Secure Socket Layer),是以安全为目标的HTTP通道,简单讲是HTTP的安全版.即HTTP下加入 ...

  10. kubernetes进阶(05)kubernetes的命令

    在Kubernetes中,Node.Pod.Replication Controller.Service等概念都可以看作一种资源对象,通过Kubernetes提供的Kubectl工具或者API调用进行 ...