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

该博客介绍了如何利用后缀数组和二分查找算法来分析连续技,通过实例展示了如何找出武功掌法中最长且频繁出现的连续动作序列,并给出其出现次数。对于输入的多个掌法序列,程序会输出最长连续技的长度及其出现的次数,当不存在符合条件的连续技时,输出-1。
摘要由CSDN通过智能技术生成

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: zchao1995@gmail.com 
    > 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;
}
whu-river数据集是一个用于河流水质监测和预测的数据集。该数据集由华中科技大学收集和创建,旨在帮助研究人员和决策者更好地理解和管理河流的水质问题。 该数据集包含了来自中国不同河流的水质监测数据,包括河水的化学指标、生物学指标、水文学指标和气象指标等。这些指标可以帮助评估河流的污染程度和水质状况,并提供对河流水质变化的预测。数据集中的每个样本都包含了特定时间点对应的水质指标数值。 使用whu-river数据集,研究人员可以进行水质分析、趋势预测和模型构建等工作。例如,可以通过分析数据集中的化学指标,了解河流中各种物质的含量和分布情况,进而评估河流的水质状况。同时,还可以利用数据集中的生物学指标,评估河流的生态系统健康度。 此外,whu-river数据集还可以用于建立水质预测模型。通过对历史数据的分析和挖掘,可以发现不同指标之间的关联性和规律,从而构建出准确预测水质的模型。这样的模型可以对未来水质状况进行预测,为决策者提供科学依据,制定更为有效的水资源管理和保护措施。 总之,whu-river数据集是一个宝贵的资源,对于研究河流水质和水资源管理具有重要意义。它提供了丰富的水质监测数据,可用于水质分析、趋势预测和模型构建等工作,为保护水环境和维护水质做出贡献。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值