HDU 5558 Alice's Classified Message(后缀数组+二分+rmq(+线段树?))
题意
大概就是给你一个串,对于每个\(i\),在\([1,i-1]\)中找到一个\(j\),使得\(lcp(i,j)\)最长,若有多个最大\(j\)选最小,求\(j\)和这个\(lcp\)长度
思路
首先我们需要知道对于每个\(i\),能与下标小于\(i\)开头的前缀构成的最大\(lcp\)是多少
这个可以在最外层枚举\(i\)的过程中维护一个\(set\),这样在插入当前的\(rk[i]\)的时候能\(O(logn)\)得到这个最长的\(lcp\)
然后根据这个值二分出\(rk[i]\)向左右能扩展的最远的地方\([L,R]\),使得\(lcp\)依然为这个值
\([L,R]\)内的最小的\(sa\)就是答案
脑子抽了用了个线段树维护已经存在过的下标,没出现的用\(inf\)表示,但实际上无论是inf还是本身它都不可能是真正的答案,所以这边还是rmq
不过复杂度不影响,少了点常数..
代码
下面注释提供了一些自己拍出来的数据
谁让我不会自动机呢
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<stack>
#include<queue>
#include<deque>
#include<set>
#include<vector>
#include<map>
#include<functional>
#define fst first
#define sc second
#define pb push_back
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,root<<1
#define rson mid+1,r,root<<1|1
#define lc root<<1
#define rc root<<1|1
using namespace std;
typedef double db;
typedef long double ldb;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PI;
typedef pair<ll,ll> PLL;
typedef pair<ll,int> PIL;
const db eps = 1e-6;
const int mod = 1e9+7;
const int maxn = 3e5+100;
const int maxm = 2e6+100;
const int inf = 0x3f3f3f3f;
const db pi = acos(-1.0);
int t;
int n,m;
char s[maxn];
int sa[maxn],rk[maxn],height[maxn];
int y[maxn],x[maxn],c[maxn];
void getSa(){
for(int i=1;i<=n;i++)++c[x[i]=s[i]];
for(int i=2;i<=m;i++)c[i]+=c[i-1];
for(int i=n;i>=1;i--){
//printf("sa[%d]=%d\n",c[x[i]],i);
sa[c[x[i]]--]=i;
}
for(int k=1;k<=n;k<<=1){
//printf("^^^^%d\n",k);
int num = 0;
for(int i=n-k+1;i<=n;i++)y[++num]=i;
for(int i=1;i<=n;i++)if(sa[i]>k)y[++num]=sa[i]-k;
for(int i=1;i<=m;i++)c[i]=0;
for(int i=1;i<=n;i++)++c[x[i]];
for(int i=2;i<=m;i++)c[i]+=c[i-1];
for(int i=n;i>=1;i--)sa[c[x[y[i]]]--]=y[i],y[i]=0;
swap(x,y);
x[sa[1]]=1;
num=1;
for(int i = 1; i <= n; i++){
//printf("y[%d]=%d\n",i,y[i]);
}
for(int i=2;i<=n;i++){
x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k])?num:++num;
/*printf("num==%d\n",num);
printf("x[%d]=(y[%d]==y[%d]&&y[%d]==y[%d])?num:++num\n",sa[i],sa[i],sa[i-1],sa[i]+k,sa[i-1]+k);*/
}
if(num==n)break;
m=num;
}
}
void getHeight(){
int k=0;
for(int i=1; i<=n; ++i)rk[sa[i]]=i;
for(int i=1; i<=n; ++i){
if(rk[i]==1) continue;
if(k)--k;
int j=sa[rk[i]-1];
while(j+k<=n&&i+k<=n&&s[i+k]==s[j+k])++k;
height[rk[i]]=k;
}
}
set<int>S;
int d[maxn][23];
void init(){
for(int i = 1; i <= n; i++)d[i][0]=height[i];
for(int j = 1; (1<<j)<=n; j++){
for(int i = 1; i+(1<<j)-1<=n; i++){
d[i][j]=min(d[i][j-1],d[i+(1<<(j-1))][j-1]);
}
}
}
int rmq(int l, int r){
if (l > r) swap (l, r);
int k = 0;
while((1<<(k+1))<=r-l+1)k++;
return min(d[l][k],d[r-(1<<k)+1][k]);
}
int mi[maxn<<2];
void build(int l, int r, int root){
int mid = l+r>>1;
if(l==r){mi[root]=inf;return;}
build(lson);build(rson);
mi[root]=min(mi[lc],mi[rc]);
return;
}
void update(int x, int y, int l, int r, int root){
int mid = l+r>>1;
if(l==r){mi[root]=y;return;}
if(x<=mid)update(x,y,lson);
else update(x,y,rson);
mi[root]=min(mi[lc],mi[rc]);
return;
}
int ask(int x, int y, int l, int r, int root){
int mid=l+r>>1;
if(x<=l&&r<=y)return mi[root];
int ans = inf;
if(x<=mid)ans=min(ans,ask(x,y,lson));
if(y>mid)ans=min(ans,ask(x,y,rson));
return ans;
}
int ntot;
PI Ans[maxn];
int main(){
scanf("%d", &t);
int ncase = 0;
while(t--){
scanf("%s",s+1);
n = strlen(s+1);
m=122;
for(int i = 0; i <= m; i++)c[i]=0;
//for(int i = 0; i <= 2*n+2; i++)x[i]=y[i]=0;
getSa();
/*
for(int i = 1; i <= n; i++){
printf(" %d %d %d\n",i,x[i],y[i]);
}*/
getHeight();
printf("Case #%d:\n",++ncase);
S.clear();
init();
build(1,n,1);
for(int i = 1; i <= n; ){
int mxLen = 0;
set<int>::iterator it = S.lower_bound(rk[i]);
if(it!=S.end()){
mxLen=max(mxLen,rmq(rk[i]+1,(*it)));
}
if(it!=S.begin()){
it--;
mxLen=max(mxLen,rmq((*it)+1,rk[i]));
}
S.insert(rk[i]);
if(!mxLen){
printf("-1 %d\n",s[i]);
update(rk[i],i,1,n,1);
i++;
continue;
}
int L = -1,R=-1;
int l=1,r=rk[i]-1;
while(l<=r){
int mid=l+r>>1;
if(rmq(mid+1,rk[i])>=mxLen){
L=mid;r=mid-1;
}
else l=mid+1;
}
l=rk[i]+1,r=n;
while(l<=r){
int mid=l+r>>1;
int tmp = rmq(rk[i]+1,mid);
if(rmq(rk[i]+1,mid)>=mxLen){
R=mid;l=mid+1;
}
else r=mid-1;
}
int to = i+mxLen;
if(L==-1)L=inf;
else L=ask(L,rk[i],1,n,1);
if(R==-1)R=inf;
else R=ask(rk[i],R,1,n,1);
int ans=min(L,R);
printf("%d %d\n",mxLen,ans-1);
while(i<to){
S.insert(rk[i]);
update(rk[i],i,1,n,1);
i++;
if(i>n)break;
}
}
}
return 0;
}
/*
1
zlikumoyyokhaokil
stuiftqhxcksooesore
txtvacgqjtvxheeq
xnghnjzwjlpiyjlwf
qjvxxmhhbhbqingby
fsutshglayemjmshy
fxoxaepodcfvxaisda
44444
sad
a
d
sa
2
psdkfdsokfpslb
svvvcxxxxaaaa
1
svvvcxxxxaaaa
dsfwfqqqfqq
svvvcxxxxaaaa
2
aaaaabbbbbaaabbc
aaaaaa
uiuiuiuiiuiui
sad
a
d
sa
sdddqqdqqqerr
psdkfdsokfpslbl
svvvcxxxxaaaa
dsfwfqqqfqq
dsjdfssdjf
sasasaosos
doodfofsdoos
osdosofifi
soodfofofo
asoidosgkoogo
sifioobibninni
sidibvhbjjj
jjajjajjajajaj
siisisiii
iiiiiiii
a
sdiibibi
vv
sddddds
*/
HDU 5558 Alice's Classified Message(后缀数组+二分+rmq(+线段树?))的更多相关文章
- HUID 5558 Alice's Classified Message 后缀数组+单调栈+二分
http://acm.hdu.edu.cn/showproblem.php?pid=5558 对于每个后缀suffix(i),想要在前面i - 1个suffix中找到一个pos,使得LCP最大.这样做 ...
- (HDU 5558) 2015ACM/ICPC亚洲区合肥站---Alice's Classified Message(后缀数组)
题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=5558 Problem Description Alice wants to send a classi ...
- POJ 2774 Long Long Message (后缀数组+二分)
题目大意:求两个字符串的最长公共子串长度 把两个串接在一起,中间放一个#,然后求出height 接下来还是老套路,二分出一个答案ans,然后去验证,如果有连续几个位置的h[i]>=ans,且存在 ...
- Luogu3732 [HAOI2017] 供给侧改革 【后缀数组】【线段树】【乱搞】
题目分析: 这道题我是乱搞的,因为他说$01$串是随机的. 那么我们可以猜测能够让LCP变大的地方很少.求出后缀数组之后可能让LCP变大的地方就等价于从大到小往height里动态加点同时维护这个点左右 ...
- Luogu5289 十二省联考2019字符串问题(后缀数组+拓扑排序+线段树/主席树/KDTree)
先考虑80分做法,即满足A串长度均不小于B串,容易发现每个B串对应的所有A串在后缀数组上都是一段连续区间,线段树优化连边然后判环求最长链即可.场上就写了这个. 100分也没有什么本质区别,没有A串长度 ...
- BZOJ 4556: [Tjoi2016&Heoi2016]字符串(后缀数组 + 二分答案 + 主席树 + ST表 or 后缀数组 + 暴力)
题意 一个长为 \(n\) 的字符串 \(s\),和 \(m\) 个询问.每次询问有 \(4\) 个参数分别为 \(a,b,c,d\). 要你告诉它 \(s[a...b]\) 中的所有子串 和 \(s ...
- 【BZOJ3277/3473】串/字符串 后缀数组+二分+RMQ+双指针
[BZOJ3277]串 Description 字符串是oi界常考的问题.现在给定你n个字符串,询问每个字符串有多少子串(不包括空串)是所有n个字符串中至少k个字符串的子串(注意包括本身). Inpu ...
- HDU 1166 敌兵布阵 (数状数组,或线段树)
题意:... 析:可以直接用数状数组进行模拟,也可以用线段树. 代码如下: #pragma comment(linker, "/STACK:1024000000,1024000000&quo ...
- BZOJ 2780: [Spoj]8093 Sevenk Love Oimaster( 后缀数组 + 二分 + RMQ + 树状数组 )
全部串起来做SA, 在按字典序排序的后缀中, 包含每个询问串必定是1段连续的区间, 对每个询问串s二分+RMQ求出包含s的区间. 然后就是求区间的不同的数的个数(经典问题), sort queries ...
随机推荐
- 27.python中excel处理库openpyxl使用详解
openpyxl是一个第三方库,可以处理xlsx格式的Excel文件.pip install openpyxl安装. 读取Excel文件 需要导入相关函数 ? 1 2 3 from openpyxl ...
- 查找2-n之间素数的个数
题目描述 查找2-n之间素数的个数.n为用户输入值.素数:一个大于1的正整数,如果除了1和它本身以外,不能被其他正整数整除,就叫素数.如2,3,5,7,11,13,17…. 输入 整数n 输出 2-n ...
- 【转】推荐!国外程序员整理的Java资源大全
构建 这里搜集了用来构建应用程序的工具. Apache Maven:Maven使用声明进行构建并进行依赖管理,偏向于使用约定而不是配置进行构建.Maven优于Apache Ant.后者采用了一种过程化 ...
- 【转】DB2数据库编目的概念以及对其的正确解析
此文章主要向大家描述的是DB2数据库编目的概念以及对DB2数据库编目的概念的正确理解,在DB2中编目(catalog)这个单词看似很难理解,我自己当初在学习DB2数据库的时候也常常被这个编目搞的很不明 ...
- springboot2 整合redis
1.添加依赖 <dependency> <groupId>org.springframework.boot</groupId> <artifactId> ...
- cogs 2632. [HZOI 2016] 数列操作d
2632. [HZOI 2016] 数列操作d ★★★ 输入文件:segment.in 输出文件:segment.out 简单对比时间限制:3 s 内存限制:512 MB [题目描述] ...
- Java多线程的创建(一)
方法一:继承Thread类实现 1.创建一个类A,并继承Thread类 2.重写A的run()方法 3.创建A的实例对象b,即创建了线程对象 4.使用b调用start()方法:启动线程(会自动调用ru ...
- 20191102Java课堂记录
1. import javax.swing.*; class AboutException { public static void main(String[] a) { int i=1, j=0, ...
- Lobooi个人作业:阅读与准备作业
这个作业属于哪个课程 https://edu.cnblogs.com/campus/xnsy/GeographicInformationScience 这个作业要求在哪里 https://edu.cn ...
- python,finally的应用
脚本执行过程中可能因为被测试的环境有改变导致中间某一部分无法继续执行下去 可以在最后一行加上finally来执行最后一句脚本 比如 最后执行退出 表示 无论中间过程失败还是成功,最终都会执行退出操作 ...