結果

問題 No.778 クリスマスツリー
ユーザー betrue12betrue12
提出日時 2019-06-22 17:22:13
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 154 ms / 2,000 ms
コード長 1,960 bytes
コンパイル時間 1,855 ms
コンパイル使用メモリ 173,212 KB
実行使用メモリ 28,132 KB
最終ジャッジ日時 2024-04-22 03:08:14
合計ジャッジ時間 3,694 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
8,272 KB
testcase_01 AC 4 ms
8,256 KB
testcase_02 AC 4 ms
8,104 KB
testcase_03 AC 4 ms
8,336 KB
testcase_04 AC 4 ms
8,132 KB
testcase_05 AC 4 ms
8,200 KB
testcase_06 AC 101 ms
28,056 KB
testcase_07 AC 50 ms
15,028 KB
testcase_08 AC 154 ms
21,172 KB
testcase_09 AC 143 ms
17,476 KB
testcase_10 AC 140 ms
17,484 KB
testcase_11 AC 140 ms
17,344 KB
testcase_12 AC 131 ms
17,504 KB
testcase_13 AC 88 ms
17,444 KB
testcase_14 AC 101 ms
28,132 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

struct EulerTour{
    // オイラーツアー
    vector<int> tour;
    // 各頂点番号に対する、tourの最左/最右の位置
    vector<int> L, R;

    EulerTour(){}
    EulerTour(int N, vector<int> G[], int root=0){ initialize(N, G, root); }

    void dfs(int i, int p, vector<int> G[]){
        L[i] = tour.size();
        tour.push_back(i);
        for(int j : G[i]) if(j != p){
            dfs(j, i, G);
            tour.push_back(i);
        }
        R[i] = tour.size()-1;
    }
    void initialize(int N, vector<int> G[], int root){
        L.resize(N);
        R.resize(N);
        dfs(root, -1, G);
    }
};

template<typename T>
struct BIT {
    int n;
    vector<T> dat;

    BIT(int n=0){
        initialize(n);
    }

    void initialize(int nin){
        n = nin;
        dat.resize(n, 0);
    }

    T sum(int i){
        T s = 0;
        while(i >= 0){
            s += dat[i];
            i = (i & (i+1)) - 1;
        }
        return s;
    }

    T sum_between(int i, int j){
        return sum(j) - sum(i-1);
    }

    void plus(int i, T x){
        while(i < n){
            dat[i] += x;
            i |= i+1;
        }
    }

    // a[0]+...+a[ret] >= x
    int lower_bound(T x){
        int ret = -1;
        int k = 1;
        while(2*k <= n) k <<= 1;
        for( ;k>0; k>>=1){
            if(ret+k < n && dat[ret+k] < x){
                x -= dat[ret+k];
                ret += k;
            }
        }
        return ret + 1;
    }
};

vector<int> child[200000];

int main(){
    int N;
    cin >> N;
    for(int i=1; i<N; i++){
        int a;
        cin >> a;
        child[a].push_back(i);
    }
    
    EulerTour et(N, child, 0);
    int len = et.tour.size();
    BIT<int64_t> bit(len);

    int64_t ans = 0;
    for(int i=N-1; i>=0; i--){
        ans += bit.sum_between(et.L[i], et.R[i]);
        bit.plus(et.L[i], 1);
    }
    cout << ans << endl;
    return 0;
}
0