Cannot sample enough valid points. (more)

\[0 \lt c_p \land 0 \lt c_n\]
\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]
\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}
double f(double c_p, double c_n, double t, double s) {
        double r5294958 = 1.0;
        double r5294959 = s;
        double r5294960 = -r5294959;
        double r5294961 = exp(r5294960);
        double r5294962 = r5294958 + r5294961;
        double r5294963 = r5294958 / r5294962;
        double r5294964 = c_p;
        double r5294965 = pow(r5294963, r5294964);
        double r5294966 = r5294958 - r5294963;
        double r5294967 = c_n;
        double r5294968 = pow(r5294966, r5294967);
        double r5294969 = r5294965 * r5294968;
        double r5294970 = t;
        double r5294971 = -r5294970;
        double r5294972 = exp(r5294971);
        double r5294973 = r5294958 + r5294972;
        double r5294974 = r5294958 / r5294973;
        double r5294975 = pow(r5294974, r5294964);
        double r5294976 = r5294958 - r5294974;
        double r5294977 = pow(r5294976, r5294967);
        double r5294978 = r5294975 * r5294977;
        double r5294979 = r5294969 / r5294978;
        return r5294979;
}