博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
hihocoder1415 重复旋律3
阅读量:4661 次
发布时间:2019-06-09

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

思路:

扫一遍height 判一下即可

//By SiriusRen#include 
#include
#include
using namespace std;#define N 200050int n,cntA[N],cntB[N],A[N],B[N],tsa[N],sa[N],rk[N],ht[N],ans;char s[N],b[N];void SA(){ for(int i=1;i<=n;i++)cntA[s[i]]++; for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1]; for(int i=n;i;i--)sa[cntA[s[i]]--]=i; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]); for(int l=1;rk[sa[n]]
<<=1){ memset(cntA,0,sizeof(cntA)); memset(cntB,0,sizeof(cntB)); for(int i=1;i<=n;i++) cntA[A[i]=rk[i]]++, cntB[B[i]=i+l<=n?rk[i+l]:0]++; for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1]; for(int i=n;i;i--)tsa[cntB[B[i]]--]=i; for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i]; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]); } for(int i=1,j=0;i<=n;i++){ j=j?j-1:0; while(s[i+j]==s[sa[rk[i]-1]+j])j++; ht[rk[i]]=j; }}int main(){ scanf("%s%s",s+1,b+1); int r=strlen(s+1),m=strlen(b+1); for(int i=1;i<=m;i++)s[r+i+1]=b[i]; s[r+1]='#';n=r+m+1,SA(); for(int i=1;i<=n;i++){ int maxx=max(sa[i],sa[i-1]),minn=min(sa[i],sa[i-1]); if(maxx>r&&minn<=r)ans=max(ans,ht[i]); } printf("%d\n",ans);}

转载于:https://www.cnblogs.com/SiriusRen/p/6532070.html

你可能感兴趣的文章
java远程调用中出现的问题(主要是在不同电脑之间出现的问题)
查看>>
PHP正则表达式
查看>>
移除字符串中的汉字
查看>>
AC日记—— codevs 1031 质数环(搜索)
查看>>
一个restframework快速实例
查看>>
easyui Combotree根据用户输入显示对应的tree值
查看>>
C++之命名空间(End Chapter)
查看>>
获取http请求的响应状态
查看>>
【MFC两种视频图像採集方法】DirectShow与Opencv
查看>>
tensorflow中的batch_normalization实现
查看>>
20145219 《信息安全系统设计基础》第05周学习总结
查看>>
C#中隐藏(new)和方法重写(override)和重载(overload)的区别
查看>>
NET Core-学习笔记(三)
查看>>
解决 Unable to load native-hadoop library for your platform
查看>>
eclipse的debug
查看>>
SqlServer 由于未在SqlServer的此实例上安装复制组件解决方法
查看>>
SQLServer (多表左连接,指定操作某几行)
查看>>
在腾讯云服务器上安装JDK+Tomcat并启动tomcat
查看>>
远程连接mysql数据库注意点记录
查看>>
UOJ #206. 【APIO2016】Gap【交互题】
查看>>