結果

問題 No.189 SUPER HAPPY DAY
ユーザー rpy3cpp
提出日時 2015-08-04 17:36:19
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 356 ms / 5,000 ms
コード長 1,113 bytes
コンパイル時間 85 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 19,456 KB
最終ジャッジ日時 2024-07-18 01:15:32
合計ジャッジ時間 4,458 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

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