結果

問題 No.847 Divisors of Power
ユーザー fujitanaoki0718fujitanaoki0718
提出日時 2019-07-05 23:08:54
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 746 bytes
コンパイル時間 245 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 17,824 KB
最終ジャッジ日時 2024-04-16 08:03:27
合計ジャッジ時間 7,334 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
17,824 KB
testcase_01 AC 86 ms
10,624 KB
testcase_02 AC 300 ms
10,624 KB
testcase_03 AC 38 ms
10,624 KB
testcase_04 AC 1,834 ms
10,752 KB
testcase_05 AC 480 ms
10,624 KB
testcase_06 AC 32 ms
10,624 KB
testcase_07 AC 846 ms
10,624 KB
testcase_08 TLE -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import collections
N,K,M = [int(x) for x in input().split()]
def prime_factors(n):
    i = 2
    factors = []
    while i * i <= n:
        if n % i:
            i += 1
        else:
            n //= i
            factors.append(i)
    if n > 1:
        factors.append(n)
    return factors

ps = sorted([1,*prime_factors(N)])
maxcount = {k:v*K for k,v in collections.Counter(ps).items()}

solve = set([1])

def minimum_next(solve):
    m=2**64
    for s in solve:
        for p in ps:
            if collections.Counter(prime_factors(p*s))[p]>maxcount[p]:
                continue
            if p*s not in solve:
                m=min(m,p*s)
    solve.add(m)
while(max(solve)<=M):
    minimum_next(solve)
print(len([x for x in solve if x<=M]))
0