HDU 2459 Maximum repetition substring
题目:Maximum repetition substring
链接:http://acm.hdu.edu.cn/showproblem.php?pid=2459
题意:给你一个字符串,求连续重复出现次数最多的子串(不重叠),如果有多个,输出字典序最小的那个。别如aacdabcdab,a连续着重复出现2次,cdab连续着出现2次,但aa的字典序小,所以输出aa。
思路:
枚举+后缀树组+rmq
枚举重复串的长度,后缀数组预处理出height数组,rmq预处理height,实现O(1)查询。
给几组测试数据:
1. aabababa,答案为ababab,不是bababa
2. atbctbctb,答案是bctbct,不是ctbctb
3. aabcabcabcab,答案是abcabcabc,不是cabcabcab
#include<stdio.h>
#include<string.h>
#include<math.h>
#define N 200020
char s1[];
int ws[N],wv[N];
int sa[N],r[N],wx[N],wy[N];
int height[N];
bool cmp(int *r,int a,int b,int l)
{
return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(int *r,int n,int m)
{
//注意,这里的n必须比原始数组大小大1
int *x=wx,*y=wy;
for(int i=;i<m;i++) ws[i]=;
for(int i=;i<n;i++) ws[x[i]=r[i]]++;
for(int i=;i<m;i++) ws[i]+=ws[i-];
for(int i=n-;i>=;i--) sa[--ws[x[i]]]=i;
//这里的x[i] 表示下标i的第一关键字排名
int i,j,p,*t;
for(j=,p=;p<n;j*=,m=p)
{
for(p=,i=n-j;i<n;i++) y[p++]=i;
for(i=;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
//此时的y[i] 表示第二关键字排第i的下标是y[i]
for(i=;i<n;i++) wv[i]=x[y[i]];
for(i=;i<m;i++) ws[i]=;
for(i=;i<n;i++) ws[wv[i]]++;
for(i=;i<m;i++) ws[i]+=ws[i-];
for(i=n-;i>=;i--) sa[--ws[wv[i]]]=y[i];
for(t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
x[sa[i]]=cmp(y,sa[i],sa[i-],j)?p-:p++;
}
for(int i=;i<n;i++)
{
r[sa[i]]=i;
}
}
void calHeight(int n)
{
int h=;
for(int i=;i<n;i++)
{
if(r[i]==) h=;
else
{
int k=sa[r[i]-];
if(--h<) h=;
while(s1[k+h]==s1[i+h]) h++;
}
height[r[i]]=h;
}
}
int f[][];
int min(int x,int y)
{
return x<y?x:y;
}
int max(int x,int y)
{
return x<y?y:x;
}
void rmq(int n)
{
for(int j=;j<n;j++) f[j][]=height[j];
for(int i=;i<;i++)
{
for(int j=;j<n;j++)
{
if(j+(<<i)- < n)
{
f[j][i]=min(f[j][i-],f[j+(<<i-)][i-]);
}
}
}
}
int look(int x,int y)
{
x=r[x];
y=r[y];
if(x>y) x^=y^=x^=y;
x++;
int k=(int)log2((double)(y-x+));
return min(f[x][k],f[y-(<<k)+][k]);
} bool cmp(int sl,int l)
{
if(r[sl]>r[l]) return ;
return ;
} void solve(int n)
{
int maxt=,sl=,sr=;
for(int i=;i<=n/;i++)
{
for(int j=;j+i<n;j+=i)
{
int l=j;
int r=j+i;
int lcp=look(l,r);
int d=lcp/i+;
int t=j-(i-lcp%i);
for(int k=l-;k>=&&k+i>j&&s1[k]==s1[k+i];k--)
{
if(k==t)
{
d++;
l=k;
}
else if(cmp(l,k))
{
l=k;
}
}
if(d>maxt)
{
sl=l;
sr=l+d*i-;
maxt=d;
}
else if(d==maxt&&cmp(sl,l))
{
sl=l;
sr=l+d*i-;
}
}
}
for(int i=sl;i<=sr;i++)
{
printf("%c",s1[i]);
}
printf("\n");
} int main()
{
int cas=;
while(scanf("%s",s1)!=EOF)
{
if(s1[]=='#') break;
int len=strlen(s1);
for(int i=;i<len;i++)
r[i]=s1[i]-'a'+;
r[len++]=;
da(r,len,);
calHeight(len);
rmq(len);
printf("Case %d: ",cas++);
solve(len-); }
return ;
}
HDU 2459 Maximum repetition substring的更多相关文章
- 【HDOJ】2459 Maximum repetition substring
后缀数组+RMQ. /* 2459 */ #include <iostream> #include <sstream> #include <string> #inc ...
- POJ3693 Maximum repetition substring [后缀数组 ST表]
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9458 Acc ...
- Maximum repetition substring 后缀数组
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 7578 Acc ...
- Maximum repetition substring (poj3693 后缀数组求重复次数最多的连续重复子串)
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 6328 Acc ...
- POJ 3693 Maximum repetition substring(最多重复次数的子串)
Maximum repetition substring Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 10461 Ac ...
- POJ3693 Maximum repetition substring —— 后缀数组 重复次数最多的连续重复子串
题目链接:https://vjudge.net/problem/POJ-3693 Maximum repetition substring Time Limit: 1000MS Memory Li ...
- POJ3693 Maximum repetition substring 后缀数组
POJ - 3693 Maximum repetition substring 题意 输入一个串,求重复次数最多的连续重复字串,如果有次数相同的,则输出字典序最小的 Sample input ccab ...
- Maximum repetition substring(POJ - 3693)(sa(后缀数组)+st表)
The repetition number of a string is defined as the maximum number \(R\) such that the string can be ...
- POJ 3693 Maximum repetition substring(后缀数组)
Description The repetition number of a string is defined as the maximum number R such that the strin ...
随机推荐
- 【Git】Git pull 强制覆盖本地文件
git fetch --all git reset --hard origin/master git pull 备注: git fetch 只是下载远程的库的内容,不做任何的合并 git reset ...
- [docker][win10]安装的坑
右键这个小图标,先signin,注意这里是ID 不是邮箱 image.png 可能starting 时候就报错说 “Containers feature is not enabled” 或者 ...
- UVA11400-Lighting System Design(动态规划基础)
Problem UVA11400-Lighting System Design Accept: 654 Submit: 4654Time Limit: 3000 mSec Problem Descr ...
- Hive窗口函数
参考地址:https://cwiki.apache.org/confluence/display/Hive/LanguageManual+WindowingAndAnalytics 环境准备: CRE ...
- LeetCode练习2 两数相加
问题如下: 给出两个 非空 的链表用来表示两个非负的整数.其中,它们各自的位数是按照 逆序 的方式存储的,并且它们的每个节点只能存储 一位 数字. 如果,我们将这两个数相加起来,则会返回一个新的链表来 ...
- flask_sqlalchemy中根据聚合分组后的结果进行排序,根据日期(datetime)按天分组
from sqlalchemy import func, desc # 根据聚合查询总收入,按总收入逆序 s= db.session.query(TpOrders.room_type_id, (fun ...
- 数组升序排序的方法Arrays.sort();的应用
package com.Summer_0421.cn; import java.util.Arrays; /** * @author Summer * 数组升序排序的方法Arrays.sort();应 ...
- centos 7 安装mqtt 修改用户名和密码
我先新买个Centos 的系统 咱登录呢就用这个软件,主要是方便,可以少安装一些东西 根据自己的系统选择,上面的是32位的. 输入 root 回车 输入密码然后回车 输入的时候什么也不显示 输入 c ...
- Codeforces Round #545 (Div. 1) Solution
人生第一场Div. 1 结果因为想D想太久不晓得Floyd判环法.C不会拆点.E想了个奇奇怪怪的set+堆+一堆乱七八糟的标记的贼难写的做法滚粗了qwq靠手速上分qwqqq A. Skyscraper ...
- Spark SQL在100TB上的自适应执行实践(转载)
Spark SQL是Apache Spark最广泛使用的一个组件,它提供了非常友好的接口来分布式处理结构化数据,在很多应用领域都有成功的生产实践,但是在超大规模集群和数据集上,Spark SQL仍然遇 ...