#include #include #include #include #include #include #include #include #include #include #include #include #include #include #define endl '\n' #define ALL(a) (a).begin(),(a).end() #define SZ(a) int((a).size()) #define FOR(i,a,b) for(int i=(a);i<(b);++i) #define RFOR(i,a,b) for (int i=(b)-1;i>=(a);i--) #define REP(i,n) FOR(i,0,n) #define RREP(i,n) for (int i=(n)-1;i>=0;i--) #define DEBUG(x) cout<<#x<<": "< P; typedef long long int LL; typedef pair LP; int yakusu[1000000+1][2]; int main() { ios::sync_with_stdio(false); cin.tie(0); int N; long double p; cin>>N>>p; FOR(i,2,N+1){ if(yakusu[i][0]==0){ yakusu[i][0]=i; yakusu[i][1]=1; for(int j=2;j*i<=N;j++){ if(yakusu[i*j][0]==0){ yakusu[i*j][0]=i; yakusu[i*j][1]=j; } } } } long double ans=0; FOR(i,2,N+1){ //DEBUG(i); int d=1; int po=i; int pre=yakusu[po][0]; int co=1; po=yakusu[po][1]; while(po!=1){ /*DEBUG(po); DEBUG(pre); DEBUG(co);*/ if(pre==yakusu[po][0]){ co++; }else{ d*=(co+1); co=1; pre=yakusu[po][0]; } po=yakusu[po][1]; } //DEBUG(co); d*=(co+1); ans+=pow(1-p,d-2); //cout<