BZOJ2251 [2010Beijing Wc]外星联络 后缀数组 + Height数组
Description
小 P 在看过电影《超时空接触》(Contact)之后被深深的打动,决心致力于寻
找外星人的事业。于是,他每天晚上都爬在屋顶上试图用自己的收音机收听外星
人发来的信息。虽然他收听到的仅仅是一些噪声,但是他还是按照这些噪声的高
低电平将接收到的信号改写为由 0 和 1 构成的串, 并坚信外星人的信息就隐藏在
其中。他认为,外星人发来的信息一定会在他接受到的 01 串中重复出现,所以
他希望找到他接受到的 01 串中所有重复出现次数大于 1 的子串。但是他收到的
信号串实在是太长了,于是,他希望你能编一个程序来帮助他。
Input
输入文件的第一行是一个整数N ,代表小 P 接收到的信号串的长度。
输入文件第二行包含一个长度为N 的 01 串,代表小 P 接收到的信号串。
Output
输出文件的每一行包含一个出现次数大于1 的子串所出现的次数。输出的顺
序按对应的子串的字典序排列。
往往,为了不算重,我们让一个排名为 $i$ 的后缀只计算 $height_{i}+1$ 到 $n$ 的贡献.
这道题也是如此.
依次枚举每一个后缀,并枚举长度为 $[height_{i}+1,maxlen]$ 的子串.
尽可能向下延申,如果出现次数 $>$ 2 则直接输出答案即可.
#include <bits/stdc++.h>
#define setIO(s) freopen(s".in", "r", stdin)
#define maxn 3000000
using namespace std;
char str[maxn];
namespace SA
{
int n, m;
int rk[maxn], sa[maxn], tp[maxn], tax[maxn], height[maxn];
void qsort()
{
for(int i = 0; i <= m ; ++i) tax[i] = 0;
for(int i = 1; i <= n ; ++i) ++tax[rk[i]];
for(int i = 1; i <= m ; ++i) tax[i] += tax[i - 1];
for(int i = n; i >= 1 ; --i) sa[tax[rk[tp[i]]] -- ] = tp[i];
}
void build()
{
for(int i = 1; i <= n ; ++i) rk[i] = str[i], tp[i] = i;
qsort();
for(int k = 1; k <= n ; k <<= 1)
{
int p = 0;
for(int i = n - k + 1; i <= n ; ++i) tp[++p] = i;
for(int i = 1; i <= n ; ++i) if(sa[i] > k) tp[++p] = sa[i] - k;
qsort(), swap(rk, tp), rk[sa[1]] = p = 1;
for(int i = 2; i <= n ; ++i)
{
rk[sa[i]] = (tp[sa[i - 1]] == tp[sa[i]] && tp[sa[i - 1] + k] == tp[sa[i] + k]) ? p : ++p;
}
if(p == n) break;
m = p;
}
}
void get_height()
{
int k = 0;
for(int i = 1; i <= n ; ++i) rk[sa[i]] = i;
for(int i = 1; i <= n ; ++i)
{
if(k) --k;
int j = sa[rk[i] - 1];
while(str[j + k] == str[i + k]) ++k;
height[rk[i]] = k;
}
}
};
int main()
{
//setIO("input");
scanf("%d%s",&SA::n,str + 1), SA::m = 122, SA::build(), SA::get_height();
for(int i = 1; i < SA::n ; ++i)
{
for(int j = SA::height[i] + 1; SA::sa[i] + j - 1 <= SA::n ; ++j)
{
int k = i;
while(SA::height[k + 1] >= j) ++k;
if(k > i) printf("%d\n", k - i + 1);
}
}
return 0;
}
BZOJ2251 [2010Beijing Wc]外星联络 后缀数组 + Height数组的更多相关文章
- [bzoj2251][2010Beijing Wc]外星联络——后缀数组+暴力求解
Brief Description 找到 01 串中所有重复出现次数大于 1 的子串.并按字典序输出他们的出现次数. Algorithm Design 求出后缀数组之后,枚举每一个后缀,对于每个后缀从 ...
- 【BZOJ2251】[2010Beijing Wc]外星联络 后缀数组
[BZOJ2251][2010Beijing Wc]外星联络 Description 小 P 在看过电影<超时空接触>(Contact)之后被深深的打动,决心致力于寻找外星人的事业.于是, ...
- bzoj 2251: [2010Beijing Wc]外星联络 后缀数组
2251: [2010Beijing Wc]外星联络 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 424 Solved: 232[Submit][ ...
- [bzoj2251][2010BeiJing Wc]外星联络_后缀数组
外星联络 bzoj-2251 2010-BeiJing Wc 题目大意:题目链接. 注释:略. 想法: 这咋做啊????一看数据范围才$3\cdot 10^3$. 建立后缀数组. 所以我们将所有后缀排 ...
- bzoj2251 [2010Beijing Wc]外星联络
因为n很小,所以对于串s的每一个后缀,都把其加入字典树中,并且经过一个字典树节点,该节点权值就+1. 输出时因为要字典序最小,所以字典树先走0分叉,再走1分叉,如果节点权值大于等于2就输出 代码 #i ...
- BZOJ_2251_[2010Beijing Wc]外星联络_后缀数组
BZOJ_2251_[2010Beijing Wc]外星联络_后缀数组 Description 小 P 在看过电影<超时空接触>(Contact)之后被深深的打动,决心致力于寻 找外星人的 ...
- 2251: [2010Beijing Wc]外星联络
2251: [2010Beijing Wc]外星联络 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 801 Solved: 481[Submit][ ...
- BZOJ 2251: [2010Beijing Wc]外星联络
2251: [2010Beijing Wc]外星联络 Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 795 Solved: 477[Submit][ ...
- 2251. [2010Beijing Wc]外星联络【后缀数组】
Description 小 P 在看过电影<超时空接触>(Contact)之后被深深的打动,决心致力于寻 找外星人的事业.于是,他每天晚上都爬在屋顶上试图用自己的收音机收听外星 人发来的信 ...
随机推荐
- Void 参数
在C程序中如果在声明函数的时候如果没有任何参数那么需要将参数定义为void以此来限定此函数不可传递任何参数,如果不进行限定让参数表默认为空其意义是可以传递任何参数,这个问题的由来实际上是由于要兼容早期 ...
- 优秀软件project师必备的7大特性
不是每个程序猿都能成为优秀的软件project师. 在过去的6年时间里,我在Ooyala.Quora和now Quip这3个创业公司面试过许很多多挺有发展潜力的"种子选手".他们都 ...
- [CSS3] Make a horizontal-scrolling Menu
Our web app on desktop may use elements that use great deal of width. On many occasions we can’t sim ...
- 571B. Minimization(Codeforces Round #317)
B. Minimization time limit per test 2 seconds memory limit per test 256 megabytes input standard inp ...
- C++ 句柄类的原理以及设计
句柄类存在的意义是为了弥补将派生类对象赋给基类对象时发生的切片效应.比如以下的程序: multimap<Base> basket; Base base; Derived derive; b ...
- mysql-5.5 for linux源码安装
mysql-5.5 for linux源码安装 1.使用Yum安装依赖软件包 # yum install -y gcc gcc-c++ gcc-g77 autoconf automake bison ...
- python lmdb demo 这接口和BDB一样恶心啊!
import lmdb lmdb_img_name = "test.lmdb" env = lmdb.open(lmdb_img_name, map_size=1e6) with ...
- 最小割板子题——[USACO5.4]奶牛的电信
今天邱神给我们讲了图论,还讲了一下网络流算法.自己找了一个洛谷板子题. 题目描述 农夫约翰的奶牛们喜欢通过电邮保持联系,于是她们建立了一个奶牛电脑网络,以便互相交流.这些机器用如下的方式发送电邮:如果 ...
- pandas把多个sheet读进一个DataFrame
#!/usr/bin/python import pandas as pd import collections df_dict = pd.read_excel('c:\data\machines.x ...
- openStack 主机流量计运行状态 随笔记录
root@ruiy-controller:~# ifconfigeth0 Link encap:Ethernet HWaddr 0c:c4:7a:0d:97:2c ine ...