题目: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的更多相关文章

  1. 【HDOJ】2459 Maximum repetition substring

    后缀数组+RMQ. /* 2459 */ #include <iostream> #include <sstream> #include <string> #inc ...

  2. POJ3693 Maximum repetition substring [后缀数组 ST表]

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 9458   Acc ...

  3. Maximum repetition substring 后缀数组

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 7578   Acc ...

  4. Maximum repetition substring (poj3693 后缀数组求重复次数最多的连续重复子串)

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 6328   Acc ...

  5. POJ 3693 Maximum repetition substring(最多重复次数的子串)

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 10461   Ac ...

  6. POJ3693 Maximum repetition substring —— 后缀数组 重复次数最多的连续重复子串

    题目链接:https://vjudge.net/problem/POJ-3693 Maximum repetition substring Time Limit: 1000MS   Memory Li ...

  7. POJ3693 Maximum repetition substring 后缀数组

    POJ - 3693 Maximum repetition substring 题意 输入一个串,求重复次数最多的连续重复字串,如果有次数相同的,则输出字典序最小的 Sample input ccab ...

  8. 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 ...

  9. POJ 3693 Maximum repetition substring(后缀数组)

    Description The repetition number of a string is defined as the maximum number R such that the strin ...

随机推荐

  1. 【Linux基础】文件处理实例

    1.文件拆分 //每4000行拆分一个文件 epms_t_ep_fx_stl_xy_20190129.dat 2.行处理 //查找第二列为711611且第三列为711100记录,打印行号和整行数据 a ...

  2. Spring Boot 菜鸟教程 application.properties 常用配置

    SPRING CONFIG (ConfigFileApplicationListener) spring.config.name 配置文件名称,默认为application spring.config ...

  3. docker学习笔记(二)-dockerfile docker仓库

    dckerfile2这样写才对,不然打印的是hello $name 例子 登陆https://hub.docker.com 上面是提交到公共仓库 创建本地私有仓库 开启另外一个虚机,将重新构建上传镜像 ...

  4. UVA11925-Generating Permutations(贪心)

    Problem UVA11925-Generating Permutations Accept: 214  Submit: 1429Time Limit: 1000 mSec Problem Desc ...

  5. UVA11694-Gokigen Naname(DFS进阶)

    Problem UVA11694-Gokigen Naname Accept: 76   Submit: 586Time Limit: 10000 mSec Problem Description I ...

  6. 转://Oracle 复合压缩索引场景及性能对比

    摘要:今天为什么提到这个话题,出于一个偶然,一个同事在优化新开发的系统时向我请教如何添加复合压缩索引的问题.我总结了一下,问题有三. 第一:需不需要压缩 第二:对第几列压缩 第三:性能对比,选出最优 ...

  7. P1316 丢瓶盖(二分+贪心)

    思路:都在注解里 #include<iostream> #include<algorithm> using namespace std; ; int a[maxn], n, m ...

  8. 5linux引导流程解析

    课程大纲  Linux引导流程  Linux运行级别  Linux启动服务管理  GRUB配置与应用  启动故障分析与解决     常用固件设置 安全设置 可引导介质列表 可引导介质搜索顺序 ...

  9. 2017-2018-2 20155314《网络对抗技术》Exp7 网络欺诈防范

    2017-2018-2 20155314<网络对抗技术>Exp7 网络欺诈防范 目录 实验目标 实验内容 实验环境 基础问题回答 预备知识 实验步骤 1 利用setoolkit建立冒名网站 ...

  10. randi( )函数--MATLAB

    randi()函数生成均匀分布的伪随机整数,范围为imin--imax,如果没指定imin,则默认为1. r = randi(imax,n):生成n*n的矩阵 r = randi(imax,m,n): ...