重复旋律3

时间限制:5000ms
单点时限:1000ms
内存限制:256MB

描述

小Hi平时的一大兴趣爱好就是演奏钢琴。我们知道一个音乐旋律被表示为长度为 N 的数构成的数列。小Hi在练习过很多曲子以后发现很多作品中的旋律有共同的部分。

旋律是一段连续的数列,如果同一段旋律在作品A和作品B中同时出现过,这段旋律就是A和B共同的部分,比如在abab 在 bababab 和 cabacababc 中都出现过。小Hi想知道两部作品的共同旋律最长是多少?

输入

共两行。一行一个仅包含小写字母的字符串。字符串长度不超过 100000。

输出

一行一个整数,表示答案。

样例输入
abcdefg
abacabca
样例输出
3
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=;
int Rank[maxn],cntA[maxn],cntB[maxn],A[maxn],B[maxn];
int sa[maxn],tsa[maxn],ht[maxn];
int ans,N,ch[maxn];
char str1[maxn],str2[maxn];
void solve()
{
int i,j;
for(i=;i<=N;i++) cntA[i]=;
for(i=;i<=N;i++) cntA[ch[i]]++;
for(i=;i<=N;i++) cntA[i]+=cntA[i-];
for(i=N;i>=;i--) sa[cntA[ch[i]]--]=i;
Rank[sa[]]=;
for(i=;i<=N;i++)
{
Rank[sa[i]]=Rank[sa[i-]];
if(ch[sa[i]]!=ch[sa[i-]]) Rank[sa[i]]++;
}
for(int L=;Rank[sa[N]]<N;L<<=)
{
for(i=;i<=N;i++) cntA[i]=cntB[i]=;
for(i=;i<=N;i++) cntA[A[i]=Rank[i]]++;
for(i=;i<=N;i++) cntB[B[i]=(i+L<=N)?Rank[i+L]:]++;
for(i=;i<=N;i++) cntA[i]+=cntA[i-];
for(i=;i<=N;i++) cntB[i]+=cntB[i-];
for(i=N;i>=;i--) tsa[cntB[B[i]]--]=i;
for(i=N;i>=;i--) sa[cntA[A[tsa[i]]]--]=tsa[i];
Rank[sa[]]=;
for(i=;i<=N;i++)
{
Rank[sa[i]]=Rank[sa[i-]];
if(A[sa[i]]!=A[sa[i-]]||B[sa[i]]!=B[sa[i-]]) Rank[sa[i]]++;
}
}
for (i=,j=;i<=N;i++)
{
if(j) j --;
while (ch[i+j]==ch[sa[Rank[i]-]+j]) j++;
ht[Rank[i]]=j;
}
}
int main()
{
scanf("%s",str1+);
scanf("%s",str2+);
int L1=strlen(str1+);
int L2=strlen(str2+);
for(int i=;i<=L1;i++) ch[i]=str1[i]-;
ch[L1+]=;
for(int i=;i<=L2;i++) ch[i+L1+]=str2[i]-;
N=L1+L2+;
solve();
for(int i = ; i <= N; i++)
{
if((sa[i]<=L1)!=(sa[i-]<=L1))
ans = max(ans, ht[i]);
}
printf("%d\n",ans);
return ;
}

没有找到原因,开始用字符串一直错,改成数组就AC了。也不知道为什么。过几天再看看吧。

最新文章

  1. Windows 10 部署Enterprise Solution 5.5
  2. 取消IE默认下载工具为迅雷
  3. GDI+ 操作TIFF ccitt t.6 压缩
  4. 4809 江哥的dp题c
  5. MyISAM与InnoDB两者之间怎么选择
  6. 初探YAML
  7. IT人的自我导向型学习:学习的4个层次
  8. 再看JavaScript线程
  9. centos 正确 安装 jdk
  10. Gson源码分析之Json结构抽象和注解使用
  11. [翻译] C++ STL容器参考手册 (总册)
  12. ES6的语法
  13. gulp详细入门
  14. beta冲刺1-咸鱼
  15. 《Ext JS 4.2 实战》可以买了
  16. 设置redis服务开机自启动
  17. 【Python】【BugList13】req = requests.get(url=target)报错: (Caused by SSLError(SSLError(1, &#39;[SSL: CERTIFICATE_VERIFY_FAILED] certificate verify failed (_ssl.c:777)&#39;)
  18. Linux时间戳转换成BCD码(转载)
  19. 高级组件——文件选择器JFileChooser
  20. Python全栈开发记录_第一篇(循环练习及杂碎的知识点)

热门文章

  1. 016_笼统概述MapReduce执行流程结合wordcount程序
  2. bootstrap table 复选框使用
  3. $.queue() 与 $.dequeue() -- 队列
  4. cocos2dx打飞机项目笔记一:项目结构介绍
  5. FreeMarker缓存处理
  6. Excel下载打不开
  7. 深入Struts2的过滤器FilterDispatcher--中文乱码及字符编码过滤器
  8. HTTP与HTTPS有什么区别?
  9. input ajax自动补全
  10. spring security在spring mvc的action中获取登录人信息