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.

A378524 Dirichlet inverse of A333794.

This page as a plain text file.
%I A378524 #6 Dec 01 2024 10:05:18
%S A378524 1,-3,-6,2,-12,23,-20,0,14,47,-36,-27,-40,79,102,0,-48,-81,-64,-57,
%T A378524 174,143,-96,10,68,159,-24,-97,-112,-517,-116,0,314,191,360,170,-128,
%U A378524 255,350,22,-144,-885,-176,-177,-400,383,-240,0,218,-393,414,-197,-216,211,656,38,566,447,-284,947,-232,463,-696,0,724
%N A378524 Dirichlet inverse of A333794.
%H A378524 Antti Karttunen, <a href="/A378524/b378524.txt">Table of n, a(n) for n = 1..20000</a>
%F A378524 a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, d<n} A333794(n/d) * a(d).
%o A378524 (PARI)
%o A378524 A052126(n) = if(1==n,n,(n/vecmax(factor(n)[, 1])));
%o A378524 A333794(n) = if(1==n,n,n + A333794(n-A052126(n)));
%o A378524 memoA378524 = Map();
%o A378524 A378524(n) = if(1==n,1,my(v); if(mapisdefined(memoA378524,n,&v), v, v = -sumdiv(n,d,if(d<n,A333794(n/d)*A378524(d),0)); mapput(memoA378524,n,v); (v)));
%Y A378524 Cf. A052126, A333794.
%Y A378524 Cf. also A378523.
%K A378524 sign
%O A378524 1,2
%A A378524 _Antti Karttunen_, Dec 01 2024