結果

問題 No.2046 Ans Mod? Mod Ans!
ユーザー titiatitia
提出日時 2022-08-21 15:52:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 625 ms / 4,000 ms
コード長 1,371 bytes
コンパイル時間 229 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 110,464 KB
最終ジャッジ日時 2024-12-16 08:35:37
合計ジャッジ時間 6,417 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
input = sys.stdin.readline

N=int(input())
A=sorted(map(int,input().split()))

LEN=500000

BIT=[0]*(LEN+1) # 1-indexedなtree. 配列BITの長さはLEN+1にしていることに注意。

def update(v,w): # index vにwを加える
    while v<=LEN:
        BIT[v]+=w
        v+=(v&(-v)) # v&(-v)で、最も下の立っているビット. 自分を含む大きなノードへ. たとえばv=3→v=4

def getvalue(v): # [1,v]の区間の和を求める
    ANS=0
    while v!=0:
        ANS+=BIT[v]
        v-=(v&(-v)) # 自分より小さい自分の和を構成するノードへ. たとえばv=14→v=12へ
    return ANS

def bisect_on_BIT(x): # [1,ind]の和がはじめてx以上になるindexを探す

    if x<=0:
        return 0
    
    ANS=0
    h=1<<(LEN.bit_length()-1) # LEN以下の最小の2ベキ
    while h>0:
        if ANS+h<=LEN and BIT[ANS+h]<x:
            x-=BIT[ANS+h]
            ANS+=h
        h//=2

    return ANS+1 # LENまでの和がx未満のとき, LEN+1を返すことに注意

for a in A:
    update(a,1)

SUM=sum(A)
ANS=0
for i in range(N):
    a=A[i]

    ANS+=a*(N-i-1)
    SUM-=a
    ANS-=SUM

    score=0

    for j in range(1,4*10**5):
        if j*a>2*10**5:
            break
        score+=(getvalue((j+1)*a-1)-getvalue(j*a-1))*j
        if j==1:
            score-=1

    ANS+=score*a

print(ANS)
        
    
0