結果

問題 No.1233 割り切れない気持ち
コンテスト
ユーザー iby
提出日時 2020-11-28 16:43:18
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,395 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 256 ms
コンパイル使用メモリ 82,032 KB
実行使用メモリ 115,364 KB
最終ジャッジ日時 2024-09-12 23:14:06
合計ジャッジ時間 15,746 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 29 WA * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

class SegTree():
    # 1-indexed
    def __init__(self, N, function, basement):
        self.n = N
        self.K = (self.n-1).bit_length()
        self.f = function
        self.b = basement
        self.seg = [basement]*(2**(self.K+1)+1)

    def update(self, k, value):
        X = 1 << (self.K)
        k += X
        self.seg[k] += value
        while k:
            k = k >> 1
            self.seg[k] = self.f(self.seg[k << 1], self.seg[(k << 1) | 1])

    def query(self, L, R):
        num = 1 << (self.K)
        L += num
        R += num
        vL = self.b
        vR = self.b
        while L < R:
            if L & 1:
                vL = self.f(vL, self.seg[L])
                L += 1
            if R & 1:
                R -= 1
                vR = self.f(self.seg[R], vR)
            L >>= 1
            R >>= 1
        return self.f(vL, vR)

def f(a,b):
    return a + b

n = int(input())
a = list(map(int, input().split()))

cash = [-1]*(200010)
segki = SegTree(200010, f, 0)
for v in a:
    segki.update(v, 1)

def get_quotient_sum(div):
    if cash[div] >= 0:
        return cash[div]
    res = 0
    for i in range(200001):
        if i*div >= 200000:
            break
        else:
            res += i*segki.query(div*i, min(div*(i + 1), 200001))
    cash[div] = res
    return res


ans = n*sum(a)
for i, v in enumerate(a):
    ans -= v*get_quotient_sum(v)
print(ans)
0