Pagini recente » Cod sursa (job #111830) | Cod sursa (job #1438021) | Cod sursa (job #1861444) | Cod sursa (job #1222845) | Cod sursa (job #153399)
Cod sursa(job #153399)
/*
* Nu ne intereseaza numerele pare
* pentru ca 2 este singurul numar prim par
* si il tratam separat
* Pentru a tine minte ca pe pozitia i nu este un numar prim,
* marcam in vector bitul [i / 2]
* Daca in vector este marcat bitul i, atunci inseamna ca
* numarul 2*i+1 nu este prim
*/
#include <cstdio>
const int NMAX = 1 << 17;
unsigned char V[NMAX];
int main(void) {
freopen("ciur.in", "rt", stdin);
freopen("ciur.out", "wt", stdout);
int N, i, j, i2, nrprime = 0;
scanf(" %d", &N);
for (i = 3; i * i <= N; i += 2) {
if (V[i >> 4] & (1 << ((i >> 1) & 7))) continue;
++nrprime;
i2 = i + i;
for (j = i * i; j <= N; j += i2)
V[j >> 4] |= 1 << ((j >> 1) & 7);
}
for (; i <= N; i += 2) {
if (V[i >> 4] & (1 << ((i >> 1) & 7))) continue;
++nrprime;
}
printf("%d\n", nrprime + 1);
return 0;
}