結果

問題 No.298 話の伝達
ユーザー te-sh
提出日時 2017-06-09 15:27:44
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 539 ms / 5,000 ms
コード長 1,016 bytes
コンパイル時間 2,285 ms
コンパイル使用メモリ 154,088 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 19:51:06
合計ジャッジ時間 5,858 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #

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

// allowable-error: 10 ** -6

void main()
{
  auto rd1 = readln.split.to!(size_t[]), n = rd1[0], m = rd1[1];

  auto p = new real[][](n, n);
  foreach (_; 0..m) {
    auto rd2 = readln.split, u = rd2[0].to!size_t, v = rd2[1].to!size_t, q = rd2[2].to!real / 100;
    p[v][u] = q;
  }

  auto s = real(0);

  foreach (i; 0..(1 << (n-2))) {
    auto j = (i << 1).bitSet(0).bitSet(n-1);
    auto a = real(1);

    foreach (k; 1..n) {
      auto r = real(1);
      foreach (l; 0..n)
        if (j.bitTest(l) && !p[k][l].isNaN)
          r *= 1 - p[k][l];

      a *= j.bitTest(k) ? 1-r : r;
    }

    s += a;
  }

  writeln(s);
}

pragma(inline) {
  pure bool bitTest(T)(T n, size_t i) { return (n & (T(1) << i)) != 0; }
  pure T bitSet(T)(T n, size_t i) { return n | (T(1) << i); }
  pure T bitReset(T)(T n, size_t i) { return n & ~(T(1) << i); }
  pure T bitComp(T)(T n, size_t i) { return n ^ (T(1) << i); }
}
0