poj 1743

题意：给出一个串，求两个不相交的长度相等的子串，使得对应位置的差相等。

题解：首先，假设串是a[0],a[1],...a[n]，先作差，即a[1]=a[1]-a[0],a[2]=a[2]-a[1],a[i]=a[i]-a[i-1]，然后我们求a[1],a[2],a[3],...a[n]的两个不相交的公共子串。

先将h数组求出来，然后二分子串长度，。。。。

#include <cstdio>
#include <cstring>
#include <iostream>
#define maxn 20010
using namespace std;

int n;
int aa[maxn], sa[maxn], rk[maxn], ht[maxn], vv[maxn];

void expand( int k, int sa[maxn], int rk[maxn], int tsa[maxn], int trk[maxn] ) {
    for( int i=1; i<=n; i++ ) vv[rk[sa[i]]]=i;
    for( int i=n; i>=1; i-- ) if( sa[i]>k ) tsa[vv[rk[sa[i]-k]]--]=sa[i]-k;
    for( int i=n-k+1; i<=n; i++ ) tsa[vv[rk[i]]--]=i;
    for( int i=1; i<=n; i++ ) trk[tsa[i]]=trk[tsa[i-1]]+(rk[tsa[i]]!=rk[tsa[i-1]]||rk[tsa[i]+k]!=rk[tsa[i-1]+k]);
}
void calcht() {
    int k=0;
    ht[sa[1]]=0;
    for( int i=1; i<=n; i++ ) {
        if( rk[i]==1 ) continue;
        int j=sa[rk[i]-1];
        while( aa[i+k]==aa[j+k] ) k++;
        ht[i] = k;
        if( k>0 ) k--;
    }
}
void suffix() {
    static int tsa[maxn], trk[maxn];
    for( int i=1; i<=180; i++ ) vv[i]=0;
    for( int i=1; i<=n; i++ ) vv[aa[i]]++;
    for( int i=1; i<=180; i++ ) vv[i]+=vv[i-1];
    for( int i=1; i<=n; i++ ) sa[vv[aa[i]]--]=i;
    for( int i=1; i<=n; i++ ) rk[sa[i]]=rk[sa[i-1]]+(aa[sa[i]]!=aa[sa[i-1]]);
    int *pa=sa, *pb=rk, *pc=tsa, *pd=trk;
    for( int k=1; k<n; k<<=1,swap(pa,pc),swap(pb,pd) )
        expand(k,pa,pb,pc,pd);
    if( pa!=sa ) memcpy( sa, tsa, sizeof(sa) ), memcpy( rk, trk, sizeof(rk) );
    calcht();
}
bool ok( int len ) {    //    len>=1
    int top, minp, maxp;
    top = minp = maxp = sa[1];
    for( int i=2; i<=n; i++ ) {
        if( ht[sa[i]]<len ) {
            if( maxp-minp>=len ) return true;
            top = minp = maxp = sa[i];
        } else {
            if( minp>sa[i] ) minp=sa[i];
            if( maxp<sa[i] ) maxp=sa[i];
        }
    }
    if( maxp-minp>=len ) return true;
    return false;
}
int main() {
    while( scanf("%d",&n)==1 ) {
        if( n==0 ) break;
        n--;
        for( int i=0; i<=n; i++ )
            scanf( "%d", aa+i );
        for( int i=n; i>=1; i-- )
            aa[i]=aa[i]-aa[i-1]+88;
        suffix();
        int lf=0, rg=n/2;
        while( lf<rg ) {
            int mid=(lf+rg+1)>>1;
            if( ok(mid) ) lf=mid;
            else rg=mid-1;
        }
        lf++;
        printf( "%d
", lf<5 ? 0 : lf  );
    }
}