結果

問題 No.1975 Zigzag Sequence
ユーザー 👑 KazunKazun
提出日時 2022-06-10 22:38:46
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 442 ms / 2,000 ms
コード長 4,140 bytes
コンパイル時間 213 ms
コンパイル使用メモリ 81,472 KB
実行使用メモリ 127,684 KB
最終ジャッジ日時 2023-10-21 05:29:16
合計ジャッジ時間 10,183 ms
ジャッジサーバーID
（参考情報） 		judge10 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
53,296 KB
testcase_01 AC 39 ms
53,296 KB
testcase_02 AC 38 ms
53,296 KB
testcase_03 AC 126 ms
78,356 KB
testcase_04 AC 99 ms
76,956 KB
testcase_05 AC 358 ms
97,888 KB
testcase_06 AC 262 ms
94,456 KB
testcase_07 AC 115 ms
78,304 KB
testcase_08 AC 244 ms
88,328 KB
testcase_09 AC 379 ms
108,960 KB
testcase_10 AC 257 ms
91,648 KB
testcase_11 AC 309 ms
95,540 KB
testcase_12 AC 122 ms
79,388 KB
testcase_13 AC 37 ms
53,296 KB
testcase_14 AC 37 ms
53,296 KB
testcase_15 AC 37 ms
53,296 KB
testcase_16 AC 37 ms
53,296 KB
testcase_17 AC 37 ms
53,296 KB
testcase_18 AC 255 ms
95,976 KB
testcase_19 AC 250 ms
95,976 KB
testcase_20 AC 293 ms
96,820 KB
testcase_21 AC 298 ms
96,820 KB
testcase_22 AC 322 ms
127,684 KB
testcase_23 AC 323 ms
127,684 KB
testcase_24 AC 384 ms
124,144 KB
testcase_25 AC 391 ms
124,144 KB
testcase_26 AC 356 ms
96,268 KB
testcase_27 AC 442 ms
113,692 KB
testcase_28 AC 345 ms
123,880 KB
testcase_29 AC 361 ms
123,880 KB
testcase_30 AC 368 ms
123,772 KB
testcase_31 AC 371 ms
123,764 KB
testcase_32 AC 297 ms
95,976 KB
testcase_33 AC 299 ms
95,976 KB
testcase_34 AC 343 ms
99,772 KB
testcase_35 AC 343 ms
99,772 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

class Binary_Indexed_Tree():
    def __init__(self, L, calc, unit, inv, index=1):
        """ calc を演算とする N 項の Binary Indexed Tree を作成
        calc: 演算 (2変数関数, 可換群)
        unit: 群 calc の単位元 (x+e=e+x=xを満たすe)
        inv : 群 calc の逆元 (1変数関数, x+inv(x)=inv(x)+x=e をみたす inv(x))
        """
        self.calc=calc
        self.unit=unit
        self.inv=inv
        self.index=index

        N=len(L)
        d=max(1,(N-1).bit_length())
        k=2**d

        X=[None]+[unit]*k

        self.num=k
        self.depth=d

        if L:
            for i in range(len(L)):
                p=i+1
                while p<=k:
                    X[p]=self.calc(X[p],L[i])
                    p+=p&(-p)
        self.data=X

    def index_number(self, k, index=1):
        """ 第 k 要素の値を出力する.
        k    : 数列の要素
        index: 先頭の要素の番号
        """
        return self.sum(k,k,index)

    def add(self, k, x, index=1):
        """ 第 k 要素に x を加え, 更新を行う.
        k    : 数列の要素
        x    : 加える値
        index: 先頭の要素の番号
        """
        data=self.data; calc=self.calc
        p=k+(1-index)
        while p<=self.num:
            data[p]=calc(self.data[p],x)
            p+=p&(-p)

    def update(self, k, x, index=1):
        """ 第 k 要素を x に変え, 更新を行う.
        k: 数列の要素
        x: 更新後の値
        """

        a=self.index_number(k,index)
        y=self.calc(self.inv(a),x)

        self.add(k,y,index)

    def sum(self, From, To, index=1):
        """ 第 From 要素から第 To 要素までの総和を求める.
        ※From!=1を使うならば, 群でなくてはならない.
        From : 始まり
        To   : 終わり
        index: 先頭の要素の番号
        """
        alpha=max(1,From+(1-index))
        beta=min(self.num,To+(1-index))

        if alpha>beta:
            return self.unit
        elif alpha==1:
            return self.__section(beta)
        else:
            return self.calc(self.inv(self.__section(alpha-1)),self.__section(beta))

    def __section(self,x):
        """ B[1]+...+B[x] を求める. """
        data=self.data; calc=self.calc
        S=self.unit
        while x>0:
            S=calc(data[x],S)
            x-=x&(-x)
        return S

    def all_sum(self):
        return self.data[-1]

    def binary_search(self, cond, index=1):
        """ cond(B[1]+...+B[k]) を満たす最小の k を返す.

        cond: 単調増加

        ※ cond(unit)=True の場合の返り値は index-1
        ※ cond(B[1]+...+B[k]) なる k が存在しない場合の返り値は self.num+index
        """

        if cond(self.unit):
            return index-1

        j=0
        r=self.num
        t=r
        data=self.data; calc=self.calc
        alpha=self.unit

        for _ in range(self.depth+1):
            if j+t<=self.num:
                beta=calc(alpha,data[j+t])
                if not cond(beta):
                    alpha=beta
                    j+=t
            t>>=1

        return j+index

    def __getitem__(self,index):
        if isinstance(index,int):
            return self.index_number(index,self.index)
        else:
            return [self.index_number(t,self.index) for t in index]

    def __setitem__(self,index,val):
        self.update(index,val,self.index)

#==================================================
from operator import add,neg

N=int(input())
A=list(map(int,input().split()))
A_ind={a:i for i,a in enumerate(sorted(set(A)))}

A=[A_ind[a] for a in A]

Mod=10**9+7

TWO=[0]*(N+1); TWO[0]=1
for i in range(1,N+1):
    TWO[i]=2*TWO[i-1]
    TWO[i]%=Mod

L=Binary_Indexed_Tree([0]*N, add,0,neg,0)
R=Binary_Indexed_Tree([0]*N, add,0,neg,0)

for i,a in enumerate(A):
    R.add(a, TWO[N-1-i],0)

X=0
for i in range(N):
    P1=L.sum(0,A[i]-1,0)
    Q1=R.sum(0,A[i]-1,0)

    P2=L.sum(A[i]+1,N-1,0)
    Q2=R.sum(A[i]+1,N-1,0)

    X+=P1*Q1+P2*Q2; X%=Mod

    L.add(A[i], TWO[i], 0)
    R.add(A[i], -TWO[N-1-i], 0)

print(X)
0