Pagini recente » Cod sursa (job #578674) | Cod sursa (job #204394) | Cod sursa (job #1091654) | Cod sursa (job #2329585) | Cod sursa (job #2229063)
#include <stdio.h>
int p[2000000];
int sieve(int n)
{
int i, j, nr = (n > 1) + (n > 2);
for(i = 5; i * i <= n; i += 4)
{
if((p[i >> 5] & (1 << (i & 31))) == 0)
{
for(j = 5 * i; j <= n; j += 4 * i)
{
p[j >> 5] |= (1 << (j & 31));
if((j += 2 * i) <= n) p[j >> 5] |= (1 << (j & 31));
}
}
i += 2;
if(i * i <= n)
{
if((p[i >> 5] & (1 << (i & 31))) == 0)
{
for(j = 5 * i; j <= n; j += 4 * i)
{
p[j >> 5] |= (1 << (j & 31));
if((j += 2 * i) <= n) p[j >> 5] |= (1 << (j & 31));
}
}
}
}
for(i = 5; i <= n; i += 4)
{
if((p[i >> 5] & (1 << (i & 31))) == 0)
{
++nr;
}
i += 2;
if(i <= n)
{
if((p[i >> 5] & (1 << (i & 31))) == 0)
{
++nr;
}
}
}
return nr;
}
int main()
{
freopen("ciur.in", "r", stdin);
freopen("ciur.out", "w", stdout);
int N; scanf("%d", &N);
printf("%d ", sieve(N));
return 0;
}