N, p = raw_input().strip().split()
N = int(N)
p = 1.0 - float(p)
dp = [0] * (N+1)
for i in xrange(2, N+1):
    for j in xrange(2*i, N+1, i):
        dp[j] += 1
ans = sum(pow(p, x) for x in dp[2:])
print('%.7f' % ans)