【后缀数组+RMQ】 SPOJ 687_后缀数组 rmq-程序员宅基地

技术标签: 算法  字符串  

题意:
问给出的串中连续重复的子串的最长的长度

思路:
先将后缀数组和高度数组先求出来,分别为sa和lcp数组。
然后用RMQ的ST算法将从位置i和j开始的后缀的最长公共前缀求出来,这个还好理解把,对与排名第一的后缀,它的高度数组为与排名第二的后缀的最长公共前缀,依次类推,比如到第四个,前四个分别为5,3,0,2;那么第一个与第三个的公共前缀就为0,以第二个为媒介,第二个与第一个的最长为5,第二个与第三个为0,那么第一个与第三个为0,与3后面的全为0,所以ST求最小即可,然后就是枚举要重复的串的长度。

#include <stdio.h>  
#include <stdlib.h>  
#include <string.h>  
#include <iostream>  
#include <algorithm>  
using namespace std;  
typedef long long ll;  
const int inf=0x3f3f3f3f;  
const int MAXN=50010;  
int wa[MAXN],wb[MAXN],wv[MAXN],ww[MAXN];  
int sa[MAXN],lcp[MAXN],Rank[MAXN],rank1[MAXN],dp[MAXN][20];  
char str[MAXN];  
int nn;  
inline bool cmp(int *r,int a,int b,int len){  
    return r[a]==r[b]&&r[a+len]==r[b+len];  
}  
 void construct_sa(int n,int m){  
     int i,j,p,*x=wa,*y=wb,*t;n++;  
     for(i=0;i<m;i++) ww[i]=0;  
     for(i=0;i<n;i++) ww[x[i]=str[i]]++;  
     for(i=1;i<m;i++) ww[i]+=ww[i-1];  
     for(i=n-1;i>=0;i--) sa[--ww[x[i]]]=i;  
     for(j=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<m;i++) ww[i]=0;  
         for(i=0;i<n;i++) ww[wv[i]=x[y[i]]]++;  
         for(i=1;i<m;i++) ww[i]+=ww[i-1];  
         for(i=n-1;i>=0;i--) sa[--ww[wv[i]]]=y[i];  
         for(t=x,x=y,y=t,x[sa[0]]=0,p=i=1;i<n;i++)  
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;  
     }  
}  
void construct_lcp(int n){  
    for(int i=0;i<=n;i++) rank1[sa[i]]=i;  
    int h=0;  
    lcp[0]=0;  
    for(int i=0;i<n;i++){  
        int j=sa[rank1[i]-1];  
        if(h>0) h--;  
        for(;j+h<n&&i+h<n;h++) if(str[i+h]!=str[j+h]) break;  
        lcp[rank1[i]-1]=h;  
    }  
}  
void RMQ_init(int n){  
    for(int i=1;i<=n;i++)  
        dp[i][0]=lcp[i-1];  
    for(int i=1;(1<<i)<=n;i++){  
        for(int j=1;j+(1<<i)-1<=n;j++){  
            dp[j][i]=min(dp[j][i-1],dp[j+(1<<(i-1))][i-1]);  
        }  
    }  
}  
int RMQ(int le,int ri){  
    le=rank1[le];ri=rank1[ri];//找到各自位置的排名  
    if(le>ri) swap(le,ri);le++;//为什么这里加1,因为dp中dp[i][0]的值为lcp[i-1],dp的i从1开始,而lcp的从0开始,在这里补回来  
                                //而为什么ri不补回来呢,因为lcp代表的意义是第i大的与第i+1大的最长前缀,我们不需要lcp[i]的值  
    int k=0;  
    while((1<<(k+1))<=ri-le+1) k++;  
    int ans2=min(dp[le][k],dp[ri-(1<<k)+1][k]);  
    return ans2;//返回排名从le到ri的最小值为从i和j开始的最长公共前缀  
}  
int main(){  
    int T;  
    scanf("%d",&T);  
    while(T--){  
        scanf("%d",&nn);  
        for(int i=0;i<nn;i++){  
            getchar();  
            scanf("%c",&str[i]);  
        }  
        construct_sa(nn,100);  
        construct_lcp(nn);  
        RMQ_init(nn);  
        int ans=0,sum;  
        for(int len=1;len<nn;len++){
   //枚举重复的串的长度  
            for(int i=0;i+len<nn;i+=len){
   //这里为什么直接加len就行,是因为  
                int t=RMQ(i,i+len);//找i和i+len的最长前缀  
                sum=t/len+1;//前缀能拼成的长度为len的个数加上i到i+len的一个  
                int pos=i-(len-t%len);//这个串的长度如果还有富裕,那么看看前面是否添加几个字符可以将重复次数在+1;  
                if(pos>=0&&t%len!=0) if(RMQ(pos,pos+len)>=(len-t%len)) sum++;//这个判断就是你缺几个字符,如果从pos开始把这几个  
                                                                            //字符填上了,那么就可以在组成一个  
                if(sum>ans) ans=sum;  
            }  
        }  
        printf("%d\n",ans);  
//        for(int i=0;i<=nn;i++) cout<<sa[i]<<" ";cout<<endl;  
//        for(int i=0;i<nn;i++) cout<<lcp[i]<<" ";cout<<endl;  
//        for(int i=0;i<nn;i++) cout<<rank1[i]<<" ";cout<<endl;  
    }  
    return 0;  
}  
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/Starry_starry_night/article/details/77341815

智能推荐

使用nginx解决浏览器跨域问题_nginx不停的xhr-程序员宅基地

文章浏览阅读1k次。通过使用ajax方法跨域请求是浏览器所不允许的,浏览器出于安全考虑是禁止的。警告信息如下:不过jQuery对跨域问题也有解决方案,使用jsonp的方式解决,方法如下:$.ajax({ async:false, url: 'http://www.mysite.com/demo.do', // 跨域URL ty..._nginx不停的xhr

在 Oracle 中配置 extproc 以访问 ST_Geometry-程序员宅基地

文章浏览阅读2k次。关于在 Oracle 中配置 extproc 以访问 ST_Geometry,也就是我们所说的 使用空间SQL 的方法,官方文档链接如下。http://desktop.arcgis.com/zh-cn/arcmap/latest/manage-data/gdbs-in-oracle/configure-oracle-extproc.htm其实简单总结一下,主要就分为以下几个步骤。..._extproc

Linux C++ gbk转为utf-8_linux c++ gbk->utf8-程序员宅基地

文章浏览阅读1.5w次。linux下没有上面的两个函数,需要使用函数 mbstowcs和wcstombsmbstowcs将多字节编码转换为宽字节编码wcstombs将宽字节编码转换为多字节编码这两个函数,转换过程中受到系统编码类型的影响,需要通过设置来设定转换前和转换后的编码类型。通过函数setlocale进行系统编码的设置。linux下输入命名locale -a查看系统支持的编码_linux c++ gbk->utf8

IMP-00009: 导出文件异常结束-程序员宅基地

文章浏览阅读750次。今天准备从生产库向测试库进行数据导入,结果在imp导入的时候遇到“ IMP-00009:导出文件异常结束” 错误,google一下,发现可能有如下原因导致imp的数据太大,没有写buffer和commit两个数据库字符集不同从低版本exp的dmp文件,向高版本imp导出的dmp文件出错传输dmp文件时,文件损坏解决办法:imp时指定..._imp-00009导出文件异常结束

python程序员需要深入掌握的技能_Python用数据说明程序员需要掌握的技能-程序员宅基地

文章浏览阅读143次。当下是一个大数据的时代,各个行业都离不开数据的支持。因此,网络爬虫就应运而生。网络爬虫当下最为火热的是Python,Python开发爬虫相对简单,而且功能库相当完善,力压众多开发语言。本次教程我们爬取前程无忧的招聘信息来分析Python程序员需要掌握那些编程技术。首先在谷歌浏览器打开前程无忧的首页,按F12打开浏览器的开发者工具。浏览器开发者工具是用于捕捉网站的请求信息,通过分析请求信息可以了解请..._初级python程序员能力要求

Spring @Service生成bean名称的规则(当类的名字是以两个或以上的大写字母开头的话,bean的名字会与类名保持一致)_@service beanname-程序员宅基地

文章浏览阅读7.6k次,点赞2次,收藏6次。@Service标注的bean,类名:ABDemoService查看源码后发现,原来是经过一个特殊处理:当类的名字是以两个或以上的大写字母开头的话,bean的名字会与类名保持一致public class AnnotationBeanNameGenerator implements BeanNameGenerator { private static final String C..._@service beanname

随便推点

二叉树的各种创建方法_二叉树的建立-程序员宅基地

文章浏览阅读6.9w次,点赞73次,收藏463次。1.前序创建#include&lt;stdio.h&gt;#include&lt;string.h&gt;#include&lt;stdlib.h&gt;#include&lt;malloc.h&gt;#include&lt;iostream&gt;#include&lt;stack&gt;#include&lt;queue&gt;using namespace std;typed_二叉树的建立

解决asp.net导出excel时中文文件名乱码_asp.net utf8 导出中文字符乱码-程序员宅基地

文章浏览阅读7.1k次。在Asp.net上使用Excel导出功能,如果文件名出现中文,便会以乱码视之。 解决方法: fileName = HttpUtility.UrlEncode(fileName, System.Text.Encoding.UTF8);_asp.net utf8 导出中文字符乱码

笔记-编译原理-实验一-词法分析器设计_对pl/0作以下修改扩充。增加单词-程序员宅基地

文章浏览阅读2.1k次,点赞4次,收藏23次。第一次实验 词法分析实验报告设计思想词法分析的主要任务是根据文法的词汇表以及对应约定的编码进行一定的识别,找出文件中所有的合法的单词,并给出一定的信息作为最后的结果,用于后续语法分析程序的使用;本实验针对 PL/0 语言 的文法、词汇表编写一个词法分析程序,对于每个单词根据词汇表输出: (单词种类, 单词的值) 二元对。词汇表:种别编码单词符号助记符0beginb..._对pl/0作以下修改扩充。增加单词

android adb shell 权限,android adb shell权限被拒绝-程序员宅基地

文章浏览阅读773次。我在使用adb.exe时遇到了麻烦.我想使用与bash相同的adb.exe shell提示符,所以我决定更改默认的bash二进制文件(当然二进制文件是交叉编译的,一切都很完美)更改bash二进制文件遵循以下顺序> adb remount> adb push bash / system / bin /> adb shell> cd / system / bin> chm..._adb shell mv 权限

投影仪-相机标定_相机-投影仪标定-程序员宅基地

文章浏览阅读6.8k次,点赞12次,收藏125次。1. 单目相机标定引言相机标定已经研究多年,标定的算法可以分为基于摄影测量的标定和自标定。其中,应用最为广泛的还是张正友标定法。这是一种简单灵活、高鲁棒性、低成本的相机标定算法。仅需要一台相机和一块平面标定板构建相机标定系统,在标定过程中,相机拍摄多个角度下(至少两个角度,推荐10~20个角度)的标定板图像(相机和标定板都可以移动),即可对相机的内外参数进行标定。下面介绍张氏标定法(以下也这么称呼)的原理。原理相机模型和单应矩阵相机标定,就是对相机的内外参数进行计算的过程,从而得到物体到图像的投影_相机-投影仪标定

Wayland架构、渲染、硬件支持-程序员宅基地

文章浏览阅读2.2k次。文章目录Wayland 架构Wayland 渲染Wayland的 硬件支持简 述: 翻译一篇关于和 wayland 有关的技术文章, 其英文标题为Wayland Architecture .Wayland 架构若是想要更好的理解 Wayland 架构及其与 X (X11 or X Window System) 结构;一种很好的方法是将事件从输入设备就开始跟踪, 查看期间所有的屏幕上出现的变化。这就是我们现在对 X 的理解。 内核是从一个输入设备中获取一个事件,并通过 evdev 输入_wayland

推荐文章

热门文章

相关标签