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.

Previous Showing 101-110 of 382 results. Next

A368951 Number of connected labeled graphs with n edges and n vertices and with loops allowed.

Original entry on oeis.org

1, 1, 2, 10, 79, 847, 11436, 185944, 3533720, 76826061, 1880107840, 51139278646, 1530376944768, 49965900317755, 1767387701671424, 67325805434672100, 2747849045156064256, 119626103584870552921, 5533218319763109888000, 270982462739224265922466
Offset: 0

Views

Author

Andrew Howroyd, Jan 10 2024

Keywords

Comments

Exponential transform appears to be A333331. - Gus Wiseman, Feb 12 2024

Examples

			From _Gus Wiseman_, Feb 12 2024: (Start)
The a(0) = 1 through a(3) = 10 loop-graphs:
  {}  {11}  {11,12}  {11,12,13}
            {22,12}  {11,12,23}
                     {11,13,23}
                     {22,12,13}
                     {22,12,23}
                     {22,13,23}
                     {33,12,13}
                     {33,12,23}
                     {33,13,23}
                     {12,13,23}
(End)
		

Crossrefs

This is the connected covering case of A014068.
The case without loops is A057500, covering case of A370317.
Allowing any number of edges gives A062740, connected case of A322661.
This is the connected case of A368597.
The unlabeled version is A368983, connected case of A368984.
For at most n edges we have A369197.
A000085 counts set partitions into singletons or pairs.
A006129 counts covering graphs, connected A001187.

Programs

  • Maple
    egf:= (L-> 1-L/2-log(1+L)/2-L^2/4)(LambertW(-x)):
    a:= n-> n!*coeff(series(egf, x, n+1), x, n):
    seq(a(n), n=0..25);  # Alois P. Heinz, Jan 10 2024
  • PARI
    seq(n)={my(t=-lambertw(-x + O(x*x^n))); Vec(serlaplace(-log(1-t)/2 + t/2 - t^2/4 + 1))}

Formula

a(n) = A000169(n) + A057500(n) for n > 0.
E.g.f.: 1 - log(1-T(x))/2 + T(x)/2 - T(x)^2/4 where T(x) = -LambertW(-x) is the e.g.f. of A000169.
From Peter Luschny, Jan 10 2024: (Start)
a(n) = (exp(n)*Gamma(n + 1, n) - (n - 1)*n^(n - 1))/(2*n) for n > 0.
a(n) = (1/2)*(A063170(n)/n - A053506(n)) for n > 0. (End)

A037205 a(n) = (n+1)^n - 1.

Original entry on oeis.org

0, 1, 8, 63, 624, 7775, 117648, 2097151, 43046720, 999999999, 25937424600, 743008370687, 23298085122480, 793714773254143, 29192926025390624, 1152921504606846975, 48661191875666868480, 2185911559738696531967, 104127350297911241532840, 5242879999999999999999999, 278218429446951548637196400, 15519448971100888972574851071
Offset: 0

Views

Author

Keywords

Comments

For n >= 1, a(n) = order of Fibonacci group F(n+1,n).
The terms, written in base n+1, are n digits of value n. For example, a(4) = 624 = 4444 in base 5. - Marc Morgenegg, Nov 30 2016
For n >= 1, in a square grid of side n, this is the number of ways to populate the grid with 1 X 1 blocks (with at least one block) so that no block falls under the effect of gravity. - Paolo Xausa, Apr 12 2021
For n > 1, (n-1)^2 | a(n). - David A. Corneth, Dec 15 2022

References

  • D. L. Johnson, Presentation of Groups, Cambridge, 1976, p. 182.
  • Richard M. Thomas, The Fibonacci groups revisited, in Groups - St. Andrews 1989, Vol. 2, 445-454, London Math. Soc. Lecture Note Ser., 160, Cambridge Univ. Press, Cambridge, 1991.

Crossrefs

A diagonal of A202624.

Programs

  • Magma
    [(n + 1)^n - 1: n in [0..25]]; // G. C. Greubel, Nov 10 2017
  • Mathematica
    Table[(n + 1)^n - 1, {n, 0, 21}] (* or *)
    Table[If[n < 1, Length@ #, FromDigits[#, n + 1]] &@ ConstantArray[n, n], {n, 0, 21}] (* Michael De Vlieger, Nov 30 2016 *)
  • PARI
    for(n=0,25, print1((n + 1)^n - 1, ", ")) \\ G. C. Greubel, Nov 10 2017
    

Formula

a(n) = A000169(n+1) - 1 = A060072(n+1)*(n-1) = A060073(n+1)*(n-1)^2.
E.g.f.: 1/(exp(LambertW(-x)) - x) - exp(x). - Ilya Gutkovskiy, Nov 30 2016
E.g.f.: -exp(x) - 1/(x + x/LambertW(-x)). - Vaclav Kotesovec, Dec 05 2016
a(n) = Sum_{k=1..n} binomial(n,k)*n^k [from Paolo Xausa's comment]. - Joerg Arndt, Apr 12 2021

Extensions

Revised by N. J. A. Sloane, Dec 30 2011

A060313 Number of homeomorphically irreducible rooted trees (also known as series-reduced rooted trees, or rooted trees without nodes of degree 2) on n labeled nodes.

Original entry on oeis.org

1, 2, 0, 16, 25, 576, 2989, 51584, 512649, 8927200, 130956001, 2533847328, 48008533885, 1059817074512, 24196291364925, 609350187214336, 16135860325700881, 459434230368302016, 13788624945433889593, 439102289933675933600, 14705223056221892676741
Offset: 1

Views

Author

Vladeta Jovovic, Mar 27 2001

Keywords

Examples

			From _Gus Wiseman_, Jan 22 2020: (Start)
The a(1) = 1 through a(4) = 16 trees (in the format root[branches], empty column shown as dot) are:
  1  1[2]  .  1[2,3,4]
     2[1]     1[2[3,4]]
              1[3[2,4]]
              1[4[2,3]]
              2[1,3,4]
              2[1[3,4]]
              2[3[1,4]]
              2[4[1,3]]
              3[1,2,4]
              3[1[2,4]]
              3[2[1,4]]
              3[4[1,2]]
              4[1,2,3]
              4[1[2,3]]
              4[2[1,3]]
              4[3[1,2]]
(End)
		

References

  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.

Crossrefs

The unlabeled unrooted version is A000014.
The unrooted version is A005512.
The unlabeled version is A001679 or A059123.
The lone-child-avoiding version is A060356.
Labeled rooted trees are A000169.

Programs

  • Magma
    [1] cat [n*Factorial(n-2)*(&+[(-1)^k*Binomial(n,k)*(n-k)^(n-k-2)/Factorial(n-k-2): k in [0..n-2]]): n in [2..20]]; // G. C. Greubel, Mar 07 2020
    
  • Maple
    seq( `if`(n=1, 1, n*(n-2)!*add((-1)^k*binomial(n, k)*(n-k)^(n-k-2)/(n-k-2)!, k=0..n-2)), n=1..20); # G. C. Greubel, Mar 07 2020
  • Mathematica
    f[n_] := If[n < 2, 1, n(n - 2)!Sum[(-1)^k*Binomial[n, k](n - k)^(n - 2 - k)/(n - 2 - k)!, {k, 0, n - 2}]]; Table[ f[n], {n, 19}] (* Robert G. Wilson v, Feb 12 2005 *)
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    lrt[set_]:=If[Length[set]==0,{},Join@@Table[Apply[root,#]&/@Join@@Table[Tuples[lrt/@stn],{stn,sps[DeleteCases[set,root]]}],{root,set}]];
    Table[Length[Select[lrt[Range[n]],Length[#]!=2&&FreeQ[Z@@#,Integer[]]&]],{n,6}] (* Gus Wiseman, Jan 22 2020 *)
  • Sage
    [1]+[n*factorial(n-2)*sum((-1)^k*binomial(n,k)*(n-k)^(n-k-2)/factorial( n-k-2) for k in (0..n-2)) for n in (2..20)] # G. C. Greubel, Mar 07 2020

Formula

a(n) = n*(n-2)!*Sum_{k=0..n-2} (-1)^k*binomial(n, k)*(n-k)^(n-k-2)/(n-k-2)!, n>1.
E.g.f.: x*(exp( - LambertW(-x/(1+x))) - (LambertW(-x/(1+x))/2 )^2).
a(n) ~ n^(n-1) * (1-exp(-1))^(n+1/2). - Vaclav Kotesovec, Oct 05 2013
E.g.f.: -(1+x)*LambertW(-x/(1+x)) - (x/2)*LambertW(-x/(1+x))^2. - G. C. Greubel, Mar 07 2020

A368924 Triangle read by rows where T(n,k) is the number of labeled loop-graphs on n vertices with k loops and n-k non-loops such that it is possible to choose a different vertex from each edge.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 1, 9, 6, 1, 15, 68, 48, 12, 1, 222, 720, 510, 150, 20, 1, 3670, 9738, 6825, 2180, 360, 30, 1, 68820, 159628, 110334, 36960, 6895, 735, 42, 1, 1456875, 3067320, 2090760, 721560, 145530, 17976, 1344, 56, 1, 34506640, 67512798, 45422928, 15989232, 3402756, 463680, 40908, 2268, 72, 1
Offset: 0

Views

Author

Gus Wiseman, Jan 10 2024

Keywords

Comments

The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			Triangle begins:
      1
      0      1
      0      2      1
      1      9      6      1
     15     68     48     12      1
    222    720    510    150     20      1
   3670   9738   6825   2180    360     30      1
  68820 159628 110334  36960   6895    735     42      1
Row n = 3 counts the following loop-graphs:
  {{1,2},{1,3},{2,3}}  {{1},{1,2},{1,3}}  {{1},{2},{1,3}}  {{1},{2},{3}}
                       {{1},{1,2},{2,3}}  {{1},{2},{2,3}}
                       {{1},{1,3},{2,3}}  {{1},{3},{1,2}}
                       {{2},{1,2},{1,3}}  {{1},{3},{2,3}}
                       {{2},{1,2},{2,3}}  {{2},{3},{1,2}}
                       {{2},{1,3},{2,3}}  {{2},{3},{1,3}}
                       {{3},{1,2},{1,3}}
                       {{3},{1,2},{2,3}}
                       {{3},{1,3},{2,3}}
		

Crossrefs

Column k = n-1 is A002378.
The case of a unique choice is A061356, row sums A000272.
Column k = 0 is A137916, unlabeled version A137917.
Row sums appear to be A333331.
The complement has row sums A368596, covering case A368730.
The unlabeled version is A368926.
Without the choice condition we have A368928, A116508, A367863, A368597.
A000085, A100861, A111924 count set partitions into singletons or pairs.
A006125 counts graphs, unlabeled A000088.
A006129 counts covering graphs, unlabeled A002494.
A014068 counts loop-graphs, unlabeled A000666.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n],{1,2}],{n}], Count[#,{_}]==k&&Length[Select[Tuples[#], UnsameQ@@#&]]!=0&]],{n,0,5},{k,0,n}]
  • PARI
    T(n)={my(t=-lambertw(-x + O(x*x^n))); [Vecrev(p) | p <- Vec(serlaplace(exp(-log(1-t)/2 - t/2 + t*y - t^2/4)))]}
    { my(A=T(8)); for(i=1, #A, print(A[i])) } \\ Andrew Howroyd, Jan 14 2024

Formula

E.g.f.: A(x,y) = exp(-log(1-T(x))/2 - T(x)/2 + y*T(x) - T(x)^2/4) where T(x) = -LambertW(-x) is the e.g.f. of A000169. - Andrew Howroyd, Jan 14 2024

Extensions

a(36) onwards from Andrew Howroyd, Jan 14 2024

A008789 a(n) = n^(n+3).

Original entry on oeis.org

0, 1, 32, 729, 16384, 390625, 10077696, 282475249, 8589934592, 282429536481, 10000000000000, 379749833583241, 15407021574586368, 665416609183179841, 30491346729331195904, 1477891880035400390625, 75557863725914323419136
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • GAP
    List([0..20], n-> n^(n+3)); # G. C. Greubel, Sep 11 2019
  • Magma
    [n^(n+3): n in [0..20]]; // Vincenzo Librandi, Jun 11 2013
    
  • Maple
    printlevel := -1; a := [0]; T := x->-LambertW(-x); f := series((T(x)*(1+8*T(x)+6*(T(x))^2)/(1-T(x))^7),x,24); for m from 1 to 23 do a := [op(a),op(2*m-1,f)*m! ] od; print(a); # Len Smiley, Nov 19 2001
  • Mathematica
    Table[n^(n+3),{n,0,20}](* Vladimir Joseph Stephan Orlovsky, Dec 26 2010 *)
  • PARI
    vector(20, n, (n-1)^(n+2)) \\ G. C. Greubel, Sep 11 2019
    
  • Sage
    [n^(n+3) for n in (0..20)] # G. C. Greubel, Sep 11 2019
    

Formula

E.g.f.(x): T*(1 +8*T +6*T^2)*(1-T)^(-7); where T=T(x) is Euler's tree function (see A000169). - Len Smiley, Nov 19 2001
See A008517 and A134991 for similar e.g.f.s and diagonals of A048993. - Tom Copeland, Oct 03 2011
E.g.f.: d^3/dx^3 {x^3/(T(x)^3*(1-T(x)))}, where T(x) = Sum_{n>=1} n^(n-1)*x^n/n! is the tree function of A000169. - Peter Bala, Aug 05 2012
a(n) = n*A008788(n). - R. J. Mathar, Oct 31 2015

A055858 Coefficient triangle for certain polynomials.

Original entry on oeis.org

1, 1, 2, 4, 9, 6, 27, 64, 48, 36, 256, 625, 500, 400, 320, 3125, 7776, 6480, 5400, 4500, 3750, 46656, 117649, 100842, 86436, 74088, 63504, 54432, 823543, 2097152, 1835008, 1605632, 1404928, 1229312, 1075648, 941192, 16777216, 43046721
Offset: 0

Views

Author

Wolfdieter Lang, Jun 20 2000

Keywords

Comments

The coefficients of the partner polynomials are found in triangle A055864.

Examples

			{1}; {1,2}; {4,9,6}; {27,64,48,36}; ...
Fourth row polynomial (n=3): p(3,x) = 27 + 64*x + 48*x^2 + 36*x^3.
		

Crossrefs

Column sequences are A000312(n), n >= 1, A055860 (A000169), A055861 (A053506), A055862-3 for m=0..4, row sums: A045531(n+1)= |A039621(n+1, 2)|, n >= 0.

Programs

  • Mathematica
    a[n_, m_] /; n < m = 0; a[0, 0] = 1; a[n_, 0] := n^n; a[n_, m_] := n^(m-1)*(n+1)^(n-m+1); Table[a[n, m], {n, 0, 8}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jun 20 2013 *)

Formula

a(n, m)=0 if n < m; a(0, 0)=1, a(n, 0) = n^n, n >= 1, a(n, m) = n^(m-1)*(n+1)^(n-m+1), n >= m >= 1;
E.g.f. for column m: A(m, x); A(0, x) = 1/(1+W(-x)); A(1, x) = -1 - (d/dx)W(-x) = -1-W(-x)/((1+W(-x))*x); A(2, x) = A(1, x)-int(A(1, x), x)/x-(1/x+x); recursion: A(m, x) = A(m-1, x)-int(A(m-1, x), x)/x-((m-1)^(m-1))*(x^(m-1))/(m-1)!, m >= 3; W(x) principal branch of Lambert's function.

A076113 a(n) = n^(n*(n-1)/2).

Original entry on oeis.org

1, 1, 2, 27, 4096, 9765625, 470184984576, 558545864083284007, 19342813113834066795298816, 22528399544939174411840147874772641, 1000000000000000000000000000000000000000000000, 1890591424712781041871514584574319778449301246603238034051
Offset: 0

Views

Author

Amarnath Murthy, Oct 09 2002

Keywords

Comments

Number of labeled commutative idempotent groupoids with n elements. [edited by Michel Marcus, Jul 10 2025]
Product of terms in n-th row of A076112.

Crossrefs

Programs

  • PARI
    a(n) = n^(n*(n-1)/2); \\ Joerg Arndt, Nov 04 2013

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Apr 20 2003
a(0)=1 prepended by Alois P. Heinz, Jun 30 2022

A128433 Triangle, read by rows, T(n,k) = numerator of the maximum of the k-th Bernstein polynomial of degree n; denominator is A128434.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 4, 4, 1, 1, 27, 3, 27, 1, 1, 256, 216, 216, 256, 1, 1, 3125, 80, 5, 80, 3125, 1, 1, 46656, 37500, 34560, 34560, 37500, 46656, 1, 1, 823543, 5103, 590625, 35, 590625, 5103, 823543, 1, 1, 16777216, 13176688, 1792, 11200000, 11200000, 1792, 13176688, 16777216, 1
Offset: 0

Views

Author

Reinhard Zumkeller, Mar 03 2007

Keywords

Examples

			Triangle begins as:
  1;
  1,        1;
  1,        1,        1;
  1,        4,        4,      1;
  1,       27,        3,     27,        1;
  1,      256,      216,    216,      256,        1;
  1,     3125,       80,      5,       80,     3125,     1;
  1,    46656,    37500,  34560,    34560,    37500, 46656,        1;
  1,   823543,     5103, 590625,       35,   590625,  5103,   823543,        1;
  1, 16777216, 13176688,   1792, 11200000, 11200000,  1792, 13176688, 16777216, 1;
		

Crossrefs

Programs

  • Mathematica
    B[n_, k_]:= If[k==0 || k==n, 1, Binomial[n, k]*k^k*(n-k)^(n-k)/n^n];
    T[n_, k_]= Numerator[B[n, k]];
    Table[T[n, k], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Jul 19 2021 *)
  • Sage
    def B(n,k): return 1 if (k==0 or k==n) else binomial(n, k)*k^k*(n-k)^(n-k)/n^n
    def T(n,k): return numerator(B(n,k))
    flatten([[T(n,k) for k in (0..n)] for n in (0..10)]) # G. C. Greubel, Jul 19 2021

Formula

T(n,k)/A128434(n,k) = Binomial(n,k) * k^k * (n-k)^(n-k) / n^n.
For n>0: Sum_{k=0..n} T(n,k)/A128434(n,k) = A090878(n)/A036505(n-1).
T(n,n-k) = T(n,k).
T(n,0) = 1.
for n>0: T(n,1)/A128434(n,1) = A000312(n-1)/A000169(n).

A128434 Triangle, read by rows, T(n,k) = denominator of the maximum of the k-th Bernstein polynomial of degree n; numerator is A128433.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 9, 9, 1, 1, 64, 8, 64, 1, 1, 625, 625, 625, 625, 1, 1, 7776, 243, 16, 243, 7776, 1, 1, 117649, 117649, 117649, 117649, 117649, 117649, 1, 1, 2097152, 16384, 2097152, 128, 2097152, 16384, 2097152, 1, 1, 43046721, 43046721, 6561, 43046721, 43046721, 6561, 43046721, 43046721, 1
Offset: 0

Views

Author

Reinhard Zumkeller, Mar 03 2007

Keywords

Examples

			Triangle begins as:
  1;
  1,       1;
  1,       2,      1;
  1,       9,      9,       1;
  1,      64,      8,      64,      1;
  1,     625,    625,     625,    625,       1;
  1,    7776     243,      16,    243,    7776,      1;
  1,  117649, 117649,  117649, 117649,  117649, 117649,       1;
  1, 2097152,  16384, 2097152,    128, 2097152,  16384, 2097152, 1;
		

Crossrefs

Programs

  • Mathematica
    B[n_, k_]:= If[k==0 || k==n, 1, Binomial[n, k]*k^k*(n-k)^(n-k)/n^n];
    T[n_, k_]= Denominator[B[n, k]];
    Table[T[n, k], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Jul 19 2021 *)
  • Sage
    def B(n,k): return 1 if (k==0 or k==n) else binomial(n, k)*k^k*(n-k)^(n-k)/n^n
    def T(n,k): return denominator(B(n,k))
    flatten([[T(n,k) for k in (0..n)] for n in (0..10)]) # G. C. Greubel, Jul 19 2021

Formula

A128433(n,k)/T(n,k) = binomial(n,k) * k^k * (n-k)^(n-k) / n^n.
For n>0: Sum_{k=0..n} A128433(n,k)/T(n,k) = A090878(n)/A036505(n-1);
T(n, n-k) = T(n,k).
T(n, 0) = T(n, n) = 1.
for n>0: A128433(n,1)/T(n,1) = A000312(n-1)/A000169(n).

A302583 a(n) = ((n + 1)^n - (n - 1)^n)/2.

Original entry on oeis.org

0, 1, 4, 28, 272, 3376, 51012, 908608, 18640960, 432891136, 11225320100, 321504185344, 10079828372880, 343360783937536, 12627774819845668, 498676704524517376, 21046391759976988928, 945381827279671853056, 45032132922921758270916, 2267322327322331161821184
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 10 2018

Keywords

Crossrefs

Programs

  • Mathematica
    Table[((n + 1)^n - (n - 1)^n)/2, {n, 0, 19}]
    nmax = 19; CoefficientList[Series[(x^2 - LambertW[-x]^2)/(2 x LambertW[-x] (1 + LambertW[-x])), {x, 0, nmax}], x] Range[0, nmax]!
    Table[n! SeriesCoefficient[Exp[n x] Sinh[x], {x, 0, n}], {n, 0, 19}]

Formula

E.g.f.: (x^2 - LambertW(-x)^2)/(2*x*LambertW(-x)*(1 + LambertW(-x))).
a(n) = n! * [x^n] exp(n*x)*sinh(x).
Previous Showing 101-110 of 382 results. Next