結果

問題 No.301 サイコロで確率問題 (1)
コンテスト
ユーザー rickytheta
提出日時 2015-11-13 23:42:25
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 957 ms / 1,000 ms
コード長 1,124 bytes
コンパイル時間 1,601 ms
コンパイル使用メモリ 158,716 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-13 15:39:46
合計ジャッジ時間 3,867 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef complex<double> P;
typedef pair<int,int> pii;
#define REP(i,n) for(ll i=0;i<n;++i)
#define REPR(i,n) for(ll i=1;i<n;++i)
#define FOR(i,a,b) for(ll i=a;i<b;++i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define MOD (ll)(1e9+7)
#define ADD(a,b) a=((a)+(b))%MOD
#define FIX(a) ((a)%MOD+MOD)%MOD

ll n;
double dp[1010];

bool check(double d){
  fill(dp,dp+1010,0);
  REP(i,n){
    REP(j,6){
      if(j>i) dp[i+1]+=d;
      if(i>=j) dp[i+1]+=dp[i-j];
    }
    dp[i+1]=dp[i+1]/6.0+1;
  }
  return dp[n]>d;
}

int main(){
  ll t;
  cin >> t;
  cout.precision(20);
  cout.setf(ios::fixed);
  while(t--){
    cin >> n;
    if(n<100){
      double l=0,r=1e9;
      REP(i,10000){
        double m=(l+r)/2.0;
        if(check(m))l=m;
        else r=m;
      }
      cout << (l+r)/2 << endl;
    }else{
      double res = (double)(n+1)+2.0/3.0;
      cout << res << endl;
    }
  }
  return 0;
}
0