結果

問題 No.76 回数の期待値で練習
コンテスト
ユーザー atetubou
提出日時 2014-11-24 00:31:25
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 1,489 ms / 5,000 ms
コード長 2,659 bytes
コンパイル時間 1,300 ms
コンパイル使用メモリ 156,064 KB
実行使用メモリ 59,520 KB
最終ジャッジ日時 2025-01-02 21:00:56
合計ジャッジ時間 3,242 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <iostream>
#include <vector>
#include <queue>
#include <string>
#include <cstring>
#include <sstream>
#include <set>
#include <map>
#include <algorithm>
#include <cstdio>
#include <cstdlib>

#define rep(i,n) for(int i=0;i<(int)(n);++i)
#define F first
#define S second
#define SZ(a) (int)((a).size())
#define sz(a) SZ(a)
#define pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define ALL(a) (a).begin(),(a).end()
using namespace std;

typedef long long ll;
typedef pair<ll,ll> PI;
typedef unsigned long long ull;

#define PR(...) do{cerr << "line : " << __LINE__ << endl; pr(#__VA_ARGS__, __VA_ARGS__);}while(0);
template<class T>
void pr(const string& name, T t){
  cerr << name << ": " << t << endl;
}

template<typename T, typename ... Types>
void pr(const string& names, T t, Types ... rest) {
  auto comma_pos = names.find(',');
  cerr << names.substr(0, comma_pos) << ": " << t << ", ";
  
  auto next_name_pos = names.find_first_not_of(" \t\n", comma_pos + 1);
  pr(string(names, next_name_pos), rest ...);
}

template<class T,class U> ostream& operator<< (ostream& o, const pair<T,U>& v){return o << "(" << v.F << ", " << v.S << ")";}
template<class T> ostream& operator<< (ostream& o, const vector<T>& v){o << "{";rep(i,SZ(v)) o << (i?", ":"") << v[i];return o << "}";}


#define endl '\n'
int k;
double memo[1100000];
bool vis[1100000];

double pro[7];

double rec(int cur){
  if(cur >= k) return 0;
  if(vis[cur]) return memo[cur];
  vis[cur] = 1;
  double& ret = memo[cur] = 1;
  for(int i = 1; i <= 6; ++i)
    ret += rec(cur+i) * pro[i];
  return ret;
}

/*
  1*(1-x) + x*2 = 1.0833333333333333
  1*(1-x1-x2) + x1*x1*3 + x1 * (1-x1) * 2 + x2 * (1-x1) * 2 = 1.2569444444444444
*/

int main(int argc, char *argv[])
{
  //rep(i,6) pro[i+1] = 1./6;
  
  double prd[] = {1.0000000000000000,1.0833333333333333,1.2569444444444444,1.5353009259259260,1.6915991512345676,2.0513639724794235};
  rep(i,5){
    double sum = 0;
    rep(j,i) sum += pro[j+1];
    double low = 0, up = 1-sum;
    //cout << sum << endl;
    k=i+2;
    rep(k,100){
      double mid = (low+up)/2;
      pro[i+1] = mid;
      pro[i+2] = 1-sum-mid;
      memset(vis,0,sizeof(vis));
      //cout << rec(0) << endl;
      if(rec(0) < prd[i+1]) low = mid;
      else up = mid;
    }
    //cout << low << endl;
    pro[i+1] = low;
    //cout << rec(0) << endl;
  }
  //return 0;  
  
  int t;
  cin >> t;
  rep(i,t){
    //int k;
    cin >> k;
    while(k>1e6);
    memset(vis,0,sizeof(vis));
    printf("%.18f\n",rec(0));
  }
  
  return 0;
}
0