N,P=map(float,input().split())
M=int(N)+1
l=[1.0]*M
p=1-P
for i in range(2,M):
	for j in range(i+i,M,i):l[j]*=p
print(sum(l[2:]))