#include using namespace std; struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_; #define FOR(i, begin, end) for(int i=(begin);i<(end);i++) #define REP(i, n) FOR(i,0,n) #define IFOR(i, begin, end) for(int i=(end)-1;i>=(begin);i--) #define IREP(i, n) IFOR(i,0,n) #define Sort(v) sort(v.begin(), v.end()) #define Reverse(v) reverse(v.begin(), v.end()) #define all(v) v.begin(),v.end() #define SZ(v) ((int)v.size()) #define Lower_bound(v, x) distance(v.begin(), lower_bound(v.begin(), v.end(), x)) #define Upper_bound(v, x) distance(v.begin(), upper_bound(v.begin(), v.end(), x)) #define Max(a, b) a = max(a, b) #define Min(a, b) a = min(a, b) #define bit(n) (1LL<<(n)) #define bit_exist(x, n) ((x >> n) & 1) #define Ans(f, y, n) if(f) cout << y << endl; else cout << n << endl; #define debug(x) cout << #x << "=" << x << endl; #define vdebug(v) cout << #v << "=" << endl; REP(i, v.size()){ cout << v[i] << ","; } cout << endl; #define mdebug(m) cout << #m << "=" << endl; REP(i, m.size()){ REP(j, m[i].size()){ cout << m[i][j] << ","; } cout << endl;} #define pb push_back #define f first #define s second #define int long long #define INF 1000000000000000000 using vec = vector; using mat = vector; using Pii = pair; using PiP = pair; using PPi = pair; using bools = vector; using pairs = vector; template void readv(vector &a){ REP(i, a.size()) cin >> a[i]; } void readv_m1(vector &a){ REP(i, a.size()){cin >> a[i]; a[i]--;} } //int dx[4] = {1,0,-1,0}; //int dy[4] = {0,1,0,-1}; int mod = 1000000007; //int mod = 998244353; #define Add(x, y) x = (x + (y)) % mod #define Mult(x, y) x = (x * (y)) % mod signed main(){ int N; cin >> N; double p, q; cin >> p >> q; int M = 10000; mat dp(M + 1, vec(2 * (N + 2), 0.0)); dp[0][1] = 1.0; REP(i, M){ dp[i + 1][0] = dp[i][0]; dp[i + 1][N + 2] = dp[i][N + 2]; dp[i + 1][N + 1] = dp[i][N + 1]; dp[i + 1][N + 1 + N + 2] = dp[i][N + 1 + N + 2]; FOR(j, 1, N + 1){ dp[i + 1][j - 1 + N + 2] += p * dp[i][j]; dp[i + 1][j + 1] += q * dp[i][j]; dp[i + 1][j - 1 + N + 2] += q * dp[i][j + N + 2]; dp[i + 1][j + 1] += p * dp[i][j + N + 2]; } } //mdebug(dp); cout << dp[M][N + 2] << endl; return 0; }