cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A022586 Expansion of Product_{m>=1} (1+x^m)^21.

This page as a plain text file.
%I A022586 #22 Sep 08 2022 08:44:46
%S A022586 1,21,231,1792,11067,58002,268093,1120899,4315269,15497986,52441347,
%T A022586 168487473,517184185,1524390777,4332440454,11914441196,31798680774,
%U A022586 82574231187,209091601271,517272712845,1252351944165,2971700764941,6920411525727,15835150526244,35640093688017
%N A022586 Expansion of Product_{m>=1} (1+x^m)^21.
%H A022586 Seiichi Manyama, <a href="/A022586/b022586.txt">Table of n, a(n) for n = 0..10000</a>
%F A022586 a(n) ~ 7^(1/4) * exp(Pi * sqrt(7*n)) / (4096 * n^(3/4)). - _Vaclav Kotesovec_, Mar 05 2015
%F A022586 a(0) = 1, a(n) = (21/n)*Sum_{k=1..n} A000593(k)*a(n-k) for n > 0. - _Seiichi Manyama_, Apr 04 2017
%t A022586 nmax=50; CoefficientList[Series[Product[(1+q^m)^21,{m,1,nmax}],{q,0,nmax}],q] (* _Vaclav Kotesovec_, Mar 05 2015 *)
%o A022586 (PARI) m=50; q='q+O('q^m); Vec(prod(n=1,m,(1+q^n)^21)) \\ _G. C. Greubel_, Feb 25 2018
%o A022586 (Magma) Coefficients(&*[(1+x^m)^21:m in [1..40]])[1..40] where x is PolynomialRing(Integers()).1; // _G. C. Greubel_, Feb 25 2018
%Y A022586 Column k=21 of A286335.
%K A022586 nonn
%O A022586 0,2
%A A022586 _N. J. A. Sloane_