#include using namespace std; #include #include #include #include #define FOR(i,a,b) for((i)=(a);i<(b);i++) typedef long long ll; typedef unsigned long long ull; typedef long double lb; /* ここからが本編 */ /* prime_rekkyo(prime,n)は  nまでの素数をまで2,3,…を格納していきます */ void prime_rekkyo(int prime[],int n) { int i,j,k; int a = 0; prime[a++] = 2; for(i=3;i<=n;i+=2) { k=0; for(j=3;j*j<=i;j+=2) { if(i%j==0) { k=1; break; } } if(k==0) prime[a++] = i; } } int main(void) { int i,j,k; int n; double p; double e[100000]; double sum = 0; int prime[10000]; prime_rekkyo(prime,100200); /* 100193は素数 */ cin >> n >> p; FOR(i,2,n+1) e[i] = 1.0; for(i=2; i <= n/2+1; i++) { for(j=2; i*j <= n; j++) { e[i*j] *= (1 - p); } } FOR(i,2,n+1) sum += e[i]; printf("%.14g\n",sum); return 0; }