The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Inactive » comp.soft-sys.math.mathematica

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Prime numbers and primality tests
Replies: 4   Last Post: Jan 22, 2013 11:18 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]

Posts: 18
Registered: 11/27/06
Re: Prime numbers and primality tests
Posted: Jan 22, 2013 11:18 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

The original post had this comment:
"It is interesting to note, however, that the Lucas pseudo prime
method of primality testing apparently does not have the handy
"feature" of the Miller-Rabin test, namely, the provable, and bounded
low probability of a wrong answer, from whence an estimate of
primality for any number can be made without finding a counter

I just found this interesting article:
"The Rabin-Monier Theorem for Lucas Pseudoprimes", by F. Arnault, in
Mathematics of Computation, April, 1997.
You can download it from

Theorem 1.3 in that article shows that the probability for each strong
Lucas pseudoprime test is 4/15 (analogous to the probability 1/4 for a
Miller-Rabin test).

The key points, though, are these: First, if N is a strong
pseudoprime base 2, it is more likely than the average number about
that size to also be a strong pseudoprime base 3 (and other bases).
So, if you make lists of strong pseudoprimes base 2 and base 3, those
lists will overlap, meaning that there are composite numbers that are
strong pseudoprimes bases 2 and 3 (and 5, and 7, ...).

Second, if you use any of the algorithms in the "Lucas Pseudoprime"
paper at
to make a list of strong Lucas pseudoprimes, then this list has no
known overlap with the list of strong pseudoprimes base 2. It has
been checked that any number up to at least 10^16 that passes both
types of these strong pseudoprime tests is a prime.

This is why many prime-testing functions (e.g, PrimeQ) do at least one
of each type of test, instead of just doing standard strong
pseudoprime tests to a bunch of different bases 2, 3, ... .

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.