【POJ1743】Musical Themes 乐曲主题 后缀数组、 (也可以用hash+二分做的~)
2015-01-04 20:35
323 查看
题意:
Time Limit: 1 Sec Memory Limit: 128 MB
Submit: 42 Solved: 15
[Submit][Status][Web
Board]
我们用N(1 <= N <=5000)个音符的序列来表示一首乐曲,每个音符都是1..88范围内的整数,每个数表示钢琴上的一个键。很不幸这种表示旋律的方法忽略了音符的时值,但这项编程任务是关于音高的,与时值无关。
许多作曲家围绕一个重复出现的“主题”来构建乐曲。在我们的乐曲表示法中,“主题”是整个音符序列的一个子序列,它需要满足如下条件:
长度至少为5个音符
在乐曲中重复出现(可能经过转调,见下)
重复出现的同一主题不能重叠
“转调”的意思是主题序列中每个音符都被加上或减去了同一个整数值。
给定一段乐曲,计算其中最长主题的长度(即音符数)。
输出文件的第一行包含整数N。下面的每一行(最后一行可能除外)包含20个整数,表示音符序列。最后一行可能少于20个音符。
输出文件应只含一个整数,即最长主题的长度。如果乐曲中没有主题,那么输出0。
3025 27 30 34 39 45 52 60 69 79 69 60 52 45 39 34 30 26 22 1882 78 74 70 66 67 64 60 65 80
5
样例中这个长度为5的主题是输入文件中第一行的最后5个音符和第二行开头5个音符
USACO
[Submit][Status][Web
Board]
Powered by LLQ. ©2012-2014 长春吉大附中实验学校.
题解:
啊,写个后缀数组,然后二分check时根据height分组,使每一段都满足height足够长,扫一遍,如果里面离的最远两个没重合(sa的差足够大),就return 1。。。
代码:
1829: Musical Themes 乐曲主题
Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 42 Solved: 15
[Submit][Status][Web
Board]
Description
我们用N(1 <= N <=5000)个音符的序列来表示一首乐曲,每个音符都是1..88范围内的整数,每个数表示钢琴上的一个键。很不幸这种表示旋律的方法忽略了音符的时值,但这项编程任务是关于音高的,与时值无关。许多作曲家围绕一个重复出现的“主题”来构建乐曲。在我们的乐曲表示法中,“主题”是整个音符序列的一个子序列,它需要满足如下条件:
长度至少为5个音符
在乐曲中重复出现(可能经过转调,见下)
重复出现的同一主题不能重叠
“转调”的意思是主题序列中每个音符都被加上或减去了同一个整数值。
给定一段乐曲,计算其中最长主题的长度(即音符数)。
Input
输出文件的第一行包含整数N。下面的每一行(最后一行可能除外)包含20个整数,表示音符序列。最后一行可能少于20个音符。
Output
输出文件应只含一个整数,即最长主题的长度。如果乐曲中没有主题,那么输出0。
Sample Input
3025 27 30 34 39 45 52 60 69 79 69 60 52 45 39 34 30 26 22 1882 78 74 70 66 67 64 60 65 80
Sample Output
5
HINT
样例中这个长度为5的主题是输入文件中第一行的最后5个音符和第二行开头5个音符
Source
USACO
Best Solution
C | C++ | Pascal |
---|---|---|
NULL | 80643.wyfcyx (0ms,1100KB,2168B) | 21804.Ares (72ms,292KB,561B) |
Board]
Powered by LLQ. ©2012-2014 长春吉大附中实验学校.
题解:
啊,写个后缀数组,然后二分check时根据height分组,使每一段都满足height足够长,扫一遍,如果里面离的最远两个没重合(sa的差足够大),就return 1。。。
代码:
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define N 21000 using namespace std; int s ; int sa ,rank ,h ,n,m,len; int cnt ,val ,stk ,_val ,top; bool issame(int a,int b,int hl) { return val[a]==val[b]&& ((a+hl>=len&&b+hl>=len)||(a+hl<len&&b+hl<len&&val[a+hl]==val[b+hl])); } void SA(int lim) { int i,j,k,hl; for(i=0;i<lim;i++)cnt[i]=0; for(i=0;i<len;i++)cnt[val[i]=s[i]]++; for(i=1;i<lim;i++)cnt[i]+=cnt[i-1]; for(i=len-1;i>=0;i--)sa[--cnt[val[i]]]=i; for(k=1;;k++) { top=0,hl=1<<(k-1); for(i=0;i<len;i++)if(sa[i]+hl>=len)stk[top++]=sa[i]; for(i=0;i<len;i++)if(sa[i]>=hl)stk[top++]=sa[i]-hl; for(i=0;i<lim;i++)cnt[i]=0; for(i=0;i<len;i++)cnt[val[i]]++; for(i=1;i<lim;i++)cnt[i]+=cnt[i-1]; for(i=len-1;i>=0;i--)sa[--cnt[val[stk[i]]]]=stk[i]; for(lim=i=0;i<len;lim++) { for(j=i;j<len-1&&issame(sa[j],sa[j+1],hl);j++); for(;i<=j;i++)_val[sa[i]]=lim; } for(i=0;i<len;i++)val[i]=_val[i]; if(lim==len)break; } for(i=0;i<len;i++)rank[sa[i]]=i; for(k=i=0;i<len;i++) { if(k)k--; if(!rank[i])continue; while(s[i+k]==s[sa[rank[i]-1]+k])k++; h[rank[i]]=k; } } bool check(int mid) { int l=sa[0],r=sa[0],i; for(i=0;i<len;i++) { if(h[i]<mid) { l=r=sa[i]; continue; } l=min(l,sa[i]); r=max(r,sa[i]); if(r-l>mid)return 1; } return 0; } int main() { // freopen("test.in","r",stdin); int i,j,k; scanf("%d",&len); scanf("%d",&s[0]); for(i=1;i<len;i++)scanf("%d",&s[i]),s[i-1]=s[i]-s[i-1]+100; SA(256); int l=0,r=len,mid,ans=0; while(l<=r) { if(r-l<=3) { for(i=l;i<=r;i++)if(check(i))ans=i; break; } mid=l+r>>1; if(check(mid))l=mid; else r=mid-1; } if(ans>=4)printf("%d\n",ans+1);else puts("0"); return 0; }
相关文章推荐
- 【POJ1743】Musical Theme【后缀数组】【二分】
- bnuoj 34990(后缀数组 或 hash+二分)
- poj1743 Musical Theme(后缀数组+二分判定)
- 【POJ1743】Musical Theme(后缀数组,二分)
- poj1743(后缀数组+二分)
- poj1743(后缀数组+二分)
- poj1743(后缀数组+二分--不可重叠最长重复子串)
- 【后缀数组】【poj1743】音乐主题
- poj 1743 Musical Theme(后缀数组+二分)
- poj 3261 后缀数组 找重复出现k次的子串(子串可以重叠)
- poj 1743 Musical Theme(最长不重叠重复子串 后缀数组+二分)
- hdu 3068 最长回文 (Manacher O(n) ) (扩展KMP和后缀数组也都可以做 O(nlogn))
- poj 1743 Musical Theme (后缀数组+二分答案)
- POJ 1226 Substrings(后缀数组+二分)
- POJ1743---Musical Theme (后缀数组+二分)
- hdu 5030(2014广州网络赛)后缀数组+二分!
- HDU 4898 The Revenge of the Princess’ Knight(后缀数组+二分+暴力)(2014 Multi-University Training Contest 4)
- poj 3294 Life Forms(后缀数组+二分)
- POJ1743 Musical Theme,后缀数组,最长重复不重叠子串,二分检索答案
- POJ 3261 Milk Patterns(后缀数组+二分答案+离散化)