#include using namespace std; double po(double p,int q){ double ans=1; for(int i=0;i>n>>p; vector divisor(n+1,0); for(int i=2;i<=n;i++){ for(int j=2;i*j<=n;j++) divisor[i*j]++; } double ans=0; for(int i=2;i<=n;i++) ans+=po(1-p,divisor[i]); printf("%.10lf\n",ans); return 0; }