#include<bits/stdc++.h>

using namespace std;

float logC(int n, int k) {
    return lgamma(n + 1) - lgamma(k + 1) - lgamma(n - k + 1);
}

float negative_binomial_pmf(int k, int r, float p) {
    if (k < r) return 0.0;
    float log_pmf = logC(k - 1, r - 1) + r * log(p) + (k - r) * log(1 - p);
    return exp(log_pmf);
}

float negative_binomial_cdf(int k, int r, float p) {
    float cdf = 0.0;
    for (int i = r; i <= k; i++) {
        cdf += negative_binomial_pmf(i, r, p);
    }
    return cdf;
}

int main() {
    int n, l, r;
    float p;
    cin >> n >> p >> l >> r;
    float expected = 0.0;
    for (int i = 1; i <= n; i++) {
        float prob_l = negative_binomial_cdf(l - 1, i, p);
        float prob_r = negative_binomial_cdf(r, i, p);
        expected += prob_r - prob_l;
    }
    cout <<fixed<<setprecision(10)<< expected << endl;
    return 0;
}