Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

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


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

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 ]
bobbaillie@frii.com

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
example!"

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 http://www.ams.org/journals/mcom/1997-66-218/S0025-5718-97-00836-3

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
http://mpqs.free.fr/LucasPseudoprimes.pdf
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]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.