結果

問題 No.515 典型LCP
ユーザー kimiyuki
提出日時 2017-05-09 23:22:19
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 2,547 bytes
コンパイル時間 1,322 ms
コンパイル使用メモリ 109,720 KB
実行使用メモリ 155,392 KB
最終ジャッジ日時 2024-09-14 20:07:24
合計ジャッジ時間 8,840 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other RE * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
#include <numeric>
#include <array>
#include <set>
#include <map>
#include <queue>
#include <tuple>
#include <unordered_set>
#include <unordered_map>
#include <functional>
#include <cassert>
#define repeat(i,n) for (int i = 0; (i) < int(n); ++(i))
#define whole(f,x,...) ([&](decltype((x)) whole) { return (f)(begin(whole), end(whole), ## __VA_ARGS__); })(x)
using ll = long long;
using namespace std;
int compute_lcp(string const & a, string const & b) {
    int i = 0;
    while (i < a.length() and i < b.length() and a[i] == b[i]) ++ i;
    return i;
}
template <class Monoid>
struct sparse_table {
    typedef typename Monoid::type T;
    vector<vector<T> > table;
    Monoid mon;
    sparse_table(vector<T> const & init, Monoid const & a_mon = Monoid())
            : mon(a_mon) {
        int n = init.size();
        int log_n = sqrt(n) + 1;
        table.resize(log_n, vector<T>(n, mon.unit));
        table[0] = init;
        for (int k = 0; k < log_n-1; ++ k) {
            for (int i = 0; i < n; ++ i) {
                table[k+1][i] = mon.append(table[k][i], i + (1<<k) < n ? table[k][i + (1<<k)] : mon.unit);
            }
        }
    }
    T operator () (int l, int r) {
        assert (0 <= l and l < r and r <= table[0].size());
        int k = log2(r - l);
        return mon.append(table[k][l], table[k][r - (1<<k)]);
    }
};
struct min_t {
    typedef int type;
    const int unit = 1e9+7;
    int append(int a, int b) { return min(a, b); }
};
int main() {
    // input
    int n; cin >> n;
    vector<string> s(n); repeat (i,n) cin >> s[i];
    int m; ll x, d; cin >> m >> x >> d;
    // generate queries
    vector<int> qi(m);
    vector<int> qj(m);
    repeat (k,m) {
        qi[k] = x / (n - 1);
        qj[k] = x % (n - 1);
        if (qi[k] > qj[k]) {
            swap(qi[k], qj[k]);
        } else {
            qj[k] += 1;
        }
        x = (x + d) % (n * (n - 1));
    }
    // construct a sparse table
    vector<int> rank(n);
    whole(iota, rank, 0);
    whole(sort, rank, [&](int i, int j) { return s[i] < s[j]; });
    vector<int> lcp_0(n-1);
    repeat (i,n-1) lcp_0[i] = compute_lcp(s[rank[i]], s[rank[i+1]]);
    sparse_table<min_t> table(lcp_0);
    vector<int> rank_of(n);
    repeat (i,n) rank_of[rank[i]] = i;
    // compute
    ll result = 0;
    repeat (q,m) {
        int l = rank_of[qi[q]];
        int r = rank_of[qj[q]];
        result += table(l, r);
    }
    // output
    cout << result << endl;
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0