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

prime = [1.0] * (N + 1)

total = 0

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

    total += prime[v]

print total