import numpy as np N,p = input().split() N = int(N); p = float(p) dp = np.zeros(N+1, np.int32) for i in range(2,N+1): dp[i:N+1:i]+=1 ans = 0 for i in range(2,N+1): ans += (1-p)**(dp[i]-1) print(ans)