結果
| 問題 |
No.1233 割り切れない気持ち
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-02-11 00:07:16 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 2,439 ms / 3,153 ms |
| コード長 | 1,992 bytes |
| コンパイル時間 | 379 ms |
| コンパイル使用メモリ | 82,456 KB |
| 実行使用メモリ | 115,928 KB |
| 最終ジャッジ日時 | 2024-09-28 17:22:21 |
| 合計ジャッジ時間 | 22,363 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 39 |
ソースコード
## https://yukicoder.me/problems/no/1233
import math
def main():
N = int(input())
A = list(map(int, input().split()))
A.sort()
cum_a_list = [0] * (N + 1)
for i in range(N):
cum_a_list[i + 1] = cum_a_list[i] + A[i]
a_map = {}
for a in A:
if a not in a_map:
a_map[a] = 0
a_map[a] += 1
a_keys = list(a_map.keys())
a_keys.sort()
sqrt_a = int(math.sqrt(a_keys[-1]))
# Aの値がmod未満のものをcount
index = 0
a0 = a_keys[index]
answer = 0
while a0 <= sqrt_a + 1:
for a in A:
m = a % a0
answer += m * a_map[a0]
if index < len(a_keys) - 1:
index += 1
a0 = a_keys[index]
else:
break
# Aの商がsqrt_a以下のものをcount
while index < len(a_keys):
max_k = a_keys[-1] // a0
v = -1
q = 0
while q <= max_k:
low = v + 1
high = len(A) - 1
if A[low] >= (q + 1) * a0:
q += 1
continue
low = v + 1
high = len(A) - 1
if high - low < 0:
break
while high - low > 1:
mid = (low + high) // 2
if A[mid] < (q + 1) * a0:
low = mid
else:
high = mid
if A[high] < (q + 1) * a0:
v0 = high
else:
v0 = low
c = v0 - v
a_sum = cum_a_list[v0 + 1] - cum_a_list[v + 1]
ans = a_sum - c * q * a0
answer += ans * a_map[a0]
v = v0
if v == len(A) - 1:
break
else:
q = A[v + 1] // a0
if index < len(a_keys) - 1:
index += 1
a0 = a_keys[index]
else:
break
print(answer)
if __name__ == "__main__":
main()