summaryrefslogtreecommitdiffhomepage
path: root/misc/examples/prime.c
diff options
context:
space:
mode:
author_Tradam <[email protected]>2023-09-08 01:29:47 +0000
committerGitHub <[email protected]>2023-09-08 01:29:47 +0000
commit3c76c7f3d5db3f9586a90d03f8fbb02d79de9acd (patch)
treeafbe4b540967223911f7c5de36559b82154f02f3 /misc/examples/prime.c
parent0841165881871ee01b782129be681209aeed2423 (diff)
parent1a72205fe05c2375cfd380dd8381a8460d9ed8d1 (diff)
downloadSTC-modified-modified.tar.gz
STC-modified-modified.zip
Merge branch 'stclib:master' into modifiedHEADmodified
Diffstat (limited to 'misc/examples/prime.c')
-rw-r--r--misc/examples/prime.c54
1 files changed, 0 insertions, 54 deletions
diff --git a/misc/examples/prime.c b/misc/examples/prime.c
deleted file mode 100644
index d0887353..00000000
--- a/misc/examples/prime.c
+++ /dev/null
@@ -1,54 +0,0 @@
-#include <stdio.h>
-#include <math.h>
-#include <time.h>
-#include <stc/cbits.h>
-#include <stc/algo/filter.h>
-#include <stc/algo/crange.h>
-
-
-cbits sieveOfEratosthenes(int64_t n)
-{
- cbits bits = cbits_with_size(n/2 + 1, true);
- int64_t q = (int64_t)sqrt((double) n) + 1;
- for (int64_t i = 3; i < q; i += 2) {
- int64_t j = i;
- for (; j < n; j += 2) {
- if (cbits_test(&bits, j>>1)) {
- i = j;
- break;
- }
- }
- for (int64_t j = i*i; j < n; j += i*2)
- cbits_reset(&bits, j>>1);
- }
- return bits;
-}
-
-int main(void)
-{
- int64_t n = 1000000000;
- printf("Computing prime numbers up to %" c_ZI "\n", n);
-
- clock_t t1 = clock();
- cbits primes = sieveOfEratosthenes(n + 1);
- int64_t np = cbits_count(&primes);
- clock_t t2 = clock();
-
- printf("Number of primes: %" c_ZI ", time: %f\n\n", np, (float)(t2 - t1) / (float)CLOCKS_PER_SEC);
- puts("Show all the primes in the range [2, 1000):");
- printf("2");
- c_forrange (i, 3, 1000, 2)
- if (cbits_test(&primes, i>>1)) printf(" %lld", i);
- puts("\n");
-
- puts("Show the last 50 primes using a temporary crange generator:");
- c_forfilter (i, crange, crange_obj(n - 1, 0, -2),
- cbits_test(&primes, *i.ref/2) &&
- c_flt_take(i, 50)
- ){
- printf("%lld ", *i.ref);
- if (c_flt_getcount(i) % 10 == 0) puts("");
- }
-
- cbits_drop(&primes);
-}