#include #include #include using namespace std; int main() { int N; cin >> N; double P; cin >> P; double res = 0.0; vector divided(N+1); for (int i=2;i<=N;i++) { res += pow(1.0-P, divided[i]); for (int j=i;j<=N;j+=i) { divided[j]++; } } cout << res << endl; }