P4094 [HEOI2016/TJOI2016]字符串 [SA + 主席树]

传送门

转换为判断，二分一个len，判断后缀起点在a -- b-mid+1 直接存不存在S，

使lcp(S,c)>= len，我们发现rank[S] 一定在c左右的一个范围内, 我们二分这个范围，用主席树判断就可以了

#include<bits/stdc++.h>
#define N 100050
using namespace std;int n,m,q;
char s[N];
int SA[N],rank[N],Height[N],y[N],c[N],tmp[N];
int st[N][20],Log[N];
struct Node{int ls,rs,val;
}t[N*30]; int rt[N],tot;void Sort(){for(int i=1;i<=m;i++) c[i] = 0;for(int i=1;i<=n;i++) c[rank[i]]++;for(int i=1;i<=m;i++) c[i] += c[i-1];for(int i=n;i>=1;i--) SA[c[rank[y[i]]]--] = y[i];
}void get_SA(){for(int i=1;i<=n;i++) rank[i] = s[i], y[i] = i;Sort();for(int k=1;k<=n;k<<=1){int ret = 0;for(int i=n-k+1;i<=n;i++) y[++ret] = i;for(int i=1;i<=n;i++) if(SA[i] > k) y[++ret] = SA[i] - k;Sort(); swap(rank,tmp); rank[SA[1]] = 1; int num = 1;for(int i=2;i<=n;i++){if(tmp[SA[i]] == tmp[SA[i-1]] && tmp[SA[i]+k] == tmp[SA[i-1]+k])rank[SA[i]] = num;else rank[SA[i]] = ++num;} m = num;}
}
void get_Height(){int k = 0;for(int i=1;i<=n;i++){if(rank[i]==1) continue;int j = SA[rank[i]-1]; if(k) k--;while(i+k<=n && j+k<=n && s[i+k]==s[j+k]) k++;Height[rank[i]] = k;}
} void get_ST(){for(int i=2;i<=n;i++) Log[i] = Log[i/2] + 1;for(int i=1;i<=n;i++) st[i][0] = Height[i];for(int i=1;(1<<i)<=n;i++)for(int j=1;j+(1<<i)-1<=n;j++)st[j][i] = min(st[j][i-1], st[j+(1<<(i-1))][i-1]);
}
int Lcp(int l,int r){l++; int x = Log[r-l+1];return min(st[l][x], st[r-(1<<x)+1][x]);
}
void Build(int &x,int l,int r){x = ++tot;if(l==r) return;int mid = (l+r) >> 1;Build(t[x].ls, l, mid);Build(t[x].rs, mid+1, r);
}
void Insert(int &x,int last,int l,int r,int pos){x = ++tot; t[x] = t[last]; t[x].val++;if(l==r) return;int mid = (l+r) >> 1;if(pos<=mid) Insert(t[x].ls, t[last].ls, l, mid, pos);else Insert(t[x].rs, t[last].rs, mid+1, r, pos);
}
int check(int a,int b,int l,int r,int L,int R){if(L<=l && r<=R) return t[b].val - t[a].val;int mid = (l+r) >> 1, ans = 0;if(L<=mid) ans += check(t[a].ls, t[b].ls, l, mid, L, R);if(R>mid) ans += check(t[a].rs, t[b].rs, mid+1, r, L, R);return ans;
}void get_LR(int st,int len,int &L,int &R){int l=1,r=rank[st];while(l<r){int mid = (l+r) >> 1;if(Lcp(mid,rank[st])>=len) r = mid;else l = mid+1;} L = l;l = rank[st], r = n;while(l<r){int mid = (l+r+1) >> 1;if(Lcp(rank[st],mid)>=len) l = mid;else r = mid-1;} R = l;
}
int calc(int a,int b,int c,int d){int l = 0, r = min(b-a, d-c) + 1;while(l<r){int mid = (l+r+1) >> 1;int L = 0, R = 0; get_LR(c,mid,L,R);if(check(rt[a-1],rt[b-mid+1],1,n,L,R)) l = mid;else r = mid-1;} return l;
}
int main(){scanf("%d%d%s",&n,&q,s+1);m = 127; get_SA(); get_Height(); get_ST();Build(rt[0],1,n); for(int i=1;i<=n;i++) Insert(rt[i], rt[i-1], 1, n, rank[i]);while(q--){int a,b,c,d;scanf("%d%d%d%d",&a,&b,&c,&d);printf("%d\n",calc(a,b,c,d));} return 0;
}

P4094 [HEOI2016/TJOI2016]字符串 [SA + 主席树]相关推荐

洛谷P4094 [HEOI2016/TJOI2016]字符串【后缀数组+主席树+st表】
时空限制 2000ms / 256MB 题目描述佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一个长为n的字符串s,和m个问题.佳媛姐姐必须正确 ...
主席树 + 后缀数组求LCP + 二分套二分 ---- P4094 [HEOI2016/TJOI2016]字符串
题目链接题目大意: 解题思路: 设我们的答案为midmidmid(注意这里有坑是[a,b][a,b][a,b]的所有子串和[c,d][c,d][c,d]这个子串的最长lcplcplcp),那么我们会 ...
洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树
题目链接后缀数组题目分析: sa[i] – 第i小的后缀的编号 rank[i] --编号为i的后缀排第几: height[i] – 第i和第i-1的最长lcp最长公共前缀: 1.二分答案,答案肯定 ...
[HEOI2016/TJOI2016]字符串（后缀数组+主席树+二分）
description 佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一个长为 n 的字符串 s,和 m 个问题.佳媛姐姐必须正确回答这 m 个问 ...
HEOI2016/TJOI2016 字符串问题
题目链接:戳我非常不好意思,因为想要排版,所以今天先只把代码贴出来,明天补题解. 40pts暴力:直接暴力匹配 #include<iostream>#include<cstring ...
[BZOJ4556][TJOI2016HEOI2016]字符串(二分答案+后缀数组+RMQ+主席树)
4556: [Tjoi2016&Heoi2016]字符串 Time Limit: 20 Sec Memory Limit: 128 MB Submit: 1360 Solved: 545 ...
洛谷P4094 - [TJOI2016]字符串
Portal Description 给出一个字符串\(s(|s|\leq10^5)\)和\(m\)次询问,每次询问子串\(s[x_1..x_2]\)的所有子串和\(s[y_1..y_2]\)的最长公 ...
BZOJ3473:字符串(后缀数组,主席树,二分,ST表)
Description 给定n个字符串,询问每个字符串有多少子串(不包括空串)是所有n个字符串中至少k个字符串的子串? Input 第一行两个整数n,k. 接下来n行每行一个字符串. Output 一 ...
[2020.11.26NOIP模拟赛]勇者的后缀【SA,RMQ,主席树,二分】
正题题目链接:https://www.luogu.com.cn/problem/U142356?contestId=37784 题目大意一个字符串,询问给出(x,l,r)(x,l,r)(x,l,r ...

P4094 [HEOI2016/TJOI2016]字符串 [SA + 主席树]

P4094 [HEOI2016/TJOI2016]字符串 [SA + 主席树]相关推荐

最新文章

热门文章