OFFSET
1,1
COMMENTS
MATHEMATICA
f1[p_, e_] := (p + 1)^e; f2[p_, e_] := p^e + 1; a1[1] = 1; a1[n_] := GCD[Times @@ f1 @@@ (f = FactorInteger[n]), Times @@ f2 @@@ f]; f3[p_, e_] := a1[p^e]; a2[n_] := Times @@ f3 @@@ FactorInteger[n]; Position[Table[a2[n] - a1[n], {n, 1, 42000}], _?(# != 0 &)] // Flatten (* Amiram Eldar, Nov 05 2021 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 05 2021
STATUS
approved