【[SCOI2012]喵星球上的点名】
好题啊
\(SA+ST\text{表}+\text{莫队}\)
我们先强行把所有的串连起来,串与串之间插入特殊字符,姓和名之间也插入特殊字符
之后跑一遍\(SA\),求出\(sa\)和\(het\)
对于所有的询问串,我们标记好他们的开头,之后我们对于排好序的后缀建一个\(st\)表,找到每个询问串往左往右最多可以扩展到哪里
扩展到哪里自然是这个区间内的\(het\)的最小值大于等于询问串的长度了,这个二分+\(st\)表就能做到
第一问就变成了区间数颜色了,这个套上一个莫队就可以了
第二问是求出一种颜色被数了多少次,我们在莫队的时候顺便打上一个时间戳,如果在一个询问里这个颜色第一次出现我们就标记好这个询问的时间,等到下一次这个颜色被完全删除的时候我们利用这个操作的时间做一个差就好了
代码
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define maxn 500005
#define re register
#define LL long long
#define lowbit(x) ((x)&(-x))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read()
{
char c=getchar();re int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
int tax[maxn],rk[maxn],sa[maxn],het[maxn],tp[maxn],f[maxn],S[maxn];
int to[maxn],c[maxn],St[maxn][19],logg[maxn],Ans[maxn],li[maxn],cnt[maxn];
int n,m,L,M,len,sz,ans;
inline void qsort()
{
for(re int i=0;i<=M;i++) tax[i]=0;
for(re int i=1;i<=L;i++) tax[rk[i]]++;
for(re int i=1;i<=M;i++) tax[i]+=tax[i-1];
for(re int i=L;i;--i) sa[tax[rk[tp[i]]]--]=tp[i];
}
struct Ask{int x,y,rk;} a[maxn];
inline void add(int x,int o) {if(!f[sa[x]]) return;if(!cnt[f[sa[x]]]) ++ans,c[f[sa[x]]]+=m-o+1;cnt[f[sa[x]]]++;}
inline void del(int x,int o) {if(!f[sa[x]]) return;cnt[f[sa[x]]]--;if(!cnt[f[sa[x]]]) --ans,c[f[sa[x]]]-=m-o+1;}
inline int query(int l,int r){int k=logg[r-l+1];return min(St[l][k],St[r-(1<<k)+1][k]);}
inline int cmp(Ask A,Ask B) { if(A.x/sz==B.x/sz) return A.y<B.y; return A.x<B.x;}
int main()
{
n=read(),m=read();
for(re int i=1;i<=n;i++)
{
len=read();
for(re int j=1;j<=len;j++) S[++L]=read()+1,f[L]=i; S[++L]=0;
len=read();
for(re int j=1;j<=len;j++) S[++L]=read()+1,f[L]=i; S[++L]=0;
}
for(re int i=1;i<=m;i++)
{
len=read();to[i]=L+1;li[i]=len;
for(re int j=1;j<=len;j++) S[++L]=read()+1; S[++L]=0;
}
for(re int i=1;i<=L;i++) M=max(M,S[i]);sz=M;
for(re int i=1;i<=L;i++) if(!S[i]) S[i]=++sz;
for(re int i=1;i<=L;i++) M=max(M,S[i]),rk[i]=S[i],tp[i]=i;
qsort();
for(re int w=1,p=0;p<L;M=p,w<<=1)
{
p=0;
for(re int i=1;i<=w;i++) tp[++p]=L-w+i;
for(re int i=1;i<=L;i++) if(sa[i]>w) tp[++p]=sa[i]-w;
qsort();
for(re int i=1;i<=L;i++) std::swap(rk[i],tp[i]);
rk[sa[1]]=p=1;
for(re int i=2;i<=L;i++)
rk[sa[i]]=(tp[sa[i-1]]==tp[sa[i]]&&tp[sa[i-1]+w]==tp[sa[i]+w])?p:++p;
}
int k=0;
for(re int i=1;i<=L;i++)
{
if(k) --k; int j=sa[rk[i]-1];
while(S[i+k]==S[j+k]) ++k; het[rk[i]]=k;
}
memset(St,20,sizeof(St));
for(re int i=1;i<=L;i++) St[i][0]=het[i];
for(re int i=2;i<=L;i++) logg[i]=1+logg[i>>1];
for(re int j=1;j<=18;j++)
for(re int i=1;i+(1<<j)-1<=L;i++)
St[i][j]=min(St[i][j-1],St[i+(1<<(j-1))][j-1]);
for(re int i=1;i<=m;i++)
{
int j=rk[to[i]];a[i].rk=i;
int l=2,r=j;
while(l<=r) {int mid=l+r>>1;if(query(mid,j)>=li[i]) a[i].x=mid-1,r=mid-1;else l=mid+1;}
if(!a[i].x) a[i].x=j;
l=j+1,r=L;
while(l<=r) {int mid=l+r>>1;if(query(j+1,mid)>=li[i]) a[i].y=mid,l=mid+1;else r=mid-1;}
if(!a[i].y) a[i].y=j;
}
sz=std::sqrt(L);
std::sort(a+1,a+m+1,cmp);
int l=0,r=0;
for(re int i=1;i<=m;i++)
{
while(l<a[i].x) del(l++,i);
while(r<a[i].y) add(++r,i);
while(l>a[i].x) add(--l,i);
while(r>a[i].y) del(r--,i);
Ans[a[i].rk]=ans;
}
for(re int i=1;i<=m;i++) printf("%d\n",Ans[i]);
for(re int i=1;i<=n;i++) printf("%d ",c[i]);
return 0;
}
【[SCOI2012]喵星球上的点名】的更多相关文章
- BZOJ 2754: [SCOI2012]喵星球上的点名
2754: [SCOI2012]喵星球上的点名 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 649 Solved: 305[Submit][Sta ...
- BZOJ2754: [SCOI2012]喵星球上的点名
2754: [SCOI2012]喵星球上的点名 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 680 Solved: 314[Submit][Sta ...
- BZOJ 2754: [SCOI2012]喵星球上的点名 [后缀数组+暴力]
2754: [SCOI2012]喵星球上的点名 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1906 Solved: 839[Submit][St ...
- BZOJ 2754: [SCOI2012]喵星球上的点名 [AC自动机+map+暴力]
2754: [SCOI2012]喵星球上的点名 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1902 Solved: 837[Submit][St ...
- P2336 [SCOI2012]喵星球上的点名(后缀自动机+莫队+dfs序)
P2336 [SCOI2012]喵星球上的点名 名字怎么存?显然是后缀自动机辣 询问点到多少个喵喵喵其实就是 查询后缀自动机上parent树的一个子树 于是我们考虑莫队 怎么树上莫队呢 我们用dfs序 ...
- 洛咕 P2336 [SCOI2012]喵星球上的点名
洛咕 P2336 [SCOI2012]喵星球上的点名 先求出SA和height,一个点名串对应的就是一段区间,还有很多个点,就转化成了 有很多个区间,很多个点集,对每个区间计算和多少个点集有交,对每个 ...
- 洛谷 P2336 [SCOI2012]喵星球上的点名 解题报告
P2336 [SCOI2012]喵星球上的点名 题目描述 a180285 幸运地被选做了地球到喵星球的留学生.他发现喵星人在上课前的点名现象非常有趣. 假设课堂上有 \(N\) 个喵星人,每个喵星人的 ...
- 【BZOJ2754】[SCOI2012]喵星球上的点名
[BZOJ2754][SCOI2012]喵星球上的点名 题面 bzoj 洛谷 题解 这题有各种神仙做法啊,什么暴力\(AC\)自动机.\(SAM\)等等五花八门 我这个蒟蒻在这里提供一种复杂度正确且常 ...
- 2754. [SCOI2012]喵星球上的点名【后缀数组】
Description a180285幸运地被选做了地球到喵星球的留学生.他发现喵星人在上课前的点名现象非常有趣. 假设课堂上有N个喵星人,每个喵星人的名字由姓和名构成.喵星球上的老师会选择M个串 ...
- SCOI2012喵星球上的点名
http://codevs.cn/problem/2403/ 2012年省队选拔赛四川 时间限制: 2 s 空间限制: 128000 KB 题目描述 Description a180285幸运 ...
随机推荐
- Java线程池详解(二)
一.前言 在总结了线程池的一些原理及实现细节之后,产出了一篇文章:Java线程池详解(一),后面的(一)是在本文出现之后加上的,而本文就成了(二).因为在写完第一篇关于java线程池的文章之后,越发觉 ...
- 小程序wxml文件引用方式
import<!--header.wxml --> <template name="item"> <text>{{text}}</text ...
- 01-oracle限定查询-20190404
关系型数据库和半结构化数据(xml文件) oracle12c:c代表云计算 PDB,CDB sql语句执行顺序: 第一步:from子句控制数据来源: 第二步:where子句使用限定符对数据行过滤: 第 ...
- c# 跨平台ide JetBrains Rider
https://www.jetbrains.com/rider/ et框架 调试hotfix用的,说是vs调试容易崩溃 破解方法 https://zhile.io/2018/08/18/jetbrai ...
- 转:JAVA线程池ThreadPoolExecutor与阻塞队列BlockingQueue
从Java5开始,Java提供了自己的线程池.每次只执行指定数量的线程,java.util.concurrent.ThreadPoolExecutor 就是这样的线程池.以下是我的学习过程. 首先是构 ...
- 使用SeaJS实现模块化JavaScript开发【转】
前言 SeaJS是一个遵循CommonJS规范的JavaScript模块加载框架,可以实现JavaScript的模块化开发及加载机制.与jQuery等JavaScript框架不同,SeaJS不会扩展封 ...
- 白话SpringCloud | 第十章:路由网关(Zuul)进阶:过滤器、异常处理
前言 简单介绍了关于Zuul的一些简单使用以及一些路由规则的简单说明.而对于一个统一网关而言,需要处理各种各类的请求,对不同的url进行拦截,或者对调用服务的异常进行二次处理等等.今天,我们就来了解下 ...
- 什么是图像 -- opencv基础
opencv基础篇--到底什么是图像 什么是图像?英语中有两个单词来形容图像,一个是picture,一个是image.这两者虽然是形容同一个东西,但却又有着区别.picture代表实而有物的真实图像: ...
- Django分页解析
分页 django中实现管理数据分页的类位于 django.core.paginator.py中 Paginator类 对列表数据进行分页处理 对象 Paginator(Post.objects.al ...
- -ms-zoom property
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...