Musical Themes

题意：转载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;
}