結果

問題 No.1069 電柱 / Pole (Hard)
ユーザー 👑 hos.lyrichos.lyric
提出日時 2020-05-29 22:24:19
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 5,412 bytes
コンパイル時間 3,622 ms
コンパイル使用メモリ 274,256 KB
実行使用メモリ 52,440 KB
最終ジャッジ日時 2023-09-04 07:59:26
合計ジャッジ時間 10,276 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 1,276 ms
48,216 KB
testcase_05 WA -
testcase_06 AC 48 ms
21,140 KB
testcase_07 AC 55 ms
20,788 KB
testcase_08 AC 24 ms
11,696 KB
testcase_09 AC 53 ms
21,488 KB
testcase_10 AC 56 ms
27,432 KB
testcase_11 AC 4 ms
5,256 KB
testcase_12 AC 5 ms
4,504 KB
testcase_13 AC 18 ms
4,856 KB
testcase_14 AC 9 ms
4,812 KB
testcase_15 AC 9 ms
5,716 KB
testcase_16 AC 23 ms
11,732 KB
testcase_17 WA -
testcase_18 AC 18 ms
7,688 KB
testcase_19 AC 38 ms
12,508 KB
testcase_20 AC 6 ms
4,380 KB
testcase_21 AC 9 ms
4,380 KB
testcase_22 AC 4 ms
4,376 KB
testcase_23 AC 23 ms
4,376 KB
testcase_24 WA -
testcase_25 AC 95 ms
38,060 KB
testcase_26 AC 81 ms
37,336 KB
testcase_27 AC 87 ms
37,828 KB
testcase_28 AC 83 ms
36,936 KB
testcase_29 AC 17 ms
11,212 KB
testcase_30 AC 4 ms
4,380 KB
testcase_31 AC 2 ms
4,380 KB
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 AC 2 ms
4,380 KB
testcase_37 AC 6 ms
4,380 KB
testcase_38 AC 6 ms
4,380 KB
testcase_39 AC 2 ms
4,380 KB
testcase_40 AC 4 ms
4,380 KB
testcase_41 AC 6 ms
4,376 KB
testcase_42 AC 6 ms
4,376 KB
testcase_43 AC 6 ms
4,380 KB
testcase_44 AC 31 ms
12,244 KB
testcase_45 AC 13 ms
7,320 KB
testcase_46 AC 22 ms
7,644 KB
testcase_47 AC 33 ms
12,988 KB
testcase_48 AC 61 ms
26,784 KB
testcase_49 AC 34 ms
16,520 KB
testcase_50 AC 4 ms
4,376 KB
testcase_51 AC 45 ms
20,828 KB
testcase_52 AC 6 ms
4,384 KB
testcase_53 AC 54 ms
26,164 KB
testcase_54 WA -
testcase_55 AC 2 ms
4,376 KB
testcase_56 AC 1 ms
4,380 KB
testcase_57 WA -
testcase_58 AC 2 ms
4,380 KB
testcase_59 WA -
testcase_60 WA -
testcase_61 AC 2 ms
4,380 KB
testcase_62 AC 1 ms
4,384 KB
testcase_63 AC 1 ms
4,376 KB
testcase_64 WA -
testcase_65 AC 1 ms
4,380 KB
testcase_66 WA -
testcase_67 AC 1 ms
4,384 KB
testcase_68 WA -
testcase_69 AC 80 ms
40,072 KB
testcase_70 AC 92 ms
50,332 KB
testcase_71 AC 91 ms
50,052 KB
testcase_72 AC 23 ms
12,776 KB
testcase_73 AC 38 ms
21,236 KB
testcase_74 AC 94 ms
51,812 KB
testcase_75 AC 93 ms
51,196 KB
testcase_76 AC 67 ms
31,432 KB
testcase_77 AC 64 ms
28,148 KB
testcase_78 AC 3 ms
4,380 KB
testcase_79 AC 19 ms
4,896 KB
testcase_80 AC 38 ms
13,696 KB
testcase_81 AC 9 ms
6,612 KB
testcase_82 AC 14 ms
4,604 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