結果
| 問題 |
No.1616 Joke
|
| コンテスト | |
| ユーザー |
hir355
|
| 提出日時 | 2021-07-22 21:23:59 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 50 ms / 2,000 ms |
| コード長 | 1,151 bytes |
| コンパイル時間 | 165 ms |
| コンパイル使用メモリ | 81,972 KB |
| 実行使用メモリ | 57,984 KB |
| 最終ジャッジ日時 | 2024-07-17 16:12:39 |
| 合計ジャッジ時間 | 1,599 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 13 |
ソースコード
class osa_k:
min_factor = []
def __init__(self, n):
self.min_factor = list(range(n + 1))
i = 2
while i * i <= n:
if self.min_factor[i] < i:
i += 1
continue
for j in range(i*i, n + 1, i):
if self.min_factor[j] == j:
self.min_factor[j] = i
i += 1
def factor(self, n):
res = []
while n > 1:
res.append(self.min_factor[n])
n //= self.min_factor[n]
return res
def divisors(self, n):
res = [1]
while n > 1:
x = self.min_factor[n]
cnt = 0
while n % x == 0:
cnt += 1
n //= x
t = []
p = 1
for i in range(cnt + 1):
for j in range(len(res)):
t.append(res[j] * p)
p *= x
res = t[:]
return res
os = osa_k(1000)
n = int(input())
d = os.divisors(n)
a = 0
b = 1
for i in range(len(d)):
a += d[i]
b *= d[i]
c = 0
for i in range(len(d)):
c += b // d[i]
print(a * b // c)
hir355