●BZOJ 3879 SvT
题链:
http://www.lydsy.com/JudgeOnline/problem.php?id=3879
题解:
后缀数组,单调栈,RMQ
其实类似 BZOJ 3238 [Ahoi2013]差异
只是不过变成了多次询问某些后缀两两之间的LCP的和。
所以对于每次询问,就先把那些后缀提出来,
按 rank排序后,求出相邻的后缀的 LCP(存储在 H 数组中),
然后再按那个题的做法做就好了(即求出 H 数组对应的 L[],R[])。
读入数据有点大,用了读入优化。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 505000
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
char S[MAXN];
int sa[MAXN],rak[MAXN],hei[MAXN],L[MAXN],R[MAXN],stm[MAXN][20],log2[MAXN];
char gc(){
//return getchar();
static char IN[MAXN];
static int bit=MAXN-5000,p=0,len=0;
if(p>=len) len=fread(IN,1,bit,stdin),IN[len]=EOF,p=0;
return IN[p++];
}
void read(int &x){
static int f;
static char ch;
x=0; f=1; ch=gc();
while(ch<'0'||'9'<ch){if(ch=='-')f=-1;ch=gc();}
while('0'<=ch&&ch<='9'){x=x*10+ch-'0';ch=gc();}
x=x*f;
}
void build(int N,int M){
static int cc[MAXN],ta[MAXN],tb[MAXN],*x,*y,h,p;
x=ta; y=tb; h=0;
for(int i=0;i<M;i++) cc[i]=0;
for(int i=0;i<N;i++) cc[x[i]=S[i]]++;
for(int i=1;i<M;i++) cc[i]+=cc[i-1];
for(int i=N-1;i>=0;i--) sa[--cc[x[i]]]=i;
for(int k=1;p=0,k<N;k<<=1){
for(int i=N-k;i<N;i++) y[p++]=i;
for(int i=0;i<N;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
for(int i=0;i<M;i++) cc[i]=0;
for(int i=0;i<N;i++) cc[x[y[i]]]++;
for(int i=1;i<M;i++) cc[i]+=cc[i-1];
for(int i=N-1;i>=0;i--) sa[--cc[x[y[i]]]]=y[i];
swap(x,y); y[N]=-1; x[sa[0]]=0; M=1;
for(int i=1;i<N;i++)
x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?M-1:M++;
if(M>=N) break;
}
for(int i=0;i<N;i++) rak[sa[i]]=i;
for(int i=0,j;i<N;i++){
if(h) h--;
if(rak[i]){
j=sa[rak[i]-1];
while(S[i+h]==S[j+h]) h++;
}
stm[rak[i]][0]=hei[rak[i]]=h;
}
for(int k=1;k<=log2[N];k++)
for(int i=(1<<k)-1;i<N;i++)
stm[i][k]=min(stm[i-(1<<(k-1))][k-1],stm[i][k-1]);
}
int LCP(int l,int r){
static int k;
if(l>r) swap(l,r); l++;
k=log2[r-l+1];
return min(stm[l+(1<<k)-1][k],stm[r][k]);
}
void preLR(int N,int *val){
static int stk[MAXN],stp[MAXN],top;
stp[top=0]=0;
for(int i=0;i<N;i++){
while(top&&stk[top]>=val[i]) top--;
L[i]=stp[top]; top++;
stk[top]=val[i]; stp[top]=i;
}
stp[top=N]=N;
for(int i=N-1;i>=0;i--){
while(top<N&&stk[top]>val[i]) top++;
R[i]=stp[top]-1; top--;
stk[top]=val[i]; stp[top]=i;
}
}
int main()
{
//filein(3879);
long long ans;
log2[1]=0; for(int i=2;i<=500000;i++) log2[i]=log2[i>>1]+1;
static int A[MAXN*10],H[MAXN*10],N,M;
read(N); read(M);
for(int i=0;i<N;i++)
do{S[i]=gc();}while(S[i]<'a'||'z'<S[i]);
S[N]=0;build(N,300);
for(int i=0,t;i<M;i++){
ans=0;
read(t); for(int j=0;j<t;j++) read(A[j]),A[j]=rak[A[j]-1];
sort(A,A+t);
t=unique(A,A+t)-A;
for(int i=1;i<t;i++) H[i]=LCP(A[i-1],A[i]);
preLR(t,H);
for(int i=1;i<t;i++)
ans+=1ll*(i-L[i])*(R[i]-i+1)*H[i];
printf("%lld\n",ans);
}
return 0;
}
●BZOJ 3879 SvT的更多相关文章
- BZOJ 3879: SvT [虚树 后缀树]
传送门 题意: 多次询问,给出一些后缀,求两两之间$LCP$之和 哈哈哈哈哈哈哈竟然$1A$了,刚才还在想如果写不好这道题下节数学就不上了,看来是上天让我上数学课啊 $Suffix\ Virtual\ ...
- bzoj 3879: SvT
Description (我并不想告诉你题目名字是什么鬼) 有一个长度为n的仅包含小写字母的字符串S,下标范围为[1,n]. 现在有若干组询问,对于每一个询问,我们给出若干个后缀(以其在S中出现的起始 ...
- BZOJ 3879: SvT 虚树 + 后缀自动机
Description (我并不想告诉你题目名字是什么鬼) 有一个长度为n的仅包含小写字母的字符串S,下标范围为[1,n]. 现在有若干组询问,对于每一个询问,我们给出若干个后缀(以其在S中出现的起始 ...
- 【BZOJ 3879】SvT
http://www.lydsy.com/JudgeOnline/problem.php?id=3879 SvT的中文是后缀虚树? 反正本蒟蒻不懂,还是$O(nlogn)$的后缀数组和单调栈维护来做, ...
- bzoj 3879 虚树
题目大意: 给一个字符串,多次询问k个后缀,求它们两两间LCP长度总和 分析: 转化为后缀树,用虚树求 注意: 后缀树中代表后缀的点都是叶子节点 题目中取模并没有卵用 #include <cst ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- Week Four
2018.12.18 1.[USACO Platinum C] 2.[Gym 102028H] 3.[BZOJ 2750] 4.[BZOJ 3238] 5.[BZOJ 4310] 6.[BZOJ 38 ...
- Week Five
2018.12.25 1.[BZOJ 4310] 2.[BZOJ 3879] 3.[BZOJ 2754] 4.[BZOJ 4698] 5.[Codeforces 914E] 6.[Codeforces ...
- HDU 3879 && BZOJ 1497:Base Station && 最大获利 (最大权闭合图)
http://acm.hdu.edu.cn/showproblem.php?pid=3879 http://www.lydsy.com/JudgeOnline/problem.php?id=1497 ...
随机推荐
- 2017-2018-1 我爱学Java 第六七周 作业
团队六七周作业 完善版需求规格说明书 制定团队编码规范 数据库设计 后端架构设计 TODOList 参考资料 完善版需求规格说明书 <需求规格说明书>初稿不足之处: 1.开发工具写错 2. ...
- Django 分类标签查找
from django.conf.urls import url from django.contrib import admin from blog.views import index,stude ...
- javascript实现小鸟飞行轨迹
javascript实现小鸟飞行轨迹 代码如下:
- php的开发的apache的配置及伪静态的应用
1.Apache之所以能够解析php代码是游览器首先发送数据到模版页面,然后模版页提交数据到php页面,然后php代码经过Apache解析过后生成结果的,所以是 在Apache的配置文件中是可以看到开 ...
- Step by Step 真正从零开始,TensorFlow详细安装入门图文教程!帮你完成那个最难的从0到1
摘要: Step by Step 真正从零开始,TensorFlow详细安装入门图文教程!帮你完成那个最难的从0到1 安装遇到问题请文末留言. 悦动智能公众号:aibbtcom AI这个概念好像突然就 ...
- restful架构风格设计准则(六)版本管理
读书笔记,原文链接:http://www.cnblogs.com/loveis715/p/4669091.html,感谢作者! 版本管理 在前面已经提到过,一个REST系统为资源所抽象出的URI实际上 ...
- css3中的动画 @keyframes animation
动画的运用比较重要.接下来我希望针对我自己学习遇到的问题,再总结一下这个属性的使用方法. 创建一个动画: @keyframes 动画名 {样式} 引用自己创建的动画: animation:动画名 时 ...
- Python/Django-Web原理(一)
Python/Django-Web原理(一) websocket webSocket协议是基于TCP的一种新的协议.WebSocket最初在HTML规范中被引用为TCP连接,作为基于TCP的套接字AP ...
- 面向对象进阶---attr家族
一. isinstance(obj,cls)和issubclass(sub,super) isinstance(obj,cls)检查obj是否是类 cls 的对象 class Foo: pass ob ...
- 前端学习之jquery/下
前端学习之jquery 一 属性操作 html(): console.log($("div").html()); $(".test").html("& ...