primesieve is a program and C/C++ library that generates primes using a highly optimized sieve of Eratosthenes implementation. It counts the primes below 10^10 in just 0.45 seconds on an Intel Core i7-6700 CPU (4 x 3.4GHz). primesieve can generate primes and prime k-tuplets up to 2^64.
- Homepage: http://primesieve.org
- Binaries: http://primesieve.org/downloads
- API: http://primesieve.org/api
primesieve generates primes using the segmented sieve of Eratosthenes with wheel factorization, this algorithm has a complexity of operations and uses space. More precisely primesieve's memory usage per thread is about bytes.
primesieve is very portable, it compiles with every C++ compiler and runs on most CPU architectures out there. The parallelization is implemented using OpenMP. The primesieve GUI application (not built by default) uses the Qt framework.
primesieve is also a library, it supports C++ and C directly and there are bindings available for a few other programming languages. The author's primecount program uses libprimesieve extensively.
Download primesieve-5.6.0.tar.gz and extract it. Then build primesieve using:
$ ./configure
$ make
$ sudo make install
If you have downloaded a zip archive from GitHub then Autotools
(a.k.a. GNU Build System) must be installed and autogen.sh
must
be executed once. To install Autotools install
GNU Autoconf,
GNU Automake and
GNU Libtool
using your operating system's package manager.
$ ./autogen.sh
$ ./configure
$ make
$ sudo make install
To enable building the example programs use:
$ ./configure --enable-examples
Open a Visual Studio Command Prompt, cd into the primesieve directory and run:
> nmake -f Makefile.msvc
To build the example programs use:
> nmake -f Makefile.msvc examples
The primesieve console application can print and count primes and prime k-tuplets and find the nth prime.
# Print the primes below 1000000
$ primesieve 1000000 --print
# Print the twin primes below 1000000
$ primesieve 1000000 --print=2
# Count the primes below 1e10 using all CPU cores
$ primesieve 1e10 --count
# Count the primes within [1e10, 2e10] using 4 threads
$ primesieve 1e10 2e10 --count --threads=4
# Print an option summary
$ primesieve --help
After having installed primesieve you can easily use it in your C++ program, below is an example. primesieve's API is documented online at http://primesieve.org/api.
#include <primesieve.hpp>
#include <iostream>
#include <vector>
int main()
{
// store the primes below 1000
std::vector<int> primes;
primesieve::generate_primes(1000, &primes);
primesieve::iterator pi;
uint64_t prime;
// iterate over the primes below 10^9
while ((prime = pi.next_prime()) < 1000000000)
std::cout << prime << std::endl;
return 0;
}
All of primesieve's functions are exposed as C API (C99 or later) via
the primesieve.h
header. You can browse primesieve's C API
online at http://primesieve.org/api.
#include <primesieve.h>
#include <stdint.h>
#include <stdio.h>
int main()
{
uint64_t start = 0;
uint64_t stop = 10000;
size_t i;
size_t size;
/* get an array with primes below 10000 */
int* primes = (int*) primesieve_generate_primes(start, stop, &size, INT_PRIMES);
for (i = 0; i < size; i++)
printf("%i\n", primes[i]);
/* deallocate primes array generated using primesieve */
primesieve_free(primes);
return 0;
}
Unix-like operating systems
$ c++ -O2 primes.cpp -lprimesieve
$ cc -O2 primes.c -lprimesieve
Note that if you have built primesieve yourself the default installation path
is /usr/local/lib
which is not part of LD_LIBRARY_PATH
on many
systems. Hence you need to export some additional environment variables:
$ export LIBRARY_PATH=/usr/local/lib:$LIBRARY_PATH
$ export LD_LIBRARY_PATH=/usr/local/lib:$LD_LIBRARY_PATH
$ export CPLUS_INCLUDE_PATH=/usr/local/include:$CPLUS_INCLUDE_PATH
$ export C_INCLUDE_PATH=/usr/local/include:$C_INCLUDE_PATH
Microsoft Visual C++ (Windows)
> cl /O2 /EHsc primes.cpp /I primesieve\include /link primesieve\primesieve.lib
primesieve supports C++ and C directly, and has bindings available for a few other languages:
Python: | primesieve-python |
Ruby: | primesieve-ruby |
Many thanks to the developers of these bindings!