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.
%I A354826 #8 Jun 09 2022 11:41:35 %S A354826 1,-2,-2,0,-2,5,-2,0,0,5,-2,-2,-2,5,5,0,-2,-2,-2,-2,5,5,-2,0,0,5,0,-2, %T A354826 -2,-17,-2,0,5,5,5,8,-2,5,5,0,-2,-17,-2,-2,-2,5,-2,0,0,-2,5,-2,-2,0,5, %U A354826 0,5,5,-2,16,-2,5,-2,0,5,-17,-2,-2,5,-17,-2,-8,-2,5,-2,-2,5,-17,-2,0,0,5,-2,16,5,5,5,0,-2,16 %N A354826 Dirichlet inverse of A238745. %H A354826 Antti Karttunen, <a href="/A354826/b354826.txt">Table of n, a(n) for n = 1..12600</a> %H A354826 Antti Karttunen, <a href="/A354826/a354826.txt">Data supplement: n, a(n) computed for n = 1..65537</a> %F A354826 a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, d<n} A238745(n/d) * a(d). %o A354826 (PARI) %o A354826 A124859(n) = { my(f=factor(n)); for(k=1, #f~, f[k, 1] = prod(j=1, f[k, 2], prime(j)); f[k, 2] = 1); factorback(f); }; \\ From A124859 %o A354826 A181819(n) = factorback(apply(e->prime(e),(factor(n)[,2]))); %o A354826 A238745(n) = A181819(A124859(n)); %o A354826 memoA354826 = Map(); %o A354826 A354826(n) = if(1==n,1,my(v); if(mapisdefined(memoA354826,n,&v), v, v = -sumdiv(n,d,if(d<n,A238745(n/d)*A354826(d),0)); mapput(memoA354826,n,v); (v))); %Y A354826 Cf. A238745. %Y A354826 Cf. also A354349, A354359. %K A354826 sign %O A354826 1,2 %A A354826 _Antti Karttunen_, Jun 09 2022