WITH FILE OUTPUT
N       Time (seconds, approximate)
10M     < 0.8
100M    9
1B      92

WITH NO OUTPUT
N       Time (seconds, approximate)
10M     < 0.3
100M    6
1B      65

#include <stdio.h>
#include <time.h>
int main(int argc, char *argv[]) {
        clock_t time = clock();
        int N = atoi(argv[1]);
        uint p[N/32+1];
        int i, j, NP = 0;
        for (i = 0; i <= N>>5; i++) p[i] = 0;
        for (i = 2; i < sqrt(N); i++)
                if (0 == (p[i>>5] & (1 << ((i-1) & 31))))
                        for (j = i<<1; j < N; j += i)
                                p[j>>5] |= 1 << ((j-1) & 31);
        for (i = 2; i < N; i+=2)
                if (0 == (p[i>>5] & (1 << ((i-1) & 31)))) {
                        NP++;
                        //printf("\n%d", i);
                }
        printf("\nPRIMES: %d\n", NP+1);
        printf("TIME: %d\n", clock() - time);
        return 0;
}

# let t = Sys.time ()
  and a = Array.make 2_097_151 true in
      for i = 1 to Array.length a - 1 do
          if a.(i) then
              for j = 2 to Array.length a / (i + 1) do
                  a.(j * (i + 1) - 1) <- false
              done
      done;
      Printf.printf "%f\n" (Sys.time () -. t);;
3.875000
- : unit = ()
#

chris@frankenstein:~/Documents$ cat primes.ml
let t = Sys.time ()
  and a = Array.make 2_097_151 true in
      for i = 1 to Array.length a - 1 do
          if a.(i) then
              for j = 2 to Array.length a / (i + 1) do
                  a.(j * (i + 1) - 1) <- false
              done
      done;
      Printf.printf "%f\n" (Sys.time () -. t);;
chris@frankenstein:~/Documents$ ocamlopt primes.ml -o primes
chris@frankenstein:~/Documents$ ./primes
0.536032
chris@frankenstein:~/Documents$

chris@frankenstein:~/Documents$ cat primes2.sml
structure ArrayUtils =
struct
    fun fromToByDo(array, startIndex, endIndex, incrementBy, func) =
        if startIndex <= endIndex then
           ( func(array, startIndex, Array.sub(array, startIndex));
             fromToByDo(array, startIndex + incrementBy, endIndex, incrementBy, func) )
        else ()
end;

let
    val t = Time.now()
    val a = Array.array(2097151, true)
in
    Array.update(a, 0, false);

    Array.appi(fn (i, x) =>
        if x then
            let
                val number = i + 1
            in
                ArrayUtils.fromToByDo(a, number * 2 - 1, Array.length a - 1, number,
                                      fn (a, i, _) => Array.update(a, i, false))
            end
        else ()) a;

    let
        val t' = Time.now()
        val diff = Time.-(t', t)
    in
        print ((IntInf.toString(Time.toMilliseconds diff)) ^ "\n") 
    end
end;

chris@frankenstein:~/Documents$ mlton -output sml_primes2 primes2.sml
chris@frankenstein:~/Documents$ ./sml_primes2
584
chris@frankenstein:~/Documents$ ./sml_primes2
584
chris@frankenstein:~/Documents$ ./sml_primes2
584
chris@frankenstein:~/Documents$