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

提出ソース

結果

問題	No.409 ダイエット
ユーザー	jell
提出日時	2018-10-23 18:52:12
言語	C++11 (gcc 11.4.0)
結果	AC
実行時間	440 ms / 2,000 ms
コード長	6,687 bytes
コンパイル時間	1,587 ms
コンパイル使用メモリ	169,604 KB
実行使用メモリ	34,792 KB
最終ジャッジ日時	2024-11-19 04:32:22
合計ジャッジ時間	13,706 ms
ジャッジサーバーID （参考情報）	judge4 / judge5

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
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	2 ms 5,248 KB
testcase_04	AC	2 ms 5,248 KB
testcase_05	AC	2 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	2 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	2 ms 5,248 KB
testcase_23	AC	2 ms 5,248 KB
testcase_24	AC	2 ms 5,248 KB
testcase_25	AC	2 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	2 ms 5,248 KB
testcase_34	AC	2 ms 5,248 KB
testcase_35	AC	4 ms 5,248 KB
testcase_36	AC	4 ms 5,248 KB
testcase_37	AC	4 ms 5,248 KB
testcase_38	AC	4 ms 5,248 KB
testcase_39	AC	4 ms 5,248 KB
testcase_40	AC	4 ms 5,248 KB
testcase_41	AC	4 ms 5,248 KB
testcase_42	AC	4 ms 5,248 KB
testcase_43	AC	4 ms 5,248 KB
testcase_44	AC	4 ms 5,248 KB
testcase_45	AC	4 ms 5,248 KB
testcase_46	AC	4 ms 5,248 KB
testcase_47	AC	4 ms 5,248 KB
testcase_48	AC	4 ms 5,248 KB
testcase_49	AC	4 ms 5,248 KB
testcase_50	AC	4 ms 5,248 KB
testcase_51	AC	4 ms 5,248 KB
testcase_52	AC	4 ms 5,248 KB
testcase_53	AC	4 ms 5,248 KB
testcase_54	AC	4 ms 5,248 KB
testcase_55	AC	182 ms 18,944 KB
testcase_56	AC	218 ms 34,560 KB
testcase_57	AC	386 ms 34,792 KB
testcase_58	AC	139 ms 12,416 KB
testcase_59	AC	226 ms 19,712 KB
testcase_60	AC	130 ms 11,904 KB
testcase_61	AC	331 ms 21,504 KB
testcase_62	AC	371 ms 34,176 KB
testcase_63	AC	320 ms 21,504 KB
testcase_64	AC	165 ms 18,816 KB
testcase_65	AC	358 ms 34,048 KB
testcase_66	AC	404 ms 34,432 KB
testcase_67	AC	276 ms 20,736 KB
testcase_68	AC	221 ms 19,712 KB
testcase_69	AC	322 ms 21,504 KB
testcase_70	AC	312 ms 20,992 KB
testcase_71	AC	157 ms 12,544 KB
testcase_72	AC	440 ms 34,560 KB
testcase_73	AC	360 ms 21,632 KB
testcase_74	AC	207 ms 19,456 KB
testcase_75	AC	357 ms 34,048 KB
testcase_76	AC	243 ms 19,968 KB
testcase_77	AC	153 ms 12,288 KB
testcase_78	AC	183 ms 19,072 KB
testcase_79	AC	140 ms 12,032 KB
testcase_80	AC	367 ms 34,304 KB
testcase_81	AC	343 ms 21,632 KB
testcase_82	AC	255 ms 20,224 KB
testcase_83	AC	286 ms 20,352 KB
testcase_84	AC	234 ms 19,584 KB
testcase_85	AC	20 ms 5,248 KB
testcase_86	AC	215 ms 19,584 KB
testcase_87	AC	310 ms 21,120 KB
testcase_88	AC	118 ms 11,776 KB
testcase_89	AC	283 ms 20,224 KB
testcase_90	AC	111 ms 11,520 KB
testcase_91	AC	300 ms 20,608 KB

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

コンパイルメッセージ

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

ソースコード

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; };

    vector<K> conv;
    vector<line> data;
    int dom = 1;    // the width of the domain [0,dom)
    int clock = 0;  // the number of total addition
    bool consec = 0;    // whether the domain is consecutive integers or not

    // for the interval [l,r)
    CHT(int l, int r) {
        consec = 1;
        while(r - l >= dom) dom <<= 1;
        data.resize(dom << 1);
        conv.resize(dom = r - l);
        for(int i = l; i < r; ++i) conv[i - l] = i;
    }

    CHT(const vector<K> &sorted) : conv(sorted)
    {
        int n = conv.size();
        while(n >= dom) dom <<= 1;
        data.resize(dom << 1);
        dom = n;
    }

    int to_idx(K x) {
        if(consec) return x - conv[0];
        return lower_bound(begin(conv),end(conv),x) - begin(conv);
    }

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

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

    bool is_lower(const line& ln, int i) { return query(conv[i]) > getval(ln,conv[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(1,N + 1);
    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

結果

テストケース

コンパイルメッセージ

ソースコード