結果

問題 No.813 ユキちゃんの冒険
コンテスト
ユーザー eQe
提出日時 2019-04-13 17:04:01
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 909 bytes
コンパイル時間 1,338 ms
コンパイル使用メモリ 159,332 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-21 08:29:57
合計ジャッジ時間 2,250 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define ft first
#define sc second
#define lb lower_bound
#define ub upper_bound
#define pb push_back
#define pt(sth) cout << sth << "\n"
#define chmax(a, b) {if(a<b) a=b;}
#define chmin(a, b) {if(a>b) a=b;}
#define moC(a, s, b) (a)=((a)s(b)+MOD)%MOD
using namespace std;
typedef long long ll;
typedef pair<ll, ll> P;
static const ll INF=1e18;
static const ll MAX=1e5+7;
static const ll MOD=1e9+7;

typedef pair<double, double> DP;

double p, q;

DP joint(DP dp1, DP dp2) {
  DP res;
  double p1=dp1.ft, q1=dp1.sc;
  double p2=dp2.ft, q2=dp2.sc;
  res.ft=p1+p2*q1*q1/(1-p1*p2);
  res.sc=q1*q2/(1-p1*p2);
  return res;
}

DP dpow(DP dp, ll n) {
  DP res=DP(p, q);
  while(n>0) {
    if(n&1) res=joint(res, dp);
    dp=joint(dp, dp);
    n>>=1;
  }
  return res;
}


int main(void) {
  ll N;
  cin >> N;
  cin >> p >> q;
  ll i, j;
  
  printf("%.10f\n", dpow(DP(p, q), N-1).ft);
}
0