import numpy as np

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

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

print(solve(N, p))