OFFSET
1,1
COMMENTS
LINKS
R. J. Mathar, Michael De Vlieger, Table of n, a(n) for n = 1..10778 (first 443 terms from R. J. Mathar)
MAPLE
isA107924 := proc(n)
if type(n, 'even') then
dvs := numtheory[divisors](n^2) ;
add(d, d=dvs)/nops(dvs) ;
if type(%, 'integer') then
true;
else
false;
end if;
else
false;
end if;
end proc:
n := 1 :
for k from 2 to 100000 do
if isA107924(k) then
printf("%d %d\n", n, k) ;
n := n+1 ;
end if;
end do: # create b-file, R. J. Mathar, Jul 28 2014
MATHEMATICA
Select[Range[2, 10000, 2], Mod[DivisorSigma[1, #^2], DivisorSigma[0, #^2]]==0&]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Zak Seidov, Jun 10 2005
STATUS
approved