# 愚直実装 N, prob = map(float, input().split()) N = int(N) survive = [1]*(N+1) survive[0]=0 survive[1]=0 for p in range(2, N+1): for j in range(p*2, N+1, p): #survive[j] *= (1-prob)*survive[p] これだとワークしない survive[j] *= (1-prob) #print('p', p, survive) ans = sum(survive) print(ans)