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;
}





版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

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

题意: 1829: Musical Themes 乐曲主题 Time Limit: 1 Sec  Memory Limit: 128 MB Submit: 42  Solved: 1...

usaco 5.1 Musical Themes(KMP或DP)

Musical Themes Brian Dean A musical melody is represented as a sequence of N (1 <= N <= 5000) note...

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

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

Fourier Analysis of Musical Tones

  • 2010-12-12 15:16
  • 859KB
  • 下载

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

From: http://www.cnblogs.com/playing/archive/2011/04/01/2002469.html       越来越多互联网企业都在Android平台上部署其...

Android API Guides---Styles and Themes

Styles and Themes 样式是指定视图或窗口的外观和格式属性的集合。样式可以指定诸如高度,填充,字体颜色,字体大小,背景颜色,以及更多的特性。样式是在从指定布局XML单独的XML资源定义...

ASP.NET2.0中Web应用程序样式切换(App_Themes)

主题是Microsoft ASP.NET 2.0 提供给开发人员的一项新增功能。          在 Asp.NET 2.0 出现以前, 给网站提供更换主题功能是一项繁重的工作,鉴于每切换一套...

【Google官方译文】Styles and Themes

说明:本文转转载于recordGrowth的翻译 http://blog.csdn.net/recordgrowth/article/details/52310974 本文为Google官方译文...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

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