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!)
A327789 a(n) is the smallest Fermat pseudoprime to base 2 such that gpf(p-1) = prime(n) for all prime factors p of a(n). 1
4369, 1387, 341, 3277, 2047, 8321, 31621, 104653, 280601, 13747, 2081713, 88357, 8902741, 741751, 665333, 680627, 2008597, 1252697, 3235699, 1293337, 513629, 8095447, 83333, 2284453, 604117, 191981609, 1530787, 13747361, 3568661, 769757, 6973063, 275887, 12854437, 16705021 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Equivalently, a(n) is the smallest composite number k such that 2^(k-1) == 1 (mod k) and gpf(p-1) = prime(n) for all prime factors p of k.

LINKS

Daniel Suteu, Table of n, a(n) for n = 1..10000

Eric Weisstein's World of Mathematics, Poulet Number

EXAMPLE

a(1) = 4369 = (2*2*2*2 + 1)(2*2*2*2*2*2*2*2 + 1).

a(2) = 1387 = (2*3*3 + 1)(2*2*2*3*3 + 1).

a(3) = 341 = (2*5 + 1)(2*3*5 + 1).

a(4) = 3277 = (2*2*7 + 1)(2*2*2*2*7 + 1).

a(5) = 2047 = (2*11 + 1)(2*2*2*11 + 1).

MATHEMATICA

pspQ[n_] := CompositeQ[n] && PowerMod[2, (n - 1), n] == 1; gpf[n_] := FactorInteger[n][[-1, 1]]; g[n_] := If[Length[(u = Union[gpf /@ (FactorInteger[n][[;; , 1]] - 1)])] == 1, u[[1]], 1]; m = 10; c = 0; k = 0; v = Table[0, {m}]; While[c < m, k++ If[! pspQ[k], Continue[]]; If[(p = g[k]) > 1, i = PrimePi[p]; If[i <= m && v[[i]] == 0, c++; v[[i]] = k]]]; v (* Amiram Eldar, Oct 08 2019 *)

PROG

(Perl) use ntheory ":all"; sub a { my $p = nth_prime(shift); for(my $k = 4; ; ++$k) { return $k if (is_pseudoprime($k, 2) and !is_prime($k) and vecall { (factor($_-1))[-1] == $p } factor($k)) } }

for my $n (1..25) { print "a($n) = ", a($n), "\n" }

CROSSREFS

Cf. A001567 (Fermat pseudoprimes to base 2), A006530 (gpf).

Sequence in context: A347603 A217345 A140937 * A043598 A043843 A043857

Adjacent sequences: A327786 A327787 A327788 * A327790 A327791 A327792

KEYWORD

nonn

AUTHOR

Daniel Suteu, Sep 25 2019

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 January 31 22:24 EST 2023. Contains 359981 sequences. (Running on oeis4.)