Lucas pseudoprime
Lucas pseudoprimes and Fibonacci pseudoprimes are composite integers that pass certain tests which all primes and very few composite numbers pass: in this case, criteria relative to some Lucas sequence.
Contents
 1 BruckmanLucas pseudoprimes
 2 BaillieWagstaffLucas pseudoprimes
 3 Lucas probable primes and pseudoprimes
 4 Strong Lucas pseudoprimes
 5 Implementing a Lucas probable prime test
 6 Comparison with the Miller–Rabin primality test
 7 Fibonacci pseudoprimes
 8 Pell pseudoprimes
 9 References
 10 External links
BruckmanLucas pseudoprimes
The Lucas pseudoprimes as defined by Bruckman are as follows.^{[1]} The Lucas numbers may be defined as:
(where n belongs to the natural numbers). For any prime p, L_{p} is congruent to 1 mod p. The converse is almost always true, that is, if n is composite, L_{n} is not congruent to 1 mod n in most cases. The exceptions are the BruckmanLucas pseudoprimes, of which the first few are n = 705, 2465, 2737, 3745, 4181, 5777, 6721, ...^{[2]}
BaillieWagstaffLucas pseudoprimes
Baillie and Wagstaff define Lucas pseudoprimes differently.^{[3]} Given integers P and Q, where P > 0 and , let U_{k}(P, Q) and V_{k}(P, Q) be the corresponding Lucas sequences.
Let n be a positive integer and let be the Jacobi symbol. We define
If n is a prime such that the greatest common divisor of n and Q (that is, GCD(n, Q)) is 1, then the following congruence condition holds:

(1)
If this equation does not hold, then n is not prime. If n is composite, then this equation usually does not hold.^{[3]} These are the key facts that make Lucas sequences useful in primality testing.
Some good references are chapter 8 of the book by Bressoud and Wagon (with Mathematica code),^{[4]} pages 142–152 of the book by Crandall and Pomerance,^{[5]} and pages 53–74 of the book by Ribenboim.^{[6]}
Lucas probable primes and pseudoprimes
A Lucas probable prime for a given (P, Q) pair is any positive integer n for which equation (1) above is true (see,^{[3]} page 1398).
A Lucas pseudoprime for a given (P, Q) pair is a positive composite integer n for which equation (1) is true (see,^{[3]} page 1391).
A Lucas probable prime test is most useful if D is chosen such that the Jacobi symbol is −1 (see pages 1401–1409 of,^{[3]} page 1024 of ,^{[7]} or pages 266–269 of ^{[4]} ). This is especially important when combining a Lucas test with a strong pseudoprime test, such as the BailliePSW primality test. Typically implementations will use a parameter selection method that ensures this condition (e.g. the Selfridge method recommended in ^{[3]} and described below).
If then equation (1) becomes

(2)
If congruence (2) is false, this constitutes a proof that n is composite.
If congruence (2) is true, then n is a Lucas probable prime. In this case, either n is prime or it is a Lucas pseudoprime. If congruence (2) is true, then n is likely to be prime (this justifies the term probable prime), but this does not prove that n is prime. As is the case with any other probabilistic primality test, if we perform additional Lucas tests with different D, P and Q, then unless one of the tests proves that n is composite, we gain more confidence that n is prime.
Examples: If P = 3, Q = −1, and D = 13, the sequence of U's is OEIS: A006190: U_{0} = 0, U_{1} = 1, U_{2} = 3, U_{3} = 10, etc.
First, let n = 19. The Jacobi symbol is −1, so δ(n) = 20, U_{20} = 6616217487 = 19·348221973 and we have
Therefore, 19 is a Lucas probable prime for this (P, Q) pair. In this case 19 is prime, so it is not a Lucas pseudoprime.
For the next example, let n = 119. We have = −1, and we can compute
However, 119 = 7·17 is not prime, so 119 is a Lucas pseudoprime for this (P, Q) pair. In fact, 119 is the smallest pseudoprime for P = 3, Q = −1.
We will see below that, in order to check equation (2) for a given n, we do not need to compute all of the first n + 1 terms in the U sequence.
Let Q = −1, the smallest Lucas pseudoprime to P = 1, 2, 3, ... are
 323, 35, 119, 9, 9, 143, 25, 33, 9, 15, 123, 35, 9, 9, 15, 129, 51, 9, 33, 15, 21, 9, 9, 49, 15, 39, 9, 35, 49, 15, 9, 9, 33, 51, 15, 9, 35, 85, 39, 9, 9, 21, 25, 51, 9, 143, 33, 119, 9, 9, 51, 33, 95, 9, 15, 301, 25, 9, 9, 15, 49, 155, 9, 399, 15, 33, 9, 9, 49, 15, 119, 9, ...
Strong Lucas pseudoprimes
Now, factor into the form where is odd.
A strong Lucas pseudoprime for a given (P, Q) pair is an odd composite number n with GCD(n, D) = 1, satisfying one of the conditions
or
for some 0 ≤ r < s; see page 1396 of.^{[3]} A strong Lucas pseudoprime is also a Lucas pseudoprime (for the same (P, Q) pair), but the converse is not necessarily true. Therefore, the strong test is a more stringent primality test than equation (1).
We can set Q = −1, then and are PFibonacci sequence and PLucas sequence, the pseudoprimes can be called strong Lucas pseudoprime in base P, for example, the least strong Lucas pseudoprime with P = 1, 2, 3, ... are 323, 169, 119, ...
An extra strong Lucas pseudoprime ^{[8]} is a strong Lucas pseudoprime for a set of parameters (P, Q) where Q = 1, satisfying one of the conditions
or
for some . An extra strong Lucas pseudoprime is also a strong Lucas pseudoprime for the same pair.
Implementing a Lucas probable prime test
Before embarking on a probable prime test, one usually verifies that n, the number to be tested for primality, is odd, is not a perfect square, and is not divisible by any small prime less than some convenient limit. Perfect squares are easy to detect using Newton's method for square roots.
We choose a Lucas sequence where the Jacobi symbol , so that δ(n) = n + 1.
Given n, one technique for choosing D is to use trial and error to find the first D in the sequence 5, −7, 9, −11, ... such that is −1. (If D and n have a prime factor in common, then ). With this sequence of D values, the average number of D values that must be tried before we encounter one whose Jacobi symbol is −1 is about 1.79; see,^{[3]} p. 1416. Once we have D, we set P = 1 and . It is a good idea to check that n has no prime factors in common with P or Q. This method of choosing D, P, and Q was suggested by John Selfridge.
(This search will never succeed if n is square, and conversely if it does succeed, that is proof that n is not square. Thus, some time can be saved by delaying testing n for squareness until after the first few search steps have all failed.)
Given D, P, and Q, there are recurrence relations that enable us to quickly compute and in steps; see Lucas sequence § Other relations. First, we can double the subscript from to in one step using the recurrence relations
 .
Next, we can increase the subscript by 1 using the recurrences
 .
(If either of these numerators is odd, we can make it be even by increasing it by n, because all of these calculations are carried out modulo n.) Observe that, for each term that we compute in the U sequence, we compute the corresponding term in the V sequence. As we proceed, we also compute the same, corresponding powers of Q.
We use the bits of the binary expansion of n + 1, starting at the leftmost bit, to determine which terms in the U sequence need to be computed. For example, if n + 1 = 44 (= 101100 in binary), then, using these bits from left to right, we compute U_{1}, U_{2}, U_{4}, U_{5}, U_{10}, U_{11}, U_{22}, and U_{44}. We also compute the samenumbered terms in the V sequence, along with Q^{1}, Q^{2}, Q^{4}, Q^{5}, Q^{10}, Q^{11}, Q^{22}, and Q^{44}.
By the end of the calculation, we will have computed U_{n+1}, V_{n+1}, and Q^{n+1}. We then check congruence (2) using our known value of U_{n+1}.
When D, P, and Q are chosen as described above, the first 10 Lucas pseudoprimes are (see page 1401 of ^{[3]}): 323, 377, 1159, 1829, 3827, 5459, 5777, 9071, 9179, and 10877 (sequence A217120 in the OEIS)
The strong versions of the Lucas test can be implemented in a similar way.
When D, P, and Q are chosen as described above, the first 10 strong Lucas pseudoprimes are: 5459, 5777, 10877, 16109, 18971, 22499, 24569, 25199, 40309, and 58519 (sequence A217255 in the OEIS)
To calculate a list of extra strong Lucas pseudoprimes, set Q = 1. Then try P = 3, 4, 5, 6, ..., until a value of is found so that the Jacobi symbol . With this method for selecting D, P, and Q, the first 10 extra strong Lucas pseudoprimes are 989, 3239, 5777, 10877, 27971, 29681, 30739, 31631, 39059, and 72389 (sequence A217719 in the OEIS)
Checking additional congruence conditions
If we have checked that congruence (2) is true, there are additional congruence conditions we can check that have almost no additional computational cost. If n happens to be composite, these additional conditions may help discover that fact.
If n is an odd prime and , then we have the following (see equation 2 on page 1392 of ^{[3]}):

(3)
Although this congruence condition is not, by definition, part of the Lucas probable prime test, it is almost free to check this condition because, as noted above, the value of V_{n+1} was computed in the process of computing U_{n+1}.
If either congruence (2) or (3) is false, this constitutes a proof that n is not prime. If both of these congruences are true, then it is even more likely that n is prime than if we had checked only congruence (2).
If Selfridge's method (above) for choosing D, P, and Q happened to set Q = −1, then we can adjust P and Q so that D and remain unchanged and P = Q = 5 (see Lucas sequenceAlgebraic relations). If we use this enhanced method for choosing P and Q, then 913 = 11·83 is the only composite less than 10^{8} for which congruence (3) is true (see page 1409 and Table 6 of;^{[3]}).
Here is another congruence condition that is true for primes and that is trivial to check.
Recall that is computed during the calculation of . It would be easy to save the previouslycomputed power of , namely, .
Next, if n is prime, then, by Euler's criterion,
 .
(Here, is the Legendre symbol; if n is prime, this is the same as the Jacobi symbol). Therefore, if n is prime, we must have
 .
The Jacobi symbol on the right side is easy to compute, so this congruence is easy to check. If this congruence does not hold, then n cannot be prime.
Additional congruence conditions that must be satisfied if n is prime are described in Section 6 of.^{[3]} If any of these conditions fails to hold, then we have proved that n is not prime.
Comparison with the Miller–Rabin primality test
k applications of the Miller–Rabin primality test declare a composite n to be probably prime with a probability at most (1/4)^{k}.
There is a similar probability estimate for the strong Lucas probable prime test.^{[9]}
Aside from two trivial exceptions (see below), the fraction of (P,Q) pairs (modulo n) that declare a composite n to be probably prime is at most (4/15).
Therefore, k applications of the strong Lucas test would declare a composite n to be probably prime with a probability at most (4/15)^{k}.
There are two trivial exceptions. One is n = 9. The other is when n = p(p+2) is the product of two twin primes. Such an n is easy to factor, because in this case, n+1 = (p+1)^{2} is a perfect square. One can quickly detect perfect squares using Newton's method for square roots.
By combining a Lucas pseudoprime test with a Fermat primality test, say, to base 2, one can obtain very powerful probabilistic tests for primality, such as the Baillie–PSW primality test.
Fibonacci pseudoprimes
As noted above, when P = 1 and Q = −1, the numbers in the U sequence are the Fibonacci numbers.
A Fibonacci pseudoprime is often (page 264 of,^{[4]} page 142 of,^{[5]} or page 127 of ^{[6]}) defined as a composite number n for which equation (1) above is true with P = 1 and Q = −1 (but n is not divisible by 5). By this definition, the first ten Fibonacci pseudoprimes are 323, 377, 1891, 3827, 4181, 5777, 6601, 6721, 8149, and 10877 (sequence A081264 in the OEIS). The references of Anderson and Jacobsen below use this definition.
If n is congruent to 2 or 3 (mod 5), then Bressoud (,^{[4]} pages 272–273) and Crandall and Pomerance (,^{[5]} page 143 and exercise 3.41 on page 168) point out that it is rare for a Fibonacci pseudoprime to also be a Fermat pseudoprime base 2. However, when n is congruent to 1 or 4 (mod 5), the opposite is true, with over 12% of Fibonacci pseudoprimes under 10^{11} also being base2 Fermat pseudoprimes.
If n is prime and GCD(n, Q) = 1, then (see equation 4 on page 1392 of ^{[3]}) we also have
This leads to an alternate definition of Fibonacci pseudoprime.^{[10]} By this definition, a Fibonacci pseudoprime is a composite number n for which equation (5) is true with P = 1 and Q = −1. Using this definition, the first ten Fibonacci pseudoprimes are 705, 2465, 2737, 3745, 4181, 5777, 6721, 10877, 13201, and 15251 (^{[6]} page 129) (sequence A005845 in the OEIS).
It has been shown that there are no even Fibonacci pseudoprimes with the second definition using equation (5).^{[11]} Using the more common first definition with equation (1), they do exist (sequence A141137 in the OEIS).
A strong Fibonacci pseudoprime may be defined as a composite number for which equation (5) holds for all P.^{[12]} It follows (^{[12]} page 460) that an odd composite integer is a strong Fibonacci pseudoprime if and only if:
 n is also a Carmichael number
 2(p_{i} + 1)  (n − 1) or 2(p_{i} + 1)  (n − p_{i}) for every prime p_{i} dividing n.
The smallest example of a strong Fibonacci pseudoprime is 443372888629441 = 17·31·41·43·89·97·167·331.
Pell pseudoprimes
A Pell pseudoprime may be defined as a composite number n for which equation (1) above is true with P = 2 and Q = −1; the sequence U_{n} then being the Pell sequence. The first pseudoprimes are then 35, 169, 385, 779, 899, 961, 1121, 1189, 2419, ...
This differs from the definition in OEIS: A099011 which may be written as:
with (P, Q) = (2, −1) again defining U_{n} as the Pell sequence. The first pseudoprimes are then 169, 385, 741, 961, 1121, 2001, 3827, 4879, 5719, 6215 ...
A third definition uses equation (5) with (P, Q) = (2, −1), leading to the pseudoprimes 169, 385, 961, 1105, 1121, 3827, 4901, 6265, 6441, 6601, 7107, 7801, 8119, ...
References
 ^ P. S. Bruckman (1994). "Lucas Pseudoprimes are odd". Fibonacci Quarterly. 32: 155–157.
 ^ Sloane, N.J.A. (ed.). "Sequence A005845 (BruckmanLucas pseudoprimes)". The OnLine Encyclopedia of Integer Sequences. OEIS Foundation.
 ^ ^{a} ^{b} ^{c} ^{d} ^{e} ^{f} ^{g} ^{h} ^{i} ^{j} ^{k} ^{l} ^{m} Robert Baillie; Samuel S. Wagstaff, Jr. (October 1980). "Lucas Pseudoprimes" (PDF). Mathematics of Computation. 35 (152): 1391–1417. doi:10.1090/S00255718198005835186. JSTOR 2006406. MR 0583518.
 ^ ^{a} ^{b} ^{c} ^{d} David Bressoud; Stan Wagon (2000). A Course in Computational Number Theory. New York: Key College Publishing in cooperation with Springer. ISBN 9781930190108.
 ^ ^{a} ^{b} ^{c} Richard E. Crandall; Carl Pomerance (2005). Prime numbers: A computational perspective (2nd ed.). SpringerVerlag. ISBN 0387252827.
 ^ ^{a} ^{b} ^{c} Paulo Ribenboim (1996). The New Book of Prime Number Records. SpringerVerlag. ISBN 0387944575.
 ^ Carl Pomerance; John L. Selfridge; Samuel S. Wagstaff, Jr. (July 1980). "The pseudoprimes to 25·10^{9}" (PDF). Mathematics of Computation. 35 (151): 1003–1026. doi:10.1090/S00255718198005728727. JSTOR 2006210.
 ^ Jon Grantham (March 2000). "Frobenius Pseudoprimes". Mathematics of Computation. 70 (234): 873–891. doi:10.1090/S0025571800011972. MR 1680879.
 ^ F. Arnault (April 1997). "The RabinMonier Theorem for Lucas Pseudoprimes". Mathematics of Computation. 66 (218): 869–881. doi:10.1090/s0025571897008363.
 ^ Di Porto, Adina; Filipponi, Piero; Montolivo, Emilio (1990). "On the generalized Fibonacci pseudoprimes". Fibonacci Quarterly. 28: 347–354. CiteSeerX 10.1.1.388.4993.
 ^ Di Porto, Adina (1993). "Nonexistence of Even Fibonacci Pseudoprimes of the First Kind". Fibonacci Quarterly. 31: 173–177. CiteSeerX 10.1.1.376.2601.
 ^ ^{a} ^{b} Müller, Winfried B.; Oswald,Alan (1993). "Generalized Fibonacci Pseudoprimes and Probable Primes". In G.E. Bergum; et al. Applications of Fibonacci Numbers. 5. Kluwer. pp. 459–464. doi:10.1007/9789401120586_45.
External links
 Anderson, Peter G. Fibonacci Pseudoprimes, their factors, and their entry points.
 Anderson, Peter G. Fibonacci Pseudoprimes under 2,217,967,487 and their factors.
 Jacobsen, Dana Pseudoprime Statistics, Tables, and Data (data for Lucas, Strong Lucas, AES Lucas, ES Lucas pseudoprimes below 10^{14}; Fibonacci and Pell pseudoprimes below 10^{12})
 Weisstein, Eric W. "Fibonacci Pseudoprime". MathWorld.
 Weisstein, Eric W. "Lucas Pseudoprime". MathWorld.
 Weisstein, Eric W. "Strong Lucas Pseudoprime". MathWorld.
 Weisstein, Eric W. "Extra Strong Lucas Pseudoprime". MathWorld.