Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

Prime numbers and primality tests

  • To: mathgroup at smc.vnet.net
  • Subject: [mg129449] Prime numbers and primality tests
  • From: johnfeth at gmail.com
  • Date: Mon, 14 Jan 2013 23:30:04 -0500 (EST)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-newout@smc.vnet.net
  • Delivered-to: mathgroup-newsend@smc.vnet.net

A straightforward way to test a prime number candidate is the Miller-Rabin test (sometimes called the Rabin-Miller test).  This well known and popular test is commonly executed 50 times on a candidate prime and has a proven probability of missing a non-prime of no more than 0.25 for each execution. Note that passing 50 Miller-Rabin tests (which is a de facto standard), the probability of non-primality is 0.25^50 ~ 7.9*10^-31, I'm satisfied that the number NextPrime gives me is "prime enough".  Mathematica uses the Miller-Rabin test, although it is not clear how many iterations are used.  As I understand it, Mathematica also the Lucas pseudo prime test on the Miller-Rabin output.

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've read that there are have been no counter-examples (viz., no non-primes that pass the the Lucas pseudo prime test) to numbers that pass the Lucas pseudo prime test, but then again, I've never found an oyster with a pearl inside.

Is the Miller-Rabin a better test that the Lucas pseudo prime test?


http://reference.wolfram.com/mathematica/tutorial/IntegerAndNumberTheoreticalFunctions.html

http://reference.wolfram.com/mathematica/tutorial/SomeNotesOnInternalImplementation.html

http://mathworld.wolfram.com/LucasPseudoprime.html



  • Prev by Date: Re: Making a palette to control dynamic variables?
  • Next by Date: Re: Export with append??
  • Previous by thread: Re: Problem with DateListPlot
  • Next by thread: Re: Prime numbers and primality tests