結果

問題 No.703 ゴミ拾い Easy
コンテスト
ユーザー te-sh
提出日時 2018-07-05 16:57:51
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 94 ms / 1,500 ms
コード長 1,252 bytes
コンパイル時間 918 ms
コンパイル使用メモリ 106,668 KB
実行使用メモリ 24,332 KB
最終ジャッジ日時 2024-06-13 01:23:03
合計ジャッジ時間 6,748 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.container, std.conv, std.math, std.range, std.typecons, std.stdio, std.string;

auto rdsp(){return readln.splitter;}
void pick(R,T)(ref R r,ref T t){t=r.front.to!T;r.popFront;}
void readV(T...)(ref T t){auto r=rdsp;foreach(ref v;t)pick(r,v);}
void readA(T)(size_t n,ref T[]t){t=new T[](n);auto r=rdsp;foreach(ref v;t)pick(r,v);}

void main()
{
  int n; readV(n);
  long[] a; readA(n, a);
  long[] x; readA(n, x);
  long[] y; readA(n, y);

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

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

  writeln(dp[n]);
}

struct ConvexHullTrickSimple(T)
{
  import std.typecons;

  alias Line = Tuple!(T, T);
  Line[] deq;
  int s = 0, t = 0;

  this(size_t n)
  {
    deq = new Line[](n);
  }

  auto add(T a, T b)
  {
    auto p = Line(a, b);
    while (s+1 < t && check(deq[t-2], deq[t-1], p)) --t;
    deq[t++] = p;
  }

  auto get(T x)
  {
    while (s+1 < t && f(deq[s], x) >= f(deq[s+1], x)) ++s;
    return f(deq[s], x);
  }

  auto check(Line p1, Line p2, Line p3)
  {
    return (p2[0]-p1[0])*(p3[1]-p2[1]) >= (p2[1]-p1[1])*(p3[0]-p2[0]);
  }

  auto f(Line p, T x)
  {
    return p[0]*x+p[1];
  }
}
0