%I #43 Jul 08 2023 16:09:21
%S 1,0,0,0,0,0,1,0,0,0,1,0,1,0,0,1,1,0,1,0,1,1,1,0,1,1,1,1,1,0,2,1,1,1,
%T 1,1,2,1,1,1,2,1,2,1,1,2,2,1,2,1,2,2,2,1,2,2,2,2,2,1,3,2,2,2,2,2,3,2,
%U 2,2,3,2,3,2,2,3,3,2,3,2,3,3,3,2,3,3,3,3,3,2,4,3,3,3,3,3,4,3,3,3,4,3,4,3,3,4
%N Molien series of binary icosahedral group.
%C Meyer's generating function h(t,G) generates the sequence of the dimensions of the spaces of G-invariant harmonic polynomials of each degree, where G is a point group on three-dimensional Euclidean space. For G=I, the icosahedral rotation group, the generating function gives rise to this sequence. See Table 1, p. 143. - _William Lionheart_, May 04 2019
%D T. A. Springer, Invariant Theory, Lecture Notes in Math., Vol. 585, Springer, p. 97.
%D S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 19.
%H G. C. Greubel, <a href="/A008651/b008651.txt">Table of n, a(n) for n = 0..1000</a>
%H Burnett Meyer, <a href="https://doi.org/10.4153/CJM-1954-016-2">On the symmetries of spherical harmonics</a>, Canadian Journal of Mathematics 6 (1954): 135-157.
%H <a href="/index/Mo#Molien">Index entries for Molien series</a>
%H <a href="/index/Rec#order_09">Index entries for linear recurrences with constant coefficients</a>, signature (-1,0,1,1,1,1,0,-1,-1).
%F G.f.: (1 +x -x^3 -x^4 -x^5 +x^7 +x^8)/((1+x)*(1-x)^2*(1+x+x^2)*(1+x+x^2+x^3+x^4)). - _R. J. Mathar_, Dec 01 2014
%F Euler transform of length 30 sequence [ 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1]. - _Michael Somos_, Dec 01 2014
%F a(n) = -a(-1-n) for all n in Z. - _Michael Somos_, Dec 01 2014
%F 0 = 1 + a(n) + 2*a(n+1) + 2*a(n+2) + a(n+3) - a(n+5) - 2*a(n+6) - 2*a(n+7) - a(n+8) for all n in Z. - _Michael Somos_, Dec 01 2014
%F G.f.: (1+x^15)/((1-x^10)*(1-x^6)) - not reduced _William Lionheart_, May 04 2019
%e The Molien series is (1+q^20+q^40)/((1-q^12)*(1-q^30)). Since every other term would be zero, we replace q^2 with x to get the sequence.
%e G.f. = 1 + x^6 + x^10 + x^12 + x^15 + x^16 + x^18 + x^20 + x^21 + x^22 + ...
%e G.f. = 1 + q^12 + q^20 + q^24 + q^30 + q^32 + q^36 + q^40 + q^42 + q^44 + ...
%p t1:=(1+x^10+x^20)/((1-x^6)*(1-x^15));
%p series(t1,x,100);
%p seriestolist(%);
%t a[ n_] := With[ {s = Boole[ n<0 ], m = If[ n<0, -1-n, n]}, (-1)^s * SeriesCoefficient[(1+x^15)/((1-x^6)*(1-x^10)), {x, 0 ,m}]]; (* _Michael Somos_, Dec 01 2014 *)
%t LinearRecurrence[{-1,0,1,1,1,1,0,-1,-1}, {1,0,0,0,0,0,1,0,0}, 100] (* _Harvey P. Dale_, May 04 2017 *)
%o (PARI) Vec(O(x^99)+(1+x^10+x^20)/((1-x^6)*(1-x^15))) \\ _M. F. Hasler_, Dec 01 2014
%o (PARI) {a(n) = my(s=n<0); if(s, n = -1-n); (-1)^s * polcoeff( (1 + x^15) / ( (1 - x^6) * (1 - x^10) ) + x * O(x^n), n)}; /* _Michael Somos_, Dec 01 2014 */
%o (PARI) {a(n) = (n\30) + [0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1][n%30 + 1]}; /* _Michael Somos_, Dec 01 2014 */
%o (Magma) I:=[1,0,0,0,0,0,1,0,0]; [n le 9 select I[n] else -Self(n-1) +Self(n-3)+Self(n-4)+Self(n-5)+Self(n-6)-Self(n-8)-Self(n-9): n in [1..100]]; // _Vincenzo Librandi_, Jun 24 2015
%o (Sage)
%o def A008651_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P((1+x^10+x^20)/((1-x^6)*(1-x^15))).list()
%o A008651_list(100) # _G. C. Greubel_, Sep 07 2019
%Y Cf. A319974 for harmonic polynomials in four variables invariant under a group.
%K nonn,easy
%O 0,31
%A _N. J. A. Sloane_