OFFSET
1,2
COMMENTS
a(13) > 10^10. - Hiroaki Yamanouchi, Mar 18 2015
PROG
(PARI) spd(n) = if( n==0, 0, sigma(n) - n); \\ A001065
sad(n) = my(k); if(n>1, k=valuation(n, 2); sigma(2*n+1)+sigma(2*n-1)+sigma(n/2^k)*2^(k+1)-6*n-2, 0); \\ A066417
isok(k) = sad(k) == spd(k); \\ Michel Marcus, Mar 30 2025
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jason Earls, Sep 06 2002
EXTENSIONS
a(12) from Hiroaki Yamanouchi, Mar 18 2015
STATUS
approved
