結果
| 問題 |
No.811 約数の個数の最大化
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2020-09-09 17:43:44 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 253 ms / 2,000 ms |
| コード長 | 958 bytes |
| コンパイル時間 | 284 ms |
| コンパイル使用メモリ | 12,800 KB |
| 実行使用メモリ | 12,032 KB |
| 最終ジャッジ日時 | 2024-12-15 22:36:19 |
| 合計ジャッジ時間 | 2,125 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 12 |
ソースコード
from collections import defaultdict
import math
def Sieve_of_Eratosthenes(maxA):
lst = [-1]*(maxA+1)
lst[0] = 0
lst[1] = 1
for i in range(2, maxA+1):
if lst[i] == -1:
for g in range(i, maxA+1, i):
if lst[g] == -1:
lst[g] = i
return lst
def factorization(n):
d = defaultdict(int)
now = n
while True:
if now == lst[now]:
if now != 1:
d[now] += 1
break
r = lst[now]
while now % r == 0:
now //= r
d[r] += 1
return d
num = 0
MAX = 0
n, k = map(int, input().split())
lst = Sieve_of_Eratosthenes(n)
for i in range(1, n):
common = sum(factorization(math.gcd(i, n)).values())
if common >= k:
v = factorization(i).values()
ans = 1
for g in v:
ans *= (g+1)
if num < ans:
num = ans
MAX = max(MAX, i)
print(MAX)