結果

問題 No.31 悪のミックスジュース
ユーザー rpy3cpprpy3cpp
提出日時 2015-09-03 18:52:06
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 1,691 bytes
コンパイル時間 174 ms
コンパイル使用メモリ 11,120 KB
実行使用メモリ 41,068 KB
最終ジャッジ日時 2023-09-26 02:34:43
合計ジャッジ時間 18,985 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 16 ms
41,068 KB
testcase_01 AC 2,509 ms
19,372 KB
testcase_02 AC 3,820 ms
25,004 KB
testcase_03 TLE -
testcase_04 AC 17 ms
8,488 KB
testcase_05 AC 17 ms
8,448 KB
testcase_06 AC 17 ms
8,400 KB
testcase_07 AC 17 ms
8,348 KB
testcase_08 AC 39 ms
8,512 KB
testcase_09 TLE -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import itertools


def read_data():
    N, V = map(int, input().split())
    Cs = list(map(int, input().split()))
    return N, V, Cs


def select_most_efficient(Ds):
    n = 0
    cost = 1
    record = None
    for k, d in enumerate(Ds, 1):
        det = k * cost - n * d
        if det > 0:
            n = k
            cost = d
            record = k - 1
    return record


def solve(N, V, Cs):
    if V <= N:
        return sum(Cs)
    if N == 1:
        return Cs[0] * V
    V -= N
    Ds = list(itertools.accumulate(Cs))
    eff = select_most_efficient(Ds)
    return calc_cost(N, V, Cs, Ds, eff) + sum(Cs)


def calc_cost(N, V, Cs, Ds, idx):
    d = Ds[idx]
    m, r = divmod(V, idx + 1)
    if r == 0:
        return d * m
    dp = get_dp(N, Ds, idx)
    # V <= istart * (idx + 1) + N * idx
    i0 = max(0, 1 + (V - N * idx - 1) // (idx + 1))
    return min(d * i + dp[V - (idx + 1) * i] for i in range(i0, m + 1))


def get_dp(N, Ds, m):
    '''コスト Ds[i], 容量 i　のアイテムについて、
    dp[k][n][v] k 番目までのアイテムを 合計 n 個使ったときに、容量を v とする最小コストの値
    dp[k][n][v] = min(dp[k - 1][n][v], dp[k][n - 1][v - i] + Ds[i])
    '''
    dp = [[float('inf')] * (N * i + 1) for i in range(m + 1)]
    dp[0][0] = 0
    for i, d in enumerate(Ds, 1):
        for n, (dp_n, prev_dp) in enumerate(zip(dp[1:], dp[:-1]), 1):
            for v in range(i, i * n + 1):
                dp_n[v] = min(dp_n[v], prev_dp[v - i] + d)
    dp_min = dp[m]
    for dp_n in dp:
        for i, dp_ni in enumerate(dp_n):
            dp_min[i] = min(dp_min[i], dp_ni)
    return dp_min

N, V, Cs = read_data()
print(solve(N, V, Cs))
0