博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
后缀数组 POJ 1743 Musical Theme
阅读量:5281 次
发布时间:2019-06-14

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

 

题意:给定n个数字,求超过5个数字的,最长的,变化相同的,不相交的重复子串

分析:男人8题中的一题!数列相邻两项做差,形成新数列,即求数列中的最长重复子串(不可相交)。

后缀数组+二分答案。假如二分得到答案L,如何知道它是可行的呢? 因为对于排序后的后缀,Lcp ( Suffix ( List [ i ] ) , Suffix ( List [ i - 1 ] ) ) 是所有与Suffix ( List [ i ] )的LCP值中最大的一个。 因为 Height [ i ] 表示的是排序后后缀数组中第i个后缀和第i-1个后缀的LCP值。 那么对于后缀数组中的一段 L - R , 若 Height [ L + 1 ] ~ Height [ R ] 全部大于等于L,那么就等价于第L到第R个后缀中任意两个后缀的LCP值都大于等于L。 那么只要取这里面相隔最远的两个后缀,若他们相距大于L,那么就是可行的。 ( 为什么不是等于L呢 ? 因为我们取的关键字是 S[i]-S[i-1] , 若相距等于L,那么两段里面的首尾相连了,是不符合条件的)

简单来说,先对height数组分段,然后看每段是否有满足题意的子串。

#include 
#include
#include
const int N = 2e4 + 5;int sa[N], rank[N], height[N];int t[N], t2[N], c[N];int a[N];void da(int *s, int n, int m = 128) { int i, p, *x = t, *y = t2; for (i=0; i
=0; --i) sa[--c[x[i]]] = i; for (int k=1; k<=n; k<<=1) { for (p=0, i=n-k; i
= k) y[p++] = sa[i] - k; for (i=0; i
=0; --i) sa[--c[x[y[i]]]] = y[i]; std::swap (x, y); p = 1; x[sa[0]] = 0; for (i=1; i
= n) break; m = p; }}void calc_height(int n) { int i, k = 0; for (i=0; i
= m) { mn = std::min (mn, std::min (sa[i], sa[i-1])); mx = std::max (mx, std::max (sa[i], sa[i-1])); if (mn + m < mx) { return true; } } else { mn = mx = sa[i]; } } return false;}int main() { while (scanf ("%d", &n) == 1) { if (!n) break; for (int i=0; i
> 1; if (check (mid)) { ans = std::max (ans, mid); left = mid + 1; } else { right = mid - 1; } } if (ans >= 4) { printf ("%d\n", ans + 1); } else { puts ("0"); } } return 0;}

  

转载于:https://www.cnblogs.com/Running-Time/p/5450033.html

你可能感兴趣的文章
MFC中theApp
查看>>
类的无参方法
查看>>
Python 开发:初识Python(记笔记)
查看>>
sqlrelay 的安装配置和应用
查看>>
idea maven指定编译参数
查看>>
python functions as objects
查看>>
【刷题】SPOJ 705 SUBST1 - New Distinct Substrings
查看>>
【C语言】输入一个整数,求它的原码,反码,补码值
查看>>
wappalyzer 上各种开源框架功能
查看>>
Linux----知识储备
查看>>
使用Docker搭建Wordpress
查看>>
Extract Fasta Sequences Sub Sets by position
查看>>
Lua 随机数生成问题
查看>>
visual studio code jupyter错误!No kernel specs found
查看>>
oracle 基础(上)
查看>>
hdu3853LOOPS(概率与期望dp)
查看>>
今天心情不太好,没什么可写的
查看>>
python --条件、循环和其他语句(一)
查看>>
Dojo树使用心得
查看>>
SQL Server 存储过程 sp_helptext的不足以及解决方案
查看>>