Cod sursa(job #547163)

Utilizator ELHoriaHoria Cretescu ELHoria Data 5 martie 2011 22:23:11
Problema Ciurul lui Eratosthenes Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.56 kb
#include<fstream>
using namespace std;

ifstream f("ciur.in");
ofstream g("ciur.out");

//(p[i>>3]&(1<<(i&7))) == 0 if 2*i + 1 is prime

char p[2000000];
int main(){
  int i, j, nr = 1,n;
  f>>n;
  for (i = 1; ((i * i) << 1) + (i << 1) <= n; i += 1) {      
    if ((p[i >> 3] & (1 << (i & 7))) == 0) {
      for (j = ((i * i) << 1) + (i << 1); (j << 1) + 1 <= n; j += (i << 1) + 1) {
        p[j >> 3] |= (1 << (j & 7));
      }
    }
  }
  for (i = 1; 2 * i + 1 <= n; ++i)  
       if ((p[i >> 3] & (1 << (i & 7))) == 0) 
           nr++;
  g<<nr;
}