标签：poj1743 musical themes 后缀数组
题意：
我们用N(1 <= N <=5000)个音符的序列来表示一首乐曲，每个音符都是1..88范围内的整数，每个数表示钢琴上的一个键。很不幸这种表示旋律的方法忽略了音符的时值，但这项编程任务是关于音高的，与时值无关。
许多作曲家围绕一个重复出现的“主题”来构建乐曲。在我们的乐曲表示法中，“主题”是整个音符序列的一个子序列，它需要满足如下条件：
“转调”的意思是主题序列中每个音符都被加上或减去了同一个整数值。
给定一段乐曲，计算其中最长主题的长度(即音符数)。
输出文件的第一行包含整数N。下面的每一行(最后一行可能除外)包含20个整数，表示音符序列。最后一行可能少于20个音符。
输出文件应只含一个整数，即最长主题的长度。如果乐曲中没有主题，那么输出0。
样例中这个长度为5的主题是输入文件中第一行的最后5个音符和第二行开头5个音符
题解：
啊，写个后缀数组，然后二分check时根据height分组，使每一段都满足height足够长，扫一遍，如果里面离的最远两个没重合（sa的差足够大），就return 1。。。
代码：
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 21000
using namespace std;
int s[N];
int sa[N],rank[N],h[N],n,m,len;
int cnt[N],val[N],stk[N],_val[N],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 Themes 乐曲主题 后缀数组、 （也可以用hash+二分做的~）
标签：poj1743 musical themes 后缀数组
原文地址：http://blog.csdn.net/vmurder/article/details/42399035
