bzoj4453
单调栈+set+后缀数组
一道奇妙的题
这道题如果对于每个询问$r$是固定的,那么就很简单了,可惜并不是
由于r会变化,那么对于两个子串$[i...r],[j...r]$,他们的大小关系随着r的变化也会变化,使得我们不能直接预处理答案
所以我们把询问离线,把每个询问按照r分类,通过考虑r的变化来完成询问
对于两个后缀$[i...r],[j...r],(i<j)$,设他们的$lcp=l$,那么当$l∈[j,j+l-1]$时,$[i...r]>[j...r]$,当$r=j+l$时,则变成$[i...r]<[j...r]$,因为$s[i+l]<s[j+l]$(假设)
我们称这样的两个后缀为$i$伴随$j$。
那么我们对答案维护一个set,里面存着一些后缀,保证后缀当前的大小单调递减,这样我们就可以每次快速询问答案。
但是由于$r$的变化每次大小关系会改变,也就是我们要把一些后缀删掉
那么我们再维护一个单调栈,这里存的是严格的后缀大小,这个单调栈是为了给每个新加入的后缀求出哪些后缀伴随他
那么如果当前的栈顶是大于新的后缀,$break$,否则标记$top$伴随$i$,并且记录当$i$,也就是$r=i+lcp$时删除栈顶
每次删除一个元素时,伴随他的元素也应该删除,所以每次$dfs$删除即可,并在$set$里删除
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
#include <vector>
using namespace std;
const int N = 1e5 + ;
int n, k, m, top;
int rnk[N], tmp[N], sa[N], ans[N], vis[N], st[N], h[N][], Log[N];
char s[N];
set<int> S;
vector<pair<int, int> > v[N];
vector<int> del[N], g[N];
bool cmp(int i, int j) {
if(rnk[i] != rnk[j]) {
return rnk[i] < rnk[j];
}
int ri = i + k <= n ? rnk[i + k] : -;
int rj = j + k <= n ? rnk[j + k] : -;
return ri < rj;
}
void dfs(int u) {
vis[u] = ;
S.erase(u);
for(int i = ; i < g[u].size(); ++i) {
if(!vis[g[u][i]]) {
dfs(g[u][i]);
}
}
}
void SA() {
for(int i = ; i <= n; ++i) {
sa[i] = i;
rnk[i] = s[i];
}
for(k = ; k <= n; k <<= ) {
sort(sa + , sa + n + , cmp);
tmp[sa[]] = ;
for(int i = ; i <= n; ++i) {
tmp[sa[i]] = tmp[sa[i - ]] + (cmp(sa[i - ], sa[i]));
}
for(int i = ; i <= n; ++i) {
rnk[i] = tmp[i];
}
}
for(int i = ; i <= n; ++i) {
rnk[sa[i]] = i;
}
int d = ;
for(int i = ; i <= n; ++i) {
if(rnk[i] <= ) {
continue;
}
int j = sa[rnk[i] - ];
if(d) {
--d;
}
for(; i + d <= n && j + d <= n; ++d) {
if(s[i + d] != s[j + d]) {
break;
}
}
h[rnk[i] - ][] = d;
}
for(int j = ; j <= ; ++j) {
for(int i = ; i + ( << j) - <= n; ++i) {
h[i][j] = min(h[i][j - ], h[i + ( << (j - ))][j - ]);
}
}
}
int rmq(int l, int r) {
if(l == r) {
return n - l + ;
}
l = rnk[l];
r = rnk[r];
if(l > r) {
swap(l, r);
}
--r;
int x = Log[r - l + ];
return min(h[l][x], h[r - ( << x) + ][x]);
}
int main() {
scanf("%s%d", s + , &m);
n = strlen(s + );
for(int i = ; i <= n; ++i) {
Log[i] = Log[i >> ] + ;
}
SA();
for(int i = ; i <= m; ++i) {
int l, r;
scanf("%d%d", &l, &r);
v[r].push_back(make_pair(l, i));
}
for(int i = ; i <= n; ++i) {
S.insert(i);
while(top) {
int lcp = rmq(st[top], i);
if(s[st[top] + lcp] > s[i + lcp]) {
break;
}
del[i + lcp].push_back(st[top]);
g[i].push_back(st[top]);
--top;
}
st[++top] = i;
for(int j = ; j < del[i].size(); ++j) {
if(!vis[del[i][j]]) {
dfs(del[i][j]);
}
}
for(int j = ; j < v[i].size(); ++j) {
pair<int, int> x = v[i][j];
ans[x.second] = *S.lower_bound(x.first);
}
}
for(int i = ; i <= m; ++i) {
printf("%d\n", ans[i]);
}
return ;
}
bzoj4453的更多相关文章
- 【BZOJ4453】cys就是要拿英魂! 后缀数组+单调栈+set
[BZOJ4453]cys就是要拿英魂! Description pps又开始dota视频直播了!一群每天被pps虐的蒟蒻决定学习pps的操作技术,他们把pps在这局放的技能记录了下来,每个技能用一个 ...
随机推荐
- 【解决】无法连接 MKS:套接字连接尝试次数太多正在放弃
https://blog.csdn.net/wjunsing/article/details/78496224 我的电脑 -> 右键 -> 管理 -> 服务和应用程序 -> 服 ...
- 如何查看selenium api文档
参考文章:https://www.cnblogs.com/yoyoketang/p/6189740.html 环境:windows + python3 + selenium3 打开cmd,执行命令:p ...
- Ejabberd作为推送服务的优化手段(转)
AVOS Cloud目前还在用Ejabberd做Android的消息推送服务.当时选择Ejabberd,是因为Ejabberd是一个发展很长时间的XMPP实现,并且基于Erlang,设想能在我们自主研 ...
- 多线程快速解压FastZipArchive介绍
本文转载至 http://blog.csdn.net/xunyn/article/details/12975937 多线程解压iosfast 在iOS项目中用到解压缩,用的是ZipArchive ...
- smarty模板里实现缓存。
smarty模板里实现缓存.分页缓存在任何里都可以用 我用了三个类 include("../init.inc.php");//模板入口类 include("../DBDA ...
- android菜鸟学习笔记15----Android Junit测试
Android中的Junit测试与Java Junit测试有所不同,不能简单的使用标注…… 假设写了一个MathUtils类,有两个静态方法: public class MathUtils { pub ...
- 查看SqlServer安装的log文件
SqlServer安装时产生的log被保存在这个目录下: "%programfiles%\Microsoft SQL Server\[SQL_VERSION]\Setup Bootstrap ...
- 我的Android进阶之旅------>Android中可替换string的使用,getString(int resId, Object... formatArgs)
官方文档如下描述: 地址:http://developer.android.com/reference/android/content/Context.html#getString%28int,%20 ...
- c语言高速推断一个数是偶数还是奇数
#include <stdio.h> int main() { int a; while(1) { printf("please input the number:\n" ...
- iOS Dev (53) 修复UIImagePickerController偷换StatusBar颜色的问题
版权声明:本文为 CSDN 博主 大锐哥(ID 为 prevention)原创文章,未经博主同意不得转载. https://blog.csdn.net/prevention/article/detai ...