# coding: utf-8
import array, bisect, collections, heapq, itertools, math, random, re, string, sys, time
sys.setrecursionlimit(10 ** 7)
inf = 10 ** 20
mod = 10 ** 9 + 7


def II(): return int(input())
def ILI(): return list(map(int, input().split()))
def IAI(LINE): return [ILI() for __ in range(LINE)]
def IDI(): return {key: value for key, value in ILI()}


def solve(N, p):
    dp = [1.0] * (N + 1)
    p_minus = 1.0 - p
    for i in range(2, N + 1):
        j = i * 2
        while j <= N:
            dp[j] *= p_minus
            j += i

    ans = sum(dp[2:])

    return ans


def main():
    N, p = input().split()
    N = int(N)
    p = float(p)
    print(solve(N, p))


if __name__ == "__main__":
    main()