結果

問題 No.31 悪のミックスジュース
ユーザー hitonanode
提出日時 2020-09-24 22:06:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,514 bytes
コンパイル時間 2,502 ms
コンパイル使用メモリ 198,696 KB
最終ジャッジ日時 2025-01-14 20:06:42
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 1 WA * 16
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
#define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i<i##_end_;i++)
#define REP(i, n) FOR(i,0,n)
template <typename T> bool chmin(T &m, const T q) { if (m > q) {m = q; return true;} else return false; }
template <typename T> istream &operator>>(istream &is, vector<T> &vec) { for (auto &v : vec) is >> v; return is; }

int main()
{
    int N, V;
    cin >> N >> V;
    vector<lint> C(N);
    cin >> C;
    vector<lint> v2c(N + 1);
    REP(i, N) v2c[i + 1] = v2c[i] + C[i];

    vector<lint> dp(N * 4 + 1, 1e18);
    dp[0] = 0;
    FOR(i, 1, dp.size())
    {
        if (i < int(v2c.size()))
        {
            chmin(dp[i], v2c[i]);
        }
        REP(j, i)
        {
            chmin(dp[i], dp[j] + dp[i - j]);
        }
    }

    auto solve = [&](auto &&solve, lint v) -> vector<lint> {
        vector<lint> ret(N, 1e18);
        if (v + N - 1 < dp.size())
        {
            copy(dp.begin() + v, dp.begin() + v + N, ret.begin());
        }
        else
        {
            int head = (v - N) / 2;
            vector<lint> half = solve(solve, head);
            REP(i, N) REP(j, N)
            {
                REP(l, N * 2)
                {
                    lint k = head * 2 + i + j + l - v;
                    if (k >= 0 and k < int(ret.size())) chmin(ret[k], half[i] + half[j] + v2c[k]);
                }
            }
        }
        return ret;
    };
    cout << v2c[N] + solve(solve, V - N)[0] << '\n';
}
0