import random, sys

strs = raw_input().split()
N = int(strs[0])
p = float(strs[1])

d = [0] * (N + 1)

for i in xrange(1, N + 1):
	for j in xrange(i, N + 1, i):
		d[j] += 1

ans = 0.0
for i in xrange(2, N + 1):
	ans += pow(1 - p, d[i] - 2)
print ans