OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
lpf[n_]:=FactorInteger[n][[1, 1]]; lpf /@ Select[Range[1, 10^4, 2], CompositeQ[#] && PowerMod[2, #-1, # ] == 1 &] (* Amiram Eldar, Sep 08 2019 *)
PROG
(PARI) terms(n) = my(i=0); forcomposite(c=1, , if(i>=n, break, if(Mod(2, c)^(c-1)==1, print1(factor(c)[1, 1], ", "); i++)))
/* Print initial 80 terms as follows */
terms(80) \\ Felix Fröhlich, Sep 08 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Dec 10 2010
EXTENSIONS
Corrected by D. S. McNeil, Dec 11 2010
STATUS
approved