結果

問題 No.1031 いたずら好きなお姉ちゃん
ユーザー 👑 hitonanodehitonanode
提出日時 2020-04-20 02:39:59
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 227 ms / 3,500 ms
コード長 4,308 bytes
コンパイル時間 2,610 ms
コンパイル使用メモリ 210,276 KB
実行使用メモリ 16,900 KB
最終ジャッジ日時 2024-04-16 00:29:29
合計ジャッジ時間 11,976 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 6 ms
5,376 KB
testcase_04 AC 6 ms
5,376 KB
testcase_05 AC 6 ms
5,376 KB
testcase_06 AC 6 ms
5,376 KB
testcase_07 AC 6 ms
5,376 KB
testcase_08 AC 6 ms
5,376 KB
testcase_09 AC 6 ms
5,376 KB
testcase_10 AC 6 ms
5,376 KB
testcase_11 AC 166 ms
14,520 KB
testcase_12 AC 165 ms
14,612 KB
testcase_13 AC 220 ms
16,900 KB
testcase_14 AC 175 ms
14,924 KB
testcase_15 AC 205 ms
15,984 KB
testcase_16 AC 170 ms
14,484 KB
testcase_17 AC 191 ms
16,104 KB
testcase_18 AC 184 ms
15,704 KB
testcase_19 AC 198 ms
16,184 KB
testcase_20 AC 211 ms
16,776 KB
testcase_21 AC 187 ms
15,588 KB
testcase_22 AC 158 ms
14,356 KB
testcase_23 AC 174 ms
14,692 KB
testcase_24 AC 200 ms
16,080 KB
testcase_25 AC 211 ms
16,632 KB
testcase_26 AC 197 ms
15,900 KB
testcase_27 AC 178 ms
15,008 KB
testcase_28 AC 222 ms
16,820 KB
testcase_29 AC 214 ms
16,696 KB
testcase_30 AC 227 ms
16,700 KB
testcase_31 AC 222 ms
16,640 KB
testcase_32 AC 216 ms
16,772 KB
testcase_33 AC 214 ms
16,692 KB
testcase_34 AC 144 ms
16,696 KB
testcase_35 AC 132 ms
16,692 KB
testcase_36 AC 133 ms
16,776 KB
testcase_37 AC 134 ms
16,776 KB
testcase_38 AC 119 ms
16,692 KB
testcase_39 AC 115 ms
16,200 KB
testcase_40 AC 116 ms
16,356 KB
testcase_41 AC 122 ms
16,688 KB
testcase_42 AC 120 ms
16,560 KB
testcase_43 AC 116 ms
16,692 KB
testcase_44 AC 115 ms
16,460 KB
testcase_45 AC 111 ms
16,172 KB
testcase_46 AC 131 ms
16,360 KB
testcase_47 AC 124 ms
16,272 KB
testcase_48 AC 129 ms
16,772 KB
testcase_49 AC 127 ms
16,772 KB
testcase_50 AC 126 ms
16,692 KB
testcase_51 AC 125 ms
16,700 KB
testcase_52 AC 126 ms
16,880 KB
testcase_53 AC 128 ms
16,640 KB
testcase_54 AC 2 ms
5,376 KB
testcase_55 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_;
#define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i<i##_end_;i++)
#define IFOR(i, begin, end) for(int i=(end)-1,i##_begin_=(begin);i>=i##_begin_;i--)
#define REP(i, n) FOR(i,0,n)
#define IREP(i, n) IFOR(i,0,n)
template<typename T> istream &operator>>(istream &is, vector<T> &vec){ for (auto &v : vec) is >> v; return is; }
#define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl;


// Binary lifting / `Doubling`
// Complexity: O(NlogN) precalculation / O(logN) per query
// <https://atcoder.jp/contests/arc060/submissions/7039451>
struct BinaryLifting
{
    int N, INVALID, lgD;
    std::vector<std::vector<int>> mat;
    BinaryLifting() : N(0), lgD(0) {}
    BinaryLifting(const std::vector<int> &vec_nxt, int INVALID = -1, int lgd = 0) : N(vec_nxt.size()), INVALID(INVALID), lgD(lgd)
    {
        while ((1 << lgD) < N) lgD++;
        mat.assign(lgD, std::vector<int>(N, INVALID));
        mat[0] = vec_nxt;
        for (int i = 0; i < N; i++) if (mat[0][i] < 0 or mat[0][i] >= N) mat[0][i] = INVALID;
        for (int d = 0; d < lgD - 1; d++) {
            for (int i = 0; i < N; i++) if (mat[d][i] != INVALID) mat[d + 1][i] = mat[d][mat[d][i]];
        }
    }
    int kth_next(int now, int k)
    {
        if (k >= (1 << lgD)) exit(8);
        for (int d = 0; k and now != INVALID; d++, k >>= 1) if (k & 1) now = mat[d][now];
        return now;
    }

    // Distance from l to [r, \infty)
    // Requirement: mat[0][i] > i for all i (monotone increasing)
    int distance(int l, int r)
    {
        if (l >= r) return 0;
        int ret = 0;
        for (int d = lgD - 1; d >= 0; d--) {
            if (mat[d][l] < r and mat[d][l] != INVALID) ret += 1 << d, l = mat[d][l];
        }
        if (mat[0][l] == INVALID or mat[0][l] >= r) return ret + 1;
        else return -1; // Unable to reach
    }
};

struct RangeMinimumQuery
{
    using T = int;
    int N;
    T defaultT = 1e9;
    int head;
    vector<T> x;
    T _get(int begin, int end, int pos, int l, int r) const
    {
        if (r <= begin or l >= end) return defaultT;
        if (l >= begin and r <= end) return x[pos];
        return min(_get(begin, end, 2 * pos + 1, l, (l + r) / 2), _get(begin, end, 2 * pos + 2, (l + r) / 2, r));
    }
    RangeMinimumQuery(int N, T defaultT) : N(N), defaultT(defaultT)
    {
        int N_tmp = 1; while (N_tmp < N) N_tmp <<= 1;
        x.assign(N_tmp*2, defaultT), head = N_tmp - 1;
    }
    void build(const vector<T> &vals)
    {
        copy(vals.begin(), vals.end(), x.begin() + head);
        for (int i = head - 1; i >= 0; i--)
        {
            x[i] = min(x[i * 2 + 1], x[i * 2 + 2]);
        }
    }
    RangeMinimumQuery(const vector<T> &vals, T defaultT) : N(vals.size()), defaultT(defaultT)
    {
        int N_tmp = 1; while (N_tmp < N) N_tmp <<= 1;
        x.assign(N_tmp*2, defaultT), head = N_tmp - 1;
        build(vals);
    }
    RangeMinimumQuery() : RangeMinimumQuery(0, T()) {}
    void update(int pos, T val)
    {
        pos += head, x[pos] = val;
        while (pos) pos = (pos - 1) / 2, x[pos] = min(x[pos * 2 + 1], x[pos * 2 + 2]);
    }
    T get(int begin, int end) const { return _get(begin, end, 0, 0, (int)x.size() / 2); }
    friend ostream &operator<<(ostream &os, const RangeMinimumQuery &s)
    {
        os << "[";
        for (int i = 0; i < s.N; i++) os << s.get(i, i + 1) << ",";
        os << "]";
        return os;
    }
};
int main()
{
    int N;
    cin >> N;
    vector<int> P(N + 1);
    REP(i, N) cin >> P[i];
    long long ret = 0;
    REP(_, 2)
    {
        vector<int> Pinv(N + 1);
        REP(i, N) Pinv[P[i]] = i;
        vector<int> nxt(N + 1, N);
        RangeMinimumQuery rmq(N + 1, N);
        IREP(i, N)
        {
            nxt[i] = rmq.get(0, P[i]);
            rmq.update(P[i], i);
        }
        BinaryLifting bl(nxt);
        set<int> ban{N};
        IFOR(v, 1, N + 1)
        {
            int i = Pinv[v];
            int hi = *ban.lower_bound(i);
            ret += bl.distance(i, hi) - 1;
            ban.insert(i);
        }
        reverse(P.begin(), P.begin() + N);
    }
    cout << ret << '\n';
}
0