N, p = map(float, raw_input().split())
N = int(N)

prime = [1.0] * (N + 1)

for v in xrange(2, N + 1):
    for i in xrange(v + v, N + 1, v):
        prime[i] *= 1 - p

total = 0
for prob in xrange(2, N + 1):
    total += prime[prob]
print total