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数组求出来，然后二分子串长度，。。。。

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 #define maxn 20010
 5 using namespace std;
 6
 7 int n;
 8 int aa[maxn], sa[maxn], rk[maxn], ht[maxn], vv[maxn];
 9
10 void expand( int k, int sa[maxn], int rk[maxn], int tsa[maxn], int trk[maxn] ) {
11     for( int i=1; i<=n; i++ ) vv[rk[sa[i]]]=i;
12     for( int i=n; i>=1; i-- ) if( sa[i]>k ) tsa[vv[rk[sa[i]-k]]--]=sa[i]-k;
13     for( int i=n-k+1; i<=n; i++ ) tsa[vv[rk[i]]--]=i;
14     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]);
15 }
16 void calcht() {
17     int k=0;
18     ht[sa[1]]=0;
19     for( int i=1; i<=n; i++ ) {
20         if( rk[i]==1 ) continue;
21         int j=sa[rk[i]-1];
22         while( aa[i+k]==aa[j+k] ) k++;
23         ht[i] = k;
24         if( k>0 ) k--;
25     }
26 }
27 void suffix() {
28     static int tsa[maxn], trk[maxn];
29     for( int i=1; i<=180; i++ ) vv[i]=0;
30     for( int i=1; i<=n; i++ ) vv[aa[i]]++;
31     for( int i=1; i<=180; i++ ) vv[i]+=vv[i-1];
32     for( int i=1; i<=n; i++ ) sa[vv[aa[i]]--]=i;
33     for( int i=1; i<=n; i++ ) rk[sa[i]]=rk[sa[i-1]]+(aa[sa[i]]!=aa[sa[i-1]]);
34     int *pa=sa, *pb=rk, *pc=tsa, *pd=trk;
35     for( int k=1; k<n; k<<=1,swap(pa,pc),swap(pb,pd) )
36         expand(k,pa,pb,pc,pd);
37     if( pa!=sa ) memcpy( sa, tsa, sizeof(sa) ), memcpy( rk, trk, sizeof(rk) );
38     calcht();
39 }
40 bool ok( int len ) {    //    len>=1
41     int top, minp, maxp;
42     top = minp = maxp = sa[1];
43     for( int i=2; i<=n; i++ ) {
44         if( ht[sa[i]]<len ) {
45             if( maxp-minp>=len ) return true;
46             top = minp = maxp = sa[i];
47         } else {
48             if( minp>sa[i] ) minp=sa[i];
49             if( maxp<sa[i] ) maxp=sa[i];
50         }
51     }
52     if( maxp-minp>=len ) return true;
53     return false;
54 }
55 int main() {
56     while( scanf("%d",&n)==1 ) {
57         if( n==0 ) break;
58         n--;
59         for( int i=0; i<=n; i++ )
60             scanf( "%d", aa+i );
61         for( int i=n; i>=1; i-- )
62             aa[i]=aa[i]-aa[i-1]+88;
63         suffix();
64         int lf=0, rg=n/2;
65         while( lf<rg ) {
66             int mid=(lf+rg+1)>>1;
67             if( ok(mid) ) lf=mid;
68             else rg=mid-1;
69         }
70         lf++;
71         printf( "%d\n", lf<5 ? 0 : lf  );
72     }
73 }