LuckCalculator.kt

// ps: Probability array of "success" for each roll
fun luckSymExact(ps: List<Double>, successes: Int): Double {
    // Check for all fail/succeed case, since this will reported wrong otherwise
    val variance = ps.sumOf { it * (1 - it) }
    if (variance == 0.0) {
        return when (successes > 0) {
            true -> 100.0
            false -> -100.0
        }
    }
    val pmf = pbPmf(ps)
    val mid = midCdf(pmf, successes)        // in [0,1]
    return (2.0 * mid - 1.0) * 100.0        // in [-100,+100]
}