def prime_sieve(n): # 素数リスト(エラトステネスの篩) # N以下の素数のリストを返す is_prime = [True for i in range(n + 1)] is_prime[0] = False is_prime[1] = False for i in range(4, n + 1, 2): is_prime[i] = False for i in range(3, int(n**0.5 + 1), 2): if is_prime[i]: for j in range(i * i, n + 1, i): is_prime[j] = False return [i for i in range(n + 1) if is_prime[i]] def get_divisors(n): # 約数列挙 # リストがソートされていないことに注意せよ。 divisors = set() for i in range(1, int(n**0.5) + 1): if n % i == 0: divisors.add(i) divisors.add(n // i) return list(divisors) N, p = input().split() N = int(N) p = float(p) primes = prime_sieve(N) primes.sort() e = 0 for i in range(2, N + 1): if i in primes: e += 1 else: e += (1 - p)**(len(get_divisors(i)) - 2) print(e)