OFFSET
1,1
COMMENTS
A bisection of A341710.
MAPLE
with(numtheory):
b:= proc(n) option remember; nops(invphi(n)) end:
g:= proc(n) option remember; `if`(n=0, 1, add(
g(n-j)*add(d*b(d), d=divisors(j)), j=1..n)/n)
end:
a:= n-> g(2*n)/2:
seq(a(n), n=1..40); # Alois P. Heinz, Feb 19 2021
MATHEMATICA
terms = 64; (* number of terms of A120963 *)
nmax = Floor[terms/2];
S[m_] := S[m] = CoefficientList[Product[1/(1 - x^EulerPhi[k]),
{k, 1, m*terms}] + O[x]^(terms+1), x];
S[m = 1];
S[m++];
While[S[m] != S[m-1], m++];
A120963 = S[m];
a[n_ /; 1 <= n <= nmax] := A120963[[2n+1]]/2;
Table[a[n], {n, 1, nmax}] (* Jean-François Alcover, May 12 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Feb 19 2021
STATUS
approved