Cod sursa(job #153333)

Utilizator ProstuStefan-Alexandru Filip Prostu Data 10 martie 2008 14:10:53
Problema Ciurul lui Eratosthenes Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.72 kb
/*
 * 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 <= N; i += 2) {
		if (V[i >> 4] & (1 << ((i >> 1) & 7))) continue;
		++nrprime;

		for (j = i + (i2 = i + i); j <= N; j += i2)
			V[j >> 4] |= 1 << ((j >> 1) & 7);
	}

	printf("%d\n", nrprime + 1);

	return 0;
}