def main():
    n,p=input().split()
    n=int(n)
    p=float(p)
    ev=[1]*(n+1)
    ans=0
    for x in range(2,n+1):
        ans+=ev[x]
        for y in range(x*2,n+1,x):
            ev[y]*=1-p
    print(ans)

main()