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 ...
随机推荐
- python 实现剪刀石头布(三局两胜)
# -*- coding:utf-8 -*- import random # best of three def finger_guess(): rule = {1:'rock', 2:'paper' ...
- 【JAVA】apachehttpclient设置http1.0短链接
HttpPost httpPost = new HttpPost(url);httpPost.setEntity(new StringEntity(text, ContentType.create(& ...
- tcl之list操作-lappend/lsearch/lsort/concat/split/join/
- String使用方法详解
标准c++中string类函数介绍 注意不是CString 之所以抛弃char*的字符串而选用C++标准程序库中的string类,是因为他和前者比较起来,不必 担心内存是否足够.字符串长度等等,而且作 ...
- BFS:HDU-1242-Rescue(带守卫的迷宫问题)(优先队列)
解题心得: 1.读清楚题意,本题的题意是有多个'r'(起点),多个r多个bfs比较最短的时间即可,但是hdoj的数据比较水,直接一个起点就行了,迷宫里有多个守卫,如果在路途中遇到守卫会多花费一个时间点 ...
- 其它- in-place
in-place 刷编程题的时候,经常遇到题目要求do in-place.所以就上网搜了相关概念,简单总结一下. in-place操作,意思是所有的操作都是”就地“操作,不允许进行移动,或者称作 ...
- laravel5.2总结--服务容器(依赖注入,控制反转)
1.依赖 我们定义两个类:class Supperman 和 class Power,现在我们要使用Supperman ,而Supperman 依赖了Power class Supperman { p ...
- FSMO
FSMO 编辑 本词条缺少信息栏.名片图,补充相关内容使词条更完整,还能快速升级,赶紧来编辑吧! FSMO是Flexible single master operation的缩写,意思就是灵活单主机操 ...
- IOS开发学习笔记005-数组
数组 数组故名思议就是一组数据的集合. int a[10];//可以存储10个整数 char c[8];//可以存储8个字符‘ 一般格式:数组类型 数组名[元素个数]: 数组元素的访问:下标,a[2 ...
- Appium Windows服务端GUI详解
Appium Windows服务端GUI各项的解释,从官方扒过来的,界面图标和最新版本有点不太一样,其他还是比较简单易懂的 原文https://github.com/appium/appium-dot ...