from math import sqrt,sin,cos,tan,ceil,radians,floor,gcd,exp,log,log10,log2,factorial
import heapq
import bisect
from collections import deque,Counter,defaultdict
from itertools import permutations,combinations
from decimal import Decimal,ROUND_HALF_UP
from functools import lru_cache
import sys
sys.setrecursionlimit(10**8)

N, p = map(float, input().split())
N = int(N)

ans = 0
memo = defaultdict(int)

for i in range(2, N+1):
    for j in range(i+i, N+1, i):  #三つ目でstep指定
        memo[j] += 1
    ans += (1-p)**memo[i]

print(ans)