OFFSET
0,1
PROG
(PARI) a(n)=if(n<1, 31*(n==0), 8*(sigma(n, 9)-if(n%4==0, 1024*sigma(n/4, 9)))) /* Michael Somos, Jul 16 2004 */
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved