login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) is the denominator of A346009(n)/A346010(n) - A001221(n)/2.
3

%I #8 Jul 02 2021 03:48:52

%S 0,1,1,6,1,1,1,4,6,1,1,6,1,1,1,10,1,6,1,6,1,1,1,4,6,1,4,6,1,1,1,3,1,1,

%T 1,3,1,1,1,4,1,1,1,6,6,1,1,10,6,6,1,6,1,4,1,4,1,1,1,6,1,1,6,14,1,1,1,

%U 6,1,1,1,12,1,1,6,6,1,1,1,10,10,1,1,6,1

%N a(n) is the denominator of A346009(n)/A346010(n) - A001221(n)/2.

%C See A346012 for details.

%H Amiram Eldar, <a href="/A346013/b346013.txt">Table of n, a(n) for n = 1..10000</a>

%t f[p_, e_] := e/(e + 1); a[1] = 0; a[n_] := Denominator[Plus @@ f @@@ (fct = FactorInteger[n]) - Length[fct]/2]; Array[a, 100]

%Y Cf. A001221, A346009, A346010, A346012 (numerators).

%K nonn,frac

%O 1,4

%A _Amiram Eldar_, Jul 01 2021