Frobenius pseudoprime
   HOME

TheInfoList



OR:

In
number theory Number theory (or arithmetic or higher arithmetic in older usage) is a branch of pure mathematics devoted primarily to the study of the integers and integer-valued functions. German mathematician Carl Friedrich Gauss (1777–1855) said, "Ma ...
, a Frobenius pseudoprime is a
pseudoprime A pseudoprime is a probable prime (an integer that shares a property common to all prime numbers) that is not actually prime. Pseudoprimes are classified according to which property of primes they satisfy. Some sources use the term pseudoprime to ...
, whose definition was inspired by the quadratic Frobenius test described by Jon Grantham in a 1998 preprint and published in 2000. Frobenius pseudoprimes can be defined with respect to
polynomial In mathematics, a polynomial is an expression consisting of indeterminates (also called variables) and coefficients, that involves only the operations of addition, subtraction, multiplication, and positive-integer powers of variables. An exampl ...
s of degree at least 2, but they have been most extensively studied in the case of
quadratic polynomial In mathematics, a quadratic polynomial is a polynomial of degree two in one or more variables. A quadratic function is the polynomial function defined by a quadratic polynomial. Before 20th century, the distinction was unclear between a polynomial ...
s.


Frobenius pseudoprimes w.r.t. quadratic polynomials

Definition of Frobenius pseudoprimes with respect to a monic quadratic polynomial x^2 - Px + Q, where the
discriminant In mathematics, the discriminant of a polynomial is a quantity that depends on the coefficients and allows deducing some properties of the roots without computing them. More precisely, it is a polynomial function of the coefficients of the orig ...
D = P^2-4Q is not a square, can be expressed in terms of
Lucas sequence In mathematics, the Lucas sequences U_n(P,Q) and V_n(P, Q) are certain constant-recursive integer sequences that satisfy the recurrence relation : x_n = P \cdot x_ - Q \cdot x_ where P and Q are fixed integers. Any sequence satisfying this r ...
s U_n(P,Q) and V_n(P,Q) as follows. A
composite number A composite number is a positive integer that can be formed by multiplying two smaller positive integers. Equivalently, it is a positive integer that has at least one divisor In mathematics, a divisor of an integer n, also called a factor ...
''n'' is a Frobenius (P,Q) pseudoprime if and only if : (1) \qquad \gcd(n,2QD)=1, : (2) \qquad U_(P,Q) \equiv 0 \pmod n, and : (3) \qquad V_(P,Q) \equiv 2Q^ \pmod, where \delta=\left(\tfrac Dn\right) is the
Jacobi symbol Jacobi symbol for various ''k'' (along top) and ''n'' (along left side). Only are shown, since due to rule (2) below any other ''k'' can be reduced modulo ''n''. Quadratic residues are highlighted in yellow — note that no entry with a ...
. When condition (2) is satisfied, condition (3) becomes equivalent to : (3') \qquad V_n(P,Q) \equiv P\pmod. Therefore, Frobenius (P,Q) pseudoprime can be equivalently defined by conditions (1-2) and (3), or by conditions (1-2) and (3′). Since conditions (2) and (3) hold for all primes which satisfy the simple condition (1), they can be used as a probable prime test. (If condition (1) fails, either the
greatest common divisor In mathematics, the greatest common divisor (GCD) of two or more integers, which are not all zero, is the largest positive integer that divides each of the integers. For two integers ''x'', ''y'', the greatest common divisor of ''x'' and ''y'' is ...
is less than , in which case it is a non-trivial factor and is composite, or the GCD equals , in which case you should try different parameters and which are not multiples of .)


Relations to other pseudoprimes

Every Frobenius (P,Q) pseudoprime is also * a
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. Baillie-Wagstaff-Lucas pseudoprimes Baill ...
with parameters (P,Q), since it is defined by conditions (1) and (2); * a Dickson pseudoprime with parameters (P,Q), since it is defined by conditions (1) and (3'); * a
Fermat pseudoprime In number theory, the Fermat pseudoprimes make up the most important class of pseudoprimes that come from Fermat's little theorem. Definition Fermat's little theorem states that if ''p'' is prime and ''a'' is coprime to ''p'', then ''a'p''− ...
base , Q, when , Q, >1. Converse of neither of these statements is true, making the Frobenius (P,Q) pseudoprimes a proper subset of each of the sets of Lucas pseudoprimes and Dickson pseudoprimes with parameters (P,Q), and Fermat pseudoprimes base , Q, when , Q, >1. Furthermore, it follows that for the same parameters (P,Q), a composite number is a Frobenius pseudoprime if and only if it is both Lucas and Dickson pseudoprime. In other words, for every fixed pair of parameters (P,Q), the set of Frobenius pseudoprimes equals the intersection of the sets of Lucas and Dickson pseudoprimes. While each Frobenius (P,Q) pseudoprime is a Lucas pseudoprime, it is not necessarily a strong Lucas pseudoprime. For example, 6721 is the first Frobenius pseudoprime for (P,Q) = (1,-1), which is not a strong Lucas pseudoprime. Every Frobenius pseudoprime to x^3-x-1 is also a restricted Perrin pseudoprime. Analogous statements hold for other cubic polynomials of the form x^3-rx^2+sx-1.


Examples

Frobenius pseudoprimes with respect to the Fibonacci polynomial x^2-x-1 are determined in terms of the
Fibonacci numbers In mathematics, the Fibonacci numbers, commonly denoted , form a sequence, the Fibonacci sequence, in which each number is the sum of the two preceding ones. The sequence commonly starts from 0 and 1, although some authors start the sequence from ...
F_n = U_n(1,-1) and
Lucas numbers The Lucas numbers or Lucas series are an integer sequence named after the mathematician François Édouard Anatole Lucas (1842–1891), who studied both that sequence and the closely related Fibonacci numbers. Lucas numbers and Fibonacci nu ...
L_n = V_n(1,-1). Such Frobenius pseudoprimes form the sequence: :4181, 5777, 6721, 10877, 13201, 15251, 34561, 51841, 64079, 64681, 67861, 68251, 75077, 90061, 96049, 97921, 100127, 113573, 118441, 146611, 161027, 162133, 163081, 186961, 197209, 219781, 231703, 252601, 254321, 257761, 268801, 272611, 283361, 302101, 303101, 330929, 399001, 430127, 433621, 438751, 489601, ... . While 323 is the first
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. Baillie-Wagstaff-Lucas pseudoprimes Baill ...
with respect to the Fibonacci polynomial x^2-x-1, the first Frobenius pseudoprime with respect to the same polynomial is 4181 (Grantham stated it as 5777 but multiple authors have noted this is incorrect and is instead the first pseudoprime with \left(\tfrac\right)=-1 for this polynomial). Another case, Frobenius pseudoprimes with respect to the quadratic polynomial x^2-3x-1 can be determined using the Lucas (3,-1) sequence and are: :119, 649, 1189, 4187, 12871, 14041, 16109, 23479, 24769, 28421, 31631, 34997, 38503, 41441, 48577, 50545, 56279, 58081, 59081, 61447, 75077, 91187, 95761, 96139, 116821, 127937, 146329, 148943, 150281, 157693, 170039, 180517, 188501, 207761, 208349, 244649, 281017, 311579, 316409, 349441, 350173, 363091, 371399, 397927, 423721, 440833, 459191, 473801, 479119, 493697, ... In this case, the first Frobenius pseudoprime with respect to the quadratic polynomial x^2-3x-1 is 119, which is also the first Lucas pseudoprime with respect to the same polynomial. Besides, \left(\tfrac\right)=-1. The quadratic polynomial x^2-3x-5, i.e. (P,Q)=(3,-5), has sparser pseudoprimes as compared to many other simple quadratics. Using the same process as above, we get the sequence: : 13333, 44801, 486157, 1615681, 3125281, 4219129, 9006401, 12589081, 13404751, 15576571, 16719781, …. Notice there are only 3 such pseudoprimes below 500000, while there are many Frobenius (1, −1) and (3, −1) pseudoprimes below 500000. Every entry in this sequence is a
Fermat pseudoprime In number theory, the Fermat pseudoprimes make up the most important class of pseudoprimes that come from Fermat's little theorem. Definition Fermat's little theorem states that if ''p'' is prime and ''a'' is coprime to ''p'', then ''a'p''− ...
to base 5 as well as a Lucas (3, −5) pseudoprime, but the converse is not true: 642001 is both a psp-5 and a Lucas (3,-5) pseudoprime, but is not a Frobenius (3, −5) pseudoprime. (note that
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. Baillie-Wagstaff-Lucas pseudoprimes Baill ...
for a pair need not to be a
Fermat pseudoprime In number theory, the Fermat pseudoprimes make up the most important class of pseudoprimes that come from Fermat's little theorem. Definition Fermat's little theorem states that if ''p'' is prime and ''a'' is coprime to ''p'', then ''a'p''− ...
for base , , , e.g. 14209 is a Lucas (1, −3) pseudoprime, but not a Fermat pseudoprime for base 3.


Strong Frobenius pseudoprimes

Strong Frobenius pseudoprimes are also defined. Details on implementation for quadratic polynomials can be found in Crandall and Pomerance.


Pseudoprimality tests

The conditions defining Frobenius pseudoprime can be used for testing a given number ''n'' for probable primality. Often such tests do not rely on fixed parameters (P,Q), but rather select them in a certain way depending on the input number ''n'' in order to decrease the proportion of
false positive A false positive is an error in binary classification in which a test result incorrectly indicates the presence of a condition (such as a disease when the disease is not present), while a false negative is the opposite error, where the test resul ...
s, i.e., composite numbers that pass the test. Sometimes such composite numbers are commonly called Frobenius pseudoprimes, although they may correspond to different parameters. Using parameter selection ideas first laid out in Baillie and Wagstaff (1980) as part of the
Baillie–PSW primality test The Baillie–PSW primality test is a probabilistic primality testing algorithm that determines whether a number is composite or is a probable prime. It is named after Robert Baillie, Carl Pomerance, John Selfridge, and Samuel Wagstaff. The Baill ...
and used by Grantham in his quadratic Frobenius test, one can create even better quadratic tests. In particular, it was shown that choosing parameters from
quadratic non-residue In number theory, an integer ''q'' is called a quadratic residue modulo ''n'' if it is congruent to a perfect square modulo ''n''; i.e., if there exists an integer ''x'' such that: :x^2\equiv q \pmod. Otherwise, ''q'' is called a quadratic ...
s modulo ''n'' (based on the
Jacobi symbol Jacobi symbol for various ''k'' (along top) and ''n'' (along left side). Only are shown, since due to rule (2) below any other ''k'' can be reduced modulo ''n''. Quadratic residues are highlighted in yellow — note that no entry with a ...
) makes far stronger tests, and is one reason for the success of the
Baillie–PSW primality test The Baillie–PSW primality test is a probabilistic primality testing algorithm that determines whether a number is composite or is a probable prime. It is named after Robert Baillie, Carl Pomerance, John Selfridge, and Samuel Wagstaff. The Baill ...
. For instance, for the parameters (''P'',2), where ''P'' is the first odd integer that satisfies \left(\tfrac\right) = -1, there are no pseudoprimes below 2^. Yet another test is proposed by Khashin. For a given non-square number ''n'', it first computes a parameter ''c'' as the smallest odd prime having Jacobi symbol \left(\tfrac\right)=-1, and then verifies the congruence: :(1 + \sqrt)^n \equiv (1 - \sqrt) \pmod n. While all prime ''n'' pass this test, a composite ''n'' passes it if and only if ''n'' is a Frobenius pseudoprime for (P,Q)=(2,1-c). Similar to the above example, Khashin notes that no pseudoprime has been found for his test. He further shows that any that exist under 260 must have a factor less than 19 or have ''c'' > 128.


Properties

The computational cost of the Frobenius pseudoprimality test with respect to quadratic polynomials is roughly three times the cost of a strong pseudoprimality test (e.g. a single round of the
Miller–Rabin primality test The Miller–Rabin primality test or Rabin–Miller primality test is a probabilistic primality test: an algorithm which determines whether a given number is likely to be prime, similar to the Fermat primality test and the Solovay–Strassen prim ...
), 1.5 times that of a Lucas pseudoprimality test, and slightly more than a
Baillie–PSW primality test The Baillie–PSW primality test is a probabilistic primality testing algorithm that determines whether a number is composite or is a probable prime. It is named after Robert Baillie, Carl Pomerance, John Selfridge, and Samuel Wagstaff. The Baill ...
. Note that the quadratic Frobenius test is stronger than the Lucas test. For example, 1763 is a Lucas pseudoprime to (''P'', ''Q'') = (3, -1) since ''U''1764(3,-1) ≡ 0 (mod 1763) (''U''(3,-1) is given in ), and it also passes the Jacobi step since \left(\tfrac\right) = -1, but it fails the Frobenius test to ''x''2 - 3''x'' - 1. This property can be clearly seen when the algorithm is formulated as shown in Crandall and Pomerance Algorithm 3.6.9 or as shown by Loebenberger, as the algorithm does a Lucas test followed by an additional check for the Frobenius condition. While the quadratic Frobenius test does not have formal error bounds beyond that of the Lucas test, it can be used as the basis for methods with much smaller error bounds. Note that these have more steps, additional requirements, and non-negligible additional computation beyond what is described on this page. It is important to note that the error bounds for these methods do not apply to the standard or strong Frobenius tests with fixed values of (P,Q) described on this page. Based on this idea of pseudoprimes, algorithms with strong worst-case error bounds can be built. The quadratic Frobenius test, using a quadratic Frobenius test plus other conditions, has a bound of \tfrac. Müller in 2001 proposed the MQFT test with bounds of essentially \tfrac. Damgård and Frandsen in 2003 proposed the EQFT with a bound of essentially \tfrac. Seysen in 2005 proposed the SQFT test with a bound of \tfrac and a SQFT3 test with a bound of \tfrac. Seysen, Martin
A Simplified Quadratic Frobenius Primality Test
2005.
Given the same computational effort, these offer better worst-case bounds than the commonly used
Miller–Rabin primality test The Miller–Rabin primality test or Rabin–Miller primality test is a probabilistic primality test: an algorithm which determines whether a given number is likely to be prime, similar to the Fermat primality test and the Solovay–Strassen prim ...
.


See also

*
Pseudoprime A pseudoprime is a probable prime (an integer that shares a property common to all prime numbers) that is not actually prime. Pseudoprimes are classified according to which property of primes they satisfy. Some sources use the term pseudoprime to ...
*
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. Baillie-Wagstaff-Lucas pseudoprimes Baill ...
*
Ferdinand Georg Frobenius Ferdinand Georg Frobenius (26 October 1849 – 3 August 1917) was a German mathematician, best known for his contributions to the theory of elliptic functions, differential equations, number theory, and to group theory. He is known for the famou ...
* Quadratic Frobenius test


References


External links

* * * Jacobsen, Dan
Pseudoprime Statistics, Tables, and Data
(data for Frobenius (1,-1) and (3,-5) pseudoprimes) {{Classes of natural numbers Pseudoprimes