poj1743 Musical Theme

神题！！

2333，要求可以给区间加减同一个数，然后我就懵逼了，，%%题解，可以用差分嘛2333

所以原来的区间长度差分之后就要减一了，

然后现在用sa搞出height数组，然后对于>=len（len为二分的答案长度-1）的区间，去里面sa的最大和最小的差，就是这两个子串相距的最远距离了，判断这个是不是满足>len就好

注意，不能==len，因为这是差分数组，==len的话就相当于重合

#include <cstdio>
#include <iostream>
#include <cstring>
#define N 100005
#define inf 0x3f3f3f3f
#define LL long long
#define eps 1e-8
using namespace std;
inline int ra()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}

int sa[2][N],rank[2][N],v[N],height[N];
int p=0,q=1,n;
int a[N];

void get_height()
{
    int k=0;
    for (int i=1; i<=n; i++)
    {
        if (rank[p][i]==1) height[1]=0;
        else
        {
            int j=sa[p][rank[p][i]-1];
            while (a[i+k]==a[j+k]) k++;
            height[rank[p][i]]=k; k?k--:k;
        }
    }
}
void cal_sa(int sa[N], int rank[N], int Sa[N], int Rank[N], int k)
{
    for (int i=1; i<=n; i++) v[rank[sa[i]]]=i;
    for (int i=n; i>=1; i--) if (sa[i]>k) Sa[v[rank[sa[i]-k]]--]=sa[i]-k;
    for (int i=n-k+1; i<=n; i++) Sa[v[rank[i]]--]=i;
    for (int i=1; i<=n; i++) Rank[Sa[i]]=Rank[Sa[i-1]]+(rank[Sa[i]]!=rank[Sa[i-1]] || rank[Sa[i]+k]!=rank[Sa[i-1]+k]);
}
void work()
{
    for (int i=1; i<=n; i++) v[a[i]]++;
    for (int i=1; i<=200; i++) v[i]+=v[i-1];
    for (int i=1; i<=n; i++) sa[p][v[a[i]]--]=i;
    for (int i=1; i<=n; i++) rank[p][sa[p][i]]=rank[p][sa[p][i-1]]+(a[sa[p][i]]!=a[sa[p][i-1]]);
    for (int k=1; k<n; k<<=1,swap(p,q)) cal_sa(sa[p],rank[p],sa[q],rank[q],k);
    get_height();
}
bool check(int len)
{
    int mn=inf,mx=-inf;
    for (int i=2; i<=n; i++)
    {
        if (height[i]>=len)
        {
            mx=max(mx,max(sa[p][i-1],sa[p][i]));
            mn=min(mn,min(sa[p][i-1],sa[p][i]));
            if (mx-mn>len) return 1;
        }
        else mn=inf,mx=-inf;
    }
    return 0;
}
void solve()
{
    int l=0,r=n>>1,ans=0;
    while (l<=r)
    {
        int mid=l+r>>1;
        if (check(mid)) ans=mid,l=mid+1; else r=mid-1;
    }
    printf("%d
",ans>=4?ans+1:0);
}

int main(int argc, char const *argv[])
{
    while (n=ra())
    {
        int last;
        memset(v,0,sizeof(v));
        for (int i=1; i<=n; i++)
        {
            int x=ra();
            if (i!=1) a[i-1]=x-last+100;
            last=x;
        }
        n--; work();
        solve();
    }
    return 0;
}