poj 3693 Maximum repetition substring 重复次数最多的连续子串
题目链接
题意
对于任意的字符串,定义它的 重复次数 为:它最多可被划分成的完全相同的子串个数。例如:ababab
的重复次数为3,ababa
的重复次数为1.
现给定一字符串,求它的一个子串,其重复次数取到最大值,且字典序取到最小值。
思路
参考 hzwer.
首先,重复次数显然至少为\(1\),所以下面只考虑重复次数\(\geq 2\)的情况。
首先枚举子串长度\(L\)。对于长度为\(L\)的子串,其重复次数至少为\(2\),意味着它的其中两个重复部分必为\(s[0],s[L],s[2L],...\)中相邻的两个。
所以只需枚举\(i\),看\(s[i*L]\)和\(s[(i+1)*L]\)向左\((le)\)向右\((ri)\)最多能匹配到多远,记总长度为\(k\),则重复次数\(=k/L+1\).
至于开头位置,则落在\(i*L-le\)和\(i*L-le+(le+ri)\%k\)之间,取字典序最小者。
对于所有重复次数相同的子串,要使其字典序取到最小值,即取\(rank\)值最小者。
上述步骤中,看最长匹配多远(即看最长公共前缀)与取rank最小值均借助\(ST\)表实现。
Code
#include <stdio.h>
#include <iostream>
#include <string.h>
#define maxn 100010
using namespace std;
typedef long long LL;
int wa[maxn], wb[maxn], wv[maxn], wt[maxn], r1[maxn], r2[maxn],
h1[maxn], h2[maxn], rk1[maxn], rk2[maxn], sa[maxn], n, kas, bin[20], Log[maxn];
struct rmqNode { int val, p; } mn[maxn][20];
struct node { int val, pos, len; };
int mn1[maxn][20], mn2[maxn][20];
char s[maxn];
void rmqInit(int n) {
Log[0] = -1; bin[0] = 1;
for (int i = 1; i < 20; ++i) bin[i] = bin[i-1] << 1;
for (int i = 1; i <= n; ++i) Log[i] = Log[i>>1] + 1;
}
bool cmp(int* r, int a, int b, int l) { return r[a] == r[b] && r[a+l] == r[b+l]; }
void init(int* r, int* rk, int* h, int n, int m) {
int* x=wa, *y=wb, *t, i, j, p;
for (i = 0; i < m; ++i) wt[i] = 0;
for (i = 0; i < n; ++i) ++wt[x[i] = r[i]];
for (i = 1; i < m; ++i) wt[i] += wt[i - 1];
for (i = n-1; i >= 0; --i) sa[--wt[x[i]]] = i;
for (j = 1, p = 1; p < n; j <<= 1, m = p) {
for (p = 0, i = n-j; i < n; ++i) y[p++] = i;
for (i = 0; i < n; ++i) if (sa[i] >= j) y[p++] = sa[i] - j;
for (i = 0; i < n; ++i) wv[i] = x[y[i]];
for (i = 0; i < m; ++i) wt[i] = 0;
for (i = 0; i < n; ++i) ++wt[wv[i]];
for (i = 1; i < m; ++i) wt[i] += wt[i - 1];
for (i = n-1; i >= 0; --i) sa[--wt[wv[i]]] = y[i];
t = x, x = y, y = t, x[sa[0]] = 0;
for (p = 1, i = 1; i < n; ++i) x[sa[i]] = cmp(y, sa[i], sa[i-1], j) ? p - 1 : p++;
}
for (i = 0; i < n; ++i) rk[sa[i]] = i;
int k = 0;
for (i = 0; i < n - 1; h[rk[i++]] = k) {
for (k = k ? --k : 0, j = sa[rk[i] - 1]; r[i+k] == r[j+k]; ++k);
}
}
void rmq1(int* a, int (*mn)[20], int n) {
for (int i = 1; i <= n; ++i) mn[i][0] = a[i];
for (int j = 1; bin[j] <= n; ++j) {
for (int i = 1; i+bin[j-1]-1 <= n; ++i) {
mn[i][j] = min(mn[i][j-1], mn[i+bin[j-1]][j-1]);
}
}
}
void rmq2(int* a, rmqNode (*mn)[20], int n) {
for (int i = 1; i <= n; ++i) mn[i][0] = {a[i-1], i};
for (int j = 1; bin[j] <= n; ++j) {
for (int i = 1; i+bin[j-1]-1 <= n; ++i) {
if (mn[i][j-1].val <= mn[i+bin[j-1]][j-1].val) mn[i][j] = mn[i][j-1];
else mn[i][j] = mn[i+bin[j-1]][j-1];
}
}
}
int query1(int (*mn)[20], int l, int r) {
int k = Log[r-l+1];
return min(mn[l][k], mn[r-bin[k]+1][k]);
}
int query2(rmqNode (*mn)[20], int l, int r) {
int k = Log[r-l+1];
return mn[l][k].val < mn[r-bin[k]+1][k].val ? mn[l][k].p : mn[r-bin[k]+1][k].p;
}
int match(int* rk, int (*mn)[20], int p, int len) {
int rk1 = rk[p], rk2 = rk[p+len];
if (rk1 > rk2) swap(rk1, rk2);
return query1(mn, rk1+1, rk2);
}
void work() {
int tot1=0, tot2=0, m=0, len=strlen(s);
for (int i = 0; i < len; ++i) m = max(r1[tot1++] = s[i], m); r1[tot1++] = 0;
for (int i = len-1; i >= 0; --i) r2[tot2++] = s[i]; r2[tot2++] = 0;
rmqInit(len);
init(r1, rk1, h1, tot1, ++m);
init(r2, rk2, h2, tot2, m);
rmq1(h1, mn1, len);
rmq1(h2, mn2, len);
rmq2(rk1, mn, len);
node ans = {1,0,0};
for (int l = 1; l <= len; ++l) {
int lim = len / l, upp;
if (ans.val > lim+1) break;
if (len % l) upp = lim; else upp = lim-1;
for (int i = 0; i < upp; ++i) {
int ri = match(rk1, mn1, i*l, l),
le = i ? match(rk2, mn2, len-(i+1)*l, l) : 0,
k = le + ri;
int cnt = k/l + 1;
if (cnt >= ans.val) {
int l1 = i*l-le, l2 = l1+k%l,
p = query2(mn, l1+1, l2+1)-1;
if (cnt > ans.val || rk1[p]<rk1[ans.pos]) ans = {cnt, p, cnt*l};
}
}
}
printf("Case %d: ", ++kas);
if (ans.val == 1) {
char ch='z'+1; for (int i = 0; s[i]; ++i) ch = min(ch, s[i]);
putchar(ch); puts("");
}
else {
s[ans.pos+ans.len] = '\0';
puts(s+ans.pos);
}
}
int main() {
while (scanf("%s", s) != EOF && s[0]!='#') work();
return 0;
}
poj 3693 Maximum repetition substring 重复次数最多的连续子串的更多相关文章
- 【POJ 3693】Maximum repetition substring 重复次数最多的连续重复子串
后缀数组的论文里的例题,论文里的题解并没有看懂,,, 求一个重复次数最多的连续重复子串,又因为要找最靠前的,所以扫的时候记录最大的重复次数为$ans$,扫完后再后从头暴力扫到尾找重复次数为$ans$的 ...
- POJ - 3693 Maximum repetition substring(重复次数最多的连续重复子串)
传送门:POJ - 3693 题意:给你一个字符串,求重复次数最多的连续重复子串,如果有一样的,取字典序小的字符串. 题解: 比较容易理解的部分就是枚举长度为L,然后看长度为L的字符串最多连续出现 ...
- POJ 3693 Maximum repetition substring(连续重复子串)
http://poj.org/problem?id=3693 题意:给定一个字符串,求重复次数最多的连续重复子串. 思路: 这道题确实是搞了很久,首先枚举连续子串的长度L,那么子串肯定包含了r[k], ...
- POJ 3693 Maximum repetition substring(后缀数组)
Description The repetition number of a string is defined as the maximum number R such that the strin ...
- POJ 3693 Maximum repetition substring(后缀数组+ST表)
[题目链接] poj.org/problem?id=3693 [题目大意] 求一个串重复次数最多的连续重复子串并输出,要求字典序最小. [题解] 考虑错位匹配,设重复部分长度为l,记s[i]和s[i+ ...
- 后缀数组 POJ 3693 Maximum repetition substring
题目链接 题意:给定一个字符串,求重复次数最多的连续重复子串. 分析:(论文上的分析)先穷举长度 L,然后求长度为 L 的子串最多能连续出现几次.首先连续出现 1 次是肯定可以的,所以这里只考虑至少 ...
- POJ 3693 Maximum repetition substring(最多重复次数的子串)
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 10461 Ac ...
- poj 3693 Maximum repetition substring (后缀数组)
其实是论文题.. 题意:求一个字符串中,能由单位串repeat得到的子串中,单位串重复次数最多的子串.若有多个重复次数相同的,输出字典序最小的那个. 解题思路:其实跟论文差不多,我看了很久没看懂,后来 ...
- POJ 3693 Maximum repetition substring ——后缀数组
重复次数最多的字串,我们可以枚举循环节的长度. 然后正反两次LCP,然后发现如果长度%L有剩余的情况时,答案是在一个区间内的. 所以需要找到区间内最小的rk值. 两个后缀数组,四个ST表,$\Thet ...
随机推荐
- 27.28. VUE学习之--事件修饰符之stop&capture&self&once实例详解
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- 树莓派下ubuntu-mate中ssh服务的安装与开机后自启动
ssh程序分为客户端程序openssh-client和服务端程序openssh-server. 如果需要ssh登陆到别的电脑,需要安装openssh-client,该程序ubuntu是默认安装的.而如 ...
- SQLite3 of python
SQLite3 of python 一.SQLite3 数据库 SQLite3 可使用 sqlite3 模块与 Python 进行集成,一般 python 2.5 以上版本默认自带了sqlite3模块 ...
- Python学习笔记(二):数据类型
一.python中的数据类型 python中的数据类型包括:整型.浮点型.布尔型.字符串类型 整型(int)和浮点型(float) Python中的整型只有int,没有short.long:浮点型fl ...
- STM32启动文件:startup_stm32f10x_hd.s等启动文件的简单描述
在官方的库文件中,分别有如下文件: startup │ │ │ ├─arm │ │ │ │ startup_stm32f10x_cl.s │ │ │ │ startup_stm32f10x_hd.s ...
- Leetcode 96. 不同的二叉搜索树
题目链接 https://leetcode.com/problems/unique-binary-search-trees/description/ 题目描述 给定一个整数 n,求以 1 ... n ...
- git---gui使用
1.登陆的命令: git config –global user.email "1455971532@qq.com" git config –global user.name &q ...
- Less Css 教程
http://www.w3cplus.com/css/less,这个东西太吊了!
- Careercup - Microsoft面试题 - 5799446021406720
2014-05-12 07:17 题目链接 原题: Given below is a tree/trie A B c D e F a<b<e<>>c<>d&l ...
- Python框架之Django学习笔记(六)
模板 上篇博文学习了动态视图,但是,视图中返回文本的方式有点特别. 也就是说,HTML被直接硬编码在 Python 代码之中. def current_datetime(request): now = ...