from collections import defaultdict as D N,P=[float(e)for e in input().split()] h=D(lambda:1) r=0 for i in range(2,int(N)+1): r+=h[i] for j in range(i*2,int(N)+1,i):h[j]*=1-P print(r)