RollingHash
(string/RollingHash.hpp)

• View this file on GitHub
• Last update: 2023-05-17 22:41:14+09:00
• Include: #include "string/RollingHash.hpp"

説明

以下、文字列の長さを $N$ とする。

コンストラクタ

roriha (string s)

s から $\bmod (2^{61}-1)$ のローリングハッシュを構築します。計算量は $\Theta(N)$ です。

get

ull get(int l, int r)

0_indexed で s[l,r) のハッシュを取得します。計算量は $\Theta(1)$ です。

get_hash

ull get_hash(string s)

s 全体のハッシュを取得します。ローリングハッシュの構築を行う必要がない場合に使います。計算量は $\Theta(N)$ です。

concat

ull concat(ull hash1, ull hash2, int len2)

string s, t について s+t のハッシュを取得します。hash1 は s のハッシュ、hash2 と len2 はそれぞれ t のハッシュおよび長さです。計算量は $\Theta(1)$ です。

使い方

int main(){
    string s; cin >> s;
    roriha hs(s);
    int l, r; cin >> l >> r;
    ull lrhash = hs.get(l,r); // [l,r), 0_indexed
    int l1, r1, l2, r2; cin >> l1 >> r1 >> l2 >> r2;
    ull hash1 = hs.get(l1,r1), hash2 = hs.get(l2,r2);
    ull hash12 = roriha::concat(hash1,hash2,r2-l2); // s[l1,r1)+s[l2,r2) のハッシュを取得する
}

Depends on

• テンプレート (template/template.hpp)

Verified with

• test/string/aoj_2444.test.cpp

Code

#pragma once

#include "../template/template.hpp"

namespace lib {

using namespace std;

struct RollingHash {
    using ull = unsigned long long;
    RollingHash(const string& s = "") {
        build(s);
    }
    ull get(int l, int r) {
        assert(0 <= l && l <= r && r <= n);
        return cal(inner_hash[r] + buf - mul(inner_hash[l], pow_base[r - l]));
    }
    static ull get_hash(const string& s) {
        int len = s.size();
        set_hash();
        extend_pow_base(len);
        ull res = 0;
        rep(i, 0, len) res = cal(mul(res, BASE) + s[i]);
        return res;
    }
    static ull concat(const ull& hash1, const ull& hash2, const int& len2) {
        return cal(cal(mul(hash1, pow_base[len2])) + hash2);
    }

  private:
    static constexpr ull MASK30 = (1ULL << 30) - 1;
    static constexpr ull MASK31 = (1ULL << 31) - 1;
    static constexpr ull MASK61 = (1ULL << 61) - 1;
    static constexpr ull MOD = (1ULL << 61) - 1;
    static constexpr ull buf = MOD * 4;
    static ull BASE;
    static vector<ull> pow_base;
    static ull mul(ull a, ull b) {
        ull au = a >> 31, ad = a & MASK31;
        ull bu = b >> 31, bd = b & MASK31;
        ull mid = ad * bu + au * bd;
        ull midu = mid >> 30, midd = mid & MASK30;
        return (au * bu * 2 + midu + (midd << 31) + ad * bd);
    }
    static ull cal(ull x) {
        ull xu = x >> 61;
        ull xd = x & MASK61;
        ull res = xu + xd;
        if (res >= MOD) res -= MOD;
        return res;
    }
    static void set_hash() {
        if (BASE == 0) BASE = (1UL << 31) + (random_device()() & MASK31);
    }
    static void extend_pow_base(int len) {
        int nlen = pow_base.size();
        if (nlen > len) return;
        pow_base.resize(len + 1);
        rep(i, nlen, len + 1) pow_base[i] = cal(mul(pow_base[i - 1], BASE));
    }
    string str;
    int n;
    vector<ull> inner_hash;
    void build(const string& s) {
        set_hash();
        str = s;
        n = s.size();
        extend_pow_base(n);
        inner_hash.resize(n + 1);
        inner_hash[0] = 0;
        rep(i, 0, n) inner_hash[i + 1] = cal(mul(inner_hash[i], BASE) + s[i]);
    }
};
ull RollingHash::BASE = 0;
vector<ull> RollingHash::pow_base = vector<ull>(1, 1);
using roriha = RollingHash;

}  // namespace lib

#line 2 "string/RollingHash.hpp"

#line 2 "template/template.hpp"

#include <bits/stdc++.h>

#define rep(i, s, n) for (int i = (int)(s); i < (int)(n); i++)
#define rrep(i, s, n) for (int i = (int)(n)-1; i >= (int)(s); i--)
#define all(v) v.begin(), v.end()

using ll = long long;
using ld = long double;
using ull = unsigned long long;

template <typename T> bool chmin(T &a, const T &b) {
    if (a <= b) return false;
    a = b;
    return true;
}
template <typename T> bool chmax(T &a, const T &b) {
    if (a >= b) return false;
    a = b;
    return true;
}

namespace lib {

using namespace std;

}  // namespace lib

// using namespace lib;
#line 4 "string/RollingHash.hpp"

namespace lib {

using namespace std;

struct RollingHash {
    using ull = unsigned long long;
    RollingHash(const string& s = "") {
        build(s);
    }
    ull get(int l, int r) {
        assert(0 <= l && l <= r && r <= n);
        return cal(inner_hash[r] + buf - mul(inner_hash[l], pow_base[r - l]));
    }
    static ull get_hash(const string& s) {
        int len = s.size();
        set_hash();
        extend_pow_base(len);
        ull res = 0;
        rep(i, 0, len) res = cal(mul(res, BASE) + s[i]);
        return res;
    }
    static ull concat(const ull& hash1, const ull& hash2, const int& len2) {
        return cal(cal(mul(hash1, pow_base[len2])) + hash2);
    }

  private:
    static constexpr ull MASK30 = (1ULL << 30) - 1;
    static constexpr ull MASK31 = (1ULL << 31) - 1;
    static constexpr ull MASK61 = (1ULL << 61) - 1;
    static constexpr ull MOD = (1ULL << 61) - 1;
    static constexpr ull buf = MOD * 4;
    static ull BASE;
    static vector<ull> pow_base;
    static ull mul(ull a, ull b) {
        ull au = a >> 31, ad = a & MASK31;
        ull bu = b >> 31, bd = b & MASK31;
        ull mid = ad * bu + au * bd;
        ull midu = mid >> 30, midd = mid & MASK30;
        return (au * bu * 2 + midu + (midd << 31) + ad * bd);
    }
    static ull cal(ull x) {
        ull xu = x >> 61;
        ull xd = x & MASK61;
        ull res = xu + xd;
        if (res >= MOD) res -= MOD;
        return res;
    }
    static void set_hash() {
        if (BASE == 0) BASE = (1UL << 31) + (random_device()() & MASK31);
    }
    static void extend_pow_base(int len) {
        int nlen = pow_base.size();
        if (nlen > len) return;
        pow_base.resize(len + 1);
        rep(i, nlen, len + 1) pow_base[i] = cal(mul(pow_base[i - 1], BASE));
    }
    string str;
    int n;
    vector<ull> inner_hash;
    void build(const string& s) {
        set_hash();
        str = s;
        n = s.size();
        extend_pow_base(n);
        inner_hash.resize(n + 1);
        inner_hash[0] = 0;
        rep(i, 0, n) inner_hash[i + 1] = cal(mul(inner_hash[i], BASE) + s[i]);
    }
};
ull RollingHash::BASE = 0;
vector<ull> RollingHash::pow_base = vector<ull>(1, 1);
using roriha = RollingHash;

}  // namespace lib

Back to top page