結果

問題 No.368 LCM of K-products
ユーザー 🍡yurahuna🍡yurahuna
提出日時 2016-04-11 22:17:14
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 569 ms / 2,000 ms
コード長 1,087 bytes
コンパイル時間 231 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,776 KB
最終ジャッジ日時 2024-09-22 14:41:35
合計ジャッジ時間 10,070 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 428 ms
11,392 KB
testcase_01 AC 510 ms
11,776 KB
testcase_02 AC 485 ms
11,392 KB
testcase_03 AC 487 ms
11,392 KB
testcase_04 AC 495 ms
11,520 KB
testcase_05 AC 569 ms
11,392 KB
testcase_06 AC 54 ms
11,264 KB
testcase_07 AC 55 ms
11,264 KB
testcase_08 AC 55 ms
11,392 KB
testcase_09 AC 56 ms
11,264 KB
testcase_10 AC 55 ms
11,264 KB
testcase_11 AC 56 ms
11,264 KB
testcase_12 AC 58 ms
11,264 KB
testcase_13 AC 313 ms
11,392 KB
testcase_14 AC 454 ms
11,264 KB
testcase_15 AC 501 ms
11,392 KB
testcase_16 AC 447 ms
11,264 KB
testcase_17 AC 396 ms
11,392 KB
testcase_18 AC 506 ms
11,520 KB
testcase_19 AC 171 ms
11,264 KB
testcase_20 AC 353 ms
11,648 KB
testcase_21 AC 127 ms
11,264 KB
testcase_22 AC 495 ms
11,392 KB
testcase_23 AC 55 ms
11,264 KB
testcase_24 AC 56 ms
11,264 KB
testcase_25 AC 57 ms
11,520 KB
testcase_26 AC 54 ms
11,264 KB
testcase_27 AC 58 ms
11,264 KB
testcase_28 AC 59 ms
11,264 KB
testcase_29 AC 57 ms
11,264 KB
testcase_30 AC 59 ms
11,520 KB
testcase_31 AC 58 ms
11,264 KB
testcase_32 AC 58 ms
11,264 KB
testcase_33 AC 233 ms
11,136 KB
testcase_34 AC 478 ms
11,264 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import defaultdict as ddict

# sqrt(10^9) < 40000未満の素数をあらかじめ列挙しておく (4203個)
prime = []
memo = [True] * 40000
i = 2
while i < 40000:
    if memo[i]:
        prime.append(i)
        j = 2 * i
        while j < 40000:
            memo[j] = False
            j += i
    i += 1

# print(len(prime))

###########################

N, K = map(int, input().split())
a = list(map(int, input().split()))

exponents = ddict(list)
for ax in a:
    # 素因数分解(作成済みの素数表primeを使う)
    cnt = ddict(int)
    for p in prime:
        while ax % p == 0:
            ax //= p
            cnt[p] += 1

    # sqrt(10^9) 以下の素数で割り切れなかったら、a_iは素数
    if ax != 1:
        cnt[ax] += 1

    for key, value in cnt.items():
        exponents[key].append(value)

# print(exponents)

mod = 10**9 + 7

ans = 1
for prime, li in exponents.items():
    li.sort(reverse = True)
    num = 0
    for i in range(min(K, len(li))):
        num += li[i]

    ans *= pow(prime, num, mod)
    ans %= mod

print(ans)
0