結果

問題 No.1069 電柱 / Pole (Hard)
ユーザー 👑 hos.lyrichos.lyric
提出日時 2020-05-29 22:28:43
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 5,418 bytes
コンパイル時間 4,527 ms
コンパイル使用メモリ 282,548 KB
実行使用メモリ 51,644 KB
最終ジャッジ日時 2024-06-22 07:14:53
合計ジャッジ時間 9,849 ms
ジャッジサーバーID
(参考情報)
judge3 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 1 ms
6,940 KB
testcase_02 AC 1 ms
6,940 KB
testcase_03 AC 1 ms
6,944 KB
testcase_04 AC 1,093 ms
47,332 KB
testcase_05 WA -
testcase_06 AC 46 ms
20,240 KB
testcase_07 AC 49 ms
20,624 KB
testcase_08 AC 21 ms
11,424 KB
testcase_09 AC 47 ms
21,424 KB
testcase_10 AC 52 ms
27,708 KB
testcase_11 AC 4 ms
6,940 KB
testcase_12 AC 3 ms
6,940 KB
testcase_13 AC 15 ms
6,940 KB
testcase_14 AC 7 ms
6,944 KB
testcase_15 AC 8 ms
6,940 KB
testcase_16 AC 22 ms
11,324 KB
testcase_17 WA -
testcase_18 AC 17 ms
7,244 KB
testcase_19 AC 33 ms
11,724 KB
testcase_20 AC 5 ms
6,940 KB
testcase_21 AC 8 ms
6,944 KB
testcase_22 AC 3 ms
6,944 KB
testcase_23 AC 21 ms
6,940 KB
testcase_24 WA -
testcase_25 AC 87 ms
37,196 KB
testcase_26 AC 77 ms
36,344 KB
testcase_27 AC 79 ms
36,692 KB
testcase_28 AC 76 ms
36,576 KB
testcase_29 AC 15 ms
10,304 KB
testcase_30 AC 2 ms
6,940 KB
testcase_31 AC 1 ms
6,940 KB
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 AC 1 ms
6,940 KB
testcase_37 AC 5 ms
6,940 KB
testcase_38 AC 4 ms
6,940 KB
testcase_39 AC 1 ms
6,940 KB
testcase_40 AC 3 ms
6,940 KB
testcase_41 AC 5 ms
6,944 KB
testcase_42 AC 5 ms
6,940 KB
testcase_43 AC 5 ms
6,940 KB
testcase_44 AC 27 ms
12,200 KB
testcase_45 AC 11 ms
6,940 KB
testcase_46 AC 19 ms
7,184 KB
testcase_47 AC 28 ms
12,964 KB
testcase_48 AC 57 ms
26,908 KB
testcase_49 AC 29 ms
16,744 KB
testcase_50 AC 3 ms
6,940 KB
testcase_51 AC 41 ms
19,840 KB
testcase_52 AC 5 ms
6,944 KB
testcase_53 AC 51 ms
25,692 KB
testcase_54 WA -
testcase_55 AC 2 ms
6,944 KB
testcase_56 AC 1 ms
6,944 KB
testcase_57 WA -
testcase_58 AC 1 ms
6,948 KB
testcase_59 WA -
testcase_60 WA -
testcase_61 AC 2 ms
6,944 KB
testcase_62 AC 1 ms
6,940 KB
testcase_63 AC 1 ms
6,944 KB
testcase_64 WA -
testcase_65 AC 1 ms
6,940 KB
testcase_66 WA -
testcase_67 AC 1 ms
6,944 KB
testcase_68 WA -
testcase_69 AC 77 ms
40,764 KB
testcase_70 AC 89 ms
49,784 KB
testcase_71 AC 87 ms
50,068 KB
testcase_72 AC 21 ms
12,584 KB
testcase_73 AC 34 ms
20,744 KB
testcase_74 AC 91 ms
51,108 KB
testcase_75 AC 89 ms
51,032 KB
testcase_76 AC 65 ms
29,720 KB
testcase_77 AC 58 ms
27,360 KB
testcase_78 AC 2 ms
6,944 KB
testcase_79 AC 18 ms
6,944 KB
testcase_80 AC 37 ms
13,092 KB
testcase_81 AC 8 ms
6,940 KB
testcase_82 AC 12 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import std.conv, std.functional, std.range, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.bitmanip, std.complex, std.container, std.math, std.mathspecial, std.numeric, 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; }

bool chmin(T)(ref T t, in T f) { if (t > f) { t = f; return true; } else { return false; } }
bool chmax(T)(ref T t, in T f) { if (t < f) { t = f; return true; } else { return false; } }

int binarySearch(alias pred, T)(in T[] as) { int lo = -1, hi = cast(int)(as.length); for (; lo + 1 < hi; ) { const mid = (lo + hi) >> 1; (unaryFun!pred(as[mid]) ? hi : lo) = mid; } return hi; }
int lowerBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a >= val)); }
int upperBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a > val)); }




void main() {
  try {
    for (; ; ) {
      const N = readInt();
      const M = readInt();
      const K = readInt();
      const S = readInt() - 1;
      const T = readInt() - 1;
      auto X = new real[N];
      auto Y = new real[N];
      foreach (u; 0 .. N) {
        X[u] = readReal();
        Y[u] = readReal();
      }
      auto A = new int[M];
      auto B = new int[M];
      foreach (i; 0 .. M) {
        A[i] = readInt() - 1;
        B[i] = readInt() - 1;
      }
      
      auto dist = new real[][](N, N);
      foreach (u; 0 .. N) foreach (v; 0 .. N) {
        dist[u][v] = sqrt((X[v] - X[u])^^2 + (Y[v] - Y[u])^^2);
      }
      auto on = new bool[N];
      on[] = true;
      auto graph = new RedBlackTree!int[N];
      foreach (u; 0 .. N) {
        graph[u] = new RedBlackTree!int;
      }
      foreach (i; 0 .. M) {
        graph[A[i]].insert(B[i]);
        graph[B[i]].insert(A[i]);
      }
      
      alias Path = Tuple!(real, "cost", int[], "us");
      Path dijkstra(int s, int t) {
        assert(on[s]);
        assert(on[t]);
        alias Entry = Tuple!(real, "c", int, "u");
        BinaryHeap!(Array!Entry, "a > b") que;
        auto ds = new real[N];
        ds[] = real.infinity;
        auto vis = new bool[N];
        auto prev = new int[N];
        prev[] = -1;
        ds[s] = 0.0L;
        que.insert(Entry(0.0L, s));
        for (; !que.empty; ) {
          const c = que.front.c;
          const u = que.front.u;
          que.removeFront;
          if (!vis[u]) {
            vis[u] = true;
            if (u == t) {
              break;
            }
            foreach (v; graph[u]) {
              if (on[v]) {
                if (!vis[v]) {
                  const cc = c + dist[u][v];
                  if (chmin(ds[v], cc)) {
                    prev[v] = u;
                    que.insert(Entry(cc, v));
                  }
                }
              }
            }
          }
        }
        Path ret;
        if (vis[t]) {
          ret.cost = ds[t];
          for (int u = t; u != -1; u = prev[u]) {
            ret.us ~= u;
          }
          ret.us.reverse;
        }
        return ret;
      }
      
      // https://en.wikipedia.org/wiki/Yen%27s_algorithm
      auto ans = new Path[K];
      ans[0] = dijkstra(S, T);
      assert(!ans[0].us.empty);
      debug {
        writeln("ans[0] = ", ans[0]);
      }
      auto bs = new RedBlackTree!Path;
      foreach (k; 1 .. K) {
        foreach (i; 0 .. cast(int)(ans[k - 1].us.length) - 1) {
          const spurNode = ans[k - 1].us[i];
          auto rootPath = ans[k - 1].us[0 .. i + 1];
          foreach (l; 0 .. k) {
            if (i + 1 <= ans[l].us.length && rootPath == ans[l].us[0 .. i + 1]) {
              const u = ans[l].us[i];
              const v = ans[l].us[i + 1];
              graph[u].removeKey(v);
              // graph[v].removeKey(u);
            }
          }
          foreach (u; rootPath) {
            if (u != S && u != spurNode) {
              on[u] = false;
            }
          }
          auto spurPath = dijkstra(spurNode, T);
          if (!spurPath.us.empty) {
            Path totalPath;
            totalPath.cost = spurPath.cost;
            foreach (j; 0 .. i) {
              totalPath.cost += dist[rootPath[j]][rootPath[j + 1]];
            }
            totalPath.us = rootPath ~ spurPath.us[1 .. $];
            bs.insert(totalPath);
          }
          foreach (l; 0 .. k) {
            if (i + 1 <= ans[l].us.length && rootPath == ans[l].us[0 .. i + 1]) {
              const u = ans[l].us[i];
              const v = ans[l].us[i + 1];
              graph[u].insert(v);
              // graph[v].insert(u);
            }
          }
          foreach (u; rootPath) {
            if (u != S && u != spurNode) {
              on[u] = true;
            }
          }
        }
        if (bs.empty) {
          break;
        }
        ans[k] = bs.front;
        bs.removeFront;
        debug {
          writefln("ans[%s] = %s", k, ans[k]);
        }
      }
      
      foreach (k; 0 .. K) {
        if (ans[k].us.empty) {
          writeln("-1");
        } else {
          writefln("%.12f", ans[k].cost);
        }
      }
    }
  } catch (EOFException e) {
  }
}
0