|
%I #16 Aug 02 2023 05:21:25
%S 0,1,2,2,2,4,4,4,4,4,6,6,6,8,8,8,8,10,10,10,10,10,10,10,10,12,14,14,
%T 14,16,16,16,16,16,18,18,18,20,20,20,20,22,22,22,22,22,22,22,22,22,24,
%U 24,24,26,26,26,26,26,28,28,28,30,30,30,30,34,34,34,34,34,34
%N a(n) = A304651(n)/4.
%H Seiichi Manyama, <a href="/A305117/b305117.txt">Table of n, a(n) for n = 0..10000</a>
%F a(n) = a(n-1) + A000089(n) for n > 0.
%t a89[n_] := a89[n] = Product[{p, e} = pe; Which[p < 3 && e == 1, 1, p == 2 && e > 1, 0, Mod[p, 4] == 1, 2, Mod[p, 4] == 3, 0, True, a89[p^e]], {pe, FactorInteger[n]}];
%t a[n_] := a[n] = If[n == 0, 0, a[n-1] + a89[n]];
%t Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Aug 02 2023 *)
%Y Cf. A000089, A176562, A304651.
%K nonn
%O 0,3
%A _Seiichi Manyama_, May 26 2018
|
