Pagini recente » Cod sursa (job #3232490) | Cod sursa (job #2618272) | Cod sursa (job #631459) | Cod sursa (job #1910622) | Cod sursa (job #1404006)
#include <fstream>
#include <iostream>
using namespace std;
//class PrimeNumbersSieve5
#define MAX 100000000/2/8+1;
char p[100000];
//(p[i>>3]&(1<<(i&7))) == 0 if 2*i + 1 is prime
int ciur(int n) {
int i, j, nr = 1;
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++;
return nr;
}
int main()
{
ifstream f("ciur.in");
ofstream g("ciur.out");
int i,n;
f >> n;
ciur(n);
for(i=1;i<=n;i++) if ((p[i>>3]&(1<<(i&7))) == 0) g << 2*i+1 << " ";
}