#include #include namespace { using ModInt [[maybe_unused]] = atcoder::modint998244353; using Num [[maybe_unused]] = long long int; using Vec [[maybe_unused]] = std::vector; using Set [[maybe_unused]] = std::set; using Mset [[maybe_unused]] = std::multiset; using Edges [[maybe_unused]] = std::vector>; template using Q [[maybe_unused]] = std::queue; template using PQ [[maybe_unused]] = std::priority_queue, std::greater>; } void solve(std::istream& is, std::ostream& os) { Num n {0}; double p; is >> n >> p; Vec ps(n+9, 0); for(Num i{2}; i<=n; ++i) { for(Num j=i+i; j<=n; j+=i) { ps.at(j) += 1; } } double total {0}; for(Num i{2}; i<=n; ++i) { const auto cnt = ps.at(i); if (ps.at(i) == 0) { total += 1; } else { total += std::pow(1.0 - p, cnt); } } os << std::setprecision(20) << total << "\n"; } int main(void) { solve(std::cin, std::cout); return 0; }