結果
問題 | No.298 話の伝達 |
ユーザー |
|
提出日時 | 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 |
ソースコード
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); } }