Muscial Theme_usaco 5.1_dp||暴力

Description

---

我们用N(1 <= N <=5000)个音符的序列来表示一首乐曲，每个音符都是1..88范围内的整数，每个数表示钢琴上的一个键。很不幸这种表示旋律的方法忽略了音符的时值，但这项编程任务是关于音高的，与时值无关。
许多作曲家围绕一个重复出现的“主题”来构建乐曲。在我们的乐曲表示法中，“主题”是整个音符序列的一个子串，它需要满足如下条件：
⒈长度至少为5个音符
⒉在乐曲中重复出现(可能经过转调，见下)
⒊重复出现的同一主题不能有公共部分。
“转调”的意思是主题序列中每个音符都被加上或减去了同一个整数值。给定一段乐曲，计算其中最长主题的长度(即音符数)。
本题时限为1秒钟!

Input

---

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

Output

---

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

Analysis

---

对于每个序列我们可以转变为他们相邻数字的差的序列，题目就变成了求最长的一段相同的子序列

做法Ⅰ

dp。我们发现这跟最长公共自序列很像，那么参考原dp方程可得
若$t[i]=t[j]$$\ \ \ \ \ f[i][j]=\max(f[i-1][j-1]+1,f[i][j])$
题目有限制不能重叠，那么就是$f[i][j]\leq j-i$
usaco的内存有毒，滚起来才不会爆

做法Ⅱ

这个就比较sb了
枚举a串起点，枚举b串起点，贪心的策略一直往后扩展，也就是$O(n^3)$的暴力
这里要加上优化，如果全部匹配都不能达到之前的最大值，那么剪。还有就是$B_{st}-A_{st}\geq ans$ 的，这样能保证尽可能更新最大值

神清气爽

Code

---

/*
ID:wjp13241
PROG:theme
LANG:C++
*/
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#define fo(i,a,b) for (int i=a;i<=b;i++)
#define dfo(i,a,b) for (int i=a;i>=b;i--)
#define fil(x,t) memset(x,t,sizeof(x))
#define min(x,y) x<y?x:y
#define max(x,y) x>y?x:y
#define ll long long
#define INF 0x7f7f7f7f
#define EPS 1e-4
#define L 351
#define N 5001
#define E 60001*2+1
using namespace std;
int t[N],f[2][N];
int main()
{
    freopen("theme.in","r",stdin);
    freopen("theme.out","w",stdout);
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    fo(i,1,n)
        cin>>t[i];
    int ans=0;
    fo(i,1,n)
        fo(j,i+1,n)
        {
            f[i&1][j]=0;
            if (t[i]-t[i-1]==t[j]-t[j-1]&&f[(i-1)&1][j-1]<j-i-1)
                f[i&1][j]=max(f[(i-1)&1][j-1]+1,f[i&1][j]);
            ans=max(ans,f[i&1][j]);
        }
    if (ans<4)
        ans=-1;
    cout<<ans+1<<endl;
    return 0;
}

/*
ID:wjp13241
PROG:theme
LANG:C++
*/
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#define fo(i,a,b) for (int i=a;i<=b;i++)
#define dfo(i,a,b) for (int i=a;i>=b;i--)
#define fil(x,t) memset(x,t,sizeof(x))
#define min(x,y) x<y?x:y
#define max(x,y) x>y?x:y
#define ll long long
#define INF 0x7f7f7f7f
#define EPS 1e-4
#define L 351
#define N 5001
#define E 60001*2+1
using namespace std;
int t[N],p[N];
int main()
{
    freopen("theme.in","r",stdin);
    freopen("theme.out","w",stdout);
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    fo(i,1,n)
        cin>>p[i];
    fo(i,1,n-1)
        t[i]=p[i+1]-p[i];
    int ans=0;
    fo(i,1,n)
    {
        fo(j,i+ans,n)
        {
            if (j+ans>=n)
                break;
            int a=i,b=j;
            while (t[a]==t[b]&&a<j-1&&b<n)
            {
                ans=max(ans,a-i+1);
                a++,b++;
            }
        }
    }
    if (ans<4)
        ans=-1;
    cout<<ans+1<<endl;
    return 0;
}