n,p=list(raw_input().split()) N,P=int(n),float(p) l=[1.0]*(N+1) q=1-P for i in xrange(2,N+1): for j in xrange(i+i,N+1,i):l[j]*=q print(sum(l[2:]))