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.

A107473 Sum of numerator and denominator of Product_{p|n, p prime} (1 - 1/p).

This page as a plain text file.
%I A107473 #12 Aug 04 2025 11:47:22
%S A107473 2,3,5,3,9,4,13,3,5,7,21,4,25,10,23,3,33,4,37,7,11,16,45,4,9,19,5,10,
%T A107473 57,19,61,3,53,25,59,4,73,28,21,7,81,9,85,16,23,34,93,4,13,7,83,19,
%U A107473 105,4,19,10,31,43,117,19,121,46,11,3,113,43,133,25,113,47,141,4,145,55,23,28
%N A107473 Sum of numerator and denominator of Product_{p|n, p prime} (1 - 1/p).
%F A107473 a(n) = A076511(n) + A076512(n). - _Michel Marcus_, Aug 16 2019
%e A107473 a(12) = 4 = 1+3 because (1 - 1/2)(1 - 1/3) = 1/3.
%p A107473 a:=proc(n) local b,ct,f; with(numtheory): b:=convert(factorset(n),list): ct:=nops(b): f:=simplify(product(1-1/b[j],j=1..ct)):numer(f)+denom(f) end: seq(a(n),n=1..100); # _Emeric Deutsch_, May 28 2005
%Y A107473 Cf. A076511, A076512.
%K A107473 nonn
%O A107473 1,1
%A A107473 _Leroy Quet_, May 27 2005
%E A107473 More terms from _Emeric Deutsch_, May 28 2005