OFFSET
1,2
COMMENTS
From Amiram Eldar, Nov 12 2024: (Start)
Equivalently, numbers k that divide A377484(k) = Product_{d|k, d>1} (d - 1).
After the first term a(1) = 1, the next odd term is a(5) = 105, the next term that is coprime to 6 is a(228) = 6545, and the next term that is coprime to 30 is a(574) = 19019. (End)
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
30 is a term because 30^2 divides A056819(30) = 5320224000.
MATHEMATICA
Select[Range[1000], Divisible[Times @@ (Rest@ Divisors[#] - 1), #] &] (* Amiram Eldar, Nov 12 2024 *)
PROG
(PARI) is(k) = if(k == 1, 1, my(d = divisors(k)); !(prod(i = 2, #d, d[i]-1) % k)); \\ Amiram Eldar, Nov 12 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Leroy Quet, Sep 06 2000
STATUS
approved