結果

問題 No.189 SUPER HAPPY DAY
ユーザー rpy3cpprpy3cpp
提出日時 2015-08-04 17:36:19
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 280 ms / 5,000 ms
コード長 1,113 bytes
コンパイル時間 75 ms
コンパイル使用メモリ 10,988 KB
実行使用メモリ 16,752 KB
最終ジャッジ日時 2023-09-25 01:04:11
合計ジャッジ時間 3,810 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 19 ms
8,172 KB
testcase_01 AC 17 ms
8,264 KB
testcase_02 AC 16 ms
8,288 KB
testcase_03 AC 16 ms
8,224 KB
testcase_04 AC 16 ms
8,272 KB
testcase_05 AC 16 ms
8,264 KB
testcase_06 AC 16 ms
8,276 KB
testcase_07 AC 16 ms
8,168 KB
testcase_08 AC 16 ms
8,292 KB
testcase_09 AC 17 ms
8,152 KB
testcase_10 AC 16 ms
8,308 KB
testcase_11 AC 16 ms
8,276 KB
testcase_12 AC 16 ms
8,152 KB
testcase_13 AC 132 ms
11,756 KB
testcase_14 AC 158 ms
11,576 KB
testcase_15 AC 239 ms
14,900 KB
testcase_16 AC 280 ms
16,444 KB
testcase_17 AC 274 ms
15,924 KB
testcase_18 AC 219 ms
14,576 KB
testcase_19 AC 233 ms
13,436 KB
testcase_20 AC 89 ms
10,520 KB
testcase_21 AC 73 ms
9,872 KB
testcase_22 AC 19 ms
8,216 KB
testcase_23 AC 204 ms
16,752 KB
testcase_24 AC 205 ms
16,752 KB
testcase_25 AC 205 ms
16,684 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def solve(M, D):
    mod = 10**9 + 9
    if M > D:
        M, D = D, M
    dp = preprocess(D, mod)
    dpM = calc(M, dp, mod)
    dpD = calc(D, dp, mod)
    n = 0
    for m, d in zip(dpM, dpD):
        n += m * d
        n %= mod
    return n


def preprocess(D, mod):
    n = len(str(D))
    dp = []
    dp_prev = [0] * (n * 9 + 1)
    dp_prev[0] = 1
    dp.append(dp_prev)
    for i in range(1, n):
        dpi = dp_prev[:]
        for j in range(1, 10):
            for k in range(0, 9 * (i - 1) + 1):
                dpi[j+k] += dp_prev[k]
        for j in range(0, 9 * i + 1):
            dpi[j] %= mod
        dp.append(dpi)
        dp_prev = dpi
    return dp


def calc(D, dp, mod):
    ds = [int(d) for d in str(D)]
    n = len(ds)
    dpD = [0] * (n * 9 + 1)
    cumd = 0
    for i, di in enumerate(ds, 1):
        dpi = dp[n - i]
        for j in range(di):
            for k in range(9*(n - i) + 1):
                dpD[cumd + j + k] += dpi[k]
        cumd += di
    dpD[cumd] += 1
    dpD[0] = 0
    return dpD


if __name__ == '__main__':
    M, D = map(int, input().split())
    print(solve(M, D))
0