cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A349453 Dirichlet inverse of A133494, 3^(n-1).

This page as a plain text file.
%I A349453 #16 Feb 24 2022 09:03:31
%S A349453 1,-3,-9,-18,-81,-189,-729,-2052,-6480,-19197,-59049,-175446,-531441,
%T A349453 -1589949,-4781511,-14335704,-43046721,-129097152,-387420489,
%U A349453 -1162141182,-3486771279,-10459998909,-31381059609,-94142073420,-282429529920,-847285420797,-2541865710960,-7625587899366,-22876792454961,-68630348286531
%N A349453 Dirichlet inverse of A133494, 3^(n-1).
%H A349453 Antti Karttunen, <a href="/A349453/b349453.txt">Table of n, a(n) for n = 1..1001</a>
%F A349453 a(1) = 1; a(n) = -Sum_{d|n, d < n} A133494(n/d) * a(d).
%F A349453 G.f. A(x) satisfies: A(x) = x - Sum_{k>=2} 3^(k-1) * A(x^k). - _Ilya Gutkovskiy_, Feb 23 2022
%t A349453 a[1] = 1; a[n_] := a[n] = -DivisorSum[n, a[#] * 3^(n/# - 1) &, # < n &]; Array[a, 30] (* _Amiram Eldar_, Nov 22 2021 *)
%o A349453 (PARI)
%o A349453 A133494(n) = max(1, 3^(n-1));
%o A349453 memoA349453 = Map();
%o A349453 A349453(n) = if(1==n,1,my(v); if(mapisdefined(memoA349453,n,&v), v, v = -sumdiv(n,d,if(d<n,A133494(n/d)*A349453(d),0)); mapput(memoA349453,n,v); (v)));
%Y A349453 Cf. A133494.
%Y A349453 Cf. also A349449, A349450, A349451, A349452, A349568.
%K A349453 sign
%O A349453 1,2
%A A349453 _Antti Karttunen_, Nov 22 2021