結果

問題 No.1975 Zigzag Sequence
ユーザー Navier_BoltzmannNavier_Boltzmann
提出日時 2022-06-11 08:23:16
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,051 ms / 2,000 ms
コード長 4,463 bytes
コンパイル時間 221 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 116,220 KB
最終ジャッジ日時 2024-09-21 16:43:49
合計ジャッジ時間 20,907 ms
ジャッジサーバーID
(参考情報)
judge5 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 43 ms
54,784 KB
testcase_01 AC 44 ms
55,256 KB
testcase_02 AC 44 ms
55,040 KB
testcase_03 AC 256 ms
81,060 KB
testcase_04 AC 206 ms
80,488 KB
testcase_05 AC 892 ms
95,236 KB
testcase_06 AC 645 ms
92,456 KB
testcase_07 AC 240 ms
80,812 KB
testcase_08 AC 627 ms
87,700 KB
testcase_09 AC 880 ms
103,140 KB
testcase_10 AC 646 ms
90,924 KB
testcase_11 AC 705 ms
95,676 KB
testcase_12 AC 278 ms
82,408 KB
testcase_13 AC 45 ms
54,784 KB
testcase_14 AC 44 ms
54,784 KB
testcase_15 AC 47 ms
54,912 KB
testcase_16 AC 45 ms
55,168 KB
testcase_17 AC 45 ms
54,912 KB
testcase_18 AC 543 ms
99,284 KB
testcase_19 AC 543 ms
99,572 KB
testcase_20 AC 668 ms
98,536 KB
testcase_21 AC 664 ms
98,932 KB
testcase_22 AC 651 ms
113,472 KB
testcase_23 AC 654 ms
114,040 KB
testcase_24 AC 792 ms
114,428 KB
testcase_25 AC 797 ms
113,840 KB
testcase_26 AC 917 ms
93,448 KB
testcase_27 AC 1,051 ms
105,636 KB
testcase_28 AC 764 ms
114,544 KB
testcase_29 AC 775 ms
114,772 KB
testcase_30 AC 793 ms
116,220 KB
testcase_31 AC 788 ms
115,060 KB
testcase_32 AC 631 ms
96,512 KB
testcase_33 AC 635 ms
95,552 KB
testcase_34 AC 922 ms
95,680 KB
testcase_35 AC 915 ms
96,016 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import bisect
from collections import *
class BIT:

    def __init__(self,len_A):
        self.N = len_A + 10
        self.bit = [0]*(len_A+10)
        
    # sum(A0 ~ Ai)
    # O(log N)
    def query(self,i):
        res = 0
        idx = i+1
        while idx:
            res += self.bit[idx]
            idx -= idx&(-idx)
        return res

    # Ai += x
    # O(log N)
    def update(self,i,x):
        idx = i+1
        while idx < self.N:
            self.bit[idx] += x
            idx += idx&(-idx)
    
    # min_i satisfying {sum(A0 ~ Ai) >= w} (Ai >= 0)
    # O(log N)
    def lower_left(self,w):
        if (w < 0):
            return -1
        x = 0
        k = 1<<(self.N.bit_length()-1)
        while k > 0:
            if x+k < self.N and self.bit[x+k] < w:
                w -= self.bit[x+k]
                x += k
            k //= 2
        return x


class OrderBIT:

    def __init__(self,all_values,sort_flag = False):
        if sort_flag:
            self.A = all_values
        else:
            self.A = sorted(all_values)
        self.B = BIT(len(all_values))
        self.num = 0
        
    def insert_val(self,x,c=1):
        k = bisect.bisect_left(self.A,x)
        self.B.update(k,c)
        self.num += c
    
    def delete_val(self,x,c=1):
        k = bisect.bisect_left(self.A,x)
        self.B.update(k,-c)
        self.num -= c
    
    # find the k-th min_val (k:0-indexed)
    def find_kth_val(self,k):
        if self.num <= k:
            ##### MINIMUM VAL #######
            return -10**9
        return self.A[self.B.lower_left(k+1)]
    
    # count the number of values lower than or equal to x
    def count_lower(self,x):
        if x < self.A[0]:
            return 0
        return self.B.query(bisect.bisect_right(self.A,x)-1)

    # min_val higher than x
    def find_higher(self,x):
        return self.find_kth_val(self.count_lower(x))
        
    # min_val lower than x
    def find_lower(self,x):
        return self.find_kth_val(self.count_lower(x)-1)
        
N = int(input())
A = list(map(int,input().split()))
mod = 10**9 + 7

def segfunc(x, y):
    return (x+y)%mod

ide_ele = 0

class SegTree:
    """
    Segment Tree
    """

    def __init__(self, init_val, segfunc, ide_ele):
        """
        初期化

        init_val: 配列の初期値
        """
        n = len(init_val)
        self.segfunc = segfunc
        self.ide_ele = ide_ele
        self.num = 1 << (n - 1).bit_length()
        self.tree = [ide_ele] * 2 * self.num
        # 配列の値を葉にセット
        for i in range(n):
            self.tree[self.num + i] = init_val[i]
        # 構築していく
        for i in range(self.num - 1, 0, -1):
            self.tree[i] = segfunc(self.tree[2 * i], self.tree[2 * i + 1])

    def update(self, k, x):
        """
        k番目の値をxに更新

        k: index(0-index)
        x: update value
        """
        k += self.num
        self.tree[k] = x
        while k > 1:
            self.tree[k >> 1] = self.segfunc(self.tree[k], self.tree[k ^ 1])
            k >>= 1

    def query(self, l, r):
        """
        [l, r)のsegfuncしたものを得る

        l: index(0-index)
        r: index(0-index)
        """
        res = self.ide_ele

        l += self.num
        r += self.num
        while l < r:
            if l & 1:
                res = self.segfunc(res, self.tree[l])
                l += 1
            if r & 1:
                res = self.segfunc(res, self.tree[r - 1])
            l >>= 1
            r >>= 1
        return res





c = Counter(A)
ans = 0
X = [(A[i],i) for i in range(N)]
X.sort()
cnd = []
U = SegTree([0 for _ in range(N)],segfunc,ide_ele)
V = SegTree([0 for _ in range(N)],segfunc,ide_ele)

for i in range(N):
    a,t = X[i]
    c[a] -= 1
    cnd.append(t)
    if c[a]==0:
        for j in cnd:
            ans += U.query(j+1,N)*V.query(0,j)
            ans %= mod
        for j in cnd:
            U.update(j,pow(2,N-1-j,mod))
            V.update(j,pow(2,j,mod))
        cnd = []

U = SegTree([0 for _ in range(N)],segfunc,ide_ele)
V = SegTree([0 for _ in range(N)],segfunc,ide_ele)
c = Counter(A)
for i in range(N-1,-1,-1):
    a,t = X[i]
    c[a] -= 1
    cnd.append(t)
    if c[a]==0:
        for j in cnd:
            ans += U.query(j+1,N)*V.query(0,j)
            ans %= mod
        for j in cnd:
            U.update(j,pow(2,N-1-j,mod))
            V.update(j,pow(2,j,mod))
        cnd = []
print(ans)
    
0