cp's OEIS Frontend

A341711 a(n) = A120963(2*n+1)/2.

%I A341711 #19 Mar 03 2023 09:31:36
%S A341711 1,5,19,59,165,419,1001,2257,4877,10133,20399,39881,76085,141877,
%T A341711 259373,465493,821813,1428725,2449573,4145249,6931259,11459483,
%U A341711 18749007,30373189,48752125,77568683,122406223,191651957,297856813,459652759,704595749,1073152385
%N A341711 a(n) = A120963(2*n+1)/2.
%C A341711 A bisection of A341710.
%H A341711 D. Weigel, R. Veysseyre, T. Phan, J. M. Effantin, and Y. Billiet, <a href="https://doi.org/10.1107/S0108767384000702">Crystallography, geometry and physics in higher dimensions. I. Point-symmetry operations</a>, Acta Cryst., A40 (1984), 323-330 (see Table 3).
%p A341711 with(numtheory):
%p A341711 b:= proc(n) option remember; nops(invphi(n)) end:
%p A341711 g:= proc(n) option remember; `if`(n=0, 1, add(
%p A341711       g(n-j)*add(d*b(d), d=divisors(j)), j=1..n)/n)
%p A341711     end:
%p A341711 a:= n-> g(2*n+1)/2:
%p A341711 seq(a(n), n=0..40);  # _Alois P. Heinz_, Feb 19 2021
%t A341711 terms = 64; (* number of terms of A120963 *)
%t A341711 nmax = Floor[terms/2] - 1;
%t A341711 S[m_] := S[m] = CoefficientList[Product[1/(1 - x^EulerPhi[k]),
%t A341711      {k, 1, m*terms}] + O[x]^(terms + 1),x];
%t A341711 S[m = 1];
%t A341711 S[m++];
%t A341711 While[S[m] != S[m - 1], m++];
%t A341711 A120963 = S[m];
%t A341711 a[n_ /; 0 <= n <= nmax] := A120963[[2 n + 2]]/2;
%t A341711 Table[a[n], {n, 0, nmax}] (* _Jean-François Alcover_, May 12 2022 *)
%Y A341711 Cf. A120963, A341710, A341712.
%K A341711 nonn
%O A341711 0,2
%A A341711 _N. J. A. Sloane_, Feb 19 2021