Pagini recente » Cod sursa (job #1376545) | Cod sursa (job #44444) | Cod sursa (job #3127613) | Cod sursa (job #2670372) | Cod sursa (job #2228860)
#include <stdio.h>
int p[(100000000 >> 6) + 1];
unsigned sieve(unsigned n)
{
register unsigned i, j, k, nr = 1 + (n > 2);
for(i = 2; ((i * i) << 1) + (i << 1) <= n - 2; i += 3)
{
if((p[i >> 5] & (1 << (i & 31))) == 0)
{
for(j = ((i * i) << 1) + (i << 1); (j << 1) + 1 <= n; j += (i << 1) + 1)
{
p[j >> 5] |= (1 << (j & 31));
}
}
k = i + 1;
if((p[k >> 5] & (1 << (k & 31))) == 0)
{
for(j = ((k * k) << 1) + (k << 1); (j << 1) + 1 <= n; j += (k << 1) + 1)
{
p[j >> 5] |= (1 << (j & 31));
}
}
}
for(i = 2; (i << 1) + 1 <= n - 2; i += 3)
{
nr += ((p[ i >> 5] & (1 << ( i & 31))) == 0) +
((p[(i + 1) >> 5] & (1 << ((i + 1) & 31))) == 0);
}
return nr;
}
int main()
{
freopen("ciur.in", "r", stdin);
freopen("ciur.out", "w", stdout);
unsigned N; scanf("%d", &N);
printf("%u ", sieve(N));
return 0;
}