洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树

题目链接
后缀数组
题目分析：
sa[i] – 第i小的后缀的编号
rank[i] --编号为i的后缀排第几；
height[i] – 第i和第i-1的最长lcp最长公共前缀；

1.二分答案，答案肯定在0—min(b-a+1, d-c+1);二分答案mid;
2.s[c–d], 满足条件的子串一定在rank[c] 开始向前后扩展的一个区间[L, R]中，（因为和某一个串排名越接近，那么与这个串的lcp(最长公共前缀)就越大）。
3.这个区间满足的条件就是这个区间任意一点i的height[i] >= mid; 那么这个区间中两两后缀的lcp一定大于等于mid。
4如果在这个区间存在i，sa[i] 编号大于等于a,并且sa[i] + mid - 1 <= b，那么这个mid就是可行的，可以继续二分下去，

5，这个区间L-R使用主席树来维持：

#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
#include<map>
using namespace std;
typedef long long lt;int read()
{int f=1,x=0;char ss=getchar();while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}return f*x;
}const int maxn=100010;
int n,m,Q;
int a,b,c,d;
int rak[maxn],sa[maxn],tp[maxn],tax[maxn];
int height[maxn],mi[maxn][20];
int lft[maxn<<5],rht[maxn<<5],rt[maxn<<5],sum[maxn<<5],cnt;
char ss[maxn];void rsort()
{for(int i=0;i<=m;++i) tax[i]=0;for(int i=1;i<=n;++i) tax[rak[i]]++;for(int i=1;i<=m;++i) tax[i]+=tax[i-1];for(int i=n;i>=1;--i) sa[tax[rak[tp[i]]]--]=tp[i];
}void ssort()
{m=127;for(int i=1;i<=n;++i)rak[i]=ss[i-1],tp[i]=i;rsort(); for(int k=1;k<=n;k<<=1){int p=0;for(int i=n-k+1;i<=n;++i) tp[++p]=i;for(int i=1;i<=n;++i) if(sa[i]>k) tp[++p]=sa[i]-k;rsort(); swap(rak,tp);rak[sa[1]]=p=1;for(int i=2;i<=n;++i)rak[sa[i]]=(tp[sa[i]]==tp[sa[i-1]]&&tp[sa[i]+k]==tp[sa[i-1]+k])?p:++p;if(p>=n)break;m=p;}
}void getH()
{int k=0;for(int i=1;i<=n;++i){if(k) k--;int j=sa[rak[i]-1];while(ss[i+k-1]==ss[j+k-1]) k++;height[rak[i]]=k;}
}void RMQ()
{for(int i=1;i<=n;++i) mi[i][0]=height[i];for(int j=1;(1<<j)<=n;++j)for(int i=1;i+(1<<j)-1<=n;++i)mi[i][j]=min(mi[i][j-1],mi[i+(1<<j-1)][j-1]);
}int qmin(int ll,int rr)
{int k=0;while((1<<k+1)<=rr-ll+1) k++;return min(mi[ll][k],mi[rr-(1<<k)+1][k]);
}int update(int pre,int ll,int rr,int x)
{int tt=++cnt; sum[tt]=sum[pre]+1;lft[tt]=lft[pre]; rht[tt]=rht[pre];int mid=ll+rr>>1;if(ll<rr){if(x<=mid) lft[tt]=update(lft[pre],ll,mid,x);else rht[tt]=update(rht[pre],mid+1,rr,x);}return tt;
}int query(int u,int v,int ll,int rr,int s,int t)
{if(ll<=s&&t<=rr) return sum[v]-sum[u];int res=0,mid=s+t>>1;if(ll<=mid) res+=query(lft[u],lft[v],ll,rr,s,mid);if(rr>mid) res+=query(rht[u],rht[v],ll,rr,mid+1,t);return res;
}int qpos1(int ll,int rr,int x)
{int L=ll,R=rr,res=rr;while(L<R){int mid=L+R>>1;if(qmin(mid+1,rr)>=x) R=mid,res=mid;else L=mid+1;}return res;
}int qpos2(int ll,int rr,int x)
{int L=ll,R=rr,res=ll;while(L<R){int mid=L+R>>1;if(qmin(ll+1,mid)>=x) L=mid+1,res=mid;else R=mid;}return res;
}int check(int x)
{   int L=qpos1(1,rak[c],x);int R=qpos2(rak[c],n,x);if(query(rt[L-1],rt[R],a,b+1-x,1,n)>0) return 1;else return 0;
}int main()
{n=read();Q=read();scanf("%s",ss);ssort(); getH(); RMQ();for(int i=1;i<=n;++i)rt[i]=update(rt[i-1],1,n,sa[i]);while(Q--){a=read();b=read();c=read();d=read();int L=1,R=min(b-a+1,d-c+1),ans=0;while(L<R){int mid=L+R>>1;if(check(mid)) L=mid+1,ans=mid;else R=mid;}if(L==R) if(check(L)) ans=L;printf("%d\n",ans);}return 0;
}

洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树相关推荐

洛谷P4094 [HEOI2016/TJOI2016]字符串【后缀数组+主席树+st表】
时空限制 2000ms / 256MB 题目描述佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了一个长为n的字符串s,和m个问题.佳媛姐姐必须正确 ...
洛谷P4559 [JSOI2018]列队【70分二分 + 主席树】
题目链接洛谷P4559 题解只会做$70$分的$O(nlog^2n)$ 如果本来就在区间内的人是不用动的,区间右边的人往区间最右的那些空位跑,区间左边的人往区间最左的那些空位跑找到这些空 ...
主席树 + 后缀数组求LCP + 二分套二分 ---- P4094 [HEOI2016/TJOI2016]字符串
题目链接题目大意: 解题思路: 设我们的答案为midmidmid(注意这里有坑是[a,b][a,b][a,b]的所有子串和[c,d][c,d][c,d]这个子串的最长lcplcplcp),那么我们会 ...
洛谷P3763 [Tjoi2017]DNA 【后缀数组】
题目链接洛谷P3763 题解后缀数组裸题在BZOJ被卡常到哭QAQ #include<algorithm> #include<iostream> #include< ...
BZOJ4553/洛谷P4093 [HEOI2016/TJOI2016]序列动态规划分治
原文链接http://www.cnblogs.com/zhouzhendong/p/8672434.html 题目传送门 - BZOJ4553 题目传送门 - 洛谷P4093 题解设$Li$表示第$ ...
P4094 [HEOI2016/TJOI2016]字符串 [SA + 主席树]
传送门转换为判断,二分一个len,判断后缀起点在a -- b-mid+1 直接存不存在S, 使lcp(S,c)>= len,我们发现rank[S] 一定在c左右的一个范围内, 我们二分这个范围 ...
[BZOJ4556][TJOI2016HEOI2016]字符串(二分答案+后缀数组+RMQ+主席树)
4556: [Tjoi2016&Heoi2016]字符串 Time Limit: 20 Sec Memory Limit: 128 MB Submit: 1360 Solved: 545 ...
洛谷 P4091 [HEOI2016/TJOI2016]求和任意模数fft
题目描述在2016年,佳媛姐姐刚刚学习了第二类斯特林数,非常开心. 现在他想计算这样一个函数的值: f(n)=∑ni=0∑ij=0S(i,j)×2j×(j!) f ( n ) = ∑ i = 0 n ...
BZOJ 2780: [Spoj]8093 Sevenk Love Oimaster( 后缀数组 + 二分 + RMQ + 树状数组 )
全部串起来做SA, 在按字典序排序的后缀中, 包含每个询问串必定是1段连续的区间, 对每个询问串s二分+RMQ求出包含s的区间. 然后就是求区间的不同的数的个数(经典问题), sort queries ...

洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树

洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树相关推荐

最新文章

热门文章

洛谷 P4094 [HEOI2016/TJOI2016]字符串 后缀数组+二分+主席树

洛谷 P4094 [HEOI2016/TJOI2016]字符串 后缀数组+二分+主席树相关推荐

最新文章

热门文章

洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树

洛谷 P4094 [HEOI2016/TJOI2016]字符串后缀数组+二分+主席树相关推荐