結果

問題 No.75 回数の期待値の問題
ユーザー te-sh
提出日時 2017-01-25 15:59:13
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 659 bytes
コンパイル時間 2,630 ms
コンパイル使用メモリ 151,852 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 06:34:45
合計ジャッジ時間 3,053 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;

void main()
{
  auto k = readln.chomp.to!int;

  auto eps = real(0.001);

  auto lower(real _, real a) {
    return calcE0(k, a) < a;
  }

  auto r = real(1);
  while (!lower(0, r))
    r *= 2;

  auto m = iota(r / 2, r, eps).assumeSorted!lower.upperBound(0);
  writeln(m.front);
}

auto calcE0(int k, real m)
{
  auto ei = new real[](k);
  foreach_reverse (i; k.iota) {
    ei[i] = 0;
    foreach (j; 1..7) {
      if (i + j > k)
        ei[i] += m;
      else if (i + j == k)
        ei[i] += 0;
      else
        ei[i] += ei[i + j];
    }
    ei[i] = ei[i] / 6 + 1;
  }
  return ei[0];
}
0