nsieve-bits
routine nsieve( bits : bitarray, m=0 )
{
bits[] = 1;
count = 0;
for( i = 2 : m ){
if( bits[i] ){
count ++;
for( k = i+i : i : m ) bits[k] = 0;
}
}
return count;
}
routine main( n = 2 )
{
m = 10000 * (1<
bits : bitarray = 1L << (m+1);
stdio.printf("Primes up to %8d %8d\n", m, nsieve(bits, m) );
m = 10000 * (1<<(n-1));
stdio.printf("Primes up to %8d %8d\n", m, nsieve(bits, m) );
m = 10000 * (1<<(n-2));
stdio.printf("Primes up to %8d %8d\n", m, nsieve(bits, m) );
}
view count 1159 times
created at 2009-02-20, 16:30 GMT
created at 2009-02-20, 16:30 GMT
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 |
fu: ... I forgot to say something about the plan for the whole new year in my previous reply. Well, besides w ... (Jan.19,01:40)
fu: ... Happy new dragon year (which will start from this sunday)! Actually, it was a busy month (I wish th ... (Jan.18,22:46)
ybabel: What's the plan for the new year ? Hello 'vry budy :- ) happy new year (when is the new year for you Fu ?) I saw you come back and comm ... (Jan.18,18:59)