Problem: Your friend recently gave you some slimes for your birthday. You have a very large amount of slimes with value 1 and 2, and you decide to invent a game using these slimes.

You initialize a row with n empty spaces. You also choose a number p to be used in the game. Then, you will perform the following steps while the last space is empty. 

  1. With probability <image>, you will choose a slime with value 1, and with probability <image>, you will choose a slime with value 2. You place the chosen slime on the last space of the board. 
  2. You will push the slime to the left as far as possible. If it encounters another slime, and they have the same value v, you will merge the slimes together to create a single slime with value v + 1. This continues on until the slime reaches the end of the board, or encounters a slime with a different value than itself. 



You have played the game a few times, but have gotten bored of it. You are now wondering, what is the expected sum of all values of the slimes on the board after you finish the game.

Input

The first line of the input will contain two integers n, p (1 ≤ n ≤ 109, 1 ≤ p < 109).

Output

Print the expected sum of all slimes on the board after the game finishes. Your answer will be considered correct if its absolute or relative error does not exceed 10 - 4.

Namely, let's assume that your answer is a and the answer of the jury is b. The checker program will consider your answer correct, if <image>.

Examples

Input

2 500000000


Output

3.562500000000000


Input

10 1


Output

64.999983360007620


Input

100 123456789


Output

269.825611298854770

Note

In the first sample, we have a board with two squares, and there is a 0.5 probability of a 1 appearing and a 0.5 probability of a 2 appearing.

Our final board states can be 1 2 with probability 0.25, 2 1 with probability 0.375, 3 2 with probability 0.1875, 3 1 with probability 0.1875. The expected value is thus (1 + 2)·0.25 + (2 + 1)·0.375 + (3 + 2)·0.1875 + (3 + 1)·0.1875 = 3.5625.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
namespace pb_ds {
namespace io {
const int MaxBuff = 1 << 15;
const int Output = 1 << 23;
char B[MaxBuff], *S = B, *T = B;
char Out[Output], *iter = Out;
inline void flush() {
  fwrite(Out, 1, iter - Out, stdout);
  iter = Out;
}
}  // namespace io
template <class Type>
inline Type Read() {
  using namespace io;
  register char ch;
  register Type ans = 0;
  register bool neg = 0;
  while (ch = ((S == T) && (T = (S = B) + fread(B, 1, MaxBuff, stdin), S == T)
                   ? 0
                   : *S++),
         (ch < '0' || ch > '9') && ch != '-')
    ;
  ch == '-' ? neg = 1 : ans = ch - '0';
  while (ch = ((S == T) && (T = (S = B) + fread(B, 1, MaxBuff, stdin), S == T)
                   ? 0
                   : *S++),
         '0' <= ch && ch <= '9')
    ans = ans * 10 + ch - '0';
  return neg ? -ans : ans;
}
template <class Type>
inline void Print(register Type x, register char ch = '\n') {
  using namespace io;
  if (!x)
    *iter++ = '0';
  else {
    if (x < 0) *iter++ = '-', x = -x;
    static int s[100];
    register int t = 0;
    while (x) s[++t] = x % 10, x /= 10;
    while (t) *iter++ = '0' + s[t--];
  }
  *iter++ = ch;
}
}  // namespace pb_ds
using namespace pb_ds;
using namespace std;
const int N = 55;
int n;
double p, ans;
double p1[N][N], p2[N][N], s[N];
double f[N][N], g[N][N], dp[N][N];
struct mat {
  double a[N][N];
  inline void init(int opt = 0) {
    memset(a, 0, sizeof(a));
    if (!opt) return;
    for (int i = 0; i <= 50; ++i) a[i][i] = 1;
  }
  inline mat operator*(const mat &u) const {
    mat res;
    res.init();
    for (int i = 0; i <= 50; ++i)
      for (int k = 0; k <= 50; ++k)
        for (int j = 0; j <= 50; ++j) res.a[i][j] += a[i][k] * u.a[k][j];
    return res;
  }
} st;
inline mat Pow(mat x, int y) {
  mat res;
  res.init(1);
  while (y) {
    if (y & 1) res = res * x;
    x = x * x;
    y >>= 1;
  }
  return res;
}
inline void init() {
  p1[1][1] = p, p1[1][2] = 1 - p;
  f[1][1] = p, f[1][2] = 1 - p, g[1][2] = 1;
  for (int i = 2; i <= 50; ++i) {
    p1[i][1] = p, p1[i][2] = 1 - p + p * p, p2[i][2] = 1;
    for (int j = 3; j <= 50; ++j)
      p1[i][j] = p1[i][j - 1] * p1[i - 1][j - 1],
      p2[i][j] = p2[i][j - 1] * p1[i - 1][j - 1];
    for (int j = 1; j <= 50; ++j)
      f[i][j] = p1[i][j] * (1 - p1[i - 1][j]),
      g[i][j] = p2[i][j] * (1 - p1[i - 1][j]);
  }
  dp[1][1] = 1, dp[1][2] = 2;
  for (int i = 2; i <= 50; ++i) {
    double sum = 0;
    dp[i][1] = 1;
    for (int j = 1; j <= 50; ++j) sum += g[i - 1][j];
    for (int j = 1; j <= 50; ++j) dp[i][1] += dp[i - 1][j] * g[i - 1][j] / sum;
    sum = 0;
    for (int j = 2; j <= 50; ++j) {
      sum += f[i - 1][j - 1];
      for (int k = 1; k < j; ++k) dp[i][j] += dp[i - 1][k] * f[i - 1][k] / sum;
      dp[i][j] += j;
    }
  }
  double sum = 0;
  st.init();
  for (int i = 1; i <= 50; ++i) sum += g[50][i];
  for (int i = 1; i <= 50; ++i) st.a[i][1] = g[50][i] / sum;
  sum = 0;
  for (int i = 2; i <= 50; ++i) {
    sum += f[50][i - 1];
    for (int j = 1; j < i; ++j) st.a[j][i] = f[50][j] / sum;
  }
  for (int i = 1; i <= 50; ++i) st.a[0][i] = i;
  st.a[0][0] = 1;
}
int main() {
  n = Read<int>(), p = Read<int>();
  p = p / (double)1000000000;
  init();
  if (n <= 50) {
    for (int i = 1; i <= 50; ++i) ans += dp[n][i] * f[n][i];
    printf("%.10lf\n", ans);
    return 0;
  }
  st = Pow(st, n - 50);
  for (int i = 1; i <= 50; ++i)
    for (int j = 1; j <= 50; ++j) s[j] += dp[50][i] * st.a[i][j];
  for (int i = 1; i <= 50; ++i) {
    s[i] += st.a[0][i];
    ans += s[i] * f[50][i];
  }
  printf("%.10lf\n", ans);
  return 0;
}