OFFSET
1,1
COMMENTS
EXAMPLE
11904 = 93*(93+31+3+1).
PROG
(PARI)
gdivodd(n)={m=n; while(m/2==m\2, m=m/2); return(m)}
{for (n=2, 2*10^8, m=gdivodd(n)*sumdiv(n, d, d*(d%2)); if(m==n, print(n)))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Antonio Roldán, May 19 2013
STATUS
approved