INPUT = raw_input().split()
N,P = [int(INPUT[0]),float(INPUT[1])]
Pmap = [0,0]+[1.0] * (N-1)
for sosu in range(2,N):
   mul=2
   while sosu*mul<=N:
     Pmap[sosu*mul]*=1-P
     mul+=1
print(sum(Pmap))