ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

「JSOI2015」串分割

2020-02-25 23:54:47  阅读:303  来源: 互联网

标签:分割 frac JSOI2015 int mid read rg sa


「JSOI2015」串分割

传送门

首先我们会有一个贪心的想法:分得越均匀越好,因为长的绝对比短的大。

那么对于最均匀的情况,也就是 \(k | n\) 的情况,我们肯定是通过枚举第一次分割的位置,然后每一段长度 \(\frac{n}{k}\) 最后取最小的。

把这个思想运用到一般情况:如果分出来两段长短不一,那么长的只会比短的那个长度多 \(1\) ,再仔细想想,所有段只会有两种不同的长度 \(\lfloor \frac{n}{k} \rfloor, \lceil \frac{n}{k} \rceil\) ,那么我们就只需要判断什么情况下取长的,什么时候取短的就好了。

于是考虑二分最大的那个段的大小,然后再枚举第一次分割的位置,对于当前我们这段,如果分成长段会超过二分值,就换成短的,不然就按长的分,是否符合条件就判一下分出来的段的长度之和的大小,如果大于 \(n\) 说明当前的二分值大了,否则就是小了。

至于一段数的大小,我们可以直接将原数字环变成链然后倍长跑一遍后缀排序,比较的时候就直接拿 \(rank_i\) 来比就是了。

参考代码:

#include <algorithm>
#include <cstdio>
#define rg register
#define file(x) freopen(x".in", "r", stdin), freopen(x".out", "w", stdout)
using namespace std;
template < class T > inline void read(T& s) {
    s = 0; int f = 0; char c = getchar();
    while ('0' > c || c > '9') f |= c == '-', c = getchar();
    while ('0' <= c && c <= '9') s = s * 10 + c - 48, c = getchar();
    s = f ? -s : s;
}

const int _ = 4e5 + 5;

int n, k, a[_], sa[_], rk[_];

inline void sufsort(int n) {
#define cmp(i, j, k) (y[i] == y[j] && y[i + k] == y[j + k])
    static int x[_], y[_], cnt[_];
    int m = 57;
    for (rg int i = 1; i <= n; ++i) ++cnt[x[i] = a[i]];
    for (rg int i = 1; i <= m; ++i) cnt[i] += cnt[i - 1];
    for (rg int i = n; i >= 1; --i) sa[cnt[x[i]]--] = i;
    for (rg int k = 1; k <= n; k <<= 1) {
        int p = 0;
        for (rg int i = n - k + 1; i <= n; ++i) y[++p] = i;
        for (rg int i = 1; i <= n; ++i) if (sa[i] > k) y[++p] = sa[i] - k;
        for (rg int i = 1; i <= m; ++i) cnt[i] = 0;
        for (rg int i = 1; i <= n; ++i) ++cnt[x[i]];
        for (rg int i = 1; i <= m; ++i) cnt[i] += cnt[i - 1];
        for (rg int i = n; i >= 1; --i) sa[cnt[x[y[i]]]--] = y[i];
        swap(x, y), x[sa[1]] = p = 1;
        for (rg int i = 2; i <= n; ++i) x[sa[i]] = cmp(sa[i], sa[i - 1], k) ? p : ++p;
        if (p == n) break ; else m = p;
    }
    for (rg int i = 1; i <= n; ++i) rk[sa[i]] = i;
}

inline bool check(int mid) {
    int x = n / k + (n % k != 0);
    for (rg int i = 1; i <= x; ++i)
        for (rg int p = i, j = 1; j <= k; ++j) {
            p += rk[p] <= mid ? x : x - 1;
            if (p >= i + n) return 1;
        }
    return 0;
}

int main() {
#ifndef ONLINE_JUDGE
    file("cpp");
#endif
    read(n), read(k);
    for (rg int i = 1; i <= n; ++i) scanf("%1d", a + i), a[n + i] = a[i];
    sufsort(n << 1);
    int l = 1, r = n << 1;
    while (l < r) {
        int mid = (l + r) >> 1;
        if (check(mid)) r = mid; else l = mid + 1;
    }
    for (rg int i = 1; i <= n; ++i)
        if (rk[i] == l) {
            for (rg int j = 1; j <= n / k + (n % k != 0); ++j) putchar(a[i + j - 1] + '0'); break ;
        }
    return 0;
}

标签:分割,frac,JSOI2015,int,mid,read,rg,sa
来源: https://www.cnblogs.com/zsbzsb/p/12364662.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有