loj 6031「雅礼集训 2017 Day1」字符串
时间:2019-09-12
本文章向大家介绍loj 6031「雅礼集训 2017 Day1」字符串,主要包括loj 6031「雅礼集训 2017 Day1」字符串使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
注意到每次询问串长度都是给定的,并且询问串长\(k*\)询问次数\(q<10^5\),所以这里面一个东西大的时候另一个东西就小,那么考虑对较小的下功夫
如果\(k\le \sqrt{n}\),那么可以\(O(k^2)\)暴力枚举询问串的每一个子串,然后在\(S\)的sam找到这个子串对应的点,算出出现次数,并且乘上在区间\([a,b]\)中这个子串询问的出现次数.找到子串对应的点为了方便,可以依次让询问串的某个前缀在sam上匹配,然后按长度从大到小枚举前缀的后缀,从匹配位置开始倍增跳父亲,直到当前点包含的串的长度区间包含当前长度;然后子串在区间\([a,b]\)询问次数可以主席树(雾).所以这部分复杂度是\(O(qk^2logn)=O(n\sqrt{n}logn)\)
如果\(k> \sqrt{n}\),那么\(q\le \sqrt{n}\),所以可以枚举\(m\)个询问区间对应的子串,分别算答案.具体实现可以参考上面,依次用前缀在sam上匹配,然后处理右端点为当前位置的询问子串贡献.这部分是\(O(qmlogn)=O(n\sqrt{n}logn)\)
然后不知道为什么我的倍增被卡了...应该是我写丑了吧qwq.然后改成暴力跳父亲就跑的飞快
#include<bits/stdc++.h>
#define LL long long
#define uLL unsigned long long
#define db double
using namespace std;
const int N=1e5+10;
int rd()
{
int x=0,w=1;char ch=0;
while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
int n,m,q,kk,lz;
char cc[N],ss[N];
struct SAM
{
int fa[N<<1][19],len[N<<1],ch[N<<1][26],sz[N<<1],la,tt;
vector<int> e[N<<1];
SAM(){la=tt=1;}
void extd(int cx)
{
int np=++tt,p=la;
len[np]=len[p]+1,sz[np]=1,la=tt;
while(!ch[p][cx]) ch[p][cx]=np,p=fa[p][0];
if(!p) fa[np][0]=1;
else
{
int q=ch[p][cx];
if(len[q]==len[p]+1) fa[np][0]=q;
else
{
int nq=++tt;
fa[nq][0]=fa[q][0],len[nq]=len[p]+1;
memcpy(ch[nq],ch[q],sizeof(int)*26);
fa[np][0]=fa[q][0]=nq;
while(ch[p][cx]==q) ch[p][cx]=nq,p=fa[p][0];
}
}
}
void dfs(int x)
{
for(int j=1;j<=lz;++j)
{
fa[x][j]=fa[fa[x][j-1]][j-1];
if(!fa[x][j]) break;
}
vector<int>::iterator it;
for(it=e[x].begin();it!=e[x].end();++it)
{
int y=*it;
dfs(y),sz[x]+=sz[y];
}
}
void inii()
{
for(int i=2;i<=tt;++i) e[fa[i][0]].push_back(i);
dfs(1);
}
}sam;
namespace ct1
{
int id[320][320],tti;
int s[N*20],ch[N*20][2],rt[N],tt;
void inst(int o1,int o2,int x)
{
int l=1,r=tti;
s[o1]=s[o2]+1;
while(l<r)
{
int mid=(l+r)>>1;
if(x<=mid)
{
ch[o1][0]=++tt,ch[o1][1]=ch[o2][1];
o1=ch[o1][0],o2=ch[o2][0];
r=mid;
}
else
{
ch[o1][0]=ch[o2][0],ch[o1][1]=++tt;
o1=ch[o1][1],o2=ch[o2][1];
l=mid+1;
}
s[o1]=s[o2]+1;
}
}
int quer(int o1,int o2,int l,int r,int lx)
{
if(!s[o1]) return 0;
if(l==r) return s[o1]-s[o2];
int mid=(l+r)>>1;
if(lx<=mid) return quer(ch[o1][0],ch[o2][0],l,mid,lx);
else return quer(ch[o1][1],ch[o2][1],mid+1,r,lx);
}
void wk()
{
for(int i=1;i<=kk;++i)
for(int j=i;j<=kk;++j)
id[i][j]=++tti;
for(int i=1;i<=m;++i)
{
int l=rd()+1,r=rd()+1;
inst(rt[i]=++tt,rt[i-1],id[l][r]);
}
while(q--)
{
scanf("%s",ss+1);
int a=rd()+1,b=rd()+1;
LL ans=0;
for(int i=1,p=1,ll=0;i<=kk;++i)
{
while(p>1&&!sam.ch[p][ss[i]-'a']) p=sam.fa[p][0],ll=sam.len[p];
if(sam.ch[p][ss[i]-'a'])
{
p=sam.ch[p][ss[i]-'a'],++ll;
int x=p,nm;
for(int j=i-ll+1;j<=i;++j)
if((nm=quer(rt[b],rt[a-1],1,tti,id[j][i])))
{
if(sam.len[sam.fa[x][0]]+1>i-j+1)
{
for(int k=lz;~k;--k)
if(sam.len[sam.fa[sam.fa[x][k]][0]]+1>i-j+1)
x=sam.fa[x][k];
x=sam.fa[x][0];
}
ans+=1ll*sam.sz[x]*nm;
}
}
else p=1,ll=0;
}
printf("%lld\n",ans);
}
}
}
namespace ct2
{
struct node
{
int l,i;
bool operator < (const node &bb) const {return l<bb.l;}
};
vector<node> qr[N];
void wk()
{
for(int i=1;i<=m;++i)
{
int l=rd()+1,r=rd()+1;
qr[r].push_back((node){l,i});
}
for(int i=1;i<=n;++i) sort(qr[i].begin(),qr[i].end());
while(q--)
{
scanf("%s",ss+1);
int a=rd()+1,b=rd()+1;
LL ans=0;
for(int i=1,p=1,ll=0;i<=kk;++i)
{
while(p>1&&!sam.ch[p][ss[i]-'a']) p=sam.fa[p][0],ll=sam.len[p];
if(sam.ch[p][ss[i]-'a'])
{
p=sam.ch[p][ss[i]-'a'],++ll;
int x=p;
vector<node>::iterator it;
for(it=qr[i].begin();it!=qr[i].end();++it)
{
int j=(*it).l,ii=(*it).i;
if(i-j+1<=ll&&ii>=a&&ii<=b)
{
/*if(sam.len[sam.fa[x][0]]+1>i-j+1)
{
for(int k=lz;~k;--k)
if(sam.len[sam.fa[sam.fa[x][k]][0]]+1>i-j+1)
x=sam.fa[x][k];
x=sam.fa[x][0];
}*/
while(sam.len[sam.fa[x][0]]+1>i-j+1) x=sam.fa[x][0];
ans+=sam.sz[x];
}
}
}
else p=1,ll=0;
}
printf("%lld\n",ans);
}
}
}
int main()
{
n=rd(),m=rd(),q=rd(),kk=rd();
scanf("%s",cc+1);
for(int i=1;i<=n;++i) sam.extd(cc[i]-'a');
lz=log2(n),sam.inii();
if(kk<=315) ct1::wk();
else ct2::wk();
return 0;
}
原文地址:https://www.cnblogs.com/smyjr/p/11512886.html
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 为PXC集群引入Mycat并构建完整的高可用集群架构
- Python3 字典
- 安装Percona Server数据库(in CentOS 8)
- Python 基础篇-正斜杠("/")和反斜杠("")的用法
- 在CentOS8下搭建PXC集群
- Python 基础篇-相对路径、绝对路径的写法
- Python3 元组
- 关于MySQL的基准测试
- Python 技术篇-操作excel,对excel进行读取和写入
- Mycat 整合 MySQL 8.x 踩坑实践
- Python 技术篇-xlwt库不新建,直接读取已存在的excel并写入
- Python3 列表
- Mycat 核心配置详解
- Python字符串
- 初识 HBase