結果

問題 No.75 回数の期待値の問題
ユーザー 👑 ナナチナナチ
提出日時 2020-07-28 21:12:42
言語 C++17(1z)
(gcc 10.2.0 + boost 1.73.0)
結果
AC  
実行時間 4 ms / 5,000 ms
コード長 1,360 Byte
コンパイル時間 1,974 ms
使用メモリ 3,564 KB
最終ジャッジ日時 2020-09-26 10:58:10
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 2 ms
3,484 KB
testcase_01 AC 1 ms
3,344 KB
testcase_02 AC 2 ms
3,412 KB
testcase_03 AC 2 ms
3,412 KB
testcase_04 AC 2 ms
3,460 KB
testcase_05 AC 2 ms
3,520 KB
testcase_06 AC 2 ms
3,408 KB
testcase_07 AC 2 ms
3,348 KB
testcase_08 AC 2 ms
3,412 KB
testcase_09 AC 2 ms
3,484 KB
testcase_10 AC 2 ms
3,520 KB
testcase_11 AC 2 ms
3,344 KB
testcase_12 AC 2 ms
3,460 KB
testcase_13 AC 2 ms
3,348 KB
testcase_14 AC 2 ms
3,480 KB
testcase_15 AC 4 ms
3,564 KB
testcase_16 AC 2 ms
3,392 KB
testcase_17 AC 2 ms
3,396 KB
testcase_18 AC 2 ms
3,412 KB
testcase_19 AC 2 ms
3,512 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef pair<ll, ll> p_ll;

template<class T>
void debug(T itr1, T itr2) { auto now = itr1; while(now<itr2) { cout << *now << " "; now++; } cout << endl; }
#define repr(i,from,to) for (int i=(int)from; i<(int)to; i++)
#define all(vec) vec.begin(), vec.end()
#define rep(i,N) repr(i,0,N)
#define per(i,N) for (int i=(int)N-1; i>=0; i--)

const ll MOD = pow(10,9)+7;
const ll LLINF = pow(2,61)-1;
const int INF = pow(2,30)-1;

vector<ll> fac;
void c_fac(int x=pow(10,6)+10) { fac.resize(x,true); rep(i,x) fac[i] = i ? (fac[i-1]*i)%MOD : 1; }
ll inv(ll a, ll m=MOD) { ll b = m, x = 1, y = 0; while (b!=0) { int d = a/b; a -= b*d; swap(a,b); x -= y*d; swap(x,y); } return (x+m)%m; }
ll nck(ll n, ll k) { return fac[n]*inv(fac[k]*fac[n-k]%MOD)%MOD; }
ll gcd(ll a, ll b) { if (a<b) swap(a,b); return b==0 ? a : gcd(b, a%b); }
ll lcm(ll a, ll b) { return a/gcd(a,b)*b; }


int main() {
  ll K; cin >> K;

  // E[i] = x[i]*E[0] + y[i]
  // E[K] =    0*E[K] + 0
  // E[0] = y[0] / (1-x[0])
  double x[K+1] = {}, y[K+1] = {};
  per(i,K) {
    repr(j,1,7) {
      if (i+j<=K) { x[i] += x[i+j]/6; y[i] += y[i+j]/6; }
      else x[i] += 1./6;
    }
    y[i] += 1;
  }
  // debug(x,x+(K+1)); debug(y,y+(K+1));

  double result = y[0] / (1-x[0]);
  cout << setprecision(10) << result << endl;
  return 0;
}
0