結果
問題 | No.2046 Ans Mod? Mod Ans! |
ユーザー | titia |
提出日時 | 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 |
ソースコード
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)