【NOI2019模拟2019.6.29】字符串(SA|SAM+主席树)
Description:
1<=n<=5e4
题解:
考虑\(f\)这个东西应该是怎样算的?
不妨建出SA,然后按height从大到小启发式合并,显然只有相邻的才可能成为最优答案。这样的只有\(O(n log n)\)个有用的串。
建SAM在fail树上启发式合并是一样的。
然后用个主席树就可以快速查询答案。
现在思考查询一个[x,y],要求f>=z怎么办?
考虑一个区间[l,r],如果a[l-1]<=max[a[l..r]]或a[r+1]<=max[a[l..r]]显然延伸会使f更大,而a不会更大。
这样从每个点开始造区间,就能造出n个区间,答案显然是这些区间中的一个。
不过有个问题,就是可能这些区间过长出界了,注意出界的话一定有一个端点是x或者y,可以二分+主席树查询。
后面看上去也要一个三维偏序,实际上不用,假设二分出来最左的y'使f[x,y']>=z,和最右的x'使f[x',y]>=z,那么只用查询l∈[x,x']或者r∈[y',y]的那些区间,这样就变成了二维偏序。
Code:
#include<bits/stdc++.h>
#define fo(i, x, y) for(int i = x, B = y; i <= B; i ++)
#define ff(i, x, y) for(int i = x, B = y; i < B; i ++)
#define fd(i, x, y) for(int i = x, B = y; i >= B; i --)
#define ll long long
#define pp printf
#define hh pp("\n")
using namespace std;
const int N = 5e4 + 5;
int n, m;
char s[N];
int a[N];
struct SA {
int rk[N], sa[N], he[N], tx[N], tp[N], m;
void rsort() {
fo(i, 1, m) tx[i] = 0;
fo(i, 1, n) tx[rk[tp[i]]] ++;
fo(i, 1, m) tx[i] += tx[i - 1];
fd(i, n, 1) sa[tx[rk[tp[i]]] --] = tp[i];
}
int cmp(int *f, int x, int y, int z) { return f[x] == f[y] && f[x + z] == f[y + z];}
void build() {
s[0] = s[n + 1] = -1; tp[n + 1] = 0;
fo(i, 1, n) rk[i] = s[i], tp[i] = i;
m = 127; rsort();
for(int w = 1, p = 0; p < n; m = p, w *= 2) {
p = 0; fo(i, n - w + 1, n) tp[++ p] = i;
fo(i, 1, n) if(sa[i] > w) tp[++ p] = sa[i] - w;
rsort();
fo(i, 1, n) tp[i] = rk[i];
rk[sa[1]] = p = 1;
fo(i, 2, n) rk[sa[i]] = cmp(tp, sa[i - 1], sa[i], w) ? p : ++ p;
}
int j, k = 0;
for(int i = 1; i <= n; he[rk[i ++]] = k)
for(k ? k -- : 0, j = sa[rk[i] - 1]; s[j + k] == s[i + k]; k ++);
// fo(i, 1, n) {
// fo(j, sa[i], n) pp("%c", s[j]);
// pp(" %d\n", he[i]);
// }
}
} suf;
struct P {
int x, y, z;
} d[N * 40]; int d0;
namespace make_d {
multiset<int> s[N];
multiset<int> :: iterator it;
int f[N], t[N];
int F(int x) { return f[x] == x ? x : (f[x] = F(f[x]));}
int cmp(int x, int y) { return suf.he[x] > suf.he[y];}
void ins(int x, multiset<int> &s, int z) {
if(*s.begin() < x) d[++ d0] = (P) {*(--s.lower_bound(x)), x, z};
if(*s.rbegin() > x) d[++ d0] = (P) {x, *s.upper_bound(x), z};
}
void build() {
fo(i, 1, n) f[i] = i, s[i].insert(suf.sa[i]);
fo(i, 2, n) t[++ t[0]] = i;
sort(t + 1, t + n, cmp);
fo(i, 1, n - 1) {
int x = t[i] - 1, y = t[i];
x = F(x), y = F(y);
if(x != y) {
if(s[x].size() > s[y].size()) swap(x, y);
for(it = s[x].begin(); it != s[x].end(); it ++) {
ins(*it, s[y], suf.he[t[i]]);
}
for(it = s[x].begin(); it != s[x].end(); it ++)
s[y].insert(*it);
s[x].clear();
f[x] = y;
}
}
}
}
int cmp_d(P a, P b) { return a.x > b.x;}
struct tree {
int l, r, x, y;
} t[N * 400];
#define i0 t[i].l
#define i1 t[i].r
int tot, pl, pr, pc, px, g[N];
void upd(int i) {
t[i].x = max(t[i0].x, t[i1].x);
t[i].y = max(t[i0].y, t[i1].y);
}
void add(int &i, int x, int y) {
t[++ tot] = t[i]; i = tot;
if(x == y) {
if(!pc) {
t[i].x = max(t[i].x, px);
} else {
t[i].y = max(t[i].y, px);
}
return;
}
int m = x + y >> 1;
if(pl <= m) add(i0, x, m); else add(i1, m + 1, y);
upd(i);
}
void ft(int i, int x, int y) {
if(y < pl || x > pr || !i) return;
if(x >= pl && y <= pr) {
// pp("%d %d %d %d %d\n", x, y, t[i].x, t[i].y, pc);
if(!pc) px = max(px, t[i].x); else px = max(px, t[i].y);
return;
}
int m = x + y >> 1;
ft(i0, x, m); ft(i1, m + 1, y);
}
void make_tree() {
t[0].x = t[0].y = -1e9;
sort(d + 1, d + d0 + 1, cmp_d);
int l = 1;
fd(i, n, 1) {
g[i] = g[i + 1];
while(l <= d0 && d[l].x >= i) {
if(d[l].z > 0) {
// pp("%d %d %d %d\n", i, d[l].x, d[l].y, d[l].z);
pl = pr = d[l].y + d[l].z - 1;
pc = 0; px = d[l].z;
add(g[i], 1, n);
pl = pr = d[l].y + d[l].z - 1;
pc = 1; px = -d[l].y + 1;
add(g[i], 1, n);
}
l ++;
}
}
// fo(i, 1, d0) pp("%d %d %d\n", d[i].x, d[i].y, d[i].z); hh;
}
int query(int x, int y) {
pl = x; pr = y; px = -1e9; pc = 0;
ft(g[x], 1, n);
int ans = px;
pl = y; pr = n; px = -1e9; pc = 1;
ft(g[x], 1, n);
ans = max(ans, y + px);
return ans;
}
namespace tr {
struct tree {
int l, r, x, y;
} t[N * 40];
#define i0 t[i].l
#define i1 t[i].r
int tot, pl, pr, pc, px, rt;
void upd(int i) {
t[i].x = min(t[i0].x, t[i1].x);
t[i].y = min(t[i0].y, t[i1].y);
}
void add(int &i, int x, int y) {
if(!i) t[++ tot] = t[0], i = tot;
if(x == y) {
if(!pc) {
t[i].x = min(t[i].x, px);
} else {
t[i].y = min(t[i].y, px);
}
return;
}
int m = x + y >> 1;
if(pl <= m) add(i0, x, m); else add(i1, m + 1, y);
upd(i);
}
void ft(int i, int x, int y) {
if(y < pl || x > pr || !i) return;
if(x >= pl && y <= pr) {
if(!pc) px = min(px, t[i].x); else px = min(px, t[i].y);
return;
}
int m = x + y >> 1;
ft(i0, x, m); ft(i1, m + 1, y);
}
void xiu(int x, int y, int z) {
pl = pr = x; px = y; pc = z;
add(rt, 1, n);
}
int find(int x, int y, int z) {
px = 1e9; pl = x, pr = y; pc = z;
ft(rt, 1, n);
return px;
}
}
int f[17][N];
int qam(int x, int y) {
int l = log2(y - x + 1);
return max(f[l][x], f[l][y - (1 << l) + 1]);
}
int l[N], r[N], v[N], ti[N], z[N], z0;
int p[N], p0;
int cmp_p(int x, int y) { return ti[x] > ti[y];}
void make_dkl() {
fo(i, 1, n) f[0][i] = a[i];
fo(j, 1, 16) fo(i, 1, n) f[j][i] = max(f[j - 1][i], f[j - 1][i + (1 << j - 1)]);
fo(i, 1, n) l[i] = 1, r[i] = n;
z0 = 0;
fo(i, 1, n) {
while(z0 > 0 && a[z[z0]] <= a[i]) z0 --;
if(z0 > 0) l[i] = z[z0] + 1;
z[++ z0] = i;
}
z0 = 0;
fd(i, n, 1) {
while(z0 > 0 && a[z[z0]] <= a[i]) z0 --;
if(z0 > 0) r[i] = z[z0] - 1;
z[++ z0] = i;
}
tr :: t[0].x = tr :: t[0].y = 1e9;
fo(i, 1, n) v[i] = qam(l[i], r[i]), ti[i] = query(l[i], r[i]), p[i] = i;
sort(p + 1, p + n + 1, cmp_p);
}
int ans[N];
struct ask {
int x, y, z, i;
int as1, as2;
} q[N];
int cmp_q(ask a, ask b) { return a.z > b.z;}
void Ask() {
fo(i, 1, m) {
int x, y, z;
scanf("%d %d %d", &q[i].x, &q[i].y, &q[i].z);
q[i].i = i;
x = q[i].x, y = q[i].y, z = q[i].z;
int as1 = -1, as2 = -1;
for(int l = x, r = y; l <= r; ) {
int m = l + r >> 1;
if(query(x, m) >= z) as1 = m, r = m - 1; else l = m + 1;
}
if(as1 == -1) ans[i] = -1;
for(int l = x, r = y; l <= r; ) {
int m = l + r >> 1;
if(query(m, y) >= z) as2 = m, l = m + 1; else r = m - 1;
}
if(as1 != -1) ans[i] = min(qam(x, as1), qam(as2, y));
q[i].as1 = as1; q[i].as2 = as2;
}
sort(q + 1, q + m + 1, cmp_q);
int w = 1;
fo(i, 1, m) {
while(w <= n && ti[p[w]] >= q[i].z) {
tr :: xiu(l[p[w]], v[p[w]], 0);
tr :: xiu(r[p[w]], v[p[w]], 1);
w ++;
}
if(q[i].as1 != -1) {
ans[q[i].i] = min(ans[q[i].i], tr :: find(q[i].x, q[i].as2, 0));
ans[q[i].i] = min(ans[q[i].i], tr :: find(q[i].as1, q[i].y, 1));
}
}
fo(i, 1, m) pp("%d\n", ans[i]);
}
int main() {
freopen("string.in", "r", stdin);
freopen("string.out", "w", stdout);
scanf("%d %d", &n, &m);
scanf("%s", s + 1);
fo(i, 1, n) scanf("%d", &a[i]);
suf.build();
make_d :: build();
make_tree();
make_dkl();
Ask();
}
【NOI2019模拟2019.6.29】字符串(SA|SAM+主席树)的更多相关文章
- [JZOJ6241]【NOI2019模拟2019.6.29】字符串【数据结构】【字符串】
Description 给出一个长为n的字符串\(S\)和一个长为n的序列\(a\) 定义一个函数\(f(l,r)\)表示子串\(S[l..r]\)的任意两个后缀的最长公共前缀的最大值. 现在有q组询 ...
- 【NOI2019模拟2019.6.29】组合数(Lucas定理、数位dp)
Description: p<=10且p是质数,n<=7,l,r<=1e18 题解: Lucas定理: \(C_{n}^m=C_{n~mod~p}^{m~mod~p}*C_{n/p} ...
- BZOJ3473:字符串(后缀数组,主席树,二分,ST表)
Description 给定n个字符串,询问每个字符串有多少子串(不包括空串)是所有n个字符串中至少k个字符串的子串? Input 第一行两个整数n,k. 接下来n行每行一个字符串. Output 一 ...
- BZOJ4556:[TJOI\HEOI2016]字符串(后缀数组,主席树,二分,ST表)
Description 佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一个长为n的字符串s,和m个问题.佳媛姐姐必须正确回答这m个问题,才能打开箱 ...
- [BZOJ4556][Tjoi2016&Heoi2016]字符串 后缀数组+主席树
4556: [Tjoi2016&Heoi2016]字符串 Time Limit: 20 Sec Memory Limit: 128 MB Description 佳媛姐姐过生日的时候,她的小 ...
- P4094 [HEOI2016/TJOI2016]字符串 后缀数组+主席树+二分答案
$ \color{#0066ff}{ 题目描述 }$ 佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一个长为n的字符串s,和m个问题.佳媛姐姐必须 ...
- [HEOI2016] 字符串 - 后缀数组,主席树,ST表,二分
[HEOI2016] 字符串 Description 给定一个字符串 \(S\), 有 \(m\) 个询问,每个询问给定参数 \((a,b,c,d)\) ,求 \(s[a..b]\) 的子串与 \(s ...
- luogu4770 [NOI2018]你的名字 (SAM+主席树)
对S建SAM,拿着T在上面跑 跑的时候不仅无法转移要跳parent,转移过去不在范围内也要跳parent(注意因为范围和长度有关,跳的时候应该把长度一点一点地缩) 这样就能得到对于T的每个前缀,它最长 ...
- bzoj 4556 [Tjoi2016&Heoi2016]字符串——后缀数组+主席树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4556 本来只要查 ht[ ] 数组上的前驱和后继就行,但有长度的限制.可以二分答案解决!然后 ...
随机推荐
- cocos2D-X 屏幕适配
{ //https://v.youku.com/v_show/id_XNTIzOTM1MDYw.html }
- 纯js的右下角弹窗
<html> <head> <title></title> <meta charset="UTF-8"> <scr ...
- 用JOptionPane类实现各种对话框
用JOptionPane类实现各种对话框 运行结果: 下面部分参考: JOptionPane类提示框的一些常用的方法 - - ITeye博客 http://847353020-qq-com.itey ...
- MySQL concat、concat_ws 和 group_concat 的用法
一.CONCAT()函数CONCAT()函数用于将多个字符串连接成一个字符串.使用数据表Info作为示例,其中SELECT id,name FROM info LIMIT 1;的返回结果为+----+ ...
- mysql存储过程、函数、触发器、
当数据库版本不允许直接使用存储过程.函数的语法时用delimiter // 将结束符改成//用完之后再写delimiter;将结束符改回来即可,调用过程.函数用call+其名字即可返回结果 delim ...
- poi提取docx中的文字和图片
package com.fry.poiDemo.dao; import java.io.File; import java.io.FileInputStream; import java.io.Fil ...
- ansible 远程执行时提示 command not found 问题
问题 最近在学习 ansible ,在使用普通用户远程执行 ip a 命令是发现提示错误:/bin/sh: ip: command not found. 原因 command not found 命令 ...
- 前端学习笔记——CSS选择器
学习css选择器之前我们先了解下css规则: CSS 规则由两个主要的部分构成:选择器,以及一条或多条声明: 选择器通常是需要改变样式的 HTML 元素. 每条声明由一个属性和一个值组成. 1.id选 ...
- shell编程:变量替换
定义变量:$ var_1="I love you, Do you love me" 输出变量:$ echo $var_1 打印结果:I love you, Do you love ...
- 使用Nuget重新安装packages.config中的组件的方法
Update-Package -ProjectName 'Ko.app.web' -Reinstall 该语句作用:按照packages.config中给出的程序组件,重新下载安装一遍.