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.

A346473 a(n) = A250469(n) - sigma(n).

This page as a plain text file.
%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