This site is a static rendering of the Trac instance that was used by R7RS-WG1 for its work on R7RS-small (PDF), which was ratified in 2013. For more information, see Home.
Source for wiki PrimesGauche version 2
author
cowan
comment
ipnr
127.11.51.1
name
PrimesGauche
readonly
0
text
{{{
#!html
<p>This module provides utilities related to prime numbers.
</p><dl>
<dt><a name="index-_002aprimes_002a"></a><code>all-primes</code> [Variable]</dt>
<dd><p>An infinite lazy sequence of primes.
</p>
<table><tr><td> </td><td><pre class="example">;; show 10 prime numbers from 100-th one.
(lseq-take (lseq-drop all-primes 100) 10)
⇒ (547 557 563 569 571 577 587 593 599 601)
</pre></td></tr></table>
</dd></dl>
<dl>
<dt><a name="index-reset_002dprimes"></a> <code>reset-primes</code></dt>
<dd><p>Once you take a prime out of <code>all-primes</code>, all smaller primes
have been calculated and remain in memory, since the
head of the sequence is held in <code>all-primes</code>. Sometimes you know
you need no more prime numbers and you wish the ones that have already been computed
to be garbage-collected. Calling <code>reset-primes</code> rebinds
<code>all-primes</code> to a fully unrealized lazy sequence, allowing the previously
realized primes to be garbage collected.
</p></dd></dl>
<dl>
<dt><a name="index-primes"></a><code>primes</code></dt>
<dd><p>Returns a fresh lazy sequence of primes. It is useful when
you need certain for a short period of time — if you don’t keep
a reference to the head of the returned sequence, it will be garbage
collected after you’re done with the primes.
(Note that the calculation of a prime number needs the
sequence of primes from the beginning,
so even if your code only keeps a reference
into the middle of the sequence, the entire sequence will be stored — you
have to release all
references in order to allow the sequence to be garbage collected.)
</p>
<p>On the other hand,
each sequence returned by <code>primes</code> is computed individually,
duplicating computation.
</p>
<p>The rule of thumb is: if you use primes repeatedly throughout
the program, just use <code>all-primes</code> and you’ll save computation.
If you need primes just once, call <code>primes</code> and abandon the result
and you’ll save space.
</p></dd></dl>
<a name="Testing-primality"></a>
<h3 class="subheading">Testing primality</h3>
<dl>
<dt><a name="index-small_002dprime_003f"></a><code>small-prime?</code><i> n</i></dt>
<dd><p>For positive integers
less than <code>small-prime-bound</code>, this procedure
determines if <var>n<i> is prime or not, quickly and deterministically.
If <var>n</var> is greater than or equal to this bound, this procedure returns <code>#f</code>.
</p>
<p>This can be used to quickly filter out known primes; it never returns
<code>#t</code> on composite numbers, but it may return <code>#f</code> on
sufficiently large prime numbers).
The Miller-Rabin test below can tell if the input is definitely composite,
but it may return <code>#t</code> on some composite numbers.
</p></dd></dl>
<dl>
<dt><a name="index-_002asmall_002dprime_002dbound_002a"></a><code>small-prime-bound</code> [Variable]</dt>
<dd><p>For all positive integers below this value
<code>small-prime?</code> can determine whether it is a prime or not.
</p></dd></dl>
<dl>
<dt><a name="index-miller_002drabin_002dprime_003f"></a><code>miller-rabin-prime?</code><i> n [ num-tests ]</i></dt>
<dd><p>Check if an exact integer <var>n</var> is a prime number, using
probabilistic Miller-Rabin algorithm (<var>n</var> must be greater than 1).
If this procedure returns <code>#f</code>,
<var>n</var> is a composite number. If this procedure returns <code>#t</code>,
<var>n</var> is <em>likely</em> a prime, but there’s a small probability
that it is a false positive.
</p>
<p>Note that if <var>n</var> is smaller than <code>small-prime-bound</code>, the algorithm is
deterministic; if it returns <code>#t</code>, <var>n</var> is certainly a prime.
</p>
<p>If <var>n</var> is greater than or equal to
<code>small-prime-bound</code>, we use a probabilistic test.
We choose a random base integer
to perform the Miller-Rabin test up to <var>num-tests</var> (7 times by default).
The probability
of returning <code>#t</code> for a composite number
is at most <code>(expt 4 (- num-tests))</code>.
</p></dd></dl>
<dl>
<dt><a name="index-bpsw_002dprime_003f"></a><code>bpsw-prime?</code><i> n</i></dt>
<dd><p>Check if an exact integer <var>n</var> is a prime number, using
the Baillie-PSW primality test
(<a href="http://www.trnicely.net/misc/bpsw.html">http://www.trnicely.net/misc/bpsw.html</a>). It is deterministic,
and is known to return a definitive answer for all numbers less than 2<super>64</super>.
For larger integers this can return <code>#t</code> on a composite number,
although no such number has been found yet. This procedure never returns <code>#f</code>
on a prime number.
</p>
<p>This test is slower than Miller-Rabin but fast enough for casual use,
so it is handy when you want a definitive answer below the above range.
</p></dd></dl>
<a name="Factorization"></a>
<h3 class="subheading">Factorization</h3>
<dl>
<dt><a name="index-naive_002dfactorize"></a><code>naive-factorize</code><i> n [ divisor-limit ]</i></dt>
<dd><p>Factorize a positive integer <var>n</var> by trying to divide it into
all primes up to <code>(sqrt n)</code>. Returns a list of prime factors,
smallest first.
</p>
<table><tr><td> </td><td><pre class="example">(naive-factorize 142857)
⇒ (3 3 3 11 13 37)
</pre></td></tr></table>
<p>Although this is a pretty naive method, this works well as far as
any of <var>n</var>’s factors are not larger than about ten million.
</p>
<table><tr><td> </td><td><pre class="example">(naive-factorize 3644357367494986671013))
⇒ (10670053 10670053 32010157)
</pre></td></tr></table>
<p>If <var>n</var> includes any larger prime factors,
the performance becomes abysmal.</p>
<p>Alternatively, providing the <var>divisor-limit</var> argument specifies
the upper bound of the prime number to be tried. If it is given,
<code>naive-factorize</code> returns a factor <var>f</var> unchanged if it can’t be
divided by any primes less than or equal to <var>divisor-limit</var>.
So, the last element of the returned list may be a composite number.
This is useful for excluding trivial factors before applying more sophisticated
factorizing algorithms.
</p>
<table><tr><td> </td><td><pre class="example">(naive-factorize 825877877739 1000)
⇒ (3 43 6402154091)
;; whereas
(naive-factorize 825877877739)
⇒ (3 43 4591 1394501)
</pre></td></tr></table>
<p>The procedure also memoizes the results on smaller values of <var>n</var> to make
things faster.
</p></dd></dl>
<dl>
<dt><a name="index-mc_002dfactorize"></a><code>mc-factorize</code><i> n</i></dt>
<dd><p>Factorize a positive integer <var>n</var> using the algorithm
described in
R. P. Brent, An improved Monte Carlo factorization algorithm, BIT 20 (1980), 176-184. <a href="http://maths-people.anu.edu.au/~brent/pub/pub051.html">http://maths-people.anu.edu.au/~brent/pub/pub051.html</a>.
</p>
<p>This one is capable of handling much larger factors than
<code>naive-factorize</code>, somewhere around 10<super>20</super> or so.
</p>
<p>Since this method is probabilistic, the execution time may vary
on the same <var>n</var>. But it will always return the definitive
results as long as every prime factor of <var>n</var> is smaller than an
implementation-specified limit. If <var>n</var> contains a prime factor greater than
the limit, the procedure may loop forever.
</p></dd></dl>
<a name="Miscellaneous"></a>
<h3 class="subheading">Miscellaneous</h3>
<dl>
<dt><a name="index-jacobi"></a><code>jacobi</code><i> a n</i></dt>
<dd><p>Calculates the Jacobi symbol <code>(<var>a</var>/<var>n</var>)</code>
(<a href="http://en.wikipedia.org/wiki/Jacobi_symbol">http://en.wikipedia.org/wiki/Jacobi_symbol</a>).
</p></dd></dl>
<dl>
<dt><a name="index-totient"></a><code>totient</code><i> n</i></dt>
<dd><p>Euler’s totient function of nonnegative integer <var>n</var>.
</p>
<p>The current implementation relies on <code>mc-factorize</code> above,
so it may take very long if <var>n</var> contains large prime factors.
</p></dd></dl>
}}}
time
2016-07-11 10:38:43
version
2