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 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
  • 628

【POJ1743】Musical Themes 乐曲主题 后缀数组、 (也可以用hash+二分做的~)

题意: 1829: Musical Themes 乐曲主题 Time Limit: 1 Sec  Memory Limit: 128 MB Submit: 42  Solved: 1...
  • Vmurder
  • Vmurder
  • 2015年01月04日 20:35
  • 851

USACO Section 5.1 Musical Themes - 题目转换以及KMP..

先写了一个N^3的裸搜试了试~~结果还过了12个点~~囧...然后就想起优化...主要思想是要利用前面已经搜索过的..尽量减少重复的判断搜索...根据这些特征自然而然的想到了KMP...但由于题目中是...
  • kk303
  • kk303
  • 2012年01月26日 04:07
  • 737

Fourier Analysis of Musical Tones

  • 2010年12月12日 15:16
  • 859KB
  • 下载

Android中的样式和主题(Styles and Themes)

样式(style),是为一个视图(View)或窗体(window)指定形状和格式的属性集。一个样式能够指定的属性包括高度、内边距、字体颜色、字体大小、背景颜色等等。样式被定义在一个XML资源文件中,跟...

风格和主题(style,themes)

当你设计你的程序的时候,你可以用风格和主题来统一格式化各种屏幕和UI元素。     * 风格(style)是一个包含一种或者多种格式化属性值的集合,你可以将其用为一个单位用在布局XML单个元素当中。...

如何正确的使用 Android 中的 themes 和 styles

如何正确的使用 Android 中的 themes 和 styles 译自:How to avoid headaches with Android themes and styles Androi...

Android应用开发中的风格和主题(style,themes)

越来越多互联网企业都在Android平台上部署其客户端,为了提升用户体验,这些客户端都做得布局合理而且美观.......Android的Style设计就是提升用户体验的关键之一。Android上的St...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Musical Themes
举报原因:
原因补充:

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