結果
問題 | No.811 約数の個数の最大化 |
ユーザー | kamojiro |
提出日時 | 2019-04-12 22:46:02 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,215 bytes |
コンパイル時間 | 253 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 11,776 KB |
最終ジャッジ日時 | 2024-09-14 20:55:46 |
合計ジャッジ時間 | 2,589 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 114 ms
11,648 KB |
testcase_01 | AC | 113 ms
11,776 KB |
testcase_02 | AC | 116 ms
11,776 KB |
testcase_03 | AC | 114 ms
11,648 KB |
testcase_04 | AC | 114 ms
11,648 KB |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | AC | 114 ms
11,776 KB |
testcase_08 | WA | - |
testcase_09 | WA | - |
testcase_10 | WA | - |
testcase_11 | WA | - |
testcase_12 | WA | - |
testcase_13 | WA | - |
testcase_14 | AC | 113 ms
11,776 KB |
ソースコード
from collections import defaultdict import heapq h = [] Prime = [1 for i in range(10**5+1)] Prime[0] = Prime[1] = 0 for i in range(2,10**5+1): if Prime[i] == 1: for j in range(2, 10**5): if i*j <= 10**5: Prime[i*j] = 0 else: break D = defaultdict( int) N, K = map( int, input().split()) A = [] now = N for i in range(2, 10**5): if now == 1: break while now%i == 0: now //= i A.append(i) D[i] += 1 def check(d): ret = 0 for f in d: ret += d[f] return(ret) ans = 1 E = defaultdict( int) while check(E) < K: for f in D: if D[f] >= 1: E[f] += 1 D[f] -= 1 if check(E) >= K: break for f in E: ans *= f**E[f] t = 2 for i in range(10**5+1): if Prime[i] == 1: if E[i] == 0: if ans*i < N: ans *= i else: break for i in range(1,17): for j in range(2, 10**5+1): if Prime[j] == 1: if E[j] <= i: if ans*j < N: ans *= j else: break if ans*2 >= N: break print( ans)