博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
后缀数组 3160 最长公共子串
阅读量:6152 次
发布时间:2019-06-21

本文共 1484 字,大约阅读时间需要 4 分钟。

/*3160 最长公共子串

 

 时间限制: 2 s
 空间限制: 128000 KB
 题目等级 : 大师 Master
 
 
 
题目描述
 Description

给出两个由小写字母组成的字符串,求它们的最长公共子串的长度。

输入描述
 Input Description

读入两个字符串

输出描述
 Output Description

输出最长公共子串的长度*/

//后缀数组的简单应用,求最长公共子串,先将两个子串连起来,求h数组中的最大值 #include
#include
#include
#define N 2000005using namespace std;int a[N],sa[2][N],rk[2][N],v[N],n,k,mid,h[N],ans;char ch[N];void jin(int sa[N],int rk[N],int SA[N],int RK[N]){ for(int i=1;i<=n;i++) v[rk[sa[i]]]=i; for(int i=n;i;i--) if(sa[i]>k) SA[v[rk[sa[i]-k]]--]=sa[i]-k; for(int i=n-k+1;i<=n;i++) SA[v[rk[i]]--]=i; for(int i=1;i<=n;i++) RK[SA[i]]=RK[SA[i-1]]+(rk[SA[i]]!=rk[SA[i-1]]||rk[SA[i]+k]!=rk[SA[i-1]+k]); return;}int main(){ scanf("%s",ch+1); mid=strlen(ch+1); ch[++mid]='z'+1; scanf("%s",ch+mid+1); n=strlen(ch+1); int p=0,q=1; for(int i=1;i<=n;i++) a[i]=ch[i]-'a'+1; for(int i=1;i<=n;i++) v[a[i]]++; for(int i=1;i<31;i++) v[i]+=v[i-1]; for(int i=1;i<=n;i++) sa[p][v[a[i]]--]=i; for(int i=1;i<=n;i++) rk[p][sa[p][i]]=rk[p][sa[p][i-1]]+(a[sa[p][i]]!=a[sa[p][i-1]]); for(k=1;k
<<=1,swap(p,q)) jin(sa[p],rk[p],sa[q],rk[q]); k=0; for(int i=1;i<=n;i++) if(rk[p][i]==1) h[1]=0; else { int j=sa[p][rk[p][i]-1]; for(;a[j+k]==a[i+k];k++); h[rk[p][i]]=k; if(k>0) k--; } for(int i=2;i<=n;i++) if((sa[p][i]>mid)^(sa[p][i-1]>mid)) ans=max(ans,h[i]); printf("%d",ans); return 0;}

转载于:https://www.cnblogs.com/xydddd/p/5153344.html

你可能感兴趣的文章
修改GRUB2背景图片
查看>>
Ajax异步
查看>>
好记性不如烂笔杆-android学习笔记<十六> switcher和gallery
查看>>
JAVA GC
查看>>
3springboot:springboot配置文件(外部配置加载顺序、自动配置原理,@Conditional)
查看>>
前端第七天
查看>>
图解SSH原理及两种登录方法
查看>>
[转载] 七龙珠第一部——第058话 魔境圣地
查看>>
【总结整理】JQuery基础学习---样式篇
查看>>
查询个人站点的文章、分类和标签查询
查看>>
基础知识:数字、字符串、列表 的类型及内置方法
查看>>
JSP的隐式对象
查看>>
JS图片跟着鼠标跑效果
查看>>
[SCOI2005][BZOJ 1084]最大子矩阵
查看>>
学习笔记之Data Visualization
查看>>
Leetcode 3. Longest Substring Without Repeating Characters
查看>>
416. Partition Equal Subset Sum
查看>>
Vue之项目搭建
查看>>
app内部H5测试点总结
查看>>
[TC13761]Mutalisk
查看>>