cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A086250 Smallest base-2 Fermat pseudoprime x that has ord(2,x) = n, or 0 if one does not exist.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 341, 2047, 0, 0, 5461, 4681, 4369, 0, 1387, 0, 13981, 42799, 15709, 8388607, 1105, 1082401, 22369621, 0, 645, 256999, 10261, 0, 16843009, 1227133513, 5726623061, 8727391, 1729, 137438953471, 91625968981, 647089, 561
Offset: 1

Views

Author

T. D. Noe, Jul 14 2003

Keywords

Comments

A base-2 Fermat pseudoprime is a composite number x such that 2^x == 2 (mod x). For such an x, ord(2,x) is the smallest positive integer m such that 2^m == 1 (mod x). For a number x to have order n, it must be a factor of 2^n-1 and not be a factor of 2^r-1 for rA086249 lists the number of pseudoprimes of order n.

Examples

			a(10) = 1 there is only 1 pseudoprime, 341 = 11*31, having order 10; that is, 2^10 = 1 mod 341.
		

Crossrefs

Cf. A001567 (base-2 pseudoprimes), A086249.

Programs

  • Mathematica
    Table[d=Divisors[2^n-1]; num=0; i=1; done=False; While[m=d[[i]]; done=!PrimeQ[m]&&PowerMod[2, m, m]==2&&MultiplicativeOrder[2, m]==n; If[done, num=m]; !done&&i
    				
  • PARI
    { a(n) = fordiv(2^n-1,d, if(d>1 && (d-1)%n==0 && !ispseudoprime(d) && znorder(Mod(2,d))==n,return(d)) ); 0 } /* Max Alekseyev, Jan 07 2015 */

A086256 Number of base-2 Fermat pseudoprimes that divide 2^n-1.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 2, 1, 4, 1, 2, 1, 1, 0, 13, 4, 5, 0, 2, 2, 1, 1, 13, 1, 1, 4, 7, 1, 11, 4, 14, 9, 4, 4, 28, 0, 12, 11, 12, 4, 2, 5, 28, 4, 26, 1, 63, 0, 1, 5, 12, 1, 29, 1, 12, 2, 44, 4, 101, 4, 11, 27, 12, 1, 26, 4, 15, 4, 11, 1, 75, 1, 11, 14, 36, 0, 40, 11
Offset: 1

Views

Author

T. D. Noe, Jul 14 2003

Keywords

Comments

A base-2 Fermat pseudoprime is a composite number x such that 2^x = 2 mod x.

Crossrefs

Cf. A001567 (base-2 pseudoprimes), A046801, A086249.

Programs

  • Mathematica
    Table[d=Divisors[2^n-1]; cnt=0; Do[m=d[[i]]; If[ !PrimeQ[m]&&PowerMod[2, m, m]==2, cnt++ ], {i, Length[d]}]; cnt, {n, 100}]

Formula

a(n) = Sum{d|n} A086249(d), the Mobius transform of A086249.

A298365 Numbers k such that there exists at least one odd pseudoprime of order k.

Original entry on oeis.org

10, 11, 14, 15, 16, 18, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 95, 96, 97
Offset: 1

Views

Author

Krzysztof Ziemak, Jan 17 2018

Keywords

Comments

A composite divisor d of M(m) := 2^m - 1 is called primitive if M(k) != 0 for any k < m.
A primitive composite divisor d of M(m) is said to have rank m, and we write rank(d)=m.
Let M(m)=2^m-1, and define D to be the set of all numbers d such that d|M(m), d==1 (mod m), and rank(d)=m. Then each element d from D is an odd pseudoprime, because if m|d-1, then M(m)|M(d-1) and thus d|M(d-1). The set D contains all composite and primitive divisors d|M(m) that have rank(d)=m and each odd pseudoprime d with rank(d)=m generates only one class [a(n)] with all pseudoprimes d, where a(n)=m, if a(n) is defined as below. See attached file with examples of pseudoprimes.

Examples

			10 is a term since 341 is an odd pseudoprime whose order is 10.
		

Crossrefs

Formula

a(n) = min{k: k>a(n-1) and M(k) has a composite divisor d and rank(d)=k and d==1 (mod k)} for n = 1,2,3,... where M(k):=2^k-1.
Showing 1-3 of 3 results.