#293955 (C++11) No.409 ダイエット

提出ソース

結果

問題	No.409 ダイエット
ユーザー	jell
提出日時	2018-10-23 08:32:16
言語	C++11 (gcc 11.4.0)
結果	AC
実行時間	567 ms / 2,000 ms
コード長	6,731 bytes
コンパイル時間	2,783 ms
コンパイル使用メモリ	162,072 KB
実行使用メモリ	57,168 KB
最終ジャッジ日時	2023-08-12 06:15:28
合計ジャッジ時間	20,259 ms
ジャッジサーバーID （参考情報）	judge12 / judge11

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	3 ms 5,488 KB
testcase_01	AC	2 ms 5,440 KB
testcase_02	AC	2 ms 5,416 KB
testcase_03	AC	2 ms 5,424 KB
testcase_04	AC	2 ms 5,444 KB
testcase_05	AC	2 ms 5,400 KB
testcase_06	AC	2 ms 5,384 KB
testcase_07	AC	2 ms 5,392 KB
testcase_08	AC	2 ms 5,388 KB
testcase_09	AC	2 ms 5,436 KB
testcase_10	AC	2 ms 5,448 KB
testcase_11	AC	2 ms 5,440 KB
testcase_12	AC	2 ms 5,436 KB
testcase_13	AC	2 ms 5,508 KB
testcase_14	AC	2 ms 5,460 KB
testcase_15	AC	2 ms 5,468 KB
testcase_16	AC	2 ms 5,384 KB
testcase_17	AC	2 ms 5,436 KB
testcase_18	AC	2 ms 5,628 KB
testcase_19	AC	2 ms 5,444 KB
testcase_20	AC	2 ms 5,488 KB
testcase_21	AC	2 ms 5,544 KB
testcase_22	AC	2 ms 5,480 KB
testcase_23	AC	2 ms 5,444 KB
testcase_24	AC	2 ms 5,488 KB
testcase_25	AC	2 ms 5,524 KB
testcase_26	AC	2 ms 5,420 KB
testcase_27	AC	2 ms 5,480 KB
testcase_28	AC	2 ms 5,600 KB
testcase_29	AC	2 ms 5,480 KB
testcase_30	AC	3 ms 5,592 KB
testcase_31	AC	2 ms 5,388 KB
testcase_32	AC	2 ms 5,420 KB
testcase_33	AC	2 ms 5,448 KB
testcase_34	AC	2 ms 5,460 KB
testcase_35	AC	4 ms 5,800 KB
testcase_36	AC	4 ms 5,924 KB
testcase_37	AC	4 ms 5,832 KB
testcase_38	AC	5 ms 5,796 KB
testcase_39	AC	4 ms 5,888 KB
testcase_40	AC	5 ms 6,040 KB
testcase_41	AC	5 ms 6,056 KB
testcase_42	AC	4 ms 5,996 KB
testcase_43	AC	5 ms 5,792 KB
testcase_44	AC	4 ms 5,844 KB
testcase_45	AC	5 ms 5,804 KB
testcase_46	AC	4 ms 5,796 KB
testcase_47	AC	5 ms 5,916 KB
testcase_48	AC	4 ms 5,864 KB
testcase_49	AC	5 ms 5,896 KB
testcase_50	AC	5 ms 6,044 KB
testcase_51	AC	4 ms 5,844 KB
testcase_52	AC	5 ms 5,808 KB
testcase_53	AC	4 ms 5,792 KB
testcase_54	AC	4 ms 5,836 KB
testcase_55	AC	232 ms 31,124 KB
testcase_56	AC	227 ms 57,080 KB
testcase_57	AC	567 ms 57,152 KB
testcase_58	AC	168 ms 18,720 KB
testcase_59	AC	292 ms 31,428 KB
testcase_60	AC	152 ms 18,612 KB
testcase_61	AC	404 ms 32,296 KB
testcase_62	AC	525 ms 57,040 KB
testcase_63	AC	406 ms 32,200 KB
testcase_64	AC	221 ms 31,004 KB
testcase_65	AC	537 ms 56,980 KB
testcase_66	AC	559 ms 57,104 KB
testcase_67	AC	345 ms 31,652 KB
testcase_68	AC	289 ms 31,372 KB
testcase_69	AC	398 ms 32,144 KB
testcase_70	AC	383 ms 31,908 KB
testcase_71	AC	178 ms 18,740 KB
testcase_72	AC	564 ms 57,168 KB
testcase_73	AC	443 ms 32,304 KB
testcase_74	AC	272 ms 31,292 KB
testcase_75	AC	532 ms 56,900 KB
testcase_76	AC	316 ms 31,536 KB
testcase_77	AC	184 ms 18,704 KB
testcase_78	AC	242 ms 31,108 KB
testcase_79	AC	162 ms 18,788 KB
testcase_80	AC	536 ms 57,000 KB
testcase_81	AC	437 ms 32,324 KB
testcase_82	AC	324 ms 31,508 KB
testcase_83	AC	344 ms 31,640 KB
testcase_84	AC	297 ms 31,284 KB
testcase_85	AC	25 ms 8,612 KB
testcase_86	AC	280 ms 31,256 KB
testcase_87	AC	381 ms 31,932 KB
testcase_88	AC	138 ms 18,624 KB
testcase_89	AC	359 ms 31,496 KB
testcase_90	AC	135 ms 18,364 KB
testcase_91	AC	381 ms 31,640 KB

権限があれば一括ダウンロードができます

コンパイルメッセージ

main.cpp: In function ‘int main()’:
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];
                ~~~~^
main.cpp:152:24: note: ‘idx’ was declared here
         int latest = 0,idx;
                        ^~~
main.cpp:161:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
         return data[idx];

ソースコード

raw source code

/*
#include <boost/functional/hash.hpp>
#include <boost/multiprecision/cpp_int.hpp>
typedef boost::multiprecision::cpp_int mlint;
*/
#include <bits/stdc++.h>
using namespace std;

#define esc(x) cout << (x) << endl,exit(0)
#define fcout(d) cout << fixed << setprecision(d)
#define each(i,v) for(auto i = begin(v),end_copy = end(v); i != end_copy; ++i)
#define reach(i,v) for(auto i = rbegin(v),rend_copy = rend(v); i != rend_copy; ++i)
#define urep(i,s,t) for(int i = (int)(s),t_copy = (int)t; i <= t_copy; ++i)
#define drep(i,s,t) for(int i = (int)(s),t_copy = (int)t; i >= t_copy; --i)
#define rep(i,n) urep(i,0,(n) - 1)
#define rep1(i,n) urep(i,1,(n))
#define rrep(i,n) drep(i,(n) - 1,0)
#define all(v) begin(v),end(v)
#define rall(v) rbegin(v),rend(v)
#define vct vector
#define prique priority_queue
#define l_bnd lower_bound
#define u_bnd upper_bound
#define rsz resize
#define ers erase
#define emp emplace
#define emf emplace_front
#define emb emplace_back
#define pof pop_front
#define pob pop_back
#define mkp make_pair
#define mkt make_tuple
#define fir first
#define sec second

typedef long long ll;
typedef double db;
typedef vct<vct<ll>> mat;
typedef pair<int,int> pii;
typedef pair<ll,int> pli;
typedef pair<db,int> pdi;
typedef map<int,int> mpii;

#ifdef Local
#define print(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << endl
clock_t starting_time = clock();
#endif

void setup() {
    cin.tie(0);
    ios::sync_with_stdio(false);

    #ifdef Local
    cout << "-- Execution At Local ---\n" << "\n<Standard Output>\n\n";
    auto print_info = []() {
        cout << "\n----------------\n";
        cout << "\nExec time : " << (double)(clock() - starting_time) / CLOCKS_PER_SEC * 1000 << " ms\n\n";
    };
    atexit(print_info);
    #endif
}

template <class T, class U> ostream& operator << (ostream& s, pair<T,U> p) { return s << p.fir << " " << p.sec; }
template <class T> ostream& operator << (ostream& s, vector<T>& v) { for(auto i = v.begin(); i != v.end(); i++) { if(i != begin(v)) s << " "; s << *i; } return s; }

template <typename T> void hash_combine(size_t &seed, T const &key) {
    hash<T> hasher;
    seed ^= hasher(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2);
}

namespace std {
    template <typename T1, typename T2> struct hash<pair<T1,T2>> {
        size_t operator()(pair<T1,T2> const &p) const {
            size_t seed(0);
            hash_combine(seed,p.first);
            hash_combine(seed,p.second);
            return seed;
        }
    };
}

bool odd(ll x) { return x & 1; }
bool even(ll x) { return !odd(x); }
bool parity(ll x, ll y) { return odd(x) ^ even(y); }
bool bit(ll n, uint8_t e) { return n & 1LL << e; }
int ilog(ll x, uint64_t b) { int ret = 0; while(x /= b) ret++; return ret; }
ll qceil(ll x, ll y) { return x > 0 ? (x - 1) / y + 1 :  x / y; }
ll gcd(ll a, ll b) { if(!b) return abs(a); if(a % b) return gcd(b, a % b); return b; }
ll lcm(ll a, ll b) { if(a & b) return a / gcd(a, b) * b; return 0; }
template <class T, class U> bool chmax(T& m, U x) { if(m < x) { m = x; return 1; } return 0; }
template <class T, class U> bool chmin(T& m, U x) { if(m > x) { m = x; return 1; } return 0; }
template <class T> T cmprs(T& v) {
    T tmp = v;
    sort(all(tmp));
    tmp.erase(unique(all(tmp)),end(tmp));
    each(i,v) *i = l_bnd(all(tmp),*i) - begin(tmp) + 1;
    return v;
}

const int dir[8][2] = { {1,0},{0,1},{-1,0},{0,-1},{1,1},{-1,1},{-1,-1},{1,-1} };
const ll inf32 = (1 << 30) - 1;
const ll inf64 = (1LL << 62) - 1;
const ll mod = 1e9 + 7;
const db eps = 1e-15;

template <class K = ll> struct CHT {
    struct line { K slop,incp; int time; };

    unordered_map<K,int> toI;
    vector<K> toK;
    vector<line> data;
    int dom = 1;    // the width of the domain [0,dom)
    int clock = 0;  // the number of total addition
    bool conv = 0;

    CHT(int n) {
        while(n >= dom) dom <<= 1;
        dom = dom * 2 + 1;
        data.resize(dom << 1);
    }

    CHT(const set<K> &S) {
        conv = 1;
        int n = S.size();
        while(n >= dom) dom <<= 1;
        dom = dom * 2 + 1;
        data.resize(dom << 1);
        toK.resize(dom = n);
        for(auto i = begin(S),j = 0; j < n; ++i,++j) {
            toK[j] = *i;
            toI[*i] = j;
        }
    }

    int conv_toI(K x) {
        if(conv) return toI[x];
        return x + dom / 2;
    }

    K conv_toK(int i) {
        if(conv) return toK[i];
        return i - dom / 2;
    }

    K getval(const line& ln, K x) { return x * ln.slop + ln.incp; }

    K query(K x) { return getval(find(conv_toI(x)),x); }

    bool is_lower(const line& ln, int i) { return query(conv_toK(i)) > getval(ln,conv_toK(i)); }

    const line& find(int i) {
        int latest = 0,idx;
        i += dom;
        while(i) {
            if(data[i].time > latest) {
                latest = data[i].time;
                idx = i;
            }
            i >>= 1;
        }
        return data[idx];
    }

    //for the interval [i,j)
    void update(const line& ln, int i, int j) {
        i += dom;
        j += dom;
        while(i < j) {
            if(i & 1) data[i++] = ln;
            if(j & 1) data[--j] = ln;
            i >>= 1, j >>= 1;
        }
    }

    int intersect(const line& ln, int low, int up) {
        while(abs(low - up) > 1) {
            int mid = (up + low) / 2;
            if(is_lower(ln,mid)) low = mid;
            else up = mid;
        }
        return low;
    }

    bool add(K s, K t) {
        line ln = {s,t,++clock};
        if(clock == 1) {
            fill(begin(data),end(data),ln);
            return true;
        }
        int ok = dom - 1;
        int ng = -1;
        while(ok - ng > 1) {
            int mid = (ok + ng) / 2;
            if(find(mid).slop > s) ng = mid;
            else ok = mid;
        }
        if(!ok) {
            if(!is_lower(ln,0)) return false;
            update(ln,0,intersect(ln,0,dom) + 1);
            return true;
        }
        if(is_lower(ln,ok)) {
            update(ln,intersect(ln,ok,0),intersect(ln,ok,dom) + 1);
            return true;
        }
        ok--;
        if(is_lower(ln,ok)) {
            update(ln,intersect(ln,ok,0),intersect(ln,ok,dom) + 1);
            return true;
        }
        return false;
    }
};

ll N,A,B;
ll D[300010];
ll dp[300010];

int main() {
    setup();

    cin >> N >> A >> B >> dp[0];
    rep1(i,N) cin >> D[i];
    CHT<ll> cht(N);
    for(ll i = 1; i <= N; i++) {
        cht.add(-B * (i - 1),dp[i - 1] + A * i + i * (i - 1) / 2 * B);
        dp[i] = D[i] - A * i + i * (i - 1) / 2 * B + cht.query(i);
    }
    ll ans = inf64;
    for(ll i = 0; i <= N; i++) {
        ans = min(ans,dp[i] + (N - i) * (N + 1 - i) / 2 * B - A * (N - i));
    }
    esc(ans);
}

yukicoder

結果

テストケース

コンパイルメッセージ

ソースコード