Cod sursa(job #232861)

Utilizator crawlerPuni Andrei Paul crawler Data 16 decembrie 2008 09:59:21
Problema Ciurul lui Eratosthenes Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 0.55 kb
#include <stdio.h>
#include <math.h>

unsigned a[2000100];

int main()
{
    freopen("ciur.in","r",stdin);
    freopen("ciur.out","w",stdout);
    
    int n,ret=0,lim;
    
    scanf("%d",&n);
    if (n>=2) ++ret;
    lim = (int)sqrt(n);
    for (unsigned i=4;i<=n;i+=2) a[i>>32u]& = 1<<(i&(1u<<32u));
    for (unsigned i=3;i<=n;i+=2) if(a[i>>32u]&(i&(1u<<32u)) == 0)
    {
        ++ret;
        if (i<=lim) for (unsigned j=i*i;j<=n;j+=2*i) a[j>>32u]& = 1<<(j&(1u<<32u));
    }
    
    printf("%d\n",ret);
        
    return 0;    
}