import numpy as np

def solve(N, p):
    ps = np.ones(N+1)
    for i in range(2, (N+1)//2):
        ind = np.arange((N+i)//i) * i
        if ps[ind[1]] > 1e-5:
            ps[ind[2:]] *= (1 - p)
    return ps[2:].sum()

N, p = input().split()
N = int(N)
p = float(p)

print(solve(N, p))