結果

問題 No.1973 Divisor Sequence
ユーザー lilictakalilictaka
提出日時 2022-06-12 02:16:34
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 685 bytes
コンパイル時間 255 ms
コンパイル使用メモリ 81,756 KB
実行使用メモリ 53,672 KB
最終ジャッジ日時 2023-10-23 18:42:42
合計ジャッジ時間 4,371 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,672 KB
testcase_01 AC 38 ms
53,672 KB
testcase_02 TLE -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
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 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import sqrt,ceil,floor
def yaku(n):#nの約数を列挙するn <= 10 ** 12(order)
    rev = []
    for i in range(1,ceil(sqrt(n))+1):
        if n % i == 0:
            rev.append(i)
            if n //i == i:
                continue
            rev.append(n//i)
    rev = list(set(rev))
    return sorted(rev)
N,M = map(int,input().split())
P = yaku(M)
Pset = set(P)
m = len(P)
dp = [1] *m 
mod = 10 ** 9+ 7
for i in range(1,N):
    nx = [0] * m
    for p in range(m):
        for q in range(m):
            if P[p] * P[q] in Pset:
                nx[q] += dp[p]
                nx[q] %= mod
    dp = nx
ans = 0
for p in range(m):
    ans += dp[p]
    ans %= mod
print(ans)
0