题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3998

题意概述:对于一个给定长度为N的字符串,求它的第K小子串是什么,T为0则表示不同位置的相同子串算作一个。T=1则表示不同位置的相同子串算作多个。N<=500000,K<=10^9.

应该是有三种做法的(当然后缀树我还没有看),于是就把这个东西当成后缀自动机的板子了(因为它很裸啊!!!)。

可以注意到当T=0的时候每走动一步的贡献是1,而T=1的时候每走动一步之后的贡献都是走到的这个状态的right集合大小。同时从同一个状态走出去得到的字符串又有很多种,不难想到可以搞个dp表示从这个点出发(包括这个点自己)可以找到的不同子串数(根据T而定)。

有个技巧是根据MAX把所有的状态排序之后就可以欢快地topo图上转移dp啦!(SAM的转移图是一张DAG)

所以还是好好理解一下SAM上的状态和转移吧,,,。。。。

 #include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<cctype>
using namespace std;
const int MAXN=; int T,K,n; char S[MAXN]; struct SAM{
static const int maxn=;
static const int sigma_sz=;
int pa[maxn<<],to[maxn<<][sigma_sz],mx[maxn<<],val[maxn<<]; long long cnt[maxn<<];
int sz,last,c[maxn],a[maxn<<];
SAM(){ last=sz=; memset(to[],,sizeof(to[])); }
int newnode(){
memset(to[++sz],,sizeof(to[sz]));
pa[sz]=mx[sz]=val[sz]=cnt[sz]=;
return sz;
}
void extend(int w){
int p=last,np=newnode(); last=np;
pa[np]=mx[np]=cnt[np]=;
mx[np]=mx[p]+,val[np]=;
while(p&&!to[p][w]) to[p][w]=np,p=pa[p];
if(!p) pa[np]=;
else{
int q=to[p][w];
if(mx[p]+==mx[q]) pa[np]=q;
else{
int nq=newnode(); mx[nq]=mx[p]+;
memcpy(to[nq],to[q],sizeof(to[nq]));
pa[nq]=pa[q];
pa[np]=pa[q]=nq;
while(p&&to[p][w]==q) to[p][w]=nq,p=pa[p];
}
}
}
void ready(){
for(int i=;i<=sz;i++) c[mx[i]]++;
for(int i=;i<=n;i++) c[i]+=c[i-];
for(int i=sz;i>;i--) a[c[mx[i]]--]=i;
for(int i=sz;i>;i--){
if(T) val[pa[a[i]]]+=val[a[i]];
else val[a[i]]=;
}
val[]=;
for(int i=sz;i>;i--){
cnt[a[i]]=val[a[i]];
for(int j=;j<;j++) cnt[a[i]]+=cnt[to[a[i]][j]];
}
}
void dfs(int i){
if(val[i]>=K) return;
K-=val[i];
for(int j=;j<;j++) if(to[i][j]){
if(cnt[to[i][j]]>=K){
putchar(j+'a'); dfs(to[i][j]);
return;
}
K-=cnt[to[i][j]];
}
}
}sam; void data_in()
{
gets(S);
scanf("%d%d",&T,&K);
}
void work()
{
n=strlen(S);
for(int i=;i<n;i++) sam.extend(S[i]-'a');
sam.ready();
if(sam.cnt[]<K) printf("%d\n",-);
else sam.dfs();
}
int main()
{
data_in();
work();
return ;
}

BZOJ 3998 TJOI2015 弦论 后缀自动机+DAG上的dp的更多相关文章

  1. BZOJ 3998: [TJOI2015]弦论 [后缀自动机 DP]

    3998: [TJOI2015]弦论 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 2152  Solved: 716[Submit][Status] ...

  2. BZOJ 3998: [TJOI2015]弦论 后缀自动机 后缀自动机求第k小子串

    http://www.lydsy.com/JudgeOnline/problem.php?id=3998 后缀自动机应用的一个模板?需要对len进行一个排序之后再统计每个出现的数量,维护的是以该字符串 ...

  3. bzoj 3998 [TJOI2015]弦论——后缀自动机

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3998 相同子串算多个的话,先求好 right ,然后求一个 sm 表示走到这个点之后有几种走 ...

  4. BZOJ 3998 [TJOI2015]弦论 ——后缀自动机

    直接构建后缀自动机. 然后. 然后只需要再后缀自动机的go树上类似二分的方法进行查找即可,实际上是“26分”. 然后遇到了处理right集合的问题,然后觉得在go和parent树上上传都是可以的,毕竟 ...

  5. BZOJ.3998.[TJOI2015]弦论(后缀自动机)

    题目链接 \(Description\) 给定字符串S,求其第K小子串.(若T=0,不同位置的相同子串算1个:否则算作多个) \(Solution\) 建SAM,处理出对于每个节点,它和它的所有后继包 ...

  6. BZOJ 3998: [TJOI2015]弦论(后缀自动机)

    传送门 解题思路 \(T=0\)时就和SP7258一样,\(T=1\)时其实也差不多,只不过要把每个点原来是\(1\)的权值改为\(Right\)集合的大小. 代码 #include<iostr ...

  7. ●BZOJ 3998 [TJOI2015]弦论

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3998题解: 后缀自动机. 当T=0时, 由于在后缀自动机上沿着trans转移,每个串都是互不 ...

  8. 【BZOJ3998】[TJOI2015]弦论 后缀自动机

    [BZOJ3998][TJOI2015]弦论 Description 对于一个给定长度为N的字符串,求它的第K小子串是什么. Input 第一行是一个仅由小写英文字母构成的字符串S 第二行为两个整数T ...

  9. [bzoj3998][TJOI2015]弦论-后缀自动机

    Brief Description 给定一个字符串, 您需要求出他的严格k小子串或非严格k小子串. Algorithm Design 考察使用后缀自动机. 首先原串建SAM, 然后如果考察每个状态代表 ...

随机推荐

  1. 菜鸟笔记 -- Chapter 6.2 类的构成

    在前面我们讲过高级开发语言大多由7种语法构成,但这是一个很空泛的概述,下,面我们仅就针对Java程序来说一下构成一个Java程序的几大部分,其中类是最小的基本元素.类是封装对象属性和行为的载体,而在J ...

  2. IOS异步获取数据并刷新界面dispatch_async的使用方法

    在ios的开发和学习中多线程编程是必须会遇到并用到的.在java中以及Android开发中,大量的后台运行,异步消息队列,基本都是运用了多线程来实现. 同样在,在ios移动开发和Android基本是很 ...

  3. c#实现的HTTP服务端

    这次在整理一个服务组件的时候,需要涉及到HTTP的请求,HTTP是应用层,建立在TCP之上的.因此,可以用TCP服务端接收HTTP请求,只需要解析请求内容.HTPP有固定的格式,大家可以直接搜索.网上 ...

  4. 对象API

    遍历对象里的每个元素 var obj ={ a:32, b:12, c :342 } for (const key of obj){ if(obj.hasOwnProperty(key)){ cons ...

  5. 小A点菜

    题目背景 uim神犇拿到了uoi的ra(镭牌)后,立刻拉着基友小A到了一家--餐馆,很低端的那种. uim指着墙上的价目表(太低级了没有菜单),说:"随便点". 题目描述 不过ui ...

  6. 快速玩转linux(3)

    Linux常用命令 软件操作命令 执行操作 命令 软件包管理器 yum 安装软件 yum install xxx 卸载软件 yum remove xxx 搜索软件 yum search xxx 清除缓 ...

  7. 【转载】在C#中主线程和子线程如何实现互相传递数据

    引用:https://blog.csdn.net/shuaihj/article/details/41316731 一.不带参数创建Thread using System; using System. ...

  8. ReentrantLock详解

    ReentrantLock概述 ReentrantLock是Lock接口的实现类,可以手动的对某一段进行加锁.ReentrantLock可重入锁,具有可重入性,并且支持可中断锁.其内部对锁的控制有两种 ...

  9. Hive优化之谓词下推

    Hive优化之谓词下推 解释 Hive谓词下推(Predicate pushdown) 关系型数据库借鉴而来,关系型数据中谓词下推到外部数据库用以减少数据传输 基本思想:尽可能早的处理表达式 属于逻辑 ...

  10. Python学习:If 语句与 While 语句

    If 语句 用以检查条件:如果条件为真(True),将运行这一块的语句(称作 if-block 或 if 块)    则将运行另一块语句(称作 else-block 或 else 块),其中 else ...