poj1743Musical Theme - 后缀数组

2024-05-12 00:18

本文主要是介绍poj1743Musical Theme - 后缀数组,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接:http://poj.org/problem?id=1743

题目大意:有N(1<=N<=20000)个音符的序列来表示一首乐曲,每个音符都是1..88范围内的整数,现在要找一个重复的子串,它需要满足如下条件:1.长度至少为5个音符。

2.在乐曲中重复出现(就是出现过至少两次)。(可能经过转调,“转调”的意思是主题序列中每个音符都被加上或减去了同一个整数值) 

3.重复出现的同一主题不能有公共部分。


题解:这道题是罗大神论文里的一道题,求不重复的最长子串。跑完倍增之后,二分答案,用h[]分组来验证。嗯是基本应用之一。可能有个难点就是对“转调”的处理,其实很容易可以发现同加同减的话差值还是不变,所以把原串变为两两间的差值再做sa

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define maxn 20100int sa[maxn],rk[maxn],s[maxn],n;
int Rsort[maxn],y[maxn],wr[maxn],h[maxn];
int mymin(int x,int y){return (x<y)?x:y;}
int mymax(int x,int y){return (x>y)?x:y;}
bool cmp(int x,int y,int ln){return (wr[x]==wr[y])&&(wr[x+ln]==wr[y+ln]);}
void get_h()
{int i,k=0;for (i=1;i<=n;i++){int j=sa[rk[i]-1];if (k>0) k--;while (s[i+k]==s[j+k]) k++;h[rk[i]]=k;}
}
void get_sa()
{int i,k,p,ln,m=200;for (i=1;i<=n;i++) rk[i]=s[i];for (i=0;i<=m;i++) Rsort[i]=0;for (i=1;i<=n;i++) Rsort[rk[i]]++;for (i=1;i<=m;i++) Rsort[i]+=Rsort[i-1];for (i=n;i>=1;i--) sa[Rsort[rk[i]]--]=i;ln=1;p=0;while (p<n){for (k=0,i=n-ln+1;i<=n;i++) y[++k]=i;for (i=1;i<=n;i++) if (sa[i]-ln>0) y[++k]=sa[i]-ln;for (i=1;i<=n;i++) wr[i]=rk[y[i]];for (i=0;i<=m;i++) Rsort[i]=0;for (i=1;i<=n;i++) Rsort[wr[i]]++;for (i=1;i<=m;i++) Rsort[i]+=Rsort[i-1];for (i=n;i>=1;i--) sa[Rsort[wr[i]]--]=y[i];memcpy(wr,rk,sizeof(wr));p=1;rk[sa[1]]=1;for (i=2;i<=n;i++){if (!cmp(sa[i],sa[i-1],ln)) p++;rk[sa[i]]=p;}m=p;ln*=2;}sa[0]=s[0]=0;
}
int main()
{//freopen("a.in","r",stdin);//freopen("a.out","w",stdout);int i,L,R,mid,mn,mx,ret;while (1){scanf("%d",&n);if (n==0) break;for (i=1;i<=n;i++) scanf("%d",&s[i]);for (i=n;i>=1;i--) s[i]=(s[i]-s[i-1])+88;get_sa();get_h();//求sa[]、h[]L=1;R=n;ret=0;while (L<R)//二分{mid=(L+R)>>1;bool bk=0;mn=n+1;mx=0;for (i=1;i<=n;i++)if (h[i]<mid) mx=mn=sa[i];else{mx=mymax(mx,sa[i]);mn=mymin(mn,sa[i]);if (mx-mn>mid) {bk=1;break;}//要求不重叠}if (bk) {L=mid+1;ret=mymax(ret,mid);}else R=mid;}ret++;//因为是用差值做的所以长度要+1if (ret<5) ret=0;printf("%d\n",ret);}return 0;
}


这篇关于poj1743Musical Theme - 后缀数组的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/981066

相关文章

hdu2241(二分+合并数组)

题意:判断是否存在a+b+c = x,a,b,c分别属于集合A,B,C 如果用暴力会超时,所以这里用到了数组合并,将b,c数组合并成d,d数组存的是b,c数组元素的和,然后对d数组进行二分就可以了 代码如下(附注释): #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

C语言:柔性数组

数组定义 柔性数组 err int arr[0] = {0}; // ERROR 柔性数组 // 常见struct Test{int len;char arr[1024];} // 柔性数组struct Test{int len;char arr[0];}struct Test *t;t = malloc(sizeof(Test) + 11);strcpy(t->arr,

C 语言基础之数组

文章目录 什么是数组数组变量的声明多维数组 什么是数组 数组,顾名思义,就是一组数。 假如班上有 30 个同学,让你编程统计每个人的分数,求最高分、最低分、平均分等。如果不知道数组,你只能这样写代码: int ZhangSan_score = 95;int LiSi_score = 90;......int LiuDong_score = 100;int Zhou

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

C# double[] 和Matlab数组MWArray[]转换

C# double[] 转换成MWArray[], 直接赋值就行             MWNumericArray[] ma = new MWNumericArray[4];             double[] dT = new double[] { 0 };             double[] dT1 = new double[] { 0,2 };

PHP7扩展开发之数组处理

前言 这次,我们将演示如何在PHP扩展中如何对数组进行处理。要实现的PHP代码如下: <?phpfunction array_concat ($arr, $prefix) {foreach($arr as $key => $val) {if (isset($prefix[$key]) && is_string($val) && is_string($prefix[$key])) {$arr[

Go 数组赋值问题

package mainimport "fmt"type Student struct {Name stringAge int}func main() {data := make(map[string]*Student)list := []Student{{Name:"a",Age:1},{Name:"b",Age:2},{Name:"c",Age:3},}// 错误 都指向了最后一个v// a

码蹄集部分题目(2024OJ赛9.4-9.8;线段树+树状数组)

1🐋🐋配对最小值(王者;树状数组) 时间限制:1秒 占用内存:64M 🐟题目思路 MT3065 配对最小值_哔哩哔哩_bilibili 🐟代码 #include<bits/stdc++.h> using namespace std;const int N=1e5+7;int a[N],b[N],c[N],n,q;struct QUERY{int l,r,id;}que

【第0006页 · 数组】寻找重复数

【前言】本文以及之后的一些题解都会陆续整理到目录中,若想了解全部题解整理,请看这里: 第0006页 · 寻找重复数         今天想讨论的一道题在 LeetCode 上评论也是颇为“不错”。有一说一,是道好题,不过我们还是得先理解了它才算真正的好题。这里我们展示一种使用二进制的做法,希望能帮到你哟! 【寻找重复数】给定一个包含 n + 1 个整数的数组 nums ,其数字都