結果

問題 No.134 走れ！サブロー君
コンテスト
ユーザー te-sh
提出日時 2017-04-13 16:27:32
言語 D
(dmd 2.111.0)
結果
AC  
実行時間 17 ms / 5,000 ms
コード長 1,607 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,434 ms
コンパイル使用メモリ 157,896 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 18:40:56
合計ジャッジ時間 3,120 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

void main()
{
  auto rd1 = readln.split.to!(int[]), s = point(rd1[0], rd1[1]);
  auto n = readln.chomp.to!size_t;
  auto pi = new point[](n), wi = new real[](n);
  foreach (i; 0..n) {
    auto rd2 = readln.split;
    pi[i] = point(rd2[0].to!int, rd2[1].to!int);
    wi[i] = rd2[2].to!real;
  }

  auto vi = new real[](1 << n);
  foreach (i; 0..(1 << n))
    vi[i] = wi.indexed(i.bitsSet).sum;

  auto dp = new real[][](1 << n, n);
  foreach (j; 0..n) {
    dp[0][j] = real(100) / 120 * manDis(pi[j], s) + wi[j];
  }

  foreach (i; 0..((1 << n) - 1)) {
    foreach (j; 0..n) {
      foreach (k; 0..n) {
        if (j == k || i.bitTest(j)) continue;
        auto ni = i.bitSet(j);
        auto nv = dp[i][j] + (vi[ni] + 100) / 120 * manDis(pi[k], pi[j]) + wi[k];
        if (dp[ni][k].isNaN)
          dp[ni][k] = nv;
        else
          dp[ni][k] = min(dp[ni][k], nv);
      }
    }
  }

  auto f = (1 << n) - 1;
  auto mi = real.max;
  foreach (i; n.iota) {
    auto ni = f.bitReset(i);
    auto nv = dp[ni][i] + (vi[f] + 100) / 120 * manDis(pi[i], s);
    mi = min(mi, nv);
  }

  writefln("%.8f", mi);
}

auto manDis(point p1, point p2)
{
  return (p1.x - p2.x).abs + (p1.y - p2.y).abs;
}

template BitOp(T) {
  bool bitTest(T n, size_t i) { return (n & (T(1) << i)) != 0; }
  T bitSet(T n, size_t i) { return n | (T(1) << i); }
  T bitReset(T n, size_t i) { return n & ~(T(1) << i); }
}

mixin BitOp!int;

struct Point(T) {
  T x, y;
}

alias Point!int point;
0