OFFSET
1,1
COMMENTS
Up to 2^64 all the 1085 terms are nonsquarefree, 2 terms have 4 prime factors: a(163) = 18362297383286473 = 3037 * 6073 * 9109 * 109297 and a(651) = 2587580959818925201 = 18121 * 36241 * 54361 * 72481, and no term have more than 4 prime factors.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1085 (terms below 2^64)
EXAMPLE
11500521553 is a term since it is a Fermat pseudoprime to both bases 2 and 3, and its proper divisors that are larger than 1 are either primes (937, 1873, 6553) or Fermat pseudoprimes to both bases 2 and 3 (1755001, 6140161, 12273769, 11500521553).
MATHEMATICA
pspQ[n_] := PrimeOmega[n] > 2 && AllTrue[Rest @ Divisors[n], PowerMod[2, # - 1, #] == 1 && PowerMod[3, # - 1, #] == 1 &]; seq = {}; Do[If[pspQ[n], AppendTo[seq, n]], {n, 1, 6*10^10}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 08 2020
STATUS
approved