Pagini recente » Cod sursa (job #771072) | Cod sursa (job #2317659) | Cod sursa (job #2382812) | Cod sursa (job #2952672) | Cod sursa (job #2228821)
#include <stdio.h>
char p[(20000000 / 8) + 10];
unsigned sieve(unsigned n)
{
register unsigned i, j, k, nr = 1 + (n > 2);
for(i = 5; i * i <= n - 2; i += 6)
{
if((p[i >> 3] & (1 << (i & 7))) == 0)
{
for(j = i * i; j <= n; j += i << 1)
{
p[j >> 3] |= (1 << (j & 7));
}
}
k = i + 2;
if((p[k >> 3] & (1 << (k & 7))) == 0)
{
for(j = k * k; j <= n; j += k << 1)
{
p[j >> 3] |= (1 << (j & 7));
}
}
}
for(i = 5; i <= n - 2; i += 6)
{
nr += ((p[ i >> 3] & (1 << ( i & 7))) == 0) +
((p[(i + 2) >> 3] & (1 << ((i + 2) & 7))) == 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;
}