import sys
input=lambda: sys.stdin.readline().rstrip()
n,p=map(float,input().split())
n=int(n)
def f(n):
  D=[]
  for i in range(1,int(n**0.5)+1):
    if n%i==0:
      D.append(i)
      if i!=n//i:
        D.append(n//i)
  return len(D)
ans=0
for i in range(2,n+1):
  ans+=pow(1-p,f(i)-2)
print(ans)