結果

問題 No.2846 Birthday Cake
ユーザー 👑 rin204rin204
提出日時 2024-08-23 22:11:39
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 648 ms / 2,000 ms
コード長 1,112 bytes
コンパイル時間 267 ms
コンパイル使用メモリ 82,344 KB
実行使用メモリ 363,604 KB
最終ジャッジ日時 2024-08-23 22:11:51
合計ジャッジ時間 10,515 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 35 ms
52,992 KB
testcase_01 AC 49 ms
65,876 KB
testcase_02 AC 543 ms
360,928 KB
testcase_03 AC 36 ms
54,232 KB
testcase_04 AC 36 ms
52,644 KB
testcase_05 AC 648 ms
362,640 KB
testcase_06 AC 507 ms
361,824 KB
testcase_07 AC 583 ms
363,604 KB
testcase_08 AC 410 ms
360,824 KB
testcase_09 AC 385 ms
360,692 KB
testcase_10 AC 370 ms
360,712 KB
testcase_11 AC 308 ms
360,544 KB
testcase_12 AC 475 ms
339,112 KB
testcase_13 AC 483 ms
340,308 KB
testcase_14 AC 38 ms
53,020 KB
testcase_15 AC 42 ms
53,720 KB
testcase_16 AC 37 ms
53,948 KB
testcase_17 AC 36 ms
53,760 KB
testcase_18 AC 35 ms
52,888 KB
testcase_19 AC 35 ms
53,396 KB
testcase_20 AC 321 ms
280,224 KB
testcase_21 AC 34 ms
52,948 KB
testcase_22 AC 492 ms
247,176 KB
testcase_23 AC 485 ms
342,408 KB
testcase_24 AC 297 ms
339,076 KB
testcase_25 AC 34 ms
52,868 KB
testcase_26 AC 175 ms
135,416 KB
testcase_27 AC 176 ms
209,892 KB
testcase_28 AC 36 ms
53,668 KB
testcase_29 AC 36 ms
52,416 KB
testcase_30 AC 42 ms
52,860 KB
testcase_31 AC 363 ms
245,228 KB
testcase_32 AC 287 ms
297,876 KB
testcase_33 AC 325 ms
338,824 KB
testcase_34 AC 229 ms
259,976 KB
testcase_35 AC 177 ms
236,952 KB
testcase_36 AC 520 ms
319,632 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import lcm

k, n = map(int, input().split())
if n == k or k == 1:
    print(1)
    exit()
nCk = [[0] * (n + 1) for _ in range(n + 1)]
nCk[0][0] = 1
for i in range(1, n + 1):
    nCk[i][0] = 1
    for j in range(1, n + 1):
        nCk[i][j] = nCk[i - 1][j - 1] + nCk[i - 1][j]

ng = [False] * 25
ng[13] = True
ng[17] = True
ng[19] = True
ng[23] = True

lc = 1
for i in range(1, n + 1):
    if ng[i]:
        continue
    lc = lcm(lc, i)

dp = [[-1] * (lc + 1) for _ in range((n + 1) * (n + 1) + n + 1)]


def dfs(t, x, tot):
    if tot == lc:
        global ans
        if x == 0:
            return 1
        return 0
    if t == 0:
        return 0

    if ng[t]:
        return dfs(t - 1, x, tot)

    if tot + (lc // t) * x > lc:
        return 0

    if dp[t * (n + 1) + x][tot] != -1:
        return dp[t * (n + 1) + x][tot]

    ret = 0
    for j in range(x + 1):
        ntot = tot + (lc // t) * j
        if ntot > lc:
            break
        ret += dfs(t - 1, x - j, ntot) * nCk[x][j]

    dp[t * (n + 1) + x][tot] = ret
    return ret


ans = dfs(n, k, 0)
if ng[k]:
    ans += 1
print(ans)
0