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.

A355829 Dirichlet inverse of A009194, the greatest common divisor of sigma(n) and n, where sigma is the sum of divisors function.

This page as a plain text file.
%I A355829 #9 Jul 20 2022 18:43:34
%S A355829 1,-1,-1,0,-1,-4,-1,0,0,0,-1,7,-1,0,-1,0,-1,8,-1,1,1,0,-1,-10,0,0,0,
%T A355829 -25,-1,10,-1,0,-1,0,1,15,-1,0,1,-8,-1,6,-1,-1,2,0,-1,16,0,2,-1,1,-1,
%U A355829 -6,1,46,1,0,-1,-9,-1,0,0,0,1,10,-1,1,-1,2,-1,-29,-1,0,4,-1,1,6,-1,16,0,0,-1,29,1,0,-1,2,-1,-8
%N A355829 Dirichlet inverse of A009194, the greatest common divisor of sigma(n) and n, where sigma is the sum of divisors function.
%H A355829 Antti Karttunen, <a href="/A355829/b355829.txt">Table of n, a(n) for n = 1..16384</a>
%H A355829 <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%F A355829 a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, d<n} A009194(n/d) * a(d).
%t A355829 s[n_] := GCD[n, DivisorSigma[1, n]]; a[1] = 1; a[n_] := - DivisorSum[n, a[#] * s[n/#] &, # < n &]; Array[a, 100] (* _Amiram Eldar_, Jul 20 2022 *)
%o A355829 (PARI)
%o A355829 A009194(n) = gcd(n, sigma(n));
%o A355829 memoA355829 = Map();
%o A355829 A355829(n) = if(1==n,1,my(v); if(mapisdefined(memoA355829,n,&v), v, v = -sumdiv(n,d,if(d<n,A009194(n/d)*A355829(d),0)); mapput(memoA355829,n,v); (v)));
%Y A355829 Cf. A000203, A009194.
%Y A355829 Cf. also A355828.
%K A355829 sign
%O A355829 1,6
%A A355829 _Antti Karttunen_, Jul 20 2022