import sys
input=lambda: sys.stdin.readline().rstrip()
n,p=map(float,input().split())
n=int(n)
D=[0]*(n+1)
for i in range(2,n+1):
  for j in range(1,n//i+1):
    D[i*j]+=1
ans=0
for i in range(2,n+1):
  ans+=pow(1-p,max(0,D[i]-1))
print(ans)