結果

問題 No.94 圏外です。(EASY)
ユーザー te-sh
提出日時 2017-05-09 11:07:08
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 15 ms / 5,000 ms
コード長 1,545 bytes
コンパイル時間 2,778 ms
コンパイル使用メモリ 164,432 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 19:04:50
合計ジャッジ時間 3,029 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.math;      // math functions

void main()
{
  auto n = readln.chomp.to!size_t;
  auto pi = n.iota.map!(_ => readln.split.to!(int[])).map!(rd => point(rd[0], rd[1])).array;

  if (n == 0) {
    writeln(1);
    return;
  }

  auto dij = new int[][](n, n);
  foreach (i; 0..n)
    foreach (j; i+1..n)
      dij[i][j] = dij[j][i] = (pi[i] - pi[j]).hypot2;

  auto uf = UnionFind!size_t(n);
  foreach (i; 0..n)
    foreach (j; i+1..n)
      if (i != j && dij[i][j] <= 100)
        uf.unite(i, j);

  auto maxD = uf.groups
      .map!(g => dij.indexed(g).map!(di => di.indexed(g)).joiner.fold!max)
      .fold!max;
  writefln("%.7f", maxD.to!real.sqrt + 2);
}

alias Point!int point;

struct UnionFind(T)
{
  import std.algorithm, std.range;

  T[] p; // parent
  const T s; // sentinel
  const T n;

  this(T n)
  {
    this.n = n;
    p = new T[](n);
    s = n + 1;
    p[] = s;
  }

  T find(T i)
  {
    if (p[i] == s) {
      return i;
    } else {
      p[i] = find(p[i]);
      return p[i];
    }
  }

  void unite(T i, T j)
  {
    auto pi = find(i), pj = find(j);
    if (pi != pj) p[pj] = pi;
  }

  bool isSame(T i, T j) { return find(i) == find(j); }

  auto groups()
  {
    auto g = new T[][](n);
    foreach (i; 0..n) g[find(i)] ~= i;
    return g.filter!(l => !l.empty);
  }
}

struct Point(T) {
  T x, y;

  point opBinary(string op)(point rhs) {
    static if (op == "-") return point(x - rhs.x, y - rhs.y);
  }

  T hypot2() { return x ^^ 2 + y ^^ 2; }
}
0