SP687 Repeats

题意

给定字符串,求重复次数最多的连续重复子串

思路

有两个方法

方法1

直接枚举单个循环节的长度 k i ki ki。

然后将字符串划分成 n / k i n/ki n/ki个区间,这样如果只用枚举每一个区间就可以做到 O ( n l o g n ) O(nlogn) O(nlogn)。

考虑到一个重复子串如果循环了 t t t次,那么有至少 t − 1 t-1 t−1个区间是被完全覆盖的。我们只需要用哈希判断相邻的两个区间是不是一样的。

假设现在连续 d d d个区间是一样的,那么我们现在的目标就是确定这个子串到底是循环了 d d d还是 d + 1 d+1 d+1。因为我们已经可以确定循环节,因此左右两边需要添加的字符都是固定的,我们只用再在相邻的两个区间上二分+哈希,就能求出这一段最多扩展的长度,然后就可以判断了。

cpp 复制代码
#include <bits/stdc++.h>
#define ll long long
using namespace std;

const int N = 5e4 + 10;
const ll mod = 998244353;

int n; char s[N];
ll val[N], hsh[N];

ll gethsh(int l, int r) {
    return (hsh[r] - hsh[l - 1] * val[r - l + 1] % mod + mod) % mod;
}

void solve() {
    cin >> n;
    char ss[2];
    val[0] = 1;
    for (int i = 1; i <= n; i++) {
        cin >> ss;
        val[i] = val[i - 1] * 357 % mod;
        hsh[i] = (hsh[i - 1] * 357 + ss[0]) % mod;
    }
    int ans = 1;
    // int ki = 3, l = 5, r;
    for (int ki = 1; ki <= n / 2; ki++) {
        for (int l = 1, r; l <= n; l = r + 1) {
            r = l + ki - 1;
            ll flag = gethsh(l, r);
            while (r + ki <= n && gethsh(r + 1, r + ki) == flag) {
                r = r + ki;
            }
            int zuo, you;
            int L = max(1, l - ki), R = l - 1, mid, pos = l;
            while (L <= R) {
                mid = (L + R) / 2;
                if (gethsh(mid, l - 1) == gethsh(mid + ki, l - 1 + ki)) {
                    pos = mid;
                    R = mid - 1;
                }
                else {
                    L = mid + 1;
                }
            }
            zuo = pos;
            L = r + 1, R = min(n, r + ki), pos = r;
            while (L <= R) {
                mid = (L + R) / 2;
                if (gethsh(r + 1, mid) == gethsh(r + 1 - ki, mid - ki)) {
                    pos = mid;
                    L = mid + 1;
                }
                else {
                    R = mid - 1;
                }
            }
            you = pos;
            ans = max(ans, (you - zuo + 1) / ki);
        }
    }
    cout << ans << endl;
}

int main() {
    // freopen("in.in", "r", stdin);
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    int T; cin >> T;
    while (T--) {
        solve();
    }
}

方法二

我采用的是方法一,还没有具体实现过方法二。方法二主要是使用了kmp求循环节的思想。

考虑一个循环子串abaabaaba。

我们可以把它划分成两个相同的子串 a b a a b a abaaba abaaba,长度 l e n = 6 len=6 len=6,而着两个串的endpos之差为 d = 3 d=3 d=3。

如果 d ∣ l e n d|len d∣len,那么就可以判定是由长度为 d d d的循环节构成。

对于sam上面的一个节点,我们先要找出最小的d,用线段树合并来维护。

然后对于在这个节点所代表的字符串长度区间里面找到最大的 d ∣ l e n d|len d∣len,然后就能统计答案了。

相关推荐
Tony沈哲1 分钟前
macOS 上为 Compose Desktop 构建跨架构图像处理 dylib:OpenCV + libraw + libheif 实践指南
opencv·算法
刘海东刘海东33 分钟前
结构型智能科技的关键可行性——信息型智能向结构型智能的转变(修改提纲)
人工智能·算法·机器学习
pumpkin845141 小时前
Rust 调用 C 函数的 FFI
c语言·算法·rust
挺菜的1 小时前
【算法刷题记录(简单题)003】统计大写字母个数(java代码实现)
java·数据结构·算法
mit6.8241 小时前
7.6 优先队列| dijkstra | hash | rust
算法
2401_858286112 小时前
125.【C语言】数据结构之归并排序递归解法
c语言·开发语言·数据结构·算法·排序算法·归并排序
guygg882 小时前
基于matlab的FIR滤波器
开发语言·算法·matlab
双叶8363 小时前
(C++)学生管理系统(正式版)(map数组的应用)(string应用)(引用)(文件储存的应用)(C++教学)(C++项目)
c语言·开发语言·数据结构·c++
ysh98883 小时前
PP-OCR:一款实用的超轻量级OCR系统
算法
遇雪长安3 小时前
差分定位技术:原理、分类与应用场景
算法·分类·数据挖掘·rtk·差分定位