n, p = map(float,input().split())
ans = 0
n = int(n)
b = [0]*(n+1)

for i in range(2,n+1):
    
    ans += pow(1-p, b[i])
    
    j = 1
    while i*j <= n:
        b[i*j] += 1
        j += 1
print(ans)