結果
問題 | No.703 ゴミ拾い Easy |
ユーザー |
|
提出日時 | 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 |
ソースコード
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]; } }