結果

問題 No.813 ユキちゃんの冒険
ユーザー ojisan_ITojisan_IT
提出日時 2019-04-12 22:17:13
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,527 bytes
コンパイル時間 1,441 ms
コンパイル使用メモリ 170,028 KB
実行使用メモリ 10,400 KB
最終ジャッジ日時 2024-06-10 13:30:05
合計ジャッジ時間 5,205 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,272 KB
testcase_01 WA -
testcase_02 AC 22 ms
6,944 KB
testcase_03 TLE -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
template<class T> using vt = vector<T>;
template<class T> using vvt = vector<vt<T>>;
template<class T> using ttt = tuple<T,T>;
using tii = tuple<int,int>;
using tiii = tuple<int,int,int>;
using vi = vector<int>;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define pb push_back
#define ALL(a) (a).begin(),(a).end()
#define FST first
#define SEC second
#define DEB cerr<<"!"<<endl
#define SHOW(a,b) cerr<<(a)<<" "<<(b)<<endl
#define DIV int(1e9+7)
const int INF = (INT_MAX/2);
const ll LLINF = (LLONG_MAX/2);
const double eps = 1e-8;
//const double PI = M_PI;  
inline ll pow(ll x,ll n,ll m){ll r=1;while(n>0){if((n&1)==1)r=r*x%m;x=x*x%m;n>>=1;}return r%m;}
inline ll lcm(ll d1, ll d2){return d1 / __gcd(d1, d2) * d2;}

/* Coding space */
std::random_device seed_gen;
std::default_random_engine engine(seed_gen());
std::uniform_real_distribution<long double> dist(0.0, 1.0);
int main(){
  int n; cin >> n;
  long double p,q; cin >> p >> q;
  long double ans = 0.0;
  long double d[1001] = {};
  d[1] = p;
  const ll X = 1000000;
  const ll iX = 1000;
  rep(i,1001) if(i!=0)d[i] *= d[i-1]*p;
  rep(i,X){
    int now = 1;
    rep(j,1000){
      long double t = dist(engine);
      //cerr << t << endl;
      if(t < p)  now--;
      else if(t < p + q)  now++;
      else  break;
      if(now == 0){ans += 1.0; break;}
      if(now >= n){break;}
      if( j == iX -1) {ans += d[now]; break;}
    }
  }
  cout <<fixed << setprecision(10) <<  ans / (1.0*X) << endl;
}
0