#312019 (D) No.409 ダイエット

提出ソース

結果

問題	No.409 ダイエット
ユーザー	👑 hos.lyric
提出日時	2019-01-22 23:00:26
言語	D (dmd 2.106.1)
結果	AC
実行時間	58 ms / 2,000 ms
コード長	2,802 bytes
コンパイル時間	1,032 ms
コンパイル使用メモリ	122,764 KB
実行使用メモリ	15,648 KB
最終ジャッジ日時	2024-06-13 03:31:34
合計ジャッジ時間	6,076 ms
ジャッジサーバーID （参考情報）	judge5 / judge2

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 6,812 KB
testcase_01	AC	1 ms 6,812 KB
testcase_02	AC	1 ms 6,944 KB
testcase_03	AC	1 ms 6,940 KB
testcase_04	AC	2 ms 6,940 KB
testcase_05	AC	1 ms 6,944 KB
testcase_06	AC	1 ms 6,940 KB
testcase_07	AC	1 ms 6,944 KB
testcase_08	AC	1 ms 6,944 KB
testcase_09	AC	1 ms 6,944 KB
testcase_10	AC	1 ms 6,944 KB
testcase_11	AC	1 ms 6,944 KB
testcase_12	AC	1 ms 6,944 KB
testcase_13	AC	1 ms 6,944 KB
testcase_14	AC	1 ms 6,944 KB
testcase_15	AC	1 ms 6,940 KB
testcase_16	AC	1 ms 6,940 KB
testcase_17	AC	2 ms 6,940 KB
testcase_18	AC	1 ms 6,944 KB
testcase_19	AC	1 ms 6,940 KB
testcase_20	AC	1 ms 6,944 KB
testcase_21	AC	2 ms 6,940 KB
testcase_22	AC	1 ms 6,944 KB
testcase_23	AC	2 ms 6,944 KB
testcase_24	AC	2 ms 6,944 KB
testcase_25	AC	1 ms 6,944 KB
testcase_26	AC	1 ms 6,940 KB
testcase_27	AC	1 ms 6,948 KB
testcase_28	AC	1 ms 6,940 KB
testcase_29	AC	1 ms 6,940 KB
testcase_30	AC	1 ms 6,944 KB
testcase_31	AC	1 ms 6,944 KB
testcase_32	AC	1 ms 6,944 KB
testcase_33	AC	1 ms 6,944 KB
testcase_34	AC	1 ms 6,940 KB
testcase_35	AC	2 ms 6,940 KB
testcase_36	AC	2 ms 6,940 KB
testcase_37	AC	2 ms 6,944 KB
testcase_38	AC	2 ms 6,944 KB
testcase_39	AC	2 ms 6,944 KB
testcase_40	AC	2 ms 6,940 KB
testcase_41	AC	2 ms 6,940 KB
testcase_42	AC	2 ms 6,940 KB
testcase_43	AC	2 ms 6,940 KB
testcase_44	AC	2 ms 6,944 KB
testcase_45	AC	2 ms 6,944 KB
testcase_46	AC	2 ms 6,944 KB
testcase_47	AC	2 ms 6,944 KB
testcase_48	AC	2 ms 6,940 KB
testcase_49	AC	2 ms 6,944 KB
testcase_50	AC	2 ms 6,940 KB
testcase_51	AC	2 ms 6,944 KB
testcase_52	AC	2 ms 6,940 KB
testcase_53	AC	2 ms 6,940 KB
testcase_54	AC	2 ms 6,940 KB
testcase_55	AC	20 ms 8,848 KB
testcase_56	AC	38 ms 15,384 KB
testcase_57	AC	42 ms 15,616 KB
testcase_58	AC	19 ms 7,676 KB
testcase_59	AC	28 ms 12,108 KB
testcase_60	AC	17 ms 7,160 KB
testcase_61	AC	39 ms 14,728 KB
testcase_62	AC	40 ms 15,036 KB
testcase_63	AC	38 ms 14,428 KB
testcase_64	AC	20 ms 8,044 KB
testcase_65	AC	42 ms 14,804 KB
testcase_66	AC	42 ms 15,548 KB
testcase_67	AC	31 ms 13,280 KB
testcase_68	AC	27 ms 12,064 KB
testcase_69	AC	38 ms 14,248 KB
testcase_70	AC	48 ms 13,864 KB
testcase_71	AC	27 ms 8,388 KB
testcase_72	AC	58 ms 15,648 KB
testcase_73	AC	55 ms 14,780 KB
testcase_74	AC	34 ms 11,692 KB
testcase_75	AC	55 ms 14,796 KB
testcase_76	AC	38 ms 12,448 KB
testcase_77	AC	27 ms 8,196 KB
testcase_78	AC	30 ms 11,164 KB
testcase_79	AC	25 ms 7,780 KB
testcase_80	AC	55 ms 15,024 KB
testcase_81	AC	54 ms 14,760 KB
testcase_82	AC	40 ms 12,636 KB
testcase_83	AC	43 ms 13,116 KB
testcase_84	AC	37 ms 11,904 KB
testcase_85	AC	5 ms 5,376 KB
testcase_86	AC	36 ms 11,712 KB
testcase_87	AC	50 ms 13,872 KB
testcase_88	AC	21 ms 6,752 KB
testcase_89	AC	42 ms 12,816 KB
testcase_90	AC	20 ms 6,504 KB
testcase_91	AC	46 ms 13,124 KB

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

ソースコード

raw source code

import std.conv, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.container, std.math, std.numeric, std.range, std.regex, std.typecons;
import core.bitop;

class EOFException : Throwable { this() { super("EOF"); } }
string[] tokens;
string readToken() { for (; tokens.empty; ) { if (stdin.eof) { throw new EOFException; } tokens = readln.split; } auto token = tokens.front; tokens.popFront; return token; }
int readInt() { return readToken.to!int; }
long readLong() { return readToken.to!long; }
real readReal() { return readToken.to!real; }

void chmin(T)(ref T t, in T f) { if (t > f) t = f; }
void chmax(T)(ref T t, in T f) { if (t < f) t = f; }

int binarySearch(T)(in T[] as, in bool delegate(T) test) { int low = -1, upp = cast(int)(as.length); for (; low + 1 < upp; ) { int mid = (low + upp) >> 1; (test(as[mid]) ? low : upp) = mid; } return upp; }
int lowerBound(T)(in T[] as, in T val) { return as.binarySearch((T a) => (a < val)); }
int upperBound(T)(in T[] as, in T val) { return as.binarySearch((T a) => (a <= val)); }


/*
  dp[i] = min{ dp[j] - A (i - j - 1) + B (i - j - 1) (i - j) / 2 + D[i] | 0 <= j < i }
        = min{ (dp[j] + A j + B j (j + 1) / 2) - (B j) i | 0 <= j < i } - A (i - 1) + B i (i - 1) / 2 + D[i]
*/

int N;
long A, B, W;
long[] D;

void main() {
  try {
    for (; ; ) {
      N = readInt();
      A = readLong();
      B = readLong();
      W = readLong();
      D = new long[N + 2];
      foreach (i; 1 .. N + 1) {
        D[i] = readLong();
      }
      
      debug {
        auto dp = new long[N + 2];
        dp[0] = W;
        foreach (i; 1 .. N + 2) {
          dp[i] = iota(0, i).map!(j => (dp[j] + A * j + B * j * (j + 1) / 2) - (B * j) * i).minElement - A * (i - 1) + B * i * (i - 1) / 2 + D[i];
        }
        writeln(dp);
      }
      
      long ans;
      auto q = new Tuple!(long, long)[N + 2];
      int qHead, qTail;
      foreach (i; 0 .. N + 2) {
        long x;
        if (i == 0) {
          x = W;
        } else {
          for (; qTail - qHead >= 2 && q[qHead][0] - q[qHead][1] * i >= q[qHead + 1][0] - q[qHead + 1][1] * i; ++qHead) {}
          x = q[qHead][0] - q[qHead][1] * i - A * (i - 1) + B * i * (i - 1) / 2 + D[i];
        }
        debug {
          writeln(i, " ", x);
        }
        if (i == N + 1) {
          ans = x;
        } else {
          const p = Tuple!(long, long)(x + A * i + B * i * (i + 1) / 2, B * i);
          // ([-1][0] - [-2][0]) / ([-1][1] - [-2][1]) >= (p[0] - [-2][0]) / (p[1] - [-2][1])
          for (; qTail - qHead >= 2 && (q[qTail - 1][0] - q[qTail - 2][0]) * (p[1] - q[qTail - 2][1]) >= (p[0] - q[qTail - 2][0]) * (q[qTail - 1][1] - q[qTail - 2][1]); --qTail) {}
          q[qTail++] = p;
        }
      }
      writeln(ans);
    }
  } catch (EOFException e) {
  }
}

yukicoder

結果

テストケース

ソースコード