WHU---1084 - 连续技 (后缀数组+二分)

Description
不管是什么武功，多少都会有一或两个连续技多次出现，这些连续技常常是发明该武功的人的习惯性动作，如果这些动作被对手分析出来了，就很容易被对手把握住先机。比如松风剑谱里面有一式叫做迎风傲骨是如下的动作：
劈 刺 削 刺 削 踢 刺 削 刺 削
很明显 刺-削 这个连续动作出现了4次，而 刺-削-刺-削 这个连续动作则出现了两次。

现在刘白宇弄到了一本魔教的掌法，想让你帮忙来分析其中最长的且出现尽量多的连续技，当然，他不好意思麻烦你太久，只想让你告诉它这个连续技有多长并且出现了多少次即可。

注意：只有一个动作或者只出现了一次的连续动作都不能被认为是连续技。
Input
输入包含多组数据，每组数据包含两行，第一行是一个数字N(1<=N<=1000)，第二行有N个数字，每个数字代表了掌法里的一个基本动作，如劈、震、拍、点等等。
Output
对于输入的每个掌法，输出其中最长的连续技的长度以及出现的次数。如果不存在这样的连续技，就输出-1。
Sample Input
1
2
10
4 1 2 1 2 5 1 2 1 2
8
1 2 1 2 1 2 1 2
Sample Output
-1
4 2
4 2
Hint
Source
Gardon - DYGG’s contest 4

后缀数组+二分能够求出最长的重复子串,然后对于每一组的后缀,选取第一个或者最后一个后缀(按位置来说的先后关系)一定是最优的,所以遍历一遍就可以知道最多重复次数了

/************************************************************************* > File Name: WHU1084.cpp > Author: ALex > Mail: [email protected] > Created Time: 2015年04月10日 星期五 16时39分10秒 ************************************************************************/

#include <functional>
#include <algorithm>
#include <iostream>
#include <fstream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <queue>
#include <stack>
#include <map>
#include <bitset>
#include <set>
#include <vector>

using namespace std;

const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
const double eps = 1e-15;
typedef long long LL;
typedef pair <int, int> PLL;

int arr[1200];
int xis[1200];
int cnt;

int search(int val)
{
    int l = 1, r = cnt, mid;
    while (l <= r)
    {
        mid = (l + r) >> 1;
        if (xis[mid] == val)
        {
            break;
        }
        else if (xis[mid] > val)
        {
            r = mid - 1;
        }
        else
        {
            l = mid + 1;
        }
    }
    return mid;
}

class SuffixArray
{
    public:
        static const int N = 1200;
        int init[N];
        int X[N];
        int Y[N];
        int Rank[N];
        int sa[N];
        int height[N];
        int buc[N];
        int LOG[N];
        int dp[N][20];
        int size;
        vector <int> task;

        void clear()
        {
            size = 0;
        }

        void insert(int n)
        {
            init[size++] = n;
        }

        bool cmp(int *r, int a, int b, int l)
        {
            return (r[a] == r[b] && r[a + l] == r[b + l]);
        }

        void getsa(int m = 256) //m一般为最大值+1
        {
            init[size] = 0;
            int l, p, *x = X, *y = Y, n = size + 1;
            for (int i = 0; i < m; ++i)
            {
                buc[i] = 0;
            }
            for (int i = 0; i < n; ++i)
            {
                ++buc[x[i] = init[i]];
            }
            for (int i = 1; i < m; ++i)
            {
                buc[i] += buc[i - 1];
            }
            for (int i = n - 1; i >= 0; --i)
            {
                sa[--buc[x[i]]] = i;
            }
            for (l = 1, p = 1; l <= n && p < n; m = p, l *= 2)
            {
                p = 0;
                for (int i = n - l; i < n; ++i)
                {
                    y[p++] = i;
                }
                for (int i = 0; i < n; ++i)
                {
                    if (sa[i] >= l)
                    {
                        y[p++] = sa[i] - l;
                    }
                }
                for (int i = 0; i < m; ++i)
                {
                    buc[i] = 0;
                }
                for (int i = 0; i < n; ++i)
                {
                    ++buc[x[y[i]]];
                }
                for (int i = 1; i < m; ++i)
                {
                    buc[i] += buc[i - 1];
                }
                for (int i = n - 1; i >= 0; --i)
                {
                    sa[--buc[x[y[i]]]] = y[i];
                }
                int i;

                for (swap(x, y), x[sa[0]] = 0, p = 1, i = 1; i < n; ++i)
                { 
                    x[sa[i]] = cmp(y, sa[i - 1], sa[i], l) ? p - 1 : p++; 
                }
            }
        }

        void getheight()
        {
            int h = 0, n = size;
            for (int i = 0; i <= n; ++i)
            {
                Rank[sa[i]] = i;
            }
            height[0] = 0;
            for (int i = 0; i < n; ++i)
            {
                if (h > 0)
                {
                    --h;
                }
                int j =sa[Rank[i] - 1];
                for (; i + h < n && j + h < n && init[i + h] == init[j + h]; ++h);
                height[Rank[i] - 1] = h;
            }
        }   

        bool check(int k)
        {
            int l = sa[1], r = sa[1];
            for (int i = 1; i < size; ++i)
            {
                if (height[i] >= k)
                {
                    l = min(l, sa[i + 1]);
                    r = max(r, sa[i + 1]);
                }
                else
                {
                    if (r - l >= k)
                    {
                        return 1;
                    }
                    l = r = sa[i + 1];
                }
            }
            return 0;
        }

        void solve()
        {
            int l = 1, r = size / 2, mid;
            int ans = -1;
            while (l <= r)
            {
                mid = (l + r) >> 1;
                if (check(mid))
                {
                    l = mid + 1;
                    ans = mid;
                }
                else
                {
                    r = mid - 1;
                }
            }
            if (ans == -1)
            {
                printf("-1\n");
                return;
            }
            task.clear();
            task.push_back(sa[1]);
            int ret = 0;
            for (int i = 1; i < size; ++i)
            {
                if (height[i] >= ans)
                {
                    task.push_back(sa[i + 1]);
                }
                else
                {
                    sort(task.begin(), task.end());
                    int tmp1 = 1, tmp2 = 1, isize = task.size(), last = task[0];
                    for (int j = 1; j < isize; ++j)
                    {
                        if (task[j] - last >= ans)
                        {
                            ++tmp1;
                            last = task[j];
                        }
                    }
                    last = task[isize - 1];
                    for (int j = isize - 2; j >= 0; --j)
                    {
                        if (last - task[j] >= ans)
                        {
                            ++tmp2;
                            last = task[j];
                        }
                    }
                    ret = max(ret, max(tmp1, tmp2));
                    task.clear();
                    task.push_back(sa[i + 1]);
                }
            }
            printf("%d %d\n", ans, ret);

        }
}SA;

int main()
{
    int n;
    while (~scanf("%d", &n))
    {
        SA.clear();
        cnt = 0;
        for (int i = 0; i < n; ++i)
        {
            scanf("%d", &arr[i]);
            xis[++cnt] = arr[i];
        }
        sort (xis + 1, xis + 1 + cnt);
        cnt = unique(xis + 1, xis + 1 + cnt) - xis - 1;
        for (int i = 0; i < n; ++i)
        {
            int val = search(arr[i]);
            SA.insert(val);
        }
        SA.getsa(cnt + 1);
        SA.getheight();
        SA.solve();
    }
    return 0;
}