OFFSET
1,1
COMMENTS
From Robert Israel, Mar 14 2017: (Start)
If x and y are coprime and a(x)=a(y)=k, then a(xy)=k as well.
If n > 1 is squarefree, then a(n^k) = k+2 for all k>=1.
Is there any n > 1 with a(n) = 2? (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
f:= proc(n) local k, s; uses numtheory;
s:= sigma(n);
for k from 2 do if sigma(n^k) mod s = 0 then return k fi
od
end proc:
map(f, [$1..200]); # Robert Israel, Mar 14 2017
MATHEMATICA
a[n_] := For[k = 2, True, k++, If[Divisible[DivisorSigma[1, n^k], DivisorSigma[1, n]], Return[k]]];
Array[a, 100] (* Jean-François Alcover, Aug 24 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Dec 01 2002
STATUS
approved