#include using namespace std; double m_pow(double x,int n){ //べき乗 double ans=1; while(n!=0){ if(n&1) ans=ans*x; x=x*x; n=n>>1; } return ans; } int main(){ int N; cin >> N; double p; cin >> p; double ans=0; for(int i=2;i<=N;i++){ vector divisor; for(int j=1;j*j<=i;j++){ if(i%j==0){ divisor.push_back(j); if(j!=i/j) divisor.push_back(i/j); } } ans+=m_pow(1-p,divisor.size()-2); } cout << setprecision(10) << ans << endl; }