Problem: You're playing a game called Osu! Here's a simplified version of it. There are n clicks in a game. For each click there are two outcomes: correct or bad. Let us denote correct as "O", bad as "X", then the whole play can be encoded as a sequence of n characters "O" and "X".

Using the play sequence you can calculate the score for the play as follows: for every maximal consecutive "O"s block, add the square of its length (the number of characters "O") to the score. For example, if your play can be encoded as "OOXOOOXXOO", then there's three maximal consecutive "O"s block "OO", "OOO", "OO", so your score will be 22 + 32 + 22 = 17. If there are no correct clicks in a play then the score for the play equals to 0.

You know that the probability to click the i-th (1 ≤ i ≤ n) click correctly is pi. In other words, the i-th character in the play sequence has pi probability to be "O", 1 - pi to be "X". You task is to calculate the expected score for your play.

Input

The first line contains an integer n (1 ≤ n ≤ 105) — the number of clicks. The second line contains n space-separated real numbers p1, p2, ..., pn (0 ≤ pi ≤ 1).

There will be at most six digits after the decimal point in the given pi.

Output

Print a single real number — the expected score for your play. Your answer will be considered correct if its absolute or relative error does not exceed 10 - 6.

Examples

Input

3
0.5 0.5 0.5


Output

2.750000000000000


Input

4
0.7 0.2 0.1 0.9


Output

2.489200000000000


Input

5
1 1 1 1 1


Output

25.000000000000000

Note

For the first example. There are 8 possible outcomes. Each has a probability of 0.125.

  * "OOO"  →  32 = 9; 
  * "OOX"  →  22 = 4; 
  * "OXO"  →  12 + 12 = 2; 
  * "OXX"  →  12 = 1; 
  * "XOO"  →  22 = 4; 
  * "XOX"  →  12 = 1; 
  * "XXO"  →  12 = 1; 
  * "XXX"  →  0. 



So the expected score is <image>

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
int diru[] = {-1, -1, -1, 0, 0, 1, 1, 1};
int dirv[] = {-1, 0, 1, -1, 1, -1, 0, 1};
using namespace std;
struct debugger {
  template <typename T>
  debugger& operator,(const T& v) {
    cerr << v << " ";
    return *this;
  }
} dbg;
template <class T>
T sq(T n) {
  return n * n;
}
template <class T>
T gcd(T a, T b) {
  return (b != 0 ? gcd<T>(b, a % b) : a);
}
template <class T>
T lcm(T a, T b) {
  return (a / gcd<T>(a, b) * b);
}
template <class T>
bool inside(T a, T b, T c) {
  return a <= b && b <= c;
}
template <class T>
void setmax(T& a, T b) {
  if (a < b) a = b;
}
template <class T>
void setmin(T& a, T b) {
  if (b < a) a = b;
}
template <class T>
T power(T N, T P) {
  return (P == 0) ? 1 : N * power(N, P - 1);
}
double dp0[100010];
double dp1[100010];
double dp2[100010];
double in[100010];
int main() {
  int n, T, t = 1, m, i, j, k;
  scanf("%d", &n);
  for (i = 0; i < n; i++) {
    cin >> in[i];
  }
  dp0[n] = 1;
  for (i = n - 1; i > -1; i--) {
    dp0[i] = in[i] * (dp0[i + 1]) + ((double)1 - in[i]);
  }
  for (i = n - 1; i > -1; i--) {
    dp1[i] = in[i] * (dp1[i + 1] + dp0[i + 1]) + ((double)1 - in[i]) * 0;
  }
  for (i = n - 1; i > -1; i--) {
    dp2[i] = in[i] * (dp2[i + 1] + (double)2 * dp1[i + 1] + dp0[i + 1]) +
             ((double)1 - in[i]) * (dp2[i + 1]);
  }
  printf("%.15lf\n", dp2[0]);
  return 0;
}