結果

問題 No.778 クリスマスツリー
ユーザー SI1000-githubSI1000-github
提出日時 2022-03-31 14:38:17
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 6,191 bytes
コンパイル時間 179 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 567,420 KB
最終ジャッジ日時 2024-11-16 21:51:35
合計ジャッジ時間 14,731 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 45 ms
54,912 KB
testcase_01 AC 43 ms
54,912 KB
testcase_02 AC 42 ms
54,784 KB
testcase_03 AC 44 ms
54,656 KB
testcase_04 AC 45 ms
55,040 KB
testcase_05 AC 47 ms
54,656 KB
testcase_06 AC 1,622 ms
500,056 KB
testcase_07 AC 884 ms
315,832 KB
testcase_08 TLE -
testcase_09 AC 1,385 ms
317,148 KB
testcase_10 AC 1,016 ms
316,240 KB
testcase_11 AC 1,084 ms
317,516 KB
testcase_12 AC 1,004 ms
316,500 KB
testcase_13 AC 740 ms
315,868 KB
testcase_14 AC 1,340 ms
494,636 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#!/usr/bin/python
# -*- coding: utf-8 -*-


# import
from sys import stdin, setrecursionlimit
setrecursionlimit(10**8)




# def

def int_map():
    return map(int,input().split())

def int_list():
    return ['Null']+list(map(int,input().split()))



# 整数のリストのリスト:N+1行のリスト(各行のリストは同じ長さでなくてよい.)
def int_mtx(N):
    x = [['Null']+list(map(int, stdin.readline().split())) for _ in range(N)]
    x.insert(0,['Null'])
    return x

# 文字のリストのリスト:N+1行のリスト(各行のリストは同じ長さでなくてよい.)
# AAAAA
# BBBBB
# CCCCC
# のような標準入力を
# '0AAAAA'
# '0BBBBB'
# '0CCCCC'
# と受け取る
def str_mtx(N):
    x = ['0'+stdin.readline()[:-1] for _ in range(N)]
    x.insert(0,'0')
    return x

# 文字のリストのリスト:N+1行のリスト(各行のリストは同じ長さでなくてよい.)
# AAAAA XXXXX
# BBBBB YYYYY
# CCCCC ZZZZZ
# のような標準入力を
# [[['0'],[AAAAA],[XXXXX]],
#  [['0'],[BBBBB],[YYYYY]],
#  [['0'],[CCCCC],[ZZZZZ]] ]
# と受け取る
def str_list(N):
    x = [['0']+list(map(str, stdin.readline().split())) for _ in range(N)]
    x.insert(0,['0'])
    return x


# 高さH+1, 幅W+1, のゼロ行列の作成
def zero_mtx(H,W):
    x = [[0]*(W+1) for i in range(H+1)]
    return x


# リストをスペースで分割する(先頭を省く)
def print_space(l):
    return print(" ".join([str(x) for x in l[1:]]))





# ゼロインデックスの場合
# 整数のリストのリスト:N行のリスト(各行のリストは同じ長さでなくてよい.)
def int_mtx_0(N):
    x = [list(map(int, stdin.readline().split())) for _ in range(N)]
    return x

# ゼロインデックスの場合
# 文字のリストのリスト:N+1行のリスト(各行のリストは同じ長さでなくてよい.)
def str_mtx_0(N):
    x = [list(stdin.readline()[:-1]) for _ in range(N)]
    return x

# ゼロインデックスの場合
# リストをスペースで分割する(先頭を省かない)
def print_space_0(l):
    return print(" ".join([str(x) for x in l]))

# ゼロインデックスの場合
# 高さH, 幅W, のゼロ行列の作成
def zero_mtx_0(H,W):
    x = [[0]*W for i in range(H)]
    return x

def int_list_0():
    return list(map(int,input().split()))




## インポート
# from collections import deque

# 順列に使う
# import itertools

# 最大公約数などに使う
# import math

# リストの要素の数をdict形式で
# import collections

# 二次元配列のコピーを作りたいとき
# a_copy = deepcopy(a)
# from copy import deepcopy





from collections import deque, defaultdict


# 絶対に有効化する
#import pypyjit
#pypyjit.set_param('max_unroll_recursion=-1')


N = int(input())
# UV については無向辺のつもりで書けば良い
cccc = int_list_0()
UV = []
for i in range(N-1):
    UV.append([i+1,cccc[i]])

root_id = 0





G = [[] for _ in range(N)]
for ui,vi in UV:
    G[ui].append(vi)
    G[vi].append(ui)

# Euler Tour Technique
S = []
FS = [0]*N
AS = defaultdict(list)
depth = [0]*N
def dfs(v, p, d):
    depth[v] = d
    FS[v] = len(S)
    AS[v].append(len(S))
    S.append(v)
    for w in G[v]:
        if w == p:
            continue
        dfs(w, v, d+1)
        AS[v].append(len(S))
        S.append(v)
dfs(root_id, -1, 0)
 
# Sparse Table
L = len(S)
lg = [0]*(L + 1)
for i in range(2, L+1):
    lg[i] = lg[i >> 1] + 1
st = [None]*(lg[L] + 1)
st0 = st[0] = S
b = 1
for i in range(lg[L]):
    st0 = st[i+1] = [p if depth[p] <= depth[q] else q for p, q in zip(st0, st0[b:])]
    b <<= 1
 
# LCA O(1)
def LCA(u, v):
    x = FS[u]; y = FS[v]
    if x > y:
        x, y = y, x
    l = lg[y - x + 1]
    px = st[l][x]; py = st[l][y - (1 << l) + 1]
    return px if depth[px] <= depth[py] else py






class SegTree:
    def __init__(self, A, SegType):
        if SegType == "min":
            self.X_f = min
            self.X_unit = float('inf')

        elif SegType == "max":
            self.X_f = max
            self.X_unit = -float('inf')

        elif SegType == "sum":
            def segfunc(x, y):
                return x+y
            self.X_f = segfunc
            self.X_unit = 0

        elif SegType == "pro":
            def segfunc(x, y):
                return x*y
            self.X_f = segfunc
            self.X_unit = 1

        elif SegType == "xor":
            def segfunc(x, y):
                return x^y
            self.X_f = segfunc
            self.X_unit = 0

        elif SegType == "gcd":
            from math import gcd
            self.X_f = gcd
            self.X_unit = 0

        elif SegType == "lcm":
            from math import gcd
            def lcm(a,b):
                return a//gcd(a,b)*b
            self.X_f = lcm
            self.X_unit = 1


        self.N = len(A)
        self.X = [self.X_unit] * (self.N + self.N)
        self.build(A)



    def leaf(self):
        return self.X[self.N:]

    def build(self, seq):
        for i, x in enumerate(seq, self.N):
            self.X[i] = x
        for i in range(self.N - 1, 0, -1):
            self.X[i] = self.X_f(self.X[i << 1], self.X[i << 1 | 1])

    def set_val(self, i, x):
        i += self.N
        self.X[i] = x
        while i > 1:
            i >>= 1
            self.X[i] = self.X_f(self.X[i << 1], self.X[i << 1 | 1])

    def add_val(self, i, x):
        i += self.N
        self.X[i] += x
        while i > 1:
            i >>= 1
            self.X[i] = self.X_f(self.X[i << 1], self.X[i << 1 | 1])

    def fold(self, L, R):
        L += self.N
        R += self.N
        vL = self.X_unit
        vR = self.X_unit
        while L < R:
            if L & 1:
                vL = self.X_f(vL, self.X[L])
                L += 1
            if R & 1:
                R -= 1
                vR = self.X_f(self.X[R], vR)
            L >>= 1
            R >>= 1
        return self.X_f(vL, vR)




A = [0]*len(S)
SegType = "sum"

seg = SegTree(A, SegType)

ans = 0
for i in reversed(range(N)):
    l = AS[i][0]
    r = AS[i][-1]
    ans += seg.fold(l,r+1)
    seg.set_val(AS[i][0],1)

print(ans)
0