OFFSET
1,1
PROG
(PARI) for(n=1, 10^8, if(sigma(n^2)%(sigma(n)+n)==0, print1(n ", "))) \\ Donovan Johnson, Nov 30 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Nov 28 2013
EXTENSIONS
a(23)-a(25) from Donovan Johnson, Nov 30 2013
STATUS
approved