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

先写了一个N^3的裸搜试了试~~结果还过了12个点~~囧...然后就想起优化...主要思想是要利用前面已经搜索过的..尽量减少重复的判断搜索...根据这些特征自然而然的想到了KMP...但由于题目中是要求两列每位对应差相等而不是说两列完全相等...所以写了一晚上搞得十分蛋疼..还是搞不出来...

再看了下样例..发现题目可以转化!!! 其实差了一个自然数的两列..若写成相邻前后之差的形式就会变成完全相同的两列..如1 2 4 5 4 5 7 1 写成相邻前后之差的形式为 1 2 1 -1 1 2 1 -6 那么看出最长公共部分是有两个1 2 ...那么这串中符合题目要求最长的主旋律是3个音节( 因为这个1,2是1 2 4 又是 4 5 7 不难理解为啥实际结果比转化后的最长结果会大一.)...

经过这么转化..题目就变成了给出一个串...找出其中最长的两段不相邻的相同子串...我最先想到的KMP就可以完美实现了!!..N^2的算法..先枚举第一个段的起点..从1~N..将这个起点到N的段构造好Fail指针( 同最本质的KMP..指向离自己最近的相同前缀尾部部分 ) .. 构造好后..从第一段起点的后面一个点扫到N...就是用后面一段在先做好的长串里找最长公共子串...但要注意题目的要求..两段不能有公共部分...并且在通过转化后..这个条件又等同了两段间至少要有一个数..不能相邻...注意了这个条件的处理...这个KMP还是很好写的..

Program:

/* ID: zzyzzy12 LANG: C++ TASK: theme */ #include<iostream> #include<istream> #include<stdio.h> #include<string.h> #include<math.h> #include<stack> #include<map> #include<algorithm> #include<queue> #define oo 2000000005 #define ll long long #define pi (atan(2)+atan(0.5))*2 using namespace std; int a[5005],n,Fail[5005]; int main() { freopen("theme.in","r",stdin); freopen("theme.out","w",stdout); int i,j,m,k,ans; bool f; scanf("%d",&n); for (i=1;i<=n;i++) scanf("%d",&a[i]); for (i=1;i<n;i++) a[i]=a[i+1]-a[i]; n--; ans=0; for (i=1;i<=n;i++) { Fail[i]=i-1; for (j=i+1;j<=n;j++) { k=Fail[j-1]; while (k>=i && a[k+1]!=a[j]) k=Fail[k]; if (a[k+1]==a[j]) Fail[j]=k+1; else Fail[j]=k; } m=0; j=i+1; while (j<=n) { f=true; while (i+m<j-m-1 && a[i+m]==a[j]) { j++; m++; f=false; } if (f) { j++; continue; } if (m+1>ans) ans=m+1; k=Fail[i+m-1]; while (k>=i && a[k+1]!=a[j]) k=Fail[k]; if (a[k+1]==a[j]) m=k-i+1; else m=0; } } if (ans<5) ans=0; printf("%d\n",ans); return 0; }