#include #include #include #include #include #include #include #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 (int i=0;i<(n);i++) #define RREP(i,n) for (int i=(n)-1;i>=0;i--) #define PB push_back #define ALL(a) (a).begin(),(a).end() #define PII pair #define VS vector #define VD vector #define VL vector #define V(type) vector #define DEBUG(x) cout<<#x<<": "< 1 ){ double piv = backup_koho[i]; while ( i + j * piv <= koho.size() - 1 ) { if ( koho[i + j * piv] >= 1 ) { koho[i + j * piv] = p; ++j; continue; } koho[i + j * piv] *= p; ++j; } // } } } else { REP(i, koho.size()){ dst.push_back(koho[i]); int j = 0; if ( koho[i] > 1 ){ double piv = backup_koho[i]; while ( i + j * piv <= koho.size() - 1 ) { koho[i + j * piv] = 0; ++j; } // } } } } return dst; } int main(int argc, const char * argv[]) { int N; double P; VD result; double result_p = 0; cin >> N >> P; result = Eratosthenes(N, P, result_p); // cout << result_p << endl; double sum = 0; FOR(i, 0, result.size()){ if ( result[i] >= 1 ){ sum++; } else { sum += result[i]; } } if ( P == 0 ){ cout << result.size() << endl; } else { cout << sum << endl; } return 0; }