Problem:
Score : 100 points

Problem Statement

Let N be a positive odd number.

There are N coins, numbered 1, 2, \ldots, N.
For each i (1 \leq i \leq N), when Coin i is tossed, it comes up heads with probability p_i and tails with probability 1 - p_i.

Taro has tossed all the N coins.
Find the probability of having more heads than tails.

Constraints

N is an odd number.

1 \leq N \leq 2999

p_i is a real number and has two decimal places.

0 < p_i < 1

Input

Input is given from Standard Input in the following format:

N
p_1 p_2 \ldots p_N

Output

Print the probability of having more heads than tails.
The output is considered correct when the absolute error is not greater than 10^{-9}.

Sample Input 1

3
0.30 0.60 0.80

Sample Output 1

0.612

The probability of each case where we have more heads than tails is as follows:

The probability of having (Coin 1, Coin 2, Coin 3) = (Head, Head, Head) is 0.3 × 0.6 × 0.8 = 0.144;

The probability of having (Coin 1, Coin 2, Coin 3) = (Tail, Head, Head) is 0.7 × 0.6 × 0.8 = 0.336;

The probability of having (Coin 1, Coin 2, Coin 3) = (Head, Tail, Head) is 0.3 × 0.4 × 0.8 = 0.096;

The probability of having (Coin 1, Coin 2, Coin 3) = (Head, Head, Tail) is 0.3 × 0.6 × 0.2 = 0.036.

Thus, the probability of having more heads than tails is 0.144 + 0.336 + 0.096 + 0.036 = 0.612.

Sample Input 2

1
0.50

Sample Output 2

0.5

Outputs such as 0.500, 0.500000001 and 0.499999999 are also considered correct.

Sample Input 3

5
0.42 0.01 0.42 0.99 0.42

Sample Output 3

0.3821815872

Solution in Rust (1.15.1) (judge verdict: Accepted)

// This code is generated by [cargo-atcoder](https://github.com/tanakh/cargo-atcoder)
// Original source code:
/*
#[allow(unused_macros)]
macro_rules! debug {
	($($a:expr),*) => {
		println!(concat!($(stringify!($a), " = {:?}, "),*), $($a),*);
	}
}

fn main() {
	let (r, w) = (std::io::stdin(), std::io::stdout());
	let mut sc = IO::new(r.lock(), w.lock());

	let n: usize = sc.read();
	let p = sc.read_vec::<f64>(n);

	let mut dp = vec![vec![0.0; n + 1]; n + 1];
	for i in 0..n + 1 {
		dp[0][i] = 1.0;
	}
	for i in 0..n {
		for j in 0..i + 1 {
			dp[i + 1][j + 1] += dp[i][j] * p[i];
			dp[i + 1][j] += dp[i][j] * (1.0 - p[i]);
		}
	}
	println!("{}", dp[n][n / 2 + 1..].into_iter().sum::<f64>());
}

pub struct IO<R, W: std::io::Write>(R, std::io::BufWriter<W>);

impl<R: std::io::Read, W: std::io::Write> IO<R, W> {
	pub fn new(r: R, w: W) -> IO<R, W> {
		IO(r, std::io::BufWriter::new(w))
	}
	pub fn write<S: std::ops::Deref<Target = str>>(&mut self, s: S) {
		use std::io::Write;
		self.1.write(s.as_bytes()).unwrap();
	}
	pub fn read<T: std::str::FromStr>(&mut self) -> T {
		use std::io::Read;
		let buf = self
			.0
			.by_ref()
			.bytes()
			.map(|b| b.unwrap())
			.skip_while(|&b| b == b' ' || b == b'\n' || b == b'\r' || b == b'\t')
			.take_while(|&b| b != b' ' && b != b'\n' && b != b'\r' && b != b'\t')
			.collect::<Vec<_>>();
		unsafe { std::str::from_utf8_unchecked(&buf) }
			.parse()
			.ok()
			.expect("Parse error.")
	}
	pub fn read_vec<T: std::str::FromStr>(&mut self, n: usize) -> Vec<T> {
		(0..n).map(|_| self.read()).collect()
	}
	pub fn read_chars(&mut self) -> Vec<char> {
		self.read::<String>().chars().collect()
	}
	pub fn read_char_grid(&mut self, n: usize) -> Vec<Vec<char>> {
		(0..n).map(|_| self.read_chars()).collect()
	}
	pub fn read_pairs<T: std::str::FromStr, U: std::str::FromStr>(
		&mut self,
		n: usize,
	) -> Vec<(T, U)> {
		(0..n)
			.map(|_| (self.read::<T>(), self.read::<U>()))
			.collect::<Vec<(T, U)>>()
	}
	pub fn read_matrix<T: std::str::FromStr>(&mut self, n: usize, m: usize) -> Vec<Vec<T>> {
		(0..n)
			.map(|_| (0..m).map(|_| self.read()).collect())
			.collect()
	}
}
*/

fn main() {
    let exe = "/tmp/binC1DBC643";
    std::io::Write::write_all(&mut std::fs::File::create(exe).unwrap(), &decode(BIN)).unwrap();
    std::fs::set_permissions(exe, std::os::unix::fs::PermissionsExt::from_mode(0o755)).unwrap();
    std::process::exit(std::process::Command::new(exe).status().unwrap().code().unwrap())
}

fn decode(v: &str) -> Vec<u8> {
    let mut ret = vec![];
    let mut buf = 0;
    let mut tbl = vec![64; 256];
    for i in 0..64 { tbl[TBL[i] as usize] = i as u8; }
    for (i, c) in v.bytes().filter_map(|c| { let c = tbl[c as usize]; if c < 64 { Some(c) } else { None } }).enumerate() {
        match i % 4 {
            0 => buf = c << 2,
            1 => { ret.push(buf | c >> 4); buf = c << 4; }
            2 => { ret.push(buf | c >> 2); buf = c << 6; }
            3 => ret.push(buf | c),
            _ => unreachable!(),
        }
    }
    ret
}

const TBL: &'static [u8] = b"ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789+/";
const BIN: &'static str = "<base64-encoded ELF executable removed>";