【LOJ】#2059. 「TJOI / HEOI2016」字符串
题解
我们冷静一下,先画一棵后缀树
然后发现我们要给c和d这一段区间在[a,b]这一段开头的串里找lcp
而lcp呢,就是c点的祖先的到根的一段,假如这个祖先的子树里有[a,b - dis[u] + 1],那么这个u就是合法的,维护每个点子树里出现过的后缀串起点可以用线段树合并
(这里的深度指后缀树上该点到根的距离)
我们先用min(b - a + 1,c - d + 1)限制一下这个祖先的深度,选择距离根深度最近且深度>= min(b - a + 1,c - d + 1),找的这个点称为p
那么我们要找最长的合法的部分,可以二分,二分的log,倍增找点的log,线段树里判断是否合法的log,一共是三个log,当然,我们冷静一下,显然有些log常数非常小嘛
但是我们再冷静一下,我们可以发现,我们要找的点可以当做c的祖先中,深度最小的,且子树中不存在[a,b - dis[u] + 1]的点u,我们可以从高位到低位枚举二进制位,按位判断让p跳还是不跳,这就是两个log了
同时这个高位到低位枚举二进制位也可用于找到p
那么这个时候难道p的父亲的深度就是答案了吗?too naive
是p的深度到p的父亲的深度都可能,我们需要再次二分,二分的左边界p的父亲的深度,右边界是min(p的深度 - 1,min(b - a + 1,d - c + 1))
复杂度最后还是\(O(n \log^2 n)\)
代码
#include <bits/stdc++.h>
#define enter putchar('\n')
#define space putchar(' ')
#define pii pair<int,int>
#define fi first
#define se second
#define MAXN 200005
#define pb push_back
#define mp make_pair
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 + c - '0';
c = getchar();
}
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) out(x / 10);
putchar('0' + x % 10);
}
int N,M;
char s[MAXN];
int id[MAXN];
struct node {
int to,next,val;
}E[MAXN * 4];
int head[MAXN],sumE,Ncnt,fa[MAXN][20],dis[MAXN],pos[MAXN];
void add(int u,int v,int c) {E[++sumE].to = v;E[sumE].next = head[u];E[sumE].val = c;head[u] = sumE;}
namespace SAM {
struct node {
node *nxt[26],*par;
int len,cnt;
}pool[MAXN * 2],*tail = pool,*root,*last,*que[MAXN * 2];
int c[MAXN];
void Init() {
root = last = tail++;
}
void build_sam(int l,int c) {
node *nowp = tail++,*p;
nowp->len = l;nowp->cnt = 1;
for(p = last ; p && !p->nxt[c]; p = p->par) {
p->nxt[c] = nowp;
}
if(!p) nowp->par = root;
else {
node *q = p->nxt[c];
if(q->len == p->len + 1) nowp->par = q;
else {
node *copyq = tail++;
*copyq = *q;
copyq->len = p->len + 1;copyq->cnt = 0;
q->par = nowp->par = copyq;
for(; p && p->nxt[c] == q ; p = p->par) {
p->nxt[c] = copyq;
}
}
}
last = nowp;
}
void build_suffix_tree() {
Ncnt = tail - pool;
for(int i = 0 ; i < Ncnt ; ++i) {
c[pool[i].len]++;
}
for(int i = 1 ; i <= N ; ++i) c[i] += c[i - 1];
for(int i = 0 ; i < Ncnt ; ++i) {
que[c[pool[i].len]--] = &pool[i];
}
for(int i = 1 ; i <= Ncnt ; ++i) {
int u = que[i] - pool + 1;
if(que[i]->par) {
int f = que[i]->par - pool + 1;
fa[u][0] = f;
add(f,u,que[i]->len - que[i]->par->len);
add(u,f,que[i]->len - que[i]->par->len);
}
if(que[i]->cnt) {id[N - que[i]->len + 1] = u;pos[u] = N - que[i]->len + 1;}
}
}
}
namespace seg_tr {
struct node {
int lc,rc;
}tr[MAXN * 80];
int rt[MAXN],Ncnt;
void Insert(int x,int &u,int L,int R,int pos) {
u = ++Ncnt;
tr[u] = tr[x];
if(L == R) return;
int mid = (L + R) >> 1;
if(pos <= mid) Insert(tr[x].lc,tr[u].lc,L,mid,pos);
else Insert(tr[x].rc,tr[u].rc,mid + 1,R,pos);
}
int Merge(int Lt,int Rt) {
if(!Lt) return Rt;
if(!Rt) return Lt;
int u = ++Ncnt;
tr[u].lc = Merge(tr[Lt].lc,tr[Rt].lc);
tr[u].rc = Merge(tr[Lt].rc,tr[Rt].rc);
return u;
}
bool Query(int u,int L,int R,int l,int r) {
if(!u) return false;
if(L == l && R == r) return 1;
int mid = (L + R) >> 1;
if(r <= mid) return Query(tr[u].lc,L,mid,l,r);
else if(l > mid) return Query(tr[u].rc,mid + 1,R,l,r);
else return Query(tr[u].lc,L,mid,l,mid) || Query(tr[u].rc,mid + 1,R,mid + 1,r);
}
}
using seg_tr::rt;
using seg_tr::Merge;
using seg_tr::Insert;
using seg_tr::Query;
void dfs(int u) {
if(pos[u]) {Insert(rt[u],rt[u],1,N,pos[u]);}
for(int i = head[u] ; i ; i = E[i].next) {
int v = E[i].to;
if(v != fa[u][0]) {
dis[v] = dis[u] + E[i].val;
dfs(v);
rt[u] = Merge(rt[u],rt[v]);
}
}
}
void Init() {
read(N);read(M);
scanf("%s",s + 1);
reverse(s + 1,s + N + 1);
SAM::Init();
for(int i = 1 ; i <= N ; ++i) {
SAM::build_sam(i,s[i] - 'a');
}
SAM::build_suffix_tree();
dfs(1);
for(int j = 1 ; j <= 17 ; ++j) {
for(int i = 1 ; i <= Ncnt ; ++i) {
fa[i][j] = fa[fa[i][j - 1]][j - 1];
}
}
}
void Solve() {
int a,b,c,d;
for(int i = 1 ; i <= M ; ++i) {
read(a);read(b);read(c);read(d);
int l = min(d - c + 1,b - a + 1);
int p = id[c];
for(int j = 17 ; j >= 0 ; --j) {
if(!fa[p][j]) continue;
if(dis[fa[p][j]] >= l) p = fa[p][j];
}
if(Query(rt[p],1,N,a,b - min(dis[p],l) + 1)) {
out(min(dis[p],l));enter;continue;
}
for(int j = 17 ; j >= 0 ; --j) {
if(!fa[p][j]) continue;
int u = fa[p][j];
if(!Query(rt[u],1,N,a,b - dis[u] + 1)) p = fa[p][j];
}
if(dis[p] == dis[fa[p][0]] + 1 || dis[fa[p][0]] >= l) {
p = fa[p][0];
out(min(dis[p],l));enter;
}
else {
int L = dis[fa[p][0]],R = min(dis[p] - 1,l);
while(L < R) {
int mid = (L + R + 1) >> 1;
if(Query(rt[p],1,N,a,b - mid + 1)) L = mid;
else R = mid - 1;
}
out(L);enter;
}
}
}
int main() {
Init();
Solve();
return 0;
}
我的代码能力真是低到要哭了!
真的哭了!
写跪了三次啊,debug到吐啊。。。
【LOJ】#2059. 「TJOI / HEOI2016」字符串的更多相关文章
- loj#2059. 「TJOI / HEOI2016」字符串 sam+线段树合并+倍增
题意:给你一个子串,m次询问,每次给你abcd,问你子串sa-b的所有子串和子串sc-d的最长公共前缀是多长 题解:首先要求两个子串的最长公共前缀就是把反过来插入变成最长公共后缀,两个节点在paren ...
- loj#2054. 「TJOI / HEOI2016」树
题目链接 loj#2054. 「TJOI / HEOI2016」树 题解 每次标记覆盖整棵字数,子树维护对于标记深度取max dfs序+线段树维护一下 代码 #include<cstdio> ...
- 「TJOI / HEOI2016」字符串
「TJOI / HEOI2016」字符串 题目描述 佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一个长为 \(n\) 的字符串 \(s\),和 ...
- loj #2055. 「TJOI / HEOI2016」排序
#2055. 「TJOI / HEOI2016」排序 题目描述 在 2016 年,佳媛姐姐喜欢上了数字序列.因而他经常研究关于序列的一些奇奇怪怪的问题,现在他在研究一个难题,需要你来帮助他. 这个 ...
- LOJ #2058「TJOI / HEOI2016」求和
不错的推柿子题 LOJ #2058 题意:求$\sum\limits_{i=0}^n\sum\limits_{j=0}^nS(i,j)·2^j·j!$其中$ S(n,m)$是第二类斯特林数 $ Sol ...
- loj2059 「TJOI / HEOI2016」字符串
字符串好难啊不会啊 #include <iostream> #include <cstdio> using namespace std; int n, m, rnk[10000 ...
- AC日记——#2057. 「TJOI / HEOI2016」游戏 LOJ
#2057. 「TJOI / HEOI2016」游戏 思路: 最大流: 代码: #include <cstdio> #include <cstring> #include &l ...
- loj2058 「TJOI / HEOI2016」求和 NTT
loj2058 「TJOI / HEOI2016」求和 NTT 链接 loj 思路 \[S(i,j)=\frac{1}{j!}\sum\limits_{k=0}^{j}(-1)^{k}C_{j}^{k ...
- AC日记——#2054. 「TJOI / HEOI2016」树
#2054. 「TJOI / HEOI2016」树 思路: 线段树: 代码: #include <cstdio> #include <cstring> #include < ...
随机推荐
- Hi35xx 通用GPIO 使用篇(板子3G电源控制脚说明)
在一个嵌入式系统中使用最多的莫过于 通用输入输出 GPIO口.看到论坛中经常有朋友问海思为什么没有提供GPIO驱动.其实不然. 在海思SDK xxx/osdrv/tools/board_tools/ ...
- jQuery技巧笔记
1.关于页面元素的引用 通过jquery的$()引用元素包括通过id.class.元素名以及元素的层级关系及dom或者xpath条件等方法,且返回的对象为jquery对象(集合对象),不能直接调用do ...
- HDU 1017 A Mathematical Curiosity 数学题
解题报告:输入两个数,n和m,求两个数a和b满足0<a<b<n,并且(a^2+b^2+m) % (a*b) =0,这样的a和b一共有多少对.注意这里的b<n,并不可以等于n. ...
- transparent 透明效果
background-color:transparent;就是把背景色设置为透明. 实际上background默认的颜色就是透明的属性.所以写和不写都是一样的 span{ width: 0; heig ...
- 让arch阻止某个软件包的升级
我更新了eclipse-java Mars版本的,感觉特别的卡,而且还有好多bug,不知道为什么,因此我去官网下载了luna版本的eclipse的安装包,不知道怎么下载的点击这里,然后安装luna版本 ...
- 对接微信支付使用HMAC-SHA256使用签名算法实现方式
最近做微信押金支付对接,很多坑,心累!这里提醒一下各位: 首先,确保自己商户号进了白名单,没有需要联系客服,否则接口是调不通的,会一直提示参数错误 其次,确保接口文档是最新的,最好去官网去看,否则可能 ...
- linux下使用indent整理代码(代码格式化)【转】
转自:https://blog.csdn.net/jiangjingui2011/article/details/7197069 常用的设置: indent -npro -kr -i8 -ts8 -s ...
- linux中core dump开启使用教程【转】
转自:http://www.111cn.net/sys/linux/67291.htm 一.什么是coredump 我们经常听到大家说到程序core掉了,需要定位解决,这里说的大部分是指对应程序由于各 ...
- Tomcat安装与优化
Tomcat安装与优化 1.安装jdk环境 最新的JDK下载地址:http://www.oracle.com/technetwork/java/javase/downloads/jdk8-downlo ...
- 002_安装第三方APP好的站点及解除安全与隐私限制
一.解除安全与隐私限制的任何来源. http://bbs.feng.com/read-htm-tid-10714286.html 接下来,我们就打开终端,然后输入以下命令: sudo spctl ...