#coding: utf-8
##yuki_144

s=raw_input().split()
N=int(s[0])
p=float(s[1])
t=[1.0 for i in xrange(10**6+1)]
res=0
for i in xrange(2,N+1):
 for j in xrange(i*2,N+1,i):
  t[j]*=(1.0-p)
 res+=t[i]

print "%.7f" % res