結果

問題 No.23 技の選択
コンテスト
ユーザー どらら
提出日時 2015-06-04 01:35:18
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,332 bytes
コンパイル時間 492 ms
コンパイル使用メモリ 80,648 KB
実行使用メモリ 11,520 KB
最終ジャッジ日時 2024-07-06 14:05:21
合計ジャッジ時間 1,867 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 13 WA * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <iostream>
#include <queue>
#include <list>
#include <map>
#include <numeric>
#include <set>
#include <sstream>
#include <string>
#include <vector>
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;

unsigned long long Comb(int n, int k){
  unsigned long long f[n+1][n+1];
  for(int i=0; i<n+1; i++) f[0][i] = f[i][0] = 1;
  for(int i=1; i<n+1; i++) for(int j=1; j<n+1; j++){
    f[i][j] = f[i-1][j]+f[i][j-1];
    if(i+j==n && j==k) return f[i][j];
  }
  return 1;
}

double calc(int B, int cnt){
  double ret = 0;
  for(int y=cnt;;){
    double sum = y * (double)Comb(y-1, cnt-1) * pow(1/3.0, y-cnt) * pow(2/3.0, cnt);
    if(sum < 10e-10) break;
    ret += sum;
    y++;
  }
  return ret;
}


int main(){
  int H, A, D;
  cin >> H >> A >> D;
  
  double ret = 1000000000.0;
  int x = 0;
  while(H-A*x>=0){
    int B = H-A*x;
    if(B==0){
      ret = min(ret, (double)x);
      break;
    }
    int cnt = B/D + ((B%D==0)?0:1);
    double y = calc(B, cnt);
    ret = min(ret, x+y);
    x++;
  }

  printf("%.6lf\n", ret);
  return 0;
}
0