結果
問題 | No.278 連続する整数の和(2) |
ユーザー |
|
提出日時 | 2015-09-05 11:25:59 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 171 ms / 2,000 ms |
コード長 | 1,026 bytes |
コンパイル時間 | 140 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-12-23 09:04:03 |
合計ジャッジ時間 | 1,867 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 18 |
ソースコード
import mathimport collectionsdef factorize(n):''' returns a list of prime factors of n.ex. factorize(24) = [2, 2, 2, 3]source: Rossetta code: prime factorization (slightly modified)http://rosettacode.org/wiki/Prime_decomposition#Python:_Using_floating_point'''step = lambda x: 1 + (x<<2) - ((x>>1)<<1)maxq = int(math.floor(math.sqrt(n)))d = 1q = n % 2 == 0 and 2 or 3while q <= maxq and n % q != 0:q = step(d)d += 1return q <= maxq and [q] + factorize(n//q) or [n]def sum_of_divisors(n):''' returns the sum of divisors of integer n.n must be a positive integer.the sum includes n itself.ex. sum_of_divisors(6) = 6 + 3 + 2 + 1 = 12'''if n == 1:return 1factors = collections.Counter(factorize(n))result = 1for p, a in factors.items():result *= (p ** (a + 1) - 1)//(p - 1)return resultN = int(input())if N & 1:print(sum_of_divisors(N))else:print(sum_of_divisors(N//2))