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.

A118930 E.g.f.: A(x) = exp( Sum_{n>=0} x^(2^n)/2^(2^n-1) ).

This page as a plain text file.
%I A118930 #19 Mar 26 2020 05:05:44
%S A118930 1,1,2,4,13,41,166,652,3494,18118,114076,681176,5016892,35377564,
%T A118930 288204008,2232198256,21124254181,191779964597,2011347229114,
%U A118930 19840403629108,231266808172181,2553719667653281,31743603728993542
%N A118930 E.g.f.: A(x) = exp( Sum_{n>=0} x^(2^n)/2^(2^n-1) ).
%C A118930 Equals invariant column vector V that satisfies matrix product A100861*V = V, where Bessel numbers A100861(n,k) = n!/[k!(n-2k)!*2^k] give the number of k-matchings of the complete graph K(n).
%C A118930 Equals Lim_{n->inf.} A144299^n, if A144299 is considered an infinite lower triangular matrix. - _Gary W. Adamson_, Dec 08 2008
%F A118930 a(n) = Sum_{k=0..[n/2]} n!/[k!*(n-2*k)!*2^k] * a(k), with a(0)=1. a(n) = Sum_{k=0..[n/2]} A100861(n,k)*a(k), with a(0)=1.
%e A118930 E.g.f. A(x) = exp( x + x^2/2 + x^4/2^3 + x^8/2^7 + x^16/2^15 +...)
%e A118930 = 1 + 1*x + 2*x^2/2! + 4*x^3/3! + 13*x^4/4! + 41*x^5/5!+ 166*x^6/6!+...
%e A118930 Using coefficients A100861(n,k) = n!/[k!(n-2k)!*2^k]:
%e A118930 a(5) = 1*a(0) +10*a(1) +15*a(2) = 1*1 +10*1 +15*2 = 41.
%e A118930 a(6) = 1*a(0) +15*a(1) +45*a(2) +15*a(3) = 1*1 +15*1 +45*2 +15*4 = 166.
%p A118930 A118930 := proc(n)
%p A118930     option remember;
%p A118930     if n<= 1 then
%p A118930         1 ;
%p A118930     else
%p A118930         n!*add(procname(k)/k!/(n-2*k)!/2^k,k=0..n/2) ;
%p A118930     end if;
%p A118930 end proc;
%p A118930 seq(A118930(n),n=0..10) ; # _R. J. Mathar_, Aug 19 2014
%t A118930 a[n_] := a[n] = If[n==0, 1, Sum[Binomial[n, 2k] (2k-1)!! a[k], {k, 0, n/2}]];
%t A118930 a /@ Range[0, 22] (* _Jean-François Alcover_, Mar 26 2020 *)
%o A118930 (PARI) {a(n)=if(n==0,1,sum(k=0,n\2,n!/(k!*(n-2*k)!*2^k)*a(k)))}
%o A118930 for(n=0,30,print1(a(n),", "))
%o A118930 (PARI) /* Defined by E.G.F.: */
%o A118930 {a(n)=n!*polcoeff( exp(sum(k=0,#binary(n),x^(2^k)/2^(2^k-1))+x*O(x^n)),n,x)}
%o A118930 for(n=0,30,print1(a(n),", "))
%Y A118930 Cf. A100861; variants: A118932, A118935.
%Y A118930 Equals row sums of triangle A152685. - _Gary W. Adamson_, Dec 10 2008
%Y A118930 Cf. A144299. - _Gary W. Adamson_, Dec 08 2008
%K A118930 nonn
%O A118930 0,3
%A A118930 _Paul D. Hanna_, May 06 2006