結果

問題 No.1069 電柱 / Pole (Hard)
ユーザー 👑 hos.lyrichos.lyric
提出日時 2020-05-29 22:34:22
言語 D
(dmd 2.106.1)
結果
AC  
実行時間 1,261 ms / 2,000 ms
コード長 5,374 bytes
コンパイル時間 3,642 ms
コンパイル使用メモリ 274,236 KB
実行使用メモリ 52,356 KB
最終ジャッジ日時 2023-09-04 08:00:34
合計ジャッジ時間 10,122 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 1,261 ms
47,728 KB
testcase_05 AC 948 ms
50,800 KB
testcase_06 AC 49 ms
21,440 KB
testcase_07 AC 55 ms
20,772 KB
testcase_08 AC 24 ms
12,244 KB
testcase_09 AC 52 ms
20,260 KB
testcase_10 AC 56 ms
27,764 KB
testcase_11 AC 5 ms
4,380 KB
testcase_12 AC 4 ms
4,376 KB
testcase_13 AC 17 ms
4,848 KB
testcase_14 AC 8 ms
4,784 KB
testcase_15 AC 10 ms
4,556 KB
testcase_16 AC 23 ms
10,920 KB
testcase_17 AC 5 ms
4,376 KB
testcase_18 AC 18 ms
7,688 KB
testcase_19 AC 39 ms
13,552 KB
testcase_20 AC 5 ms
4,380 KB
testcase_21 AC 9 ms
4,376 KB
testcase_22 AC 5 ms
4,500 KB
testcase_23 AC 22 ms
4,380 KB
testcase_24 AC 79 ms
36,856 KB
testcase_25 AC 85 ms
36,956 KB
testcase_26 AC 78 ms
36,712 KB
testcase_27 AC 81 ms
37,024 KB
testcase_28 AC 79 ms
37,920 KB
testcase_29 AC 16 ms
11,976 KB
testcase_30 AC 3 ms
4,380 KB
testcase_31 AC 2 ms
4,380 KB
testcase_32 AC 2 ms
4,376 KB
testcase_33 AC 1 ms
4,376 KB
testcase_34 AC 2 ms
4,380 KB
testcase_35 AC 1 ms
4,376 KB
testcase_36 AC 2 ms
4,376 KB
testcase_37 AC 5 ms
4,380 KB
testcase_38 AC 6 ms
4,380 KB
testcase_39 AC 2 ms
4,376 KB
testcase_40 AC 4 ms
4,380 KB
testcase_41 AC 7 ms
4,380 KB
testcase_42 AC 6 ms
4,376 KB
testcase_43 AC 6 ms
4,380 KB
testcase_44 AC 31 ms
12,240 KB
testcase_45 AC 13 ms
6,580 KB
testcase_46 AC 22 ms
7,700 KB
testcase_47 AC 34 ms
12,976 KB
testcase_48 AC 63 ms
27,372 KB
testcase_49 AC 34 ms
16,484 KB
testcase_50 AC 4 ms
4,376 KB
testcase_51 AC 47 ms
20,288 KB
testcase_52 AC 6 ms
4,380 KB
testcase_53 AC 56 ms
25,972 KB
testcase_54 AC 2 ms
4,376 KB
testcase_55 AC 2 ms
4,376 KB
testcase_56 AC 1 ms
4,380 KB
testcase_57 AC 2 ms
4,380 KB
testcase_58 AC 2 ms
4,376 KB
testcase_59 AC 2 ms
4,380 KB
testcase_60 AC 2 ms
4,376 KB
testcase_61 AC 2 ms
4,380 KB
testcase_62 AC 1 ms
4,376 KB
testcase_63 AC 2 ms
4,380 KB
testcase_64 AC 2 ms
4,376 KB
testcase_65 AC 1 ms
4,376 KB
testcase_66 AC 2 ms
4,376 KB
testcase_67 AC 2 ms
4,380 KB
testcase_68 AC 1 ms
4,504 KB
testcase_69 AC 79 ms
40,708 KB
testcase_70 AC 95 ms
50,328 KB
testcase_71 AC 93 ms
50,300 KB
testcase_72 AC 24 ms
14,036 KB
testcase_73 AC 40 ms
20,260 KB
testcase_74 AC 96 ms
52,076 KB
testcase_75 AC 96 ms
52,356 KB
testcase_76 AC 69 ms
30,508 KB
testcase_77 AC 64 ms
27,952 KB
testcase_78 AC 3 ms
4,376 KB
testcase_79 AC 19 ms
6,604 KB
testcase_80 AC 38 ms
13,036 KB
testcase_81 AC 9 ms
4,584 KB
testcase_82 AC 14 ms
4,556 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[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 != 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 != 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