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

提出ソース

結果

問題	No.409 ダイエット
ユーザー	jell
提出日時	2018-10-23 18:19:34
言語	C++11 (gcc 11.4.0)
結果	AC
実行時間	531 ms / 2,000 ms
コード長	6,854 bytes
コンパイル時間	1,667 ms
コンパイル使用メモリ	168,908 KB
実行使用メモリ	61,772 KB
最終ジャッジ日時	2024-11-19 04:28:25
合計ジャッジ時間	16,871 ms
ジャッジサーバーID （参考情報）	judge5 / judge4

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,248 KB
testcase_02	AC	2 ms 5,248 KB
testcase_03	AC	1 ms 5,248 KB
testcase_04	AC	2 ms 5,248 KB
testcase_05	AC	1 ms 5,248 KB
testcase_06	AC	2 ms 5,248 KB
testcase_07	AC	2 ms 5,248 KB
testcase_08	AC	2 ms 5,248 KB
testcase_09	AC	2 ms 5,248 KB
testcase_10	AC	2 ms 5,248 KB
testcase_11	AC	2 ms 5,248 KB
testcase_12	AC	2 ms 5,248 KB
testcase_13	AC	2 ms 5,248 KB
testcase_14	AC	1 ms 5,248 KB
testcase_15	AC	2 ms 5,248 KB
testcase_16	AC	2 ms 5,248 KB
testcase_17	AC	2 ms 5,248 KB
testcase_18	AC	2 ms 5,248 KB
testcase_19	AC	2 ms 5,248 KB
testcase_20	AC	2 ms 5,248 KB
testcase_21	AC	2 ms 5,248 KB
testcase_22	AC	1 ms 5,248 KB
testcase_23	AC	2 ms 5,248 KB
testcase_24	AC	2 ms 5,248 KB
testcase_25	AC	1 ms 5,248 KB
testcase_26	AC	2 ms 5,248 KB
testcase_27	AC	2 ms 5,248 KB
testcase_28	AC	2 ms 5,248 KB
testcase_29	AC	2 ms 5,248 KB
testcase_30	AC	2 ms 5,248 KB
testcase_31	AC	2 ms 5,248 KB
testcase_32	AC	2 ms 5,248 KB
testcase_33	AC	1 ms 5,248 KB
testcase_34	AC	2 ms 5,248 KB
testcase_35	AC	5 ms 5,248 KB
testcase_36	AC	5 ms 5,248 KB
testcase_37	AC	5 ms 5,248 KB
testcase_38	AC	5 ms 5,248 KB
testcase_39	AC	5 ms 5,248 KB
testcase_40	AC	5 ms 5,248 KB
testcase_41	AC	5 ms 5,248 KB
testcase_42	AC	5 ms 5,248 KB
testcase_43	AC	5 ms 5,248 KB
testcase_44	AC	5 ms 5,248 KB
testcase_45	AC	5 ms 5,248 KB
testcase_46	AC	5 ms 5,248 KB
testcase_47	AC	5 ms 5,248 KB
testcase_48	AC	5 ms 5,248 KB
testcase_49	AC	5 ms 5,248 KB
testcase_50	AC	5 ms 5,248 KB
testcase_51	AC	5 ms 5,248 KB
testcase_52	AC	5 ms 5,248 KB
testcase_53	AC	4 ms 5,248 KB
testcase_54	AC	5 ms 5,248 KB
testcase_55	AC	239 ms 32,508 KB
testcase_56	AC	187 ms 61,304 KB
testcase_57	AC	531 ms 61,772 KB
testcase_58	AC	197 ms 19,628 KB
testcase_59	AC	297 ms 33,492 KB
testcase_60	AC	172 ms 19,140 KB
testcase_61	AC	442 ms 35,768 KB
testcase_62	AC	506 ms 60,984 KB
testcase_63	AC	435 ms 35,864 KB
testcase_64	AC	222 ms 32,124 KB
testcase_65	AC	491 ms 60,780 KB
testcase_66	AC	520 ms 61,592 KB
testcase_67	AC	362 ms 34,808 KB
testcase_68	AC	291 ms 33,464 KB
testcase_69	AC	425 ms 35,828 KB
testcase_70	AC	402 ms 35,336 KB
testcase_71	AC	213 ms 19,680 KB
testcase_72	AC	529 ms 61,476 KB
testcase_73	AC	453 ms 36,072 KB
testcase_74	AC	272 ms 33,124 KB
testcase_75	AC	480 ms 60,780 KB
testcase_76	AC	321 ms 33,840 KB
testcase_77	AC	203 ms 19,648 KB
testcase_78	AC	248 ms 32,472 KB
testcase_79	AC	186 ms 19,436 KB
testcase_80	AC	494 ms 61,092 KB
testcase_81	AC	453 ms 36,060 KB
testcase_82	AC	331 ms 34,124 KB
testcase_83	AC	351 ms 34,632 KB
testcase_84	AC	308 ms 33,272 KB
testcase_85	AC	28 ms 7,164 KB
testcase_86	AC	284 ms 33,144 KB
testcase_87	AC	402 ms 35,344 KB
testcase_88	AC	148 ms 18,760 KB
testcase_89	AC	380 ms 34,144 KB
testcase_90	AC	147 ms 18,584 KB
testcase_91	AC	408 ms 34,636 KB

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

コンパイルメッセージ

main.cpp: In member function ‘int CHT<K>::intersect(const CHT<K>::line&, int, int) [with K = long long int]’:
main.cpp:162:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
  162 |         return data[idx];
      |                ~~~~^
main.cpp:153:24: note: ‘idx’ was declared here
  153 |         int latest = 0,idx;
      |                        ^~~
main.cpp: In function ‘int main()’:
main.cpp:162:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
  162 |         return data[idx];
      |                ~~~~^
main.cpp:153:24: note: ‘idx’ was declared here
  153 |         int latest = 0,idx;
      |                        ^~~
main.cpp:162:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
  162 |         return data[idx];
      |                ~~~~^
main.cpp:153:24: note: ‘idx’ was declared here
  153 |         int latest = 0,idx;
      |                        ^~~
main.cpp:162:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
  162 |         return data[idx];
      |                ~~~~^
main.cpp:153:24: note: ‘idx’ was declared here
  153 |         int latest = 0,idx;
      |                        ^~~
main.cpp:162:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
  162 |         return data[idx];
      |                ~~~~^
main.cpp:153:24: note: ‘idx’ was declared here
  153 |         int latest = 0,idx;
      |                        ^~~
main.cpp:162:20: warning: ‘idx’ may be used uninitialized in this function [-Wmaybe-uninitialized]
  162 |         return data[idx];
      |                ~~~~^
main.cpp:153:24: note: ‘idx’ was declared here
  153 |         int latest = 0,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 vector<K> &sorted) : toK(sorted)
    {
        conv = 1;
        int n = toK.size();
        while(n >= dom) dom <<= 1;
        dom = dom * 2 + 1;
        data.resize(dom << 1);
        dom = n;
        // for(int i = 0; i < n; ++i) {
        //     toI[S[i]] = i;
        // }
    }

    int conv_toI(K x) {
        // if(conv) return toI[x];
        if(conv) return lower_bound(begin(toK),end(toK),x) - begin(toK);
        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];
    vector<ll> s;
    for(int i = 1; i <= N; i++) s.emb(i);
    CHT<ll> cht(s);
    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

結果

テストケース

コンパイルメッセージ

ソースコード