Cod sursa(job #222980)

Utilizator alex.cepoiAlexandru Cepoi alex.cepoi Data 26 noiembrie 2008 16:12:27
Problema Ciurul lui Eratosthenes Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.48 kb
#include <cstdio>
#define nmax 2000000

int N;
char a[nmax/16+1];

int main ()
{
	freopen ("ciur.in", "r", stdin);
	freopen ("ciur.out", "w", stdout);

	scanf ("%d", &N);


	int i, nr=1;
	for (i=1; i*(i+1)<=N/2; ++i)
		if (!(a[i/8] & (1<<(i%8))))
		{
			nr++;
			for (int j=2*i+1; j*(2*i+1)<=N; j+=2)
				a[(j*(2*i+1)-1)/2 /8] |= (1<<((j*(2*i+1)-1)/2)%8);
		}

	for (; 2*i+1<=N; ++i)
		if (!(a[i/8] & (1<<(i%8)))) nr++;

	printf ("\n%d\n", nr);
	return 0;
}