OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Antti Karttunen, Scheme-program for computing this sequence.
FORMULA
a(n) = Sum_{d|n} A318622(d). - Robert Israel, Aug 30 2018
MAPLE
h:= proc(n) option remember; numtheory:-phi(n)/numtheory:-order(2, n/2^padic:-ordp(n, 2)) end proc:
f:= n -> add(h(d), d=numtheory:-divisors(n)):
map(f, [$1..100]); # Robert Israel, Aug 30 2018
MATHEMATICA
h[n_] := EulerPhi[n]/MultiplicativeOrder[2, n/2^IntegerExponent[n, 2]];
a[n_] := DivisorSum[n, h];
Array[a, 100] (* Jean-François Alcover, Aug 19 2022, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 03 2004
STATUS
approved