Topic Links:
Poj 1743 Musical Theme
Title Description:
Give a string of numbers (the number range is in [1,88]), to find a topic in this string of numbers that satisfies:
1: Topic length is greater than or equal to 5.
2: The theme repeats in the text string (or, after a twist, the subject adds or subtracts the same integer)
3: Duplicate themes cannot overlap
Problem Solving Ideas:
To reverse the repeated occurrences of the substring, the difference between the themes must be constant. You can find the difference between two adjacent numbers in the text string s, re-form a new text string s, and then look for the longest public non-overlapping prefix in the s suffix string. Rank adjacent suffix string, the public prefix must be the longest, but it is possible to overlap. We can divide the length of the subject by K, and then verify whether K is set. According to the nature of height, the more similar suffix string rank difference is smaller, then we can find I=min (Sa[rank]) in the adjacent interval of Height[rank]>=k, J=max (Sa[rank]), if j-i>= K. K is established.
1#include <cstdio>2#include <cstring>3#include <iostream>4#include <algorithm>5 using namespacestd;6typedefLong LongLL;7 Const intMAXN =20010;8 9 intSA[MAXN], RANK[MAXN], HEIGHT[MAXN];Ten intT1[MAXN], T2[MAXN], R[MAXN], C[MAXN]; One BOOLCMP (int*STR,intAintBintk) A { - returnSTR[A]==STR[B] && str[a+k]==str[b+K]; - } the voidDa (int*STR,intNintm) - { -n + +; - int*x = T1, *y =T2, I, J; + for(i=0; i<m; i++) C[i] =0; - for(i=0; i<n; i++) C[x[i]=str[i] + +; + for(i=1; i<m; i++) C[i] + = c[i-1]; A for(i=n-1; i>=0; i--) sa[--C[x[i]] =i; at for(j=1; j<=n; j*=2) - { - intp =0; - for(i=n-j; i<n; i++) y[p++] =i; - for(i=0; i<n; i++)if(Sa[i] >= j) y[p++] = Sa[i]-J; - in for(i=0; i<m; i++) C[i] =0; - for(i=0; i<n; i++) C[x[y[i]] + +; to for(i=1; i<m; i++) C[i] + = c[i-1]; + for(i=n-1; i>=0; i--) sa[--C[x[y[i]] =Y[i]; - the swap (x, y); *p =1; $x[sa[0]] =0;Panax Notoginseng for(intI=1; i<n; i++)//i is rank -X[sa[i]] = cmp (y, sa[i-1], Sa[i], j)? p1:p + +; the if(P >=N) + Break; Am =p; the } + for(i=1; i<n; i++) -Rank[sa[i]] =i; $ intK =0; $N--; - for(intI=0; i<n; i++) - { the if(k) K--; - intj = Sa[rank[i]-1];Wuyi while(Str[i+k] = = Str[j+k]) k++; theHeight[rank[i]] =K; - } Wu } - BOOLSolve (intXintN) About { $ intMa, MI; -Ma = mi = sa[1]; - for(intI=2; i<=n; i++) - { A if(Height[i]>=x && i<=N) + { theMi =min (mi, sa[i]); -Ma =Max (MA, sa[i]); $ if(Ma-mi >= x)return true; the Continue; the } the theMa = mi =Sa[i]; - } in return false; the } the intMain () About { the intN; the while(SCANF ("%d", &N), N) the { + intS, E; -scanf ("%d", &s); the for(intI=1; i<n; i++)Bayi { thescanf ("%d", &e); ther[i-1] = S-e + -; -s =e; - } theR[--n] =0; the theDa (r, N, $); the intAns =0, high = n/2, low =1, mid; - while(Low <=High ) the { theMid = (low + high)/2; the if(Solve (Mid, n))94 { theAns =mid; theLow = mid +1; the }98 Else AboutHigh = mid-1; - }101printf ("%d\n", ans<4?0: ans+1);102 }103 return 0;104}
Poj 1743 musical Theme (suffix array + dichotomy)