結果

問題 No.703 ゴミ拾い Easy
ユーザー EmKjpEmKjp
提出日時 2020-02-20 22:42:41
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 253 ms / 1,500 ms
コード長 7,961 bytes
コンパイル時間 2,469 ms
コンパイル使用メモリ 110,884 KB
実行使用メモリ 35,424 KB
最終ジャッジ日時 2024-06-10 18:48:39
合計ジャッジ時間 9,375 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 23 ms
21,688 KB
testcase_01 AC 25 ms
23,584 KB
testcase_02 AC 24 ms
19,764 KB
testcase_03 AC 25 ms
21,688 KB
testcase_04 AC 25 ms
23,536 KB
testcase_05 AC 24 ms
25,876 KB
testcase_06 AC 23 ms
23,920 KB
testcase_07 AC 23 ms
23,664 KB
testcase_08 AC 22 ms
23,588 KB
testcase_09 AC 21 ms
25,756 KB
testcase_10 AC 22 ms
25,636 KB
testcase_11 AC 21 ms
21,812 KB
testcase_12 AC 22 ms
23,448 KB
testcase_13 AC 22 ms
23,712 KB
testcase_14 AC 22 ms
23,604 KB
testcase_15 AC 23 ms
23,792 KB
testcase_16 AC 23 ms
25,904 KB
testcase_17 AC 23 ms
23,784 KB
testcase_18 AC 23 ms
23,532 KB
testcase_19 AC 23 ms
23,792 KB
testcase_20 AC 21 ms
21,424 KB
testcase_21 AC 22 ms
23,916 KB
testcase_22 AC 22 ms
23,664 KB
testcase_23 AC 25 ms
23,604 KB
testcase_24 AC 240 ms
35,032 KB
testcase_25 AC 236 ms
35,156 KB
testcase_26 AC 243 ms
33,360 KB
testcase_27 AC 245 ms
32,980 KB
testcase_28 AC 253 ms
33,108 KB
testcase_29 AC 243 ms
32,852 KB
testcase_30 AC 245 ms
32,988 KB
testcase_31 AC 244 ms
35,160 KB
testcase_32 AC 251 ms
35,028 KB
testcase_33 AC 240 ms
35,276 KB
testcase_34 AC 201 ms
35,168 KB
testcase_35 AC 195 ms
35,028 KB
testcase_36 AC 191 ms
33,060 KB
testcase_37 AC 192 ms
35,168 KB
testcase_38 AC 184 ms
35,160 KB
testcase_39 AC 184 ms
33,120 KB
testcase_40 AC 182 ms
35,044 KB
testcase_41 AC 183 ms
35,424 KB
testcase_42 AC 181 ms
33,372 KB
testcase_43 AC 184 ms
35,296 KB
testcase_44 AC 22 ms
23,592 KB
testcase_45 AC 22 ms
23,856 KB
testcase_46 AC 238 ms
32,888 KB
testcase_47 AC 231 ms
35,284 KB
testcase_48 AC 23 ms
25,524 KB
testcase_49 AC 24 ms
23,592 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;
using System.Collections;
using System.Collections.Generic;
using System.Globalization;
using System.IO;
using System.Linq;

using E = System.Linq.Enumerable;
using CHTType = System.Int64;
using System.Text;

internal partial class Solver {
    public void Run() {
        var n = ni();
        var a = ni(n);
        var x = ni(n);
        var y = ni(n);
        var dp = new long[n + 1];
        var inf = long.MaxValue >> 10;
        for (int i = 0; i <= n; i++) {
            dp[i] = inf;
        }
        dp[0] = 0;

        var cht = new ConvexHullTrick(true);

        for (int i = 0; i < n; i++) {
            //for (int j = 0; j <= i; j++) {
            //v = Math.Min(v, dp[j] + 1L * (x[j] - a[i]) * (x[j] - a[i]) + 1L * y[j] * y[j]);
            //}
            // min_j dp[j] + (x[j] - a[i]) ^ 2 + y[j] ^ 2
            // min_j dp[j] + x[j]^2 - 2 * a[i] * x[j] + a[i]^2 + y[j] ^ 2
            // min_j { dp[j] + x[j]^2 - 2 * a[i] * x[j] + y[j]^2 } + a[i]^2
            cht.Add(-2 * x[i], dp[i] + 1L * x[i] * x[i] + 1L * y[i] * y[i]);
            dp[i + 1] = cht.Get(a[i]) + 1L * a[i] * a[i];
        }

        cout.WriteLine(dp[n]);
    }

}

public class ConvexHullTrick {
    private class Line {
        public readonly CHTType Slope;
        public readonly CHTType Intercept;

        public Line(CHTType slope, CHTType intercept) {
            Slope = slope;
            Intercept = intercept;
        }

        public CHTType Evaluate(CHTType x) {
            return Slope * x + Intercept;
        }
    }

    private readonly List<Line> _envelope = new List<Line>();
    private readonly bool _isMin;

    public ConvexHullTrick(bool isMin = true) {
        _isMin = isMin;
    }

    // Add ax + b
    // <code>a</code> must be monotone decreasing
    public void Add(CHTType a, CHTType b) {
        if (!_isMin) {
            a = -a;
            b = -b;
        }
        var line = new Line(a, b);

        if (_envelope.Count > 0 && _envelope[_envelope.Count - 1].Slope < line.Slope) {
            throw new InvalidOperationException("Slope must be monotone " + (_isMin ? "decreasing" : "increasing"));
        }

        while (_envelope.Count >= 2 &&
            IsUnnecessary(_envelope[_envelope.Count - 2], _envelope[_envelope.Count - 1], line)) {
            _envelope.RemoveAt(_envelope.Count - 1);
        }
        _envelope.Add(line);
    }

    private static bool IsUnnecessary(Line l1, Line l2, Line l3) {
        return
            (l3.Intercept - l2.Intercept) * (l2.Slope - l1.Slope) >=
            (l2.Intercept - l1.Intercept) * (l3.Slope - l2.Slope);
    }

    public CHTType Get(CHTType x) {
        var best = CHTType.MaxValue;
        if (_envelope.Count == 1) {
            best = _envelope[0].Evaluate(x);
        } else {
            int lower = 0, upper = _envelope.Count - 2;
            while (lower <= upper) {
                var mid = (lower + upper) >> 1;
                var y0 = _envelope[mid].Evaluate(x);
                var y1 = _envelope[mid + 1].Evaluate(x);
                if (y0 > y1) {
                    best = y1;
                    lower = mid + 1;
                } else {
                    best = y0;
                    upper = mid - 1;
                }
            }
        }
        return _isMin ? best : -best;
    }
}

// PREWRITEN CODE BEGINS FROM HERE
internal partial class Solver : Scanner {
    public static void Main(string[] args) {
#if LOCAL
        byte[] inputBuffer = new byte[1000000];
        var inputStream = Console.OpenStandardInput(inputBuffer.Length);
        using (var reader = new StreamReader(inputStream, Console.InputEncoding, false, inputBuffer.Length)) {
            Console.SetIn(reader);
            new Solver(Console.In, Console.Out).Run();
        }
#else
        Console.SetOut(new StreamWriter(Console.OpenStandardOutput()) { AutoFlush = false });
        new Solver(Console.In, Console.Out).Run();
        Console.Out.Flush();
#endif
    }

#pragma warning disable IDE0052
    private readonly TextReader cin;
    private readonly TextWriter cout;
#pragma warning restore IDE0052

    public Solver(TextReader reader, TextWriter writer)
        : base(reader) {
        cin = reader;
        cout = writer;
    }
    public Solver(string input, TextWriter writer)
        : this(new StringReader(input), writer) {
    }

#pragma warning disable IDE1006
#pragma warning disable IDE0051
    private int ni() { return NextInt(); }
    private int[] ni(int n) { return NextIntArray(n); }
    private long nl() { return NextLong(); }
    private long[] nl(int n) { return NextLongArray(n); }
    private double nd() { return NextDouble(); }
    private double[] nd(int n) { return NextDoubleArray(n); }
    private string ns() { return Next(); }
    private string[] ns(int n) { return NextArray(n); }
#pragma warning restore IDE1006
#pragma warning restore IDE0051
}

internal static class LinqPadExtension {
    public static T Dump<T>(this T obj) {
#if LOCAL
        return LINQPad.Extensions.Dump(obj);
#else
        return obj;
#endif
    }
}

public class Scanner {
    private readonly TextReader Reader;
    private readonly CultureInfo ci = CultureInfo.InvariantCulture;

    private readonly char[] buffer = new char[2 * 1024];
    private int cursor = 0, length = 0;
    private string Token;
    private readonly StringBuilder sb = new StringBuilder(1024);

    public Scanner()
        : this(Console.In) {
    }

    public Scanner(TextReader reader) {
        Reader = reader;
    }

    public int NextInt() { return checked((int)NextLong()); }
    public long NextLong() {
        var s = Next();
        long r = 0;
        int i = 0;
        bool negative = false;
        if (s[i] == '-') {
            negative = true;
            i++;
        }
        for (; i < s.Length; i++) {
            r = r * 10 + (s[i] - '0');
#if DEBUG
            if (!char.IsDigit(s[i])) throw new FormatException();
#endif
        }
        return negative ? -r : r;
    }
    public double NextDouble() { return double.Parse(Next(), ci); }
    public string[] NextArray(int size) {
        string[] array = new string[size];
        for (int i = 0; i < size; i++) {
            array[i] = Next();
        }

        return array;
    }
    public int[] NextIntArray(int size) {
        int[] array = new int[size];
        for (int i = 0; i < size; i++) {
            array[i] = NextInt();
        }

        return array;
    }

    public long[] NextLongArray(int size) {
        long[] array = new long[size];
        for (int i = 0; i < size; i++) {
            array[i] = NextLong();
        }

        return array;
    }

    public double[] NextDoubleArray(int size) {
        double[] array = new double[size];
        for (int i = 0; i < size; i++) {
            array[i] = NextDouble();
        }

        return array;
    }

    public string Next() {
        if (Token == null) {
            if (!StockToken()) {
                throw new InvalidOperationException();
            }
        }
        var token = Token;
        Token = null;
        return token;
    }

    public bool HasNext() {
        if (Token != null) {
            return true;
        }

        return StockToken();
    }

    private bool StockToken() {
        while (true) {
            sb.Clear();
            while (true) {
                if (cursor >= length) {
                    cursor = 0;
                    if ((length = Reader.Read(buffer, 0, buffer.Length)) <= 0) {
                        break;
                    }
                }
                var c = buffer[cursor++];
                if (33 <= c && c <= 126) {
                    sb.Append(c);
                } else {
                    if (sb.Length > 0) break;
                }
            }

            if (sb.Length > 0) {
                Token = sb.ToString();
                return true;
            }

            return false;
        }
    }
}
0