题目大意:给一个字符串S,令F(x)表示S的所有长度为x的子串中,出现次数的最大值。F(1)..F(Length(S))

建出SAM, 然后求出Right, 求Right可以按拓扑序dp..Right就是某个点到结束状态的路径数, parent树上last的那一条链都是结束状态...然后用Right去更新答案..

spoj卡常数..一开始用DFS就炸了, 改用BFS就A了..

(贴一下丽洁姐的题解: 我们构造S的SAM,那么对于一个节点s,它的长度范围是[Min(s),Max(s)],同时他的出现次数是|Right(s)|。那么我们用|Right(s)|去更新F(Max(s))的值。同时最后从大到小依次用F(i)去更新F(i-1)即可。)

-----------------------------------------------------------------------------

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
 
using namespace std;
 
const int cn = 26;
const int maxn = 1000009;
 
bool vis[maxn];
int Right[maxn], ans[maxn], deg[maxn], n = 0, N;
char s[maxn];
 
struct Node {
Node *ch[cn], *fa;
int len, id;
} pool[maxn], *pt, *root, *last;
 
Node* newNode(int v) {
memset(pt->ch, 0, sizeof pt->ch);
pt->fa = 0;
pt->len = v;
pt->id = n++;
return pt++;
}
 
void SAM_init() {
pt = pool;
root = last = newNode(0);
}
 
void Extend(int c) {
Node *p = last, *np = newNode(p->len + 1);
for(; p && !p->ch[c]; p = p->fa)
p->ch[c] = np;
if(!p)
np->fa = root;
else {
Node* q = p->ch[c];
if(p->len + 1 == q->len)
np->fa = q;
else {
Node* nq = newNode(p->len + 1);
memcpy(nq->ch, q->ch, sizeof q->ch);
nq->fa = q->fa;
q->fa = np->fa = nq;
for(; p && p->ch[c] == q; p = p->fa)
p->ch[c] = nq;
}
}
last = np;
}
 
struct edge {
int to;
edge* next;
} E[maxn], *Pt = E, *head[maxn];
 
void AddEdge(int u, int v) {
deg[Pt->to = v]++; Pt->next = head[u]; head[u] = Pt++;
}
 
void SAM_build() {
scanf("%s", s);
N = strlen(s);
for(int i = 0; i < N; i++)
Extend(s[i] - 'a');
}
 
queue<Node*> q;
queue<int> Q;
 
void ADDEDGE() {
memset(deg, 0, sizeof deg);
memset(vis, 0, sizeof vis);
q.push(root);
vis[root->id] = true;
while(!q.empty()) {
Node* t = q.front(); q.pop();
for(int i = 0; i < cn; i++) if(t->ch[i]) {
AddEdge(t->ch[i]->id, t->id);
if(!vis[t->ch[i]->id]) {
q.push(t->ch[i]);
vis[t->ch[i]->id] = true;
}
}
}
}
 
void getRight() {
memset(Right, 0, sizeof Right);
for(Node* t = last; t; t = t->fa)
Right[t->id] = 1;
Q.push(last->id);
while(!Q.empty()) {
int x = Q.front(); Q.pop();
for(edge* e = head[x]; e; e = e->next) {
Right[e->to] += Right[x];
if(!--deg[e->to])
Q.push(e->to);
}
}
}
 
void update() {
memset(vis, 0, sizeof vis);
q.push(root);
vis[root->id] = true;
while(!q.empty()) {
Node* t = q.front(); q.pop();
ans[t->len] = max(ans[t->len], Right[t->id]);
for(int i = 0; i < cn; i++) if(t->ch[i] && !vis[t->ch[i]->id]) {
q.push(t->ch[i]);
vis[t->ch[i]->id] = true;
}
}
}
 
void solve() {
getRight();
update();
for(int i = N; --i; )
ans[i] = max(ans[i], ans[i + 1]);
for(int i = 1; i <= N; i++)
printf("%d\n", ans[i]);
}
 
int main() {
SAM_init();
SAM_build();
ADDEDGE();
solve();
return 0;
}

-----------------------------------------------------------------------------

SPOJ8222 Substrings( 后缀自动机 + dp )的更多相关文章

  1. SP8222 NSUBSTR - Substrings(后缀自动机+dp)

    传送门 解题思路 首先建出\(sam\),然后把\(siz\)集合通过拓扑排序算出来.对于每个点只更新它的\(maxlen\),然后再从大到小\(dp\)一次就行了.因为\(f[maxlen-1]&g ...

  2. 【bzoj3998】[TJOI2015]弦论 后缀自动机+dp

    题目描述 对于一个给定长度为N的字符串,求它的第K小子串是什么. 输入 第一行是一个仅由小写英文字母构成的字符串S 第二行为两个整数T和K,T为0则表示不同位置的相同子串算作一个.T=1则表示不同位置 ...

  3. 【SPOJ -NSUBSTR】Substrings 【后缀自动机+dp】

    题意 给出一个字符串,要你找出所有长度的子串分别的最多出现次数. 分析 我们建出后缀自动机,然后预处理出每个状态的cnt,cnt[u]指的是u这个状态的right集合大小.我们设f[len]为长度为l ...

  4. Substrings(SPOJ8222) (sam(后缀自动机))

    You are given a string \(S\) which consists of 250000 lowercase latin letters at most. We define \(F ...

  5. SPOJ8222 NSUBSTR - Substrings(后缀自动机)

    You are given a string S which consists of 250000 lowercase latin letters at most. We define F(x) as ...

  6. SPOJ NSUBSTR Substrings 后缀自动机

    人生第一道后缀自动机,总是值得纪念的嘛.. 后缀自动机学了很久很久,先是看CJL的论文,看懂了很多概念,关于right集,关于pre,关于自动机的术语,关于为什么它是线性的结点,线性的连边.许多铺垫的 ...

  7. bzoj 2806: [Ctsc2012]Cheat 后缀自动机DP

    2806: [Ctsc2012]Cheat Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 583  Solved: 330[Submit][Statu ...

  8. fjwc2019 D1T2 原样输出(后缀自动机+dp)

    #179. 「2019冬令营提高组」原样输出 暴力对每个串建后缀自动机,然后暴力枚举每个自动机的子串.可以拿到部分分. 然鹅我们可以把每个后缀自动机连起来. 我们知道,后缀自动机是用最少的点(空间)表 ...

  9. 【CF316G3】Good Substrings 后缀自动机

    [CF316G3]Good Substrings 题意:给出n个限制(p,l,r),我们称一个字符串满足一个限制当且仅当这个字符串在p中的出现次数在[l,r]之间.现在想问你S的所有本质不同的子串中, ...

随机推荐

  1. The error indicates that IIS is in 32 bit mode, while this application is a 64 b it application and thus not compatible.

    I was trying to install a new WSS v3 Sharepoint on a 64 bit Windows 2003 server today but the instal ...

  2. linux内核升级图文攻略

    Linux内核概览 Linux是一个一体化内核(monolithic kernel)系统. 设备驱动程序可以完全访问硬件. Linux内的设备驱动程序可以方便地以模块化(modularize)的形式设 ...

  3. ubuntu 12.04下安装和配置kohana 3.3.3 的方法

    一.先到官网下载3.3.3版本的压缩包到/var/www/1117/目录下(提前建好1117的目录)解压 解压好的文件有(applications\modules\system\build.xml\c ...

  4. javascrip cookie

    首先要明白一下cookie的概念.由于HTTP协议是一种无状态协议,也就是说一旦server和client的数据交换完成后,他们之间的连接就会被断开.再次交换数据的时候就须要再次建立连接.这就意味着s ...

  5. C/C++指针的指针(**p)和指针的引用(*&amp;)使用案例分析

    #include <stdio.h> #include <unistd.h> #include <fcntl.h> #include <stdlib.h> ...

  6. Java 程序中的多线程

    概述 synchronized  关键字,代表这个方法加锁,相当于不管哪一个线程(例如线程A),运行到这个方法时,都要检查有没有其它线程B(或者C. D等)正在用这个方法,有的话要等正在使用synch ...

  7. 有感PMI Exam Dev Workshop

    有幸參加了PMI协会在上海举办的PMI Exam Development Workshop活动.这是PMI协会第二次在中国举办此活动,上一次是2009年北京. 我第一次參加,感觉收获非常多. 我们知道 ...

  8. 在SSIS包中使用 Checkpoint从失败处重新启动包

    使用SSIS做ETL的过程中会遇到各种各样的错误,对于一些大数据量的Job失败以后我们不希望重新运行,因为重新运行的时间开销是非常大的,我们只希望从失败的部分开始运行,这样可以省去很多的时间. SSI ...

  9. XMLHttpRequest对象的使用

    1.首先要创建XMLHttpRequest对象,这个对象是前台与后台进行异步的重要对象,现在的浏览器有很多种,创建 XMLHttpRequest 的方法也不相同,所以为了兼容各种浏览器,在创建XMLH ...

  10. iOS 饼状图

    首先先看一下效果: 一.创建饼状图对象   创建饼状图对象用到类是PieChartView.h, 代码如下: self.pieChartView = [[PieChartView alloc] ini ...