結果

問題 No.1069 電柱 / Pole (Hard)
ユーザー 👑 hos.lyrichos.lyric
提出日時 2020-05-29 22:32:41
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 5,425 bytes
コンパイル時間 4,893 ms
コンパイル使用メモリ 282,760 KB
実行使用メモリ 51,444 KB
最終ジャッジ日時 2024-06-22 07:15:22
合計ジャッジ時間 11,375 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,944 KB
testcase_02 AC 1 ms
6,944 KB
testcase_03 AC 1 ms
6,940 KB
testcase_04 AC 1,206 ms
47,076 KB
testcase_05 WA -
testcase_06 AC 48 ms
20,028 KB
testcase_07 AC 55 ms
21,724 KB
testcase_08 AC 23 ms
10,728 KB
testcase_09 AC 51 ms
20,772 KB
testcase_10 AC 55 ms
27,440 KB
testcase_11 AC 4 ms
6,940 KB
testcase_12 AC 4 ms
6,940 KB
testcase_13 AC 17 ms
6,940 KB
testcase_14 AC 8 ms
6,940 KB
testcase_15 AC 9 ms
6,940 KB
testcase_16 AC 22 ms
11,460 KB
testcase_17 WA -
testcase_18 AC 18 ms
7,036 KB
testcase_19 AC 37 ms
12,256 KB
testcase_20 AC 6 ms
6,940 KB
testcase_21 AC 9 ms
6,940 KB
testcase_22 AC 4 ms
6,944 KB
testcase_23 AC 22 ms
6,944 KB
testcase_24 WA -
testcase_25 AC 93 ms
36,204 KB
testcase_26 AC 80 ms
36,480 KB
testcase_27 AC 84 ms
36,652 KB
testcase_28 AC 81 ms
36,236 KB
testcase_29 AC 16 ms
10,264 KB
testcase_30 AC 3 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,944 KB
testcase_38 AC 5 ms
6,944 KB
testcase_39 AC 2 ms
6,944 KB
testcase_40 AC 4 ms
6,940 KB
testcase_41 AC 6 ms
6,940 KB
testcase_42 AC 5 ms
6,944 KB
testcase_43 AC 6 ms
6,944 KB
testcase_44 AC 30 ms
13,160 KB
testcase_45 AC 12 ms
6,940 KB
testcase_46 AC 21 ms
7,096 KB
testcase_47 AC 32 ms
12,732 KB
testcase_48 AC 62 ms
26,880 KB
testcase_49 AC 33 ms
16,828 KB
testcase_50 AC 4 ms
6,940 KB
testcase_51 AC 45 ms
20,480 KB
testcase_52 AC 5 ms
6,940 KB
testcase_53 AC 54 ms
25,480 KB
testcase_54 WA -
testcase_55 AC 1 ms
6,944 KB
testcase_56 AC 1 ms
6,944 KB
testcase_57 WA -
testcase_58 AC 1 ms
6,944 KB
testcase_59 WA -
testcase_60 WA -
testcase_61 AC 2 ms
6,940 KB
testcase_62 AC 1 ms
6,944 KB
testcase_63 AC 1 ms
6,944 KB
testcase_64 WA -
testcase_65 AC 2 ms
6,944 KB
testcase_66 WA -
testcase_67 AC 1 ms
6,944 KB
testcase_68 WA -
testcase_69 AC 76 ms
40,464 KB
testcase_70 AC 93 ms
49,892 KB
testcase_71 AC 91 ms
50,852 KB
testcase_72 AC 22 ms
12,616 KB
testcase_73 AC 36 ms
20,200 KB
testcase_74 AC 93 ms
51,132 KB
testcase_75 AC 93 ms
51,444 KB
testcase_76 AC 66 ms
29,076 KB
testcase_77 AC 63 ms
27,436 KB
testcase_78 AC 2 ms
6,944 KB
testcase_79 AC 18 ms
6,940 KB
testcase_80 AC 37 ms
12,268 KB
testcase_81 AC 9 ms
6,940 KB
testcase_82 AC 13 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);
        }
      }
break;
    }
  } catch (EOFException e) {
  }
}
0