#488904 (D) No.1069 電柱 / Pole (Hard)

提出ソース

結果

問題	No.1069 電柱 / Pole (Hard)
ユーザー	👑 hos.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

権限があれば一括ダウンロードができます

ソースコード

raw source code

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) {
  }
}

yukicoder

結果

テストケース

ソースコード