本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。

本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!

题目链接:SPOJ8222

正解:后缀自动机

解题报告:

  我好菜啊,现在才学SAM…

  大概的构造就是在线的增量法,看看代码还是挺好懂的,简洁明了…

  有几个基础性质:自动机上的每个点上代表的字符串集合的右端点相同,$right$集合相同。

  而且两个点的$right$集合要么是包含关系,要么就不相交。

  每个节点的控制范围就是$(maxs_{fa},maxs_u]$,这样就不用维护自身的$mins$了,比较方便。

  儿子节点的$right$集合是父状态的子集,在$parent$树上往上走相当于砍掉了一个前缀。

  根据前面的性质,$|right_{fa}|= \sum {|right_u|}$,那么这道题就能做了…

  自下往上搞出每个点的$right$集合大小,用集合大小去更新$maxs$的$ans$就可以了。

  要用基数排序。

//It is made by ljh2000
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
#include <bitset>
using namespace std;
typedef long long LL;
typedef long double LB;
typedef complex<double> C;
const double pi = acos(-1);
const int MAXN = 500011;
int n,last,S,cnt,a[MAXN],ch[MAXN][26],len[MAXN],fa[MAXN],ans[MAXN];
int tong[MAXN],sa[MAXN],R[MAXN],f[MAXN];
char s[MAXN]; inline int getint(){
int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
} inline void add(int x){
int c=a[x]; int p=last,np=++cnt; last=np;
len[np]=x; for(;p && !ch[p][c];p=fa[p]) ch[p][c]=np;
if(!p) fa[np]=S;
else {
int q=ch[p][c];
if(len[q]==len[p]+1) fa[np]=q;
else {//不满足SAM上每个节点关于right集合的性质,拆成两个来分别接收
int nq=++cnt; len[nq]=len[p]+1;
memcpy(ch[nq],ch[q],sizeof(ch[q]));
fa[nq]=fa[q];
fa[q]=fa[np]=nq;
for(;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;//把对应的走到q的边全改成走到nq
}
}
} inline void work(){
scanf("%s",s+1); n=strlen(s+1);
S=last=++cnt; for(int i=1;i<=n;i++) a[i]=s[i]-'a';
for(int i=1;i<=n;i++) add(i);
for(int i=1,p=S;i<=n;i++) p=ch[p][a[i]],R[p]++;//主链上的right集合初值设为1
for(int i=1;i<=cnt;i++) tong[ len[i] ]++;
for(int i=1;i<=n;i++) tong[i]+=tong[i-1];
for(int i=1;i<=cnt;i++) sa[ tong[ len[i] ] -- ]=i;
for(int i=cnt;i>=1;i--) R[ fa[ sa[i] ] ]+=R[ sa[i] ];
for(int i=1;i<=cnt;i++) f[ len[i] ]=max(f[ len[i] ],R[i]);//统计每种长度子串的最大出现次数,就是用right集合去更新maxs
for(int i=n;i>=1;i--) f[i]=max(f[i+1],f[i]);
for(int i=1;i<=n;i++) printf("%d\n",f[i]);
} int main()
{
work();
return 0;
}
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。

  

SPOJ8222 NSUBSTR - Substrings的更多相关文章

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

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

  2. SPOJ8222 NSUBSTR - Substrings 后缀自动机_动态规划

    讲起来不是特别好讲.总之,如果 $dp[i+1]>=dp[i]$,故$dp[i]=max(dp[i],dp[i+1])$ Code: #include <cstdio> #inclu ...

  3. SPOJ8222/NSUBSTR:Substrings——题解

    https://www.luogu.org/problemnew/show/SP8222#sub http://www.spoj.com/problems/NSUBSTR/ 翻译来自洛谷. 你得到一个 ...

  4. 【spoj8222】 Substrings

    http://www.spoj.com/problems/NSUBSTR/ (题目链接) 题意 给出一个字符串S,令${F(x)}$表示S的所有长度为x的子串出现次数的最大值.求${F(1)..... ...

  5. 【spoj8222】Substrings

    #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #i ...

  6. ●SPOJ 8222 NSUBSTR–Substrings

    题链: http://www.spoj.com/problems/NSUBSTR/题解: 后缀自动机. 不难发现,对于自动机里面的一个状态s, 如果其允许的最大长度为maxs[s],其right集合的 ...

  7. ●SPOJ 8222 NSUBSTR - Substrings(后缀数组)

    题链: http://www.spoj.com/problems/NSUBSTR/ 题解: 同届红太阳 --WSY给出的后缀数组解法!!! 首先用倍增算法求出 sa[i],rak[i],hei[i]然 ...

  8. ●SPOJ 8222 NSUBSTR–Substrings(后缀自动机)

    题链: http://www.spoj.com/problems/NSUBSTR/ 题解: 后缀自动机的水好深啊!懂不了相关证明,带着结论把这个题做了.看来这滩深水要以后再来了. 本题要用到一个叫 R ...

  9. SPOJ 8222 NSUBSTR - Substrings

    http://www.spoj.com/problems/NSUBSTR/ 题意: F(x)定义为字符串S中所有长度为x的子串重复出现的最大次数 输出F[1]~F[len(S)] 用字符串S构建后缀自 ...

随机推荐

  1. MySql库、表权限管理

    #授权表user #该表放行的权限,针对:所有数据,所有库下所有表,以及表下的所有字段db #该表放行的权限,针对:某一数据库,该数据库下的所有表,以及表下的所有字段tables_priv #该表放行 ...

  2. docker-compose部署zk集群、kafka集群以及kafka-manager,及其遇到的问题和解决

    zk集群docker-compose.yml 1.新建网络 docker network create --driver bridge --subnet --gateway 172.23.0.1 zo ...

  3. PAT 1151 LCA in a Binary Tree[难][二叉树]

    1151 LCA in a Binary Tree (30 分) The lowest common ancestor (LCA) of two nodes U and V in a tree is ...

  4. genymotion——VT-x is not available (VERR_VMX_NO_VMX) 的解决方案

    进入虚拟机设置页面

  5. bit,byte,char,位,字节,字符 的区别

    bit,byte,char,位,字节,字符 的区别 原创文章,未经作者允许,禁止转载!!!

  6. Creating an generated Earth AVI with C++

    Creating an generated Earth AVI with C++        EarthGenerator.cpp /*    EarthGenerator.cpp An examp ...

  7. Activity 间 bundle 传递参数

    activity_main.xml <TableLayout xmlns:android="http://schemas.android.com/apk/res/android&quo ...

  8. Extjs3.4 grid中添加一列复选框

    var sm = new Ext.grid.CheckboxSelectionModel(); var cm = new Ext.grid.ColumnModel( [ sm, new Ext.gri ...

  9. 09 grep、正则表达式和sed

    作业一:整理正则表达式博客 ^ 行首$ 行尾. 除了换行符以外的任意单个字符* 前导字符的零个或多个.* 所有字符[] 字符组内的任一字符[^] 对字符组内的每个字符取反(不匹配字符组内的每个字符)^ ...

  10. Linux系统——awk命令

    awk命令不仅仅是Linux系统的命令,也是一种编程语言,用来处理数据和生成报告(Exel),处理的数据可以是一个或多个文件(标准输入和管道获取标准输入).可在命令行上编辑操作,也可以写成awk程序运 ...