login
A349452
Dirichlet inverse of A011782, 2^(n-1).
9
1, -2, -4, -4, -16, -16, -64, -104, -240, -448, -1024, -1904, -4096, -7936, -16256, -32272, -65536, -129888, -262144, -522176, -1048064, -2093056, -4194304, -8379520, -16776960, -33538048, -67106880, -134184704, -268435456, -536801024, -1073741824, -2147352224, -4294959104, -8589672448, -17179867136, -34359197184
OFFSET
1,2
LINKS
FORMULA
a(1) = 1; a(n) = -Sum_{d|n, d < n} A011782(n/d) * a(d).
G.f. A(x) satisfies: A(x) = x - Sum_{k>=2} 2^(k-1) * A(x^k). - Ilya Gutkovskiy, Feb 23 2022
MATHEMATICA
a[1] = 1; a[n_] := a[n] = -DivisorSum[n, a[#] * 2^(n/# - 1) &, # < n &]; Array[a, 36] (* Amiram Eldar, Nov 22 2021 *)
PROG
(PARI)
A011782(n) = (2^(n-1));
memoA349452 = Map();
A349452(n) = if(1==n, 1, my(v); if(mapisdefined(memoA349452, n, &v), v, v = -sumdiv(n, d, if(d<n, A011782(n/d)*A349452(d), 0)); mapput(memoA349452, n, v); (v)));
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Nov 22 2021
STATUS
approved