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-3 of 3 results.

A003011 Number of permutations of up to n kinds of objects, where each kind of object can occur at most two times.

Original entry on oeis.org

1, 3, 19, 271, 7365, 326011, 21295783, 1924223799, 229714292041, 35007742568755, 6630796801779771, 1527863209528564063, 420814980652048751629, 136526522051229388285611
Offset: 0

Views

Author

Keywords

Comments

E.g.f. A(x)=y satisfies 0=(2x^3+2x^2)y''+(-3x^3+4x-1)y'+(x^3-x^2-2x+3)y. - Michael Somos, Mar 15 2004
Number of ways to use the elements of {1,..,k}, 0<=k<=2n, once each to form a sequence of n (possibly empty) sets, each having at most 2 elements. - Bob Proctor, Apr 18 2005

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 17.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = Sum[C(n, k)*A105749(k), 0<=k<=n]
Replace "sequence" with "collection" in comment: A105748.
Replace "sets" with "lists" in comment: A082765.

Programs

  • Mathematica
    Table[nn=2n;a=1+x+x^2/2!;Total[Range[0,nn]!CoefficientList[Series[a^n,{x,0,nn}],x]],{n,0,15}]  (* Geoffrey Critzer, Dec 23 2011 *)
  • PARI
    a(n)=local(A);if(n<0,0,A=(1+x+x^2/2)^n;sum(k=0,2*n,k!*polcoeff(A,k)))

Formula

n*a(n) = (2*n^3 - n^2 + n + 1)*a(n-1) + (-3*n^3 + 4*n^2 + 2*n - 3)*a(n-2) + (n^3 - 2*n^2 - n + 2)*a(n-3).
a(n) ~ sqrt(Pi)*2^(n+1)*n^(2*n+1/2)/exp(2*n-1). - Vaclav Kotesovec, Oct 19 2013

Extensions

More terms from Vladeta Jovovic, Aug 18 2002

A099022 a(n) = Sum_{k=0..n} C(n,k)*(2*n-k)!.

Original entry on oeis.org

1, 3, 38, 1158, 65304, 5900520, 780827760, 142358474160, 34209760152960, 10478436416945280, 3984884716852972800, 1842169367191937414400, 1017403495472574045158400, 661599650478455071589606400, 500354503197888042597961267200, 435447353708763072625260119808000
Offset: 0

Views

Author

Ralf Stephan, Sep 23 2004

Keywords

Comments

Diagonal of Euler-Seidel matrix with start sequence n!.
Number of ways to use the elements of {1,..,k}, n<=k<=2n, once each to form a sequence of n lists, each having length 1 or 2. - Bob Proctor, Apr 18 2005, Jun 26 2006
Replace "lists" by "sets": A105749.

Crossrefs

Cf. A001517, A076571, A082765 (binomial transform), A105749, row sums of A328826.

Programs

  • Maple
    f:= gfun:-rectoproc({a(n)=2*n*(2*n-1)*a(n-1)+n*(n-1)*a(n-2), a(0)=1,a(1)=3},a(n),remember):
    map(f, [$0..20]); # Robert Israel, Feb 15 2017
  • Mathematica
    Table[(2k)! Hypergeometric1F1[-k, -2k, 1], {k, 0, 10}] (* Vladimir Reshetnikov, Feb 16 2011 *)
    Table[Sum[Binomial[n,k](2n-k)!,{k,0,n}],{n,0,20}] (* Harvey P. Dale, Nov 22 2021 *)
  • PARI
    for(n=0,25, print1(sum(k=0,n, binomial(n,k)*(2*n-k)!), ", ")) \\ G. C. Greubel, Dec 31 2017

Formula

T(2*n, n), where T is the triangle in A076571.
a(n) = n!*A001517(n).
A082765(n) = Sum[C(n, k)*a(k), 0<=k<=n].
a(n) = 2*n*(2*n-1)*a(n-1)+n*(n-1)*a(n-2). - Vladeta Jovovic, Sep 27 2004
a(n) = int {x = 0..inf} exp(-x)*(x + x^2)^n dx. Applying the results of Nicolaescu, Section 3.2 to this integral we obtain the asymptotic expansion a(n) ~ (2*n)!*exp(1/2)*( 1 - 1/(16*n) - 191/(6144*n^2) + O(1/n^3) ). - Peter Bala, Jul 07 2014

A105747 Number of ways to use the elements of {1,..,k}, 0<=k<=2n, once each to form a collection of n (possibly empty) lists, each of length at most 2.

Original entry on oeis.org

1, 4, 23, 216, 2937, 52108, 1136591, 29382320, 877838673, 29753600404, 1127881002535, 47278107653768, 2171286661012617, 108417864555606300, 5847857079417024031, 338841578119273846112
Offset: 0

Views

Author

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

Keywords

Examples

			a(2)=23:
{(),()},
{(),(1)},
{(),(1,2)},
{(),(2,1)},
{(1),(2)},
{(1),(2,3)},
{(1),(3,2)},
...,
{(1,4),(2,3)},
{(1,4),(3,2)},
{(4,1),(2,3)},
{(4,1),(3,2)}.
		

Crossrefs

First differences: A001517.
Replace "collection" by "sequence": A082765.
Replace "lists" by "sets": A105748.

Programs

  • Mathematica
    Table[Sum[(k+i)!/i!/(k-i)!, {k, 0, n}, {i, 0, k}], {n, 0, 20}]

Formula

a(n) = Sum_{0<=i<=k<=n} (k+i)!/i!/(k-i)!.
a(n+3) = (4*n+11)*a(n+2) - (4*n+9)*a(n+1) - a(n) - Benoit Cloitre, May 26 2006
G.f.: 1/(1-x)/Q(0), where Q(k)= 1 - x - 2*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 17 2013
a(n) ~ 2^(2*n + 1/2) * n^n / exp(n - 1/2). - Vaclav Kotesovec, May 15 2022
Showing 1-3 of 3 results.