n,p=map(float,input().split())
n=int(n)
ct=[0]*(n+1)
for i in range(1,n+1):
  for j in range(n//i):
    ct[i*(j+1)]+=1
ans=0
for i in range(2,n+1):
  ans+=(1-p)**(ct[i]-2)
print(ans)