結果

問題 No.813 ユキちゃんの冒険
ユーザー ojisan_ITojisan_IT
提出日時 2019-04-12 22:09:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,446 bytes
コンパイル時間 1,383 ms
コンパイル使用メモリ 163,944 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2025-01-02 03:21:06
合計ジャッジ時間 7,195 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3 WA * 23
権限があれば一括ダウンロードができます

ソースコード

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;
  
  const ll X = 500000;
  long double d = p;
  rep(i,99) d *= p;
  rep(i,X){
    int now = 1;
    while(1){
      long double t = dist(engine);
      //cerr << t << endl;
      if(t < p + 0.001)  now--;
      else if(t < p + q + 0.002)  now++;
      else  break;
      if(now == 0){ans += 1.0; break;}
      if(now > 100){ans += d; break;}
    }
  }
  cout <<fixed << setprecision(10) <<  ans / (1.0*X) << endl;
}
0