結果

問題 No.703 ゴミ拾い Easy
ユーザー te-shte-sh
提出日時 2020-01-29 16:21:03
言語 D
(dmd 2.106.1)
結果
AC  
実行時間 416 ms / 1,500 ms
コード長 3,212 bytes
コンパイル時間 1,353 ms
コンパイル使用メモリ 145,004 KB
実行使用メモリ 59,904 KB
最終ジャッジ日時 2023-09-04 05:18:02
合計ジャッジ時間 11,760 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 1 ms
4,384 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 1 ms
4,380 KB
testcase_07 AC 1 ms
4,376 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 2 ms
4,376 KB
testcase_10 AC 2 ms
4,380 KB
testcase_11 AC 1 ms
4,376 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 2 ms
4,376 KB
testcase_14 AC 2 ms
4,376 KB
testcase_15 AC 3 ms
4,380 KB
testcase_16 AC 2 ms
4,376 KB
testcase_17 AC 2 ms
4,380 KB
testcase_18 AC 2 ms
4,380 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 2 ms
4,376 KB
testcase_21 AC 2 ms
4,376 KB
testcase_22 AC 2 ms
4,380 KB
testcase_23 AC 2 ms
4,376 KB
testcase_24 AC 353 ms
59,904 KB
testcase_25 AC 313 ms
43,484 KB
testcase_26 AC 393 ms
43,420 KB
testcase_27 AC 344 ms
43,396 KB
testcase_28 AC 416 ms
59,848 KB
testcase_29 AC 364 ms
43,588 KB
testcase_30 AC 299 ms
43,456 KB
testcase_31 AC 326 ms
59,892 KB
testcase_32 AC 348 ms
59,812 KB
testcase_33 AC 374 ms
59,872 KB
testcase_34 AC 184 ms
28,556 KB
testcase_35 AC 180 ms
22,664 KB
testcase_36 AC 185 ms
28,288 KB
testcase_37 AC 184 ms
29,056 KB
testcase_38 AC 182 ms
23,660 KB
testcase_39 AC 181 ms
23,948 KB
testcase_40 AC 185 ms
29,384 KB
testcase_41 AC 185 ms
35,424 KB
testcase_42 AC 185 ms
28,812 KB
testcase_43 AC 184 ms
29,084 KB
testcase_44 AC 1 ms
4,380 KB
testcase_45 AC 1 ms
4,376 KB
testcase_46 AC 223 ms
30,432 KB
testcase_47 AC 216 ms
27,212 KB
testcase_48 AC 3 ms
4,380 KB
testcase_49 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// URL: https://yukicoder.me/problems/no/703

import std.algorithm, std.array, std.container, std.math, std.range, std.typecons, std.string;

version(unittest) {} else
void main()
{
  int n; io.getV(n);
  long[] a; io.getA(n, a);
  long[] x; io.getA(n, x);
  long[] y; io.getA(n, y);

  auto dp = new long[](n+1), cht = ConvexHullTrick!long();

  foreach (i; 0..n) {
    cht.add(-2*x[i], dp[i]+x[i]^^2+y[i]^^2);
    dp[i+1] = cht.query(a[i]) + a[i]^^2;
  }

  io.put(dp[n]);
}

struct ConvexHullTrick(T, alias comp = "a>=b")
{
  import std.functional;
  alias compFun = binaryFun!comp;

  void add(T a, T b)
  {
    auto line = Line(a, b);
    while (lines.length >= 2 && check(lines[$-2], lines[$-1], line))
      lines = lines[0..$-1];
    lines ~= line;
  }

  pure T query(T x)
  {
    ptrdiff_t low = -1, high = lines.length-1;
    while (high-low > 1) {
      auto mid = (high+low)/2;
      (compFun(lines[mid].f(x), lines[mid+1].f(x)) ? low : high) = mid;
    }
    return lines[high].f(x);
  }

  private
  {
    struct Line
    {
      T a, b;
      auto opCmp(Line r)
      {
        if (a < r.a) return -1;
        else if (a > r.a) return 1;
        if (b < r.b) return -1;
        else if (b > r.b) return 1;
        return 0;
      }
      auto f(T x) { return a*x+b; }
    }
    Line[] lines;

    bool check(Line l1, Line l2, Line l3)
    {
      if (l1 < l3) swap(l1, l3);
      return (l3.b-l2.b)*(l2.a-l1.a) >= (l2.b-l1.b)*(l3.a-l2.a);
    }
  }
}

auto io = IO!()();
import std.stdio;
struct IO(string floatFormat = "%.10f", string delimiter = " ", alias IN = stdin, alias OUT = stdout)
{
  import std.conv, std.format, std.meta, std.traits;

  auto getV(T...)(ref T v) { foreach (ref w; v) get(w); }
  auto getA(T)(size_t n, ref T v) if (hasAssignableElements!T)
  { v = new T(n); foreach (ref w; v) get(w); }
  auto getC(T...)(size_t n, ref T v)
  if (allSatisfy!(hasAssignableElements, T))
  {
    foreach (ref w; v) w = new typeof(w)(n);
    foreach (i; 0..n) foreach (ref w; v) get(w[i]);
  }
  auto getM(T)(size_t r, size_t c, ref T v)
  if (hasAssignableElements!T && hasAssignableElements!(ElementType!T))
  {
    v = new T(r); foreach (ref w; v) getA(c, w);
  }
  template getS(E...)
  {
    auto getS(T)(size_t n, ref T v)
    {
      v = new T(n); foreach (ref w; v) foreach (e; E) mixin("get(w."~e~");");
    }
  }

  auto put(bool flush = false, T...)(T v)
  {
    foreach (i, w; v) { putA(w); if (i < v.length-1) OUT.write(delimiter); }
    OUT.writeln;
    static if (flush) OUT.flush();
  }
  auto putB(S, T)(bool c, S t, T f) { if (c) put(t); else put(f); }
  auto putRaw(T...)(T v) { OUT.write(v); OUT.writeln; }

  private
  {
    dchar[] buf;
    auto sp = (new dchar[](0)).splitter;
    void nextLine() { IN.readln(buf); sp = buf.splitter; }
    auto get(T)(ref T v) { if (sp.empty) nextLine(); v = sp.front.to!T; sp.popFront(); }

    auto putR(T)(T v)
    {
      auto w = v;
      while (!w.empty) { putA(w.front); w.popFront(); if (!w.empty) OUT.write(delimiter); }
    }
    auto putA(T)(T v)
    {
      static if (isInputRange!T && !isSomeString!T) putR(v);
      else if (isFloatingPoint!T) OUT.write(format(floatFormat, v));
      else OUT.write(v);
    }
  }
}
0