Musical Themes

原创 2013年12月02日 22:57:37

题意:转载NOCOW翻译(http://www.nocow.cn/index.php/Translate:USACO/theme

描述

我们用N(1 <= N <=5000)个音符的序列来表示一首乐曲,每个音符都是1..88范围内的整数,每个数表示钢琴上的一个键。很不幸这种表示旋律的方法忽略了音符的时值,但这项编程任务是关于音高的,与时值无关。

许多作曲家围绕一个重复出现的“主题”来构建乐曲。在我们的乐曲表示法中,“主题”是整个音符序列的一个子串,它需要满足如下条件:

⒈长度至少为5个音符

⒉在乐曲中重复出现(可能经过转调,见下)

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

“转调”的意思是主题序列中每个音符都被加上或减去了同一个整数值。 给定一段乐曲,计算其中最长主题的长度(即音符数)。

本题时限为1秒钟!

[编辑]格式

PROGRAM NAME: theme

INPUT FORMAT

输入文件的第一行包含整数N。下面的每一行(最后一行可能除外)包含20个整数,表示音符序列。最后一行可能少于20个音符。

OUTPUT FORMAT

输出文件应只含一个整数,即最长主题的长度。如果乐曲中没有主题,那么输出0。

[编辑]SAMPLE INPUT (file theme.in)

30
25 27 30 34 39 45 52 60 69 79 69 60 52 45 39 34 30 26 22 18
82 78 74 70 66 67 64 60 65 80


[编辑]SAMPLE OUTPUT (file theme.out)

5

(这个长度为5的主题是输入文件中第一行的最后5个音符和第二行开头5个音符)


解题思路:

  1. 对序列进行差分(信号有电压漂移的情况下,也就是叠加了一个线性偏移的时候,对信号进行差分可以消除线性偏移)
  2. 求后缀数组sa[1...N-1],rank[1...N-1],height[1...N-1]
  3. 利用后缀数组求最大不重叠的重复子串

代码

/*
ID: zc.rene1
LANG: C
PROG: theme
 */

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#define MAX_N 5000
		
int N;
int diff[MAX_N];
int rank[MAX_N];
int sa[MAX_N];
int height[MAX_N];
int bucket[10][MAX_N + 1][2];

void GetInput(FILE *fin)
{
    int i;
    int prev, current, min = 100;

    fscanf(fin, "%d %d", &N, &prev);
    for (i=1; i<N; i++)
    {
	fscanf(fin, "%d", ¤t);
	diff[i] = current - prev;
	if (diff[i] < min)
	{
	    min = diff[i];
	}
	prev = current;
    }
    if (min < 0)
    {
	for (i=1; i<N; i++)
	{
	    diff[i] -= min;
	}
    }
}

void UpdateRank(int jump)
{
    int diff_tmp[2][MAX_N];
    int i, j, k, m, temp, prev;

    for (i=1; i<N; i++)
    {
	diff_tmp[0][i] = rank[i];
	diff_tmp[1][i] = i;
    }
    for (i=1; i<N; i++)
    {
	if (jump != 0)
	{
	    diff_tmp[0][i] *= 10000;
	    if (i + jump < N)
	    {
		diff_tmp[0][i] += diff_tmp[0][i + jump];
	    }
	}
    }
    /*bucket sort*/
    for (i=0; i<8; i++)
    {
	for (j=0; j<=9; j++)
	{
	    bucket[j][0][0] = 0;
	}

	for (j=1; j<N; j++)
	{
	    temp = diff_tmp[0][j];
	    for (k=0; k<i; k++)
	    {
		temp /= 10;
	    }
	    temp %= 10;
	    ++bucket[temp][0][0];
	    bucket[temp][bucket[temp][0][0]][0] = diff_tmp[0][j];
	    bucket[temp][bucket[temp][0][0]][1] = diff_tmp[1][j];
	}
	m = 0;
	for (j=0; j<=9; j++)
	{
	    for (k=1; k<=bucket[j][0][0]; k++)
	    {
		++m;
		diff_tmp[0][m] = bucket[j][k][0];
		diff_tmp[1][m] = bucket[j][k][1];
	    }
	}
    }
    /*generate rank*/
    temp = 1;
    rank[diff_tmp[1][1]] = temp;
    prev = diff_tmp[0][1];

    for (i=2; i<N; i++)
    {
	if (diff_tmp[0][i] > prev)
	{
	    temp++;
	}
	rank[diff_tmp[1][i]] = temp;
	prev = diff_tmp[0][i];
    }
}

void GenerateRank(void)
{
    int len;

    memcpy(rank, diff, MAX_N * sizeof(int));

    UpdateRank(0); 

    for (len=1; len/2 < N; len*=2)
    {
	UpdateRank(len);
    }
}

void GenerateSa(void)
{
    int i;

    for (i=1; i<N; i++)
    {
	sa[rank[i]] = i;
    }
}

void GenerateHeight(void)
{
    int i, j = 0;

    for (i=1; i<N; i++)
    {
	if (rank[i] == 1)
	{
	    continue;
	}
	for (; (i+j<N)&&(sa[rank[i]-1]+j<N)&&(diff[i+j]==diff[sa[rank[i]-1]+j]); )
	{
	    j++;
	}
	height[rank[i]] = j;
	if (j > 0)
	{
	    j--;
	}
    }
}

int Test(int ht)
{
    int i, j, min, max;

    for (i=1; i<N; i++)
    {
	if (height[i] >= ht)
	{
	    min = sa[i - 1];
	    max = sa[i - 1]; 
	    for (j=i; j<N&&height[j]>=ht; j++)
	    {
		if (sa[j] < min)
		{
		    min = sa[j];
		}
		if (sa[j] > max)
		{
		    max = sa[j];
		}
	    }
	    if (max - min > ht)
	    {
		return 1;
	    }
	    j--;
	    i = j;
	}
    }
    return 0;
}

void GetResult(FILE *fout)
{
    int left = 0, right = N - 1, middle, result = 0;

    while (left <= right)
    {
	middle = (left + right) / 2;
	if (Test(middle) == 1)
	{
	    result = middle;
	    left = middle + 1;
	}
	else
	{
	    right = middle - 1;
	}
    }
    if (result + 1 >= 5)
    {
	fprintf(fout, "%d\n", result + 1);
    }
    else
    {
	fprintf(fout, "%d\n", 0);
    }
}

int main(void)
{
    FILE *fin, *fout;

    fin = fopen("theme.in", "r");
    fout = fopen("theme.out", "w");

    GetInput(fin);
    GenerateRank();
    GenerateSa();
    GenerateHeight();
    GetResult(fout);

    return 0;
}





usaco 5.1.3 Musical Themes

/* ID: daniel.20 LANG: JAVA TASK: theme */ import java.util.*; import java.io.*; class problem2{ ...
  • dyhu083
  • dyhu083
  • 2014年06月05日 20:59
  • 264

USACO5.1.3 Musical Themes (theme)

dp[i][j]表示第一个串匹配到了i,第二个串匹配到了j的最长公共部分的长度。 只不过这题不能直接匹配,要计算后一个和前一个的差值再匹配。 dp[i][j]=max(dp[i][j],dp[i-...
  • u011265346
  • u011265346
  • 2014年12月10日 09:36
  • 459

Section 5.1 Musical Themes

/*ID: niepeng1PROG: themeLANG: C++*//*典型的动归题目,用双重循环,第一重循环theme的长度第二重循环起始点。如果两个之间相减相同就将当前长度加一,但是注意要小于...
  • china1000
  • china1000
  • 2010年03月12日 21:42
  • 296

USACO Musical Themes 解题报告

这道题的解法应该是很多的,但是这道题的时间和ko
  • thestoryofsnow
  • thestoryofsnow
  • 2014年09月30日 08:55
  • 548

USACO5.1.3 Musical Themes(theme)

寻找最长的重复序列,枚举两个序列的开头的序列的长度 采用了一个力度很大的优化: 每次枚举末端那一位时都从开始的一位加上当前最大长度开始枚举,前一个序列开头只用枚举到(N-当前的最大长度)即可。 /...
  • qq_25471829
  • qq_25471829
  • 2015年02月04日 17:42
  • 198

usaco5.1.3 Musical Themes

一 原题 Musical Themes Brian Dean A musical melody is represented as a sequence of N (1
  • max_kibble
  • max_kibble
  • 2017年05月04日 11:21
  • 130

USACO Section 5.1 Musical Themes

题意: 找到乐谱中重复出现过的旋律  旋律ru
  • u013351160
  • u013351160
  • 2014年05月28日 16:28
  • 643

C++——【USACO 5.1.3】——Musical Themes

Musical Themes Brian Dean A musical melody is represented as a sequence of N (1
  • McDonnell_Douglas
  • McDonnell_Douglas
  • 2017年10月23日 20:22
  • 34

[USACO 5.1.3]Musical Themes(后缀数组)

【题目描述】: 我们用N(1 许多作曲家围绕一个重复出现的“主题”来构建乐曲。在我们的乐曲表示法中,“主题”是整个音符序列的一个子串,它需要满足如下条件:⒈长度至少为5个音符⒉在乐曲中重复出现(可能经...
  • Skyprophet
  • Skyprophet
  • 2010年03月13日 18:27
  • 910

usaco 5.1 Musical Themes(KMP或DP)

Musical Themes Brian Dean A musical melody is represented as a sequence of N (1
  • fp_hzq
  • fp_hzq
  • 2012年09月21日 17:22
  • 644
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Musical Themes
举报原因:
原因补充:

(最多只允许输入30个字)