The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A321866 Indices of tetrahedral numbers that are Fermat pseudoprimes to base 2. 2
3457, 16705, 21169, 28297, 30577, 45481, 114601, 123121, 127297, 140977, 156601, 159337, 312841, 393121, 418177, 437977, 443017, 453601, 509737, 518017, 521137, 539401, 545161, 545617, 657841, 679297, 704161, 717817, 762121, 775057, 832801, 904801, 996601 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Numbers n such that n(n+1)(n+2)/6 is a Fermat pseudoprimes to base 2.
The corresponding tetrahedral Fermat pseudoprimes are 6891657409, 777080801185, 1581289265305, 3776730328549, 4765143438329, 15680770945781, 250856489370101, 311068284648121, 343806031110049, ...
Sierpinski asked for the existence of these numbers in 1965.
LINKS
Andrzej Rotkiewicz, On some problems of W. Sierpinski, Acta Arithmetica, Vol. 21 (1972), pp. 251-259.
EXAMPLE
3457 is in the sequence since A000292(3457) = 6891657409 is a Fermat pseudoprime to base 2.
MATHEMATICA
fermatQ[n_, k_] := CompositeQ[n] && PowerMod[k, n-1, n]==1; p[n_] := n(n+1)(n+2)/6; seq={}; Do[p1=p[n]; If[fermatQ[p1, 2], AppendTo[seq, n]], {n, 1, 1000000, 2}]; seq
PROG
(PARI) isok(n) = my(t = n*(n+1)*(n+2)/6); (t != 1) && (Mod(2, t)^t == 2); \\ Michel Marcus, Nov 20 2018
CROSSREFS
Sequence in context: A234893 A251485 A082242 * A235835 A235830 A235586
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 20 2018
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 13 16:07 EDT 2024. Contains 372522 sequences. (Running on oeis4.)