#include #define DEBUG(x) cerr << #x << "=" << x << endl #define BINARY(x) static_cast >(x) #define rep(i,n) for(int i=0;i<(int)(n);i++) #define REP(i,m,n) for (int i=m;i<(int)(n);i++) #define if_range(x, y, w, h) if (0<=(int)(x) && (int)(x)<(int)(w) && 0<=(int)(y) && (int)(y)<(int)(h)) #define is_in(x, y, w, h) (0<=(int)(x) && (int)(x)<(int)(w) && 0<=(int)(y) && (int)(y)<(int)(h)) #define ALL(a) (a).begin(),(a).end() int dx[4]={0, -1, 1, 0}, dy[4]={-1, 0, 0, 1}; using namespace std; typedef long long ll; typedef vector > VVI; const int INF = 1e9; typedef pair PII; int d[1000001]; int main() { int N; cin >> N; double P; cin >> P; for (int i=2; i<=N; i++) { for (int j=i*2; j<=N; j+=i) { d[j]++; } } double ans = 0; for (int i=2; i<=N; i++) { ans += pow(1-P, d[i]); } cout << fixed << setprecision(10); cout << ans << endl; }