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.

Showing 1-5 of 5 results.

A001517 Bessel polynomials y_n(x) (see A001498) evaluated at 2.

Original entry on oeis.org

1, 3, 19, 193, 2721, 49171, 1084483, 28245729, 848456353, 28875761731, 1098127402131, 46150226651233, 2124008553358849, 106246577894593683, 5739439214861417731, 332993721039856822081, 20651350143685984386753
Offset: 0

Views

Author

Keywords

Comments

Numerators of successive convergents to e using continued fraction 1 + 2/(1 + 1/(6 + 1/(10 + 1/(14 + 1/(18 + 1/(22 + 1/26 + ...)))))).
Number of ways to use the elements of {1,...,k}, n <= k <= 2n, once each to form a collection of n lists, each having length 1 or 2. - Bob Proctor, Apr 18 2005, Jun 26 2006

References

  • L. Euler, 1737.
  • I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products, 6th ed., Section 0.126, p. 2.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Essentially the same as A080893.
a(n) = A099022(n)/n!.
Partial sums: A105747.
Replace "lists" with "sets" in comment: A001515.

Programs

  • Maple
    A:= gfun:-rectoproc({a(n) = (4*n-2)*a(n-1) + a(n-2),a(0)=1,a(1)=3},a(n),remember):
    map(A, [$0..20]); # Robert Israel, Jul 22 2015
    f:=proc(n) option remember; if n = 0 then 1 elif n=1 then 3 else f(n-2)+(4*n-2)*f(n-1); fi; end;
    [seq(f(n), n=0..20)]; # N. J. A. Sloane, May 09 2016
    seq(simplify(KummerU(-n, -2*n, 1)), n = 0..16); # Peter Luschny, May 10 2022
  • Mathematica
    Table[(2k)! Hypergeometric1F1[-k, -2k, 1]/k!, {k, 0, 10}] (* Vladimir Reshetnikov, Feb 16 2011 *)
  • PARI
    a(n)=sum(k=0,n,(n+k)!/k!/(n-k)!)
    
  • Sage
    A001517 = lambda n: hypergeometric([-n, n+1], [], -1)
    [simplify(A001517(n)) for n in (0..16)] # Peter Luschny, Oct 17 2014

Formula

a(n) = Sum_{k=0..n} (n+k)!/(k!*(n-k)!) = (e/Pi)^(1/2) K_{n+1/2}(1/2).
D-finite with recurrence a(n) = (4*n-2)*a(n-1) + a(n-2), n >= 2.
a(n) = (1/n!)*Sum_{k=0..n} (-1)^(n+k)*binomial(n,k)*A000522(n+k). - Vladeta Jovovic, Sep 30 2006
E.g.f. (for offset 1): exp(x*c(x)), where c(x)=(1-sqrt(1-4*x))/(2*x) (cf. A000108). - Vladimir Kruchinin, Aug 10 2010
G.f.: 1/Q(0), where Q(k) = 1 - x - 2*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 17 2013
a(n) = (1/n!)*Integral_{x>=0} (x*(1 + x))^n*exp(-x) dx. Expansion of exp(x) in powers of y = x*(1 - x): exp(x) = 1 + y + 3*y^2/2! + 19*y^3/3! + 193*y^4/4! + 2721*y^5/5! + .... - Peter Bala, Dec 15 2013
a(n) = exp(1/2) / sqrt(Pi) * BesselK(n+1/2, 1/2). - Vaclav Kotesovec, Mar 15 2014
a(n) ~ 2^(2*n+1/2) * n^n / exp(n-1/2). - Vaclav Kotesovec, Mar 15 2014
a(n) = hypergeom([-n, n+1], [], -1). - Peter Luschny, Oct 17 2014
From G. C. Greubel, Aug 16 2017: (Start)
a(n) = (1/2)_{n} * 4^n * hypergeometric1f1(-n; -2*n; 1).
G.f.: (1/(1-t))*hypergeometric2f0(1, 1/2; -; 4*t/(1-t)^2). (End)
a(n) = Sum_{k=0..n} binomial(n,k)*binomial(n+k,k)*k!. - Ilya Gutkovskiy, Nov 24 2017
a(n) = KummerU(-n, -2*n, 1). - Peter Luschny, May 10 2022

Extensions

More terms from Vladeta Jovovic, Apr 03 2000
Additional comments from Michael Somos, Jul 15 2002

A105749 Number of ways to use the elements of {1,...,k}, 0 <= k <= 2n, once each to form a sequence of n sets, each having 1 or 2 elements.

Original entry on oeis.org

1, 2, 14, 222, 6384, 291720, 19445040, 1781750880, 214899027840, 33007837322880, 6290830003852800, 1456812592995513600, 402910665227270323200, 131173228963370155161600, 49656810289225281849907200, 21628258853895305337293568000, 10739534026001485514941587456000
Offset: 0

Views

Author

Robert A. Proctor (www.math.unc.edu/Faculty/rap/), Apr 18 2005

Keywords

Comments

Equivalently, number of sequences of n labeled items such that each item occurs just once or twice. - David Applegate, Dec 08 2008
Also, number of assembly trees for a certain star graph, see Vince-Bona, Theorem 4. - N. J. A. Sloane, Oct 08 2012

Examples

			a(2) = 14 = |{ ({1},{2}), ({2},{1}), ({1},{2,3}), ({2,3},{1}), ({2},{1,3}), ({1,3},{2}), ({3},{1,2}), ({1,2},{3}), ({1,2},{3,4}), ({3,4},{1,2}), ({1,3},{2,4}), ({2,4},{1,3}), ({1,4},{2,3}), ({2,3},{1,4}) }|.
		

Crossrefs

Replace "sets" by "lists": A099022.
Column n=2 of A181731.

Programs

  • Magma
    [(&+[Binomial(n,j)*Factorial(n+j)/2^j: j in [0..n]]): n in [0..30]]; // G. C. Greubel, Sep 26 2023
    
  • Maple
    a:= n-> add(binomial(n, k)*(n+k)!/2^k, k=0..n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Jul 21 2012
  • Mathematica
    f[n_]:= Sum[Binomial[n,k]*(n+k)!/2^k, {k,0,n}]; Table[f[n], {n,0,20}]
  • SageMath
    [sum(binomial(n,j)*factorial(n+j)//2^j for j in range(n+1)) for n in range(31)] # G. C. Greubel, Sep 26 2023

Formula

a(n) = Sum_{k=0..n} C(n,k) * (n+k)! / 2^k.
a(n) = n! * A001515(n).
A003011(n) = Sum_{k=0..n} C(n, k)*a(k).
a(n) = Gamma(n+1)*Hypergeometric2F0([-n, n+1], [], -1/2). - Peter Luschny, Jul 29 2014
a(n) ~ sqrt(Pi) * 2^(n + 1) * n^(2*n + 1/2) / exp(2*n - 1). - Vaclav Kotesovec, Nov 27 2017
From G. C. Greubel, Sep 26 2023: (Start)
a(n) = n*(2*n-1)*a(n-1) + n*(n-1)*a(n-2).
a(n) = e * sqrt(2/Pi) * n! * BesselK(n+1/2, 1).
a(n) = ((2*n)!/2^n) * Hypergeometric1F1(-n, -2*n, 2).
G.f.: (-2/x) * Integrate_{t=0..oo} exp(-t)/((t+1)^2 - 1 - 2/x) dt.
G.f.: e*( exp(-sqrt(1 + 2/x)) * ExpIntegralEi(-1 + sqrt(1 + 2/x)) - exp(sqrt(1 + 2/x)) * ExpIntegralEi(-1 - sqrt(1 + 2/x)) )/sqrt(x^2 + 2*x).
E.g.f.: ((1-x)/x) * Hypergeometric1F1(1, 3/2, -(1-x)^2/(2*x)).
E.g.f.: (1/(1-x))*Hypergeometric2F0([1, 1/2]; []; 2*x/(1-x)^2). (End)

Extensions

More terms from Robert G. Wilson v, Apr 23 2005

A076571 Binomial triangle based on factorials.

Original entry on oeis.org

1, 1, 2, 2, 3, 5, 6, 8, 11, 16, 24, 30, 38, 49, 65, 120, 144, 174, 212, 261, 326, 720, 840, 984, 1158, 1370, 1631, 1957, 5040, 5760, 6600, 7584, 8742, 10112, 11743, 13700, 40320, 45360, 51120, 57720, 65304, 74046, 84158, 95901, 109601
Offset: 0

Views

Author

Henry Bottomley, Oct 19 2002

Keywords

Examples

			Rows start:
    1;
    1,   2;
    2,   3,   5;
    6,   8,  11,  16;
   24,  30,  38,  49,  65;
  120, 144, 174, 212, 261, 326;
		

Crossrefs

Right hand columns include A000522, A001339, A001340, A001341, A001342.
Cf. A002627 (row sums), A099022.

Programs

  • Magma
    A076571:= func< n,k| (&+[Binomial(k,j)*Factorial(n-j): j in [0..k]]) >;
    [A076571(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 05 2023
    
  • Mathematica
    A076571[n_, k_]:= n!*Hypergeometric1F1[-k,-n,1];
    Table[A076571[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 05 2023 *)
  • SageMath
    def A076571(n,k): return sum(binomial(k,j)*factorial(n-j) for j in range(k+1))
    flatten([[A076571(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Oct 05 2023

Formula

T(n, k) = Sum_{j=0..k} binomial(k, j)*(n-j)!.
T(n, k) = T(n, k-1) + T(n-1, k-1) with T(n, 0) = n!.
T(n, n) = A000522(n).
Sum_{k=0..n} T(n, k) = A002627(n+1).
From G. C. Greubel, Oct 05 2023: (Start)
T(n, k) = n! * Hypergeometric1F1([-k], [-n], 1).
T(2*n, n) = A099022(n). (End)

A328826 Triangle read by rows: binomial(n,k)*(2*n-k)!, n>=0, 0<=k<=n.

Original entry on oeis.org

1, 2, 1, 24, 12, 2, 720, 360, 72, 6, 40320, 20160, 4320, 480, 24, 3628800, 1814400, 403200, 50400, 3600, 120, 479001600, 239500800, 54432000, 7257600, 604800, 30240, 720, 87178291200, 43589145600, 10059033600, 1397088000, 127008000, 7620480, 282240, 5040, 20922789888000
Offset: 0

Views

Author

R. J. Mathar, Oct 28 2019

Keywords

Comments

Vertex-labeled disconnected Goldstone diagrams with n vertices and k single-particle potentials.

Examples

			The triangle starts
    1;
    2     1;
   24    12     2;
  720   360    72     6;
40320 20160  4320   480    24;
		

Crossrefs

Cf. A099022 (row sums), A000142 (diagonal), A010050 (column k=0), A002674 (k=1).

Programs

  • Maple
    A328826 := proc(n,k)
            binomial(n,k)*(2*n-k)! ;
    end proc:
  • Mathematica
    Table[Binomial[n,k](2n-k)!,{n,0,10},{k,0,n}]//Flatten (* Harvey P. Dale, Feb 03 2022 *)

Formula

T(n,k)= binomial(n,k)*(2*n-k)!.
T(n,k) = A328921(n,k) + A328922(n,k). - R. J. Mathar, Nov 02 2019

A082765 Trinomial transform of the factorial numbers (A000142).

Original entry on oeis.org

1, 4, 45, 1282, 70177, 6239016, 817234189, 147950506390, 35370826189857, 10791515504716012, 4091225768720823181, 1886585105032464025674, 1039774852573506696192385, 674970732343624159361034832
Offset: 0

Views

Author

Emanuele Munarini, May 21 2003

Keywords

Comments

Number of ways to use the elements of {1,..,k}, 0<=k<=2n, once each to form a sequence of n (possibly empty) lists, each of length at most 2. - Bob Proctor, Apr 18 2005

Crossrefs

a(n) = Sum[C(n, k)*A099022(k), 0<=k<=n]
Replace "sequence" by "collection" in comment: A105747.
Replace "lists" by "sets" in comment: A003011.

Formula

a(n) = Sum[ Trinomial[n, k] k!, {k, 0, 2n} ] where Trinomial[n, k] = trinomial coefficients (A027907)
Integral_{x=0..infinity} (x^2+x+1)^n*exp(-x) dx - Gerald McGarvey, Oct 14 2006
Showing 1-5 of 5 results.