#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; using P = pair; constexpr int INF = 1001001001; constexpr int mod = 1000000007; // constexpr int mod = 998244353; template inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int n; double p; cin >> n >> p; vector dp(n + 1, 1.0); for(int i = 2; i < n; ++i){ for(int j = i + i; j <= n; j += i){ dp[j] *= 1.0 - p; } } double ans = 0.0; for(int i = 2; i <= n; ++i) ans += dp[i]; cout << fixed << setprecision(15); cout << ans << endl; return 0; }