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 A346473 #15 Nov 02 2021 10:29:12 %S A346473 0,0,1,2,1,3,3,6,12,9,1,5,3,15,11,14,1,12,3,15,23,27,5,9,18,33,25,25, %T A346473 1,15,5,30,37,45,29,14,3,51,39,27,1,27,3,45,37,63,5,17,64,54,53,55,5, %U A346473 39,19,45,65,81,1,9,5,87,51,62,35,51,3,75,79,63,1,18,5,105,61,85,47,63,3,51,84,117,5,25,25,123,95 %N A346473 a(n) = A250469(n) - sigma(n). %C A346473 The first negative term is a(120) = -3. %H A346473 Antti Karttunen, <a href="/A346473/b346473.txt">Table of n, a(n) for n = 1..10000</a> %H A346473 <a href="/index/Si#sieve">Index entries for sequences generated by sieves</a> %H A346473 <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a> %F A346473 a(n) = A250469(n) - A000203(n). %F A346473 a(n) = A286385(n) - A280692(n). %F A346473 a(A001359(n)) = 1 for all n >= 1. %t A346473 Block[{g}, g[n_] := If[n == 1, 0, PrimePi@FactorInteger[n][[1, 1]]]; Function[s, MapIndexed[Lookup[s, g[First@ #2] + 1][[#1]] - Boole[First@ #2 == 1] - DivisorSigma[1, First@ #2] &, #] &@ Map[Position[Lookup[s, g@ #], #][[1, 1]] &, Range@ 120]]@ PositionIndex@ Array[g, 10^4]] (* _Michael De Vlieger_, Oct 18 2021 *) %o A346473 (PARI) A346473(n) = A250469(n)-sigma(n); %Y A346473 Cf. A000203, A001359, A250469, A280692, A286385. %K A346473 sign,look %O A346473 1,4 %A A346473 _Antti Karttunen_, Jul 28 2021