def solve():
    inp = input().split()
    N = int(inp[0])
    p = 1 - float(inp[1])

    s = [1.0] * (N + 1)
    for x in range(2, N + 1):
        for i in range(x + x, N + 1, x):
            s[i] *= p

    s[0] = 0
    s[1] = 0

    print(sum(s))

if __name__ == '__main__':
    solve()