結果

問題 No.75 回数の期待値の問題
ユーザー どらら
提出日時 2021-01-24 02:55:02
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 995 bytes
コンパイル時間 1,787 ms
コンパイル使用メモリ 163,200 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2025-01-02 05:24:39
合計ジャッジ時間 2,607 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)
#define ALLOF(c) (c).begin(), (c).end()
typedef long long ll;
typedef unsigned long long ull;

struct ST {
  double a, b;
  ST():a(0),b(0){}
  ST(double a, double b):a(a),b(b){}
  ST operator+=(const ST& x){
    a += x.a;
    b += x.b;
    return *this;
  }
  ST operator-=(const ST& x){
    a -= x.a;
    b -= x.b;
    return *this;
  }
  ST operator/=(double x){
    a /= x;
    b /= x;
    return *this;
  }
};

int main(){
  int K;
  cin >> K;
  vector<ST> dp(K+10);
  REP(i,K+1,K+10) dp[i] = ST(1,0);

  for(int i=K-1; i>=0; i--){
    ST sum;
    rep(j,6){
      sum += dp[i+1+j];
    }
    sum /= 6;
    sum += ST(0,1);
    dp[i] = sum;
  }

  double a = dp[0].a;
  double b = dp[0].b;
  double ret = b / (1-a);

  cout << fixed << setprecision(16) << ret << endl;
  
  return 0;
}

0