結果
| 問題 |
No.847 Divisors of Power
|
| コンテスト | |
| ユーザー |
FromBooska
|
| 提出日時 | 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)
FromBooska