結果

問題 No.1247 ブロック登り
ユーザー 👑 hitonanodehitonanode
提出日時 2020-10-03 03:28:01
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 7,847 bytes
コンパイル時間 2,669 ms
コンパイル使用メモリ 221,148 KB
実行使用メモリ 645,904 KB
最終ジャッジ日時 2023-09-25 04:38:04
合計ジャッジ時間 7,989 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 AC 1 ms
4,380 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 MLE -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
using pint = pair<int, int>;
using plint = pair<lint, lint>;
struct fast_ios {
    fast_ios() { cin.tie(nullptr), ios::sync_with_stdio(false), cout << fixed << setprecision(20); };
} fast_ios_;
#define FOR(i, begin, end) for (int i = (begin), i##_end_ = (end); i < i##_end_; i++)
#define IFOR(i, begin, end) for (int i = (end)-1, i##_begin_ = (begin); i >= i##_begin_; i--)
#define REP(i, n) FOR(i, 0, n)
#define IREP(i, n) IFOR(i, 0, n)
#define ALL(x) (x).begin(), (x).end()
//
template <typename T, typename V>
void ndarray(vector<T>& vec, const V& val, int len) { vec.assign(len, val); }
template <typename T, typename V, typename... Args>
void ndarray(vector<T>& vec, const V& val, int len, Args... args)
{
    vec.resize(len), for_each(begin(vec), end(vec), [&](T& v) { ndarray(v, val, args...); });
}
template <typename T>
bool chmax(T& m, const T q) { return m < q ? (m = q, true) : false; }
template <typename T>
bool chmin(T& m, const T q) { return m > q ? (m = q, true) : false; }
template <typename T1, typename T2>
pair<T1, T2> operator+(const pair<T1, T2>& l, const pair<T1, T2>& r) { return make_pair(l.first + r.first, l.second + r.second); }
template <typename T1, typename T2>
pair<T1, T2> operator-(const pair<T1, T2>& l, const pair<T1, T2>& r) { return make_pair(l.first - r.first, l.second - r.second); }
template <typename T>
vector<T> srtunq(vector<T> vec) { return sort(vec.begin(), vec.end()), vec.erase(unique(vec.begin(), vec.end()), vec.end()), vec; }
template <typename T>
istream& operator>>(istream& is, vector<T>& vec)
{
    return for_each(begin(vec), end(vec), [&](T& v) { is >> v; }), is;
}

// output
template <typename T, typename V>
ostream& dmpseq(ostream&, const T&, const string&, const string&, const string&);
#if __cplusplus >= 201703L
template <typename... T>
ostream& operator<<(ostream& os, const tuple<T...>& tpl)
{
    return apply([&os](auto&&... args) { ((os << args << ','), ...); }, tpl), os;
}
#endif
//
template <typename T1, typename T2>
ostream& operator<<(ostream& os, const pair<T1, T2>& p) { return os << '(' << p.first << ',' << p.second << ')'; }
template <typename T>
ostream& operator<<(ostream& os, const vector<T>& x) { return dmpseq<vector<T>, T>(os, x, "[", ",", "]"); }
template <typename T>
ostream& operator<<(ostream& os, const deque<T>& x) { return dmpseq<deque<T>, T>(os, x, "deq[", ",", "]"); }
template <typename T>
ostream& operator<<(ostream& os, const set<T>& x) { return dmpseq<set<T>, T>(os, x, "{", ",", "}"); }
template <typename T, typename TH>
ostream& operator<<(ostream& os, const unordered_set<T, TH>& x) { return dmpseq<unordered_set<T, TH>, T>(os, x, "{", ",", "}"); }
template <typename T>
ostream& operator<<(ostream& os, const multiset<T>& x) { return dmpseq<multiset<T>, T>(os, x, "{", ",", "}"); }
template <typename TK, typename T>
ostream& operator<<(ostream& os, const map<TK, T>& x) { return dmpseq<map<TK, T>, pair<TK, T>>(os, x, "{", ",", "}"); }
template <typename TK, typename T, typename TH>
ostream& operator<<(ostream& os, const unordered_map<TK, T, TH>& x) { return dmpseq<unordered_map<TK, T, TH>, pair<TK, T>>(os, x, "{", ",", "}"); }
template <typename T, typename V>
ostream& dmpseq(ostream& os, const T& seq, const string& pre, const string& sp, const string& suf)
{
    return os << pre, for_each(begin(seq), end(seq), [&](V x) { os << x << sp; }), os << suf;
}
template <typename T>
void print(const vector<T>& x) { dmpseq<vector<T>, T>(cout, x, "", " ", "\n"); }
#ifdef HITONANODE_LOCAL
#define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl
#else
#define dbg(x) {}
#endif

vector<int> solve(int N, int K, const vector<int> A)
{
    vector<int> ret(N, -2e9);

    vector<vector<int>> matl(N, vector<int>(K + 1, -2e9));
    vector<vector<int>> matr(N, vector<int>(K + 1, -2e9));
    REP(i, N)
    {
        int cost = K * A[i];
        FOR(j, i + 1, N)
        {
            int h = K - (j - i);
            if (h < 0) continue;

            cost += h * A[j];
            int tmp = cost;

            chmax(matr[j][h], tmp);
            if (K - 2 * (j - i) >= 0) chmax(matl[i][K - 2 * (j - i)], tmp);

            vector<vector<int>> Begin(N), End(N);
            int l = max(i, j - h), r = j;
            assert(l <= r);
            if (abs(i - l) % 2 != K % 2) l++;
            if (abs(i - r) % 2 != K % 2) r--;
            Begin[l].emplace_back(tmp), End[r].emplace_back(tmp);

            IREP(k, i)
            {
                int h = K - (j - i) - (j - k);
                if (h < 0) break;
                tmp += h * A[k];
                chmax(matl[k][h], tmp);
                if (h - (j - k) >= 0) chmax(matr[j][h - (j - k)], tmp);

                int l = k, r = min(k + h, j);
                if (abs(i - l) % 2 != K % 2) l++;
                if (abs(i - r) % 2 != K % 2) r--;
                assert(l <= r);
                Begin[l].emplace_back(tmp), End[r].emplace_back(tmp);
            }

            vector<multiset<int>> ms(2);
            REP(i, N)
            {
                ms[i % 2].insert(Begin[i].begin(), Begin[i].end());
                if (ms[i % 2].size()) chmax(ret[i], *prev(ms[i % 2].end()));
                for (auto x : End[i])
                {
                    ms[i % 2].erase(ms[i % 2].lower_bound(x));
                }
            }
        }
    }
    REP(i, N) REP(j, N) if (abs(i - j) % 2 == K % 2)
    {
        int cost = 0;
        if (j <= i)
        {
            FOR(s, j + 1, i)
            {
                cost += A[s] * (i - s);
            }
        }
        else
        {
            FOR(s, i, j)
            {
                cost += A[s] * (s - i);
            }
        }
        int hd = abs(i - j);
        FOR(k, hd, K + 1) if (hd % 2 == k % 2)
        {
            if (j <= i) chmax(ret[i], cost + matr[j][k]);
            else chmax(ret[i], cost + matl[j][k]);
        }
    }
    return ret;
}

int main()
{
    int N, K;
    cin >> N >> K;
    vector<int> A(N);
    cin >> A;

    constexpr int INF = 2e9;
    vector dpl(N, vector(N, vector(K + 1, -INF)));
    vector dpr(N, vector(N, vector(K + 1, -INF)));
    REP(i, N) dpl[i][i][K] = dpr[i][i][K] = A[i] * K;

    IFOR(k, 1, K + 1)
    {
        REP(l, N) FOR(r, l, N) if (dpl[l][r][k] > -INF)
        {
            if (l) chmax(dpl[l - 1][r][k - 1], dpl[l][r][k] + A[l - 1] * (k - 1));
            if (r + 1 < N and k - (r - l) - 1 >= 0) chmax(dpr[l][r + 1][k - (r - l) - 1], dpl[l][r][k] + A[r + 1] * (k - (r - l) - 1));
        }
        REP(l, N) FOR(r, l, N) if (dpr[l][r][k] > -INF)
        {
            if (r + 1 < N) chmax(dpr[l][r + 1][k - 1], dpr[l][r][k] + A[r + 1] * (k - 1));
            if (l and k - (r - l) - 1 >= 0) chmax(dpl[l - 1][r][k - (r - l) - 1], dpr[l][r][k] + A[l - 1] * (k - (r - l) - 1));
        }
    }
    vector<vector<int>> Begin(N), End(N);
    REP(l, N) FOR(r, l, N) REP(k, K)
    {
        int L = l, R = min(r, l + k);
        if (k % 2) L++;
        if (k % 2 != (R - l) % 2) R--;
        if (L <= R)
        {
            Begin[L].emplace_back(dpl[l][r][k]);
            End[R].emplace_back(dpl[l][r][k]);
        }
    }
    REP(l, N) FOR(r, l, N) REP(k, K)
    {
        int R = r, L = max(l, r - k);
        if (k % 2) R--;
        if (k % 2 != (r - L) % 2) L++;
        if (L <= R)
        {
            Begin[L].emplace_back(dpr[l][r][k]);
            End[R].emplace_back(dpr[l][r][k]);
        }
    }

    vector<multiset<int>> ms(2);
    REP(i, N)
    {
        int ret = -INF;
        ms[i % 2].insert(Begin[i].begin(), Begin[i].end());
        if (ms[i % 2].size())
            chmax(ret, *prev(ms[i % 2].end()));
        for (auto x : End[i]) {
            ms[i % 2].erase(ms[i % 2].lower_bound(x));
        }
        cout << ret << endl;
    }
}
0