#include #include #include #include using namespace std; using ld = long double; int main() { int n; ld p; cin >> n >> p; vector v(100, 1); for (int i = 1; i < 100; i++)v[i] = v[i - 1] * (1 - p); ld ret = 1; for (int i = 3; i <= n; i++) { ld count = 0; for (int j = 2; j * j <= i; j++) { if (i % j == 0) { count++; if (j * j != i)count++; } } ret += v[count]; } cout << fixed << setprecision(15) << ret << endl; return 0; }