OFFSET
1,2
MATHEMATICA
f[p_, e_] := p^e - 1; uphi[n_] := Times @@ f @@@ FactorInteger[n]; uphi[1] = 1; q[n_] := Module[{m = 1, u = uphi[n]}, While[m <= u && ! Divisible[m*(m + 1), n], m++]; m == u]; Select[Range[160], q] (* Amiram Eldar, Jun 13 2022 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 13 2022
STATUS
approved