OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
PROG
(PARI) \\ here b(n) is A023136.
b(n)=sumdiv(n>>valuation(n, 2), d, eulerphi(d)/znorder(Mod(4, d)))
a(n)=sumdiv(n, d, b(d)) \\ Andrew Howroyd, Sep 22 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Gary W. Adamson, Sep 21 2007
EXTENSIONS
More terms from R. J. Mathar, Jan 19 2009
STATUS
approved
