結果
問題 |
No.847 Divisors of Power
|
ユーザー |
![]() |
提出日時 | 2023-05-25 21:31:25 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 131 ms / 2,000 ms |
コード長 | 1,091 bytes |
コンパイル時間 | 346 ms |
コンパイル使用メモリ | 82,088 KB |
実行使用メモリ | 77,524 KB |
最終ジャッジ日時 | 2024-12-24 11:53:05 |
合計ジャッジ時間 | 3,084 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 26 |
ソースコード
# N**Kは素因数分解できないが、Nを素因数分解して、乗数をK倍すればいい # それからdfsで探したらどうか def factorization(n): arr = [] temp = n for i in range(2, int(-(-n**0.5//1))+1): if temp%i==0: cnt=0 while temp%i==0: cnt+=1 temp //= i arr.append([i, cnt]) if temp!=1: arr.append([temp, 1]) if arr==[]: arr.append([n, 1]) return arr N, K, M = map(int, input().split()) N_factors = factorization(N) for i in range(len(N_factors)): if N_factors[i][0] != 1: N_factors[i][1] *= K #print(N_factors) import sys sys.setrecursionlimit(10**7) def dfs(i, num): global ans if i == len(N_factors): if num <= M: ans += 1 #print('num', num) return for j in range(N_factors[i][1]+1): calc = num*pow(N_factors[i][0], j) if calc <= M: dfs(i+1, calc) else: break ans = 0 dfs(0, 1) if N == 1: print(1) else: print(ans)