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 41-50 of 361 results. Next

A006253 Number of perfect matchings (or domino tilings) in C_4 X P_n.

Original entry on oeis.org

1, 2, 9, 32, 121, 450, 1681, 6272, 23409, 87362, 326041, 1216800, 4541161, 16947842, 63250209, 236052992, 880961761, 3287794050, 12270214441, 45793063712, 170902040409, 637815097922, 2380358351281, 8883618307200, 33154114877521, 123732841202882
Offset: 0

Views

Author

Keywords

Comments

Number of tilings of a box with sides 2 X 2 X n in R^3 by boxes of sides 2 X 1 X 1 (3-dimensional dominoes). - Frans J. Faase
The number of domino tilings in A006253, A004003, A006125 is the number of perfect matchings in the relevant graphs. There are results of Jockusch and Ciucu that if a planar graph has a rotational symmetry then the number of perfect matchings is a square or twice a square - this applies to these 3 sequences. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 12 2001
Also stacking bricks.
a(n)*(-1)^n = (1-T(n+1,-2))/3, n>=0, with Chebyshev's polynomials T(n,x) of the first kind, is the r=-2 member of the r-family of sequences S_r(n) defined in A092184 where more information can be found. - Wolfdieter Lang, Oct 18 2004
Partial sums of A217233. - Bruno Berselli, Oct 01 2012
The sequence is the case P1 = 2, P2 = -8, Q = 1 of the 3 parameter family of 4th-order linear divisibility sequences found by Williams and Guy. - Peter Bala, Apr 03 2014

Examples

			G.f. = 1 + 2*x + 9*x^2 + 32*x^3 + 121*x^4 + 450*x^5 + ... - _Michael Somos_, Mar 17 2022
		

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 360.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A002530, A004003, A006125, A217233 (first differences), A109437 (partial sums).
Column k=2 of A181206, A189650, A233308.
Cf. A100047.

Programs

  • GAP
    a:=[1,2,9];; for n in [4..30] do a[n]:=3*a[n-1]+3*a[n-2]-a[n-3]; od; a; # G. C. Greubel, Nov 16 2018
  • Magma
    m:=30; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((1-x)/(1-3*x-3*x^2+x^3))); // G. C. Greubel, Nov 16 2018
    
  • Mathematica
    CoefficientList[Series[(1 - x)/(1 - 3 x - 3 x^2 + x^3), {x, 0, 30}], x] (* Vincenzo Librandi, Oct 15 2012 *)
    RecurrenceTable[{a[1] == 1, a[2] == 2, a[n] == BitXor[1, a[n - 1]]^2/a[n - 2]}, a, {n, 30}] (* Jon Maiga, Nov 16 2018 *)
    LinearRecurrence[{3,3,-1}, {1,2,9}, 30] (* G. C. Greubel, Nov 16 2018 *)
    a[ n_] := (-1)^n * ChebyshevU[n, Sqrt[-1/2]]^2; (* Michael Somos, Mar 17 2022 *)
  • PARI
    a(n)=(sqrt(3)+2)^(n+1) \/ 6 \\ Charles R Greathouse IV, Aug 18 2016
    
  • PARI
    a(n)=([0,1,0; 0,0,1; -1,3,3]^n*[1;2;9])[1,1] \\ Charles R Greathouse IV, Aug 18 2016
    
  • PARI
    Vec((1 - x) / ((1 + x)*(1 - 4*x + x^2)) + O(x^40)) \\ Colin Barker, Dec 16 2017
    
  • PARI
    {a(n) = simplify((-1)^n * polchebyshev(n, 2, quadgen(-8)/2)^2)}; /* Michael Somos, Mar 17 2022 */
    
  • Sage
    s=((1-x)/(1-3*x-3*x^2+x^3)).series(x,30); s.coefficients(x, sparse=False) # G. C. Greubel, Nov 16 2018
    

Formula

G.f.: (1-x)/((1+x)*(1-4*x+x^2)) = (1-x)/(1-3*x-3*x^2+x^3). - Simon Plouffe in his 1992 dissertation; typo corrected by Vincenzo Librandi, Oct 15 2012
Nearest integer to (1/6)*(2+sqrt(3))^(n+1). - Don Knuth, Jul 15 1995
For n >= 4, a(n) = 3a(n-1) + 3a(n-2) - a(n-3). - Avi Peretz (njk(AT)netvision.net.il), Mar 30 2001
For n >= 3, a(n) = 4a(n-1) - a(n-2) + 2*(-1)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 14 2001
From Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 11 2001: The values are a(1) = 2 * 1^2, a(2) = 3^2, a(3) = 2 * 4^2, a(4) = 11^2, a(5) = 2 * 15^2, ... and in general for odd n a(n) is twice a square, for even n a(n) is a square. If we define b(n) by b(n) = sqrt(a(n)) for even n, b(n) = sqrt(a(n)/2) for odd n then apart from the first 2 elements b(n) is A002530(n+1).
a(n) + a(n+1) = A001835(n+2). - R. J. Mathar, Dec 06 2013
From Peter Bala, Apr 03 2014: (Start)
a(n) = |U(n,i/sqrt(2))|^2 where U(n,x) denotes the Chebyshev polynomial of the second kind.
a(n-1) = the bottom left entry of the 2 X 2 matrix T(n, M), where M is the 2 X 2 matrix [0, 2; 1, 1] and T(n,x) denotes the Chebyshev polynomial of the first kind.
See the remarks in A100047 for the general connection between Chebyshev polynomials of the first kind and 4th-order linear divisibility sequences. (End)
a(n) = (2*(-1)^n + (2-sqrt(3))^(1+n) + (2+sqrt(3))^(1+n)) / 6. - Colin Barker, Dec 16 2017
a(n) = (1 XOR a(n-1))^2/a(n-2). - Jon Maiga, Nov 16 2018
a(n) = a(-2-n) for all n in Z. - Michael Somos, Mar 17 2022
INVERT transform of sequence p(n), n > 0, where p is the number of nonreducible tilings by height of 2 X 2 X n using dicubes; p is (2, 5, 4, 4, 4, 4...). - Nicolas Bělohoubek, Jun 04 2024

A005703 Number of n-node connected graphs with at most one cycle.

Original entry on oeis.org

1, 1, 1, 2, 4, 8, 19, 44, 112, 287, 763, 2041, 5577, 15300, 42419, 118122, 330785, 929469, 2621272, 7411706, 21010378, 59682057, 169859257, 484234165, 1382567947, 3952860475, 11315775161, 32430737380, 93044797486, 267211342954, 768096496093, 2209772802169
Offset: 0

Views

Author

Keywords

Comments

a(n) is the number of pseudotrees on n nodes. - Eric W. Weisstein, Jun 11 2012
Also unlabeled connected graphs covering n vertices with at most n edges. For this definition we have a(1) = 0 and possibly a(0) = 0. - Gus Wiseman, Feb 20 2024

Examples

			From _Gus Wiseman_, Feb 20 2024: (Start)
Representatives of the a(0) = 1 through a(5) = 8 graphs:
  {}  .  {12}  {12,13}     {12,13,14}     {12,13,14,15}
               {12,13,23}  {12,13,24}     {12,13,14,25}
                           {12,13,14,23}  {12,13,24,35}
                           {12,13,24,34}  {12,13,14,15,23}
                                          {12,13,14,23,25}
                                          {12,13,14,23,45}
                                          {12,13,14,25,35}
                                          {12,13,24,35,45}
(End)
		

References

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

Crossrefs

Cf. A000055, A000081, A001429 (labeled A057500), A134964 (number of pseudoforests, labeled A133686).
The labeled version is A129271.
The connected complement is A140636, labeled A140638.
Non-connected: A368834 (labeled A367869) or A370316 (labeled A369191).
A001187 counts connected graphs, unlabeled A001349.
A006125 counts simple graphs, unlabeled A000088.
A006129 counts covering graphs, unlabeled A002494.
A062734 counts connected graphs by number of edges.

Programs

  • Mathematica
    Needs["Combinatorica`"]; nn = 20; t[x_] := Sum[a[n] x^n, {n, 1, nn}];
    a[0] = 0;
    b = Drop[Flatten[
        sol = SolveAlways[
          0 == Series[
            t[x] - x Product[1/(1 - x^i)^a[i], {i, 1, nn}], {x, 0, nn}],
          x]; Table[a[n], {n, 0, nn}] /. sol], 1];
    r[x_] := Sum[b[[n]] x^n, {n, 1, nn}]; c =
    Drop[Table[
        CoefficientList[
         Series[CycleIndex[DihedralGroup[n], s] /.
           Table[s[i] -> r[x^i], {i, 1, n}], {x, 0, nn}], x], {n, 3,
         nn}] // Total, 1];
    d[x_] := Sum[c[[n]] x^n, {n, 1, nn}]; CoefficientList[
    Series[r[x] - (r[x]^2 - r[x^2])/2 + d[x] + 1, {x, 0, nn}], x] (* Geoffrey Critzer, Nov 17 2014 *)
  • PARI
    \\ TreeGf gives gf of A000081.
    TreeGf(N)={my(A=vector(N, j, 1)); for (n=1, N-1, A[n+1] = 1/n * sum(k=1, n, sumdiv(k, d, d*A[d]) * A[n-k+1] ) ); x*Ser(A)}
    seq(n)={my(t=TreeGf(n)); my(g(e)=subst(t + O(x*x^(n\e)), x, x^e) + O(x*x^n)); Vec(1 + g(1) + (g(2) - g(1)^2)/2 + sum(k=3, n, sumdiv(k, d, eulerphi(d)*g(d)^(k/d))/k + if(k%2, g(1)*g(2)^(k\2), (g(1)^2+g(2))*g(2)^(k/2-1)/2))/2)}; \\ Andrew Howroyd and Washington Bomfim, May 15 2021

Formula

a(n) = A000055(n) + A001429(n).

Extensions

More terms from Vladeta Jovovic, Apr 19 2000 and from Michael Somos, Apr 26 2000
a(27) corrected and a(28) and a(29) computed by Washington Bomfim, May 14 2008

A140638 Number of connected graphs on n labeled nodes that contain at least two cycles.

Original entry on oeis.org

0, 0, 0, 7, 381, 21748, 1781154, 249849880, 66257728763, 34495508486976, 35641629989151608, 73354595357480683904, 301272202621204113362497, 2471648811029413368450098688, 40527680937730440155535277704046, 1328578958335783199341353852258282496
Offset: 1

Views

Author

Washington Bomfim, May 21 2008

Keywords

Comments

These are the connected graphs that are neither trees nor unicyclic.
Also connected non-choosable graphs covering n vertices, where a graph is choosable iff it is possible to choose a different vertex from each edge. The unlabeled version is A140636. The complement is counted by A129271. - Gus Wiseman, Feb 20 2024

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Dover, 2002, p. 2.

Crossrefs

The unlabeled version is A140636.
Cf. A000272 (trees), A001187 (connected graphs), A057500 (connected unicyclic graphs).
The complement is counted by A129271, unlabeled A005703.
The non-connected complement is A133686, covering A367869.
The non-connected version is A367867, unlabeled A140637.
The non-connected covering version is A367868.
A006125 counts graphs, A000088 unlabeled.
A006129 counts covering graphs, A002494 unlabeled.
A143543 counts simple labeled graphs by number of connected components.

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Subsets[Range[Length[s]],{2}],Length[Intersection@@s[[#]]]>0&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],Union@@#==Range[n]&&Length[csm[#]]<=1&&Select[Tuples[#],UnsameQ@@#&]=={}&]],{n,0,5}] (* Gus Wiseman, Feb 19 2024 *)
  • PARI
    seq(n)={my(A=O(x*x^n), t=-lambertw(-x + A)); Vec(serlaplace( log(sum(k=0, n, 2^binomial(k, 2)*x^k/k!, A)) - log(1/(1-t))/2 - t/2 + 3*t^2/4), -n)} \\ Andrew Howroyd, Jan 15 2022

Formula

a(n) = A001187(n) - A129271(n).
a(n) = A001187(n) - A000272(n) - A057500(n).

Extensions

Definition clarified by Andrew Howroyd, Jan 15 2022

A143543 Triangle read by rows: T(n,k) = number of labeled graphs on n nodes with k connected components, 1<=k<=n.

Original entry on oeis.org

1, 1, 1, 4, 3, 1, 38, 19, 6, 1, 728, 230, 55, 10, 1, 26704, 5098, 825, 125, 15, 1, 1866256, 207536, 20818, 2275, 245, 21, 1, 251548592, 15891372, 925036, 64673, 5320, 434, 28, 1, 66296291072, 2343580752, 76321756, 3102204, 169113, 11088, 714, 36, 1
Offset: 1

Views

Author

Max Alekseyev, Aug 23 2008

Keywords

Comments

The Bell transform of A001187(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 17 2016

Examples

			The triangle T(n,k) starts as:
n=1:     1;
n=2:     1,    1;
n=3:     4,    3,   1;
n=4:    38,   19,   6,   1;
n=5:   728,  230,  55,  10,  1;
n=6: 26704, 5098, 825, 125, 15, 1;
...
		

Crossrefs

Cf. A001187 (first column), A006125 (row sums), A106240 (unlabeled variant).
Cf. A125207.
T(2n,n) gives A369827.

Programs

  • Maple
    g:= proc(n) option remember; `if`(n=0, 1, 2^(n*(n-1)/2)-add(
          binomial(n, k)*2^((n-k)*(n-k-1)/2)*g(k)*k, k=1..n-1)/n)
        end:
    b:= proc(n) option remember; `if`(n=0, 1, add(expand(
          b(n-j)*binomial(n-1, j-1)*g(j)*x), j=1..n))
        end:
    T:= (n, k)-> coeff(b(n$2), x, k):
    seq(seq(T(n, k), k=1..n), n=1..10);  # Alois P. Heinz, Feb 02 2024
  • Mathematica
    a= Sum[2^Binomial[n,2] x^n/n!,{n,0,10}];
    Rest[Transpose[Table[Range[0, 10]! CoefficientList[Series[Log[a]^n/n!, {x, 0, 10}], x], {n, 1, 10}]]] // Grid (* Geoffrey Critzer, Mar 15 2011 *)
  • PARI
    T(n)={[Vecrev(p/y) | p <- Vec(serlaplace(exp(y*log(sum(k=0, n, 2^binomial(k,2)*x^k/k!, O(x*x^n))))))]}
    { foreach(T(8), row, print(row)) } \\ Andrew Howroyd, Jun 14 2025
  • Sage
    # uses[bell_matrix from A264428, A001187]
    # Adds a column 1,0,0,0, ... at the left side of the triangle.
    bell_matrix(lambda n: A001187(n+1), 9) # Peter Luschny, Jan 17 2016
    

Formula

SUM[n,k=0..oo] T(n,k) * x^n * y^k / n! = exp( y*( F(x) - 1 ) ) = ( SUM[n=0..oo] 2^binomial(n, 2)*x^n/n! )^y, where F(x) is e.g.f. of A001187.
T(n,k) = Sum_{q=0..n-1} C(n-1, q) T(q, k-1) 2^C(n-q,2) - Sum_{q=0..n-2} C(n-1, q) T(q+1, k) 2^C(n-1-q, 2) where T(0,0) = 1 and T(0,k) = 0 and T(n,0) = 0. - Marko Riedel, Feb 04 2019
Sum_{k=1..n} k * T(n,k) = A125207(n) - Alois P. Heinz, Feb 02 2024

A088957 Hyperbinomial transform of the sequence of 1's.

Original entry on oeis.org

1, 2, 6, 29, 212, 2117, 26830, 412015, 7433032, 154076201, 3608522954, 94238893883, 2715385121740, 85574061070045, 2928110179818478, 108110945014584623, 4284188833355367440, 181370804507130015569, 8169524599872649117330, 390114757072969964280163
Offset: 0

Views

Author

Paul D. Hanna, Oct 26 2003

Keywords

Comments

See A088956 for the definition of the hyperbinomial transform.
a(n) is the number of partial functions on {1,2,...,n} that are endofunctions with no cycles of length > 1. The triangle A088956 classifies these functions according to the number of undefined elements in the domain. The triangle A144289 classifies these functions according to the number of edges in their digraph representation (considering the empty function to have 1 edge). The triangle A203092 classifies these functions according to the number of connected components. - Geoffrey Critzer, Dec 29 2011
a(n) is the number of rooted subtrees (for a fixed root) in the complete graph on n+1 vertices: a(3) = 29 is the number of rooted subtrees in K_4: 1 of size 1, 3 of size 2, 9 of size 3, and 16 spanning subtrees. - Alex Chin, Jul 25 2013 [corrected by Marko Riedel, Mar 31 2019]
From Gus Wiseman, Jan 28 2024: (Start)
Also the number of labeled loop-graphs on n vertices such that it is possible to choose a different vertex from each edge in exactly one way. For example, the a(3) = 29 uniquely choosable loop-graphs (loops shown as singletons) are:
{} {1} {1,2} {1,12} {1,2,13} {1,12,13}
{2} {1,3} {1,13} {1,2,23} {1,12,23}
{3} {2,3} {2,12} {1,3,12} {1,13,23}
{2,23} {1,3,23} {2,12,13}
{3,13} {2,3,12} {2,12,23}
{3,23} {2,3,13} {2,13,23}
{1,2,3} {3,12,13}
{3,12,23}
{3,13,23}
(End)

Examples

			a(5) = 2117 = 1296 + 625 + 160 + 30 + 5 + 1 = sum of row 5 of triangle A088956.
		

Crossrefs

Cf. A088956 (triangle).
Row sums of A144289. - Alois P. Heinz, Jun 01 2009
Column k=1 of A144303. - Alois P. Heinz, Oct 30 2012
The covering case is A000272, also the case of exactly n edges.
Without the choice condition we have A006125 (shifted left).
The unlabeled version is A087803.
The choosable version is A368927, covering A369140, loopless A133686.
The non-choosable version is A369141, covering A369142, loopless A367867.

Programs

  • Haskell
    a088957 = sum . a088956_row  -- Reinhard Zumkeller, Jul 07 2013
    
  • Maple
    a:= n-> add((n-j+1)^(n-j-1)*binomial(n,j), j=0..n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Oct 30 2012
  • Mathematica
    nn = 16; t = Sum[n^(n - 1) x^n/n!, {n, 1, nn}];
    Range[0, nn]! CoefficientList[Series[Exp[x] Exp[t], {x, 0, nn}], x]  (* Geoffrey Critzer, Dec 29 2011 *)
    With[{nmax = 50}, CoefficientList[Series[-LambertW[-x]*Exp[x]/x, {x, 0, nmax}], x]*Range[0, nmax]!] (* G. C. Greubel, Nov 14 2017 *)
  • PARI
    x='x+O('x^10); Vec(serlaplace(-lambertw(-x)*exp(x)/x)) \\ G. C. Greubel, Nov 14 2017

Formula

a(n) = Sum_{k=0..n} (n-k+1)^(n-k-1)*C(n, k).
E.g.f.: A(x) = exp(x+sum(n>=1, n^(n-1)*x^n/n!)).
E.g.f.: -LambertW(-x)*exp(x)/x. - Vladeta Jovovic, Oct 27 2003
a(n) ~ exp(1+exp(-1))*n^(n-1). - Vaclav Kotesovec, Jul 08 2013
Binomial transform of A000272. - Gus Wiseman, Jan 25 2024

A327071 Number of labeled simple connected graphs with n vertices and at least one bridge, or graphs with spanning edge-connectivity 1.

Original entry on oeis.org

0, 0, 1, 3, 28, 475, 14736, 818643, 82367552, 15278576679, 5316021393280, 3519977478407687, 4487518206535452672, 11116767463976825779115, 53887635281876408097483776, 513758302006787897939587736715, 9668884580476067306398361085853696
Offset: 0

Views

Author

Gus Wiseman, Aug 24 2019

Keywords

Comments

A bridge is an edge that, if removed without removing any incident vertices, disconnects the graph. Connected graphs with no bridges are counted by A095983 (2-edge-connected graphs).
The spanning edge-connectivity of a graph is the minimum number of edges that must be removed (without removing incident vertices) to obtain a disconnected or empty graph.

Crossrefs

Column k = 1 of A327069.
The unlabeled version is A052446.
Connected graphs without bridges are A007146.
The enumeration of labeled connected graphs by number of bridges is A327072.
Connected graphs with exactly one bridge are A327073.
Graphs with non-spanning edge-connectivity 1 are A327079.
BII-numbers of set-systems with spanning edge-connectivity 1 are A327111.
Covering set-systems with spanning edge-connectivity 1 are A327145.
Graphs with spanning edge-connectivity 2 are A327146.

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    spanEdgeConn[vts_,eds_]:=Length[eds]-Max@@Length/@Select[Subsets[eds],Union@@#!=vts||Length[csm[#]]!=1&];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],spanEdgeConn[Range[n],#]==1&]],{n,0,4}]

Formula

a(1) = 0; a(n > 1) = A001187(n) - A095983(n).

A000568 Number of outcomes of unlabeled n-team round-robin tournaments.

Original entry on oeis.org

1, 1, 1, 2, 4, 12, 56, 456, 6880, 191536, 9733056, 903753248, 154108311168, 48542114686912, 28401423719122304, 31021002160355166848, 63530415842308265100288, 244912778438520759443245824, 1783398846284777975419600287232, 24605641171260376770598003978281472
Offset: 0

Views

Author

Keywords

Comments

Harary and Palmer give incorrect values for a(24) and a(25); the correct values are a(24) = 195692027657521876084316842660833482785173437775365039898624 and a(25) = 131326696677895002131450257709457767457170027052967027982788816896. - Vladeta Jovovic, Apr 08 2001
a(n) appears to be the number of even graphs with n vertices; see comment in A334335. - Pontus von Brömssen, May 05 2020 [This has been proved by Royle et al. 2023. - Pontus von Brömssen, Apr 06 2022]

References

  • R. L. Davis, Structure of dominance relations, Bull. Math. Biophys., 16 (1954), 131-140.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 157 and 523.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, pp. 126 and 245.
  • J. W. Moon, Topics on Tournaments. Holt, NY, 1968, p. 87.
  • K. B. Reid and L. W. Beineke "Tournaments", pp. 169-204 in L. W. Beineke and R. J. Wilson, editors, Selected Topics in Graph Theory, Academic Press, NY, 1978.
  • 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

Cf. A006125 for the labeled analog, A051337.
Euler transform of A334335.

Programs

  • Maple
    with(combinat):with(numtheory): for n from 1 to 30 do p:=partition(n): s:=0:for k from 1 to nops(p) do ex:=1:for i from 1 to nops(p[k]) do if p[k][i] mod 2=0 then ex:=0:break:fi:od:
    if ex=1 then q:=convert(p[k],multiset): for i from 1 to n do a(i):=0:od:for i from 1 to nops(q) do a(q[i][1]):=q[i][2]:od:
    c:=1:ord:=1:for i from 1 to n do c:=c*a(i)!*i^a(i): if a(i)<>0 then ord:=lcm(ord,i):fi:od: g:=0:for d from 1 to ord do if ord mod d=0 then g1:=0:for del from 1 to n do if d mod del=0 then g1:=g1+del*a(del):fi:od:g:=g+phi(ord/d)*g1*(g1-1):fi:od: s:=s+2^(g/ord/2)/c:fi:
    od: print(n,s); od: # Vladeta Jovovic
  • Mathematica
    permcount[v_] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
    edges[v_] := Sum[Sum[GCD[v[[i]], v[[j]]], {j, 1, i-1}], {i, 2, Length[v]}] + Sum[Quotient[v[[i]], 2], {i, 1, Length[v]}];
    oddp[v_] := (For[i = 1, i <= Length[v], i++, If[BitAnd[v[[i]], 1] == 0, Return[0]]]; 1);
    a[n_] := a[n] = (s = 0; Do[If[oddp[p] == 1, s += permcount[p]*2^edges[p]], {p, IntegerPartitions[n]}]; s/n!);
    Table[Print["a(", n, ") = ", a[n]]; a[n], {n, 0, 25}] (* Jean-François Alcover, Nov 13 2017, after Andrew Howroyd *)
  • PARI
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
    edges(v) = {sum(i=2, #v, sum(j=1, i-1, gcd(v[i],v[j]))) + sum(i=1, #v, v[i]\2)}
    oddp(v) = {for(i=1, #v, if(bitand(v[i],1)==0, return(0)));1}
    a(n) = {my(s=0); forpart(p=n, if(oddp(p), s+=permcount(p)*2^edges(p))); s/n!} \\ Andrew Howroyd, Oct 22 2017
    
  • Python
    from itertools import product
    from math import prod, factorial, gcd
    from fractions import Fraction
    from sympy.utilities.iterables import partitions
    def A000568(n): return int(sum(Fraction(1<<(sum(p[r]*p[s]*gcd(r,s) for r,s in product(p.keys(),repeat=2))-sum(p.values())>>1),prod(q**p[q]*factorial(p[q]) for q in p)) for p in partitions(n) if all(q&1 for q in p))) # Chai Wah Wu, Jul 01 2024

Formula

Davis's formula: a(n) = Sum_{j} (1/(Product (k^(j_k) (j_k)!))) * 2^{t_j},
where j runs through all partitions of n into odd parts, say with j_1 parts of size 1, j_3 parts of size 3, etc.,
and t_j = (1/2)*[ Sum_{r=1..n, s=1..n} j_r j_s gcd(r,s) - Sum_{r} j_r ].

Extensions

More terms from Vladeta Jovovic

A004251 Number of graphical partitions (degree-vectors for simple graphs with n vertices, or possible ordered row-sum vectors for a symmetric 0-1 matrix with diagonal values 0).

Original entry on oeis.org

1, 1, 2, 4, 11, 31, 102, 342, 1213, 4361, 16016, 59348, 222117, 836315, 3166852, 12042620, 45967479, 176005709, 675759564, 2600672458, 10029832754, 38753710486, 149990133774, 581393603996, 2256710139346, 8770547818956, 34125389919850, 132919443189544, 518232001761434, 2022337118015338, 7898574056034636, 30873421455729728
Offset: 0

Views

Author

Keywords

Comments

In other words, a(n) is the number of graphic sequences of length n, where a graphic sequence is a sequence of numbers which can be the degree sequence of some graph.
In the article by A. Iványi, G. Gombos, L. Lucz, and T. Matuszka, "Parallel enumeration of degree sequences of simple graphs II", in Table 4 on page 260 the values a(30) = 7898574056034638 and a(31) = 30873429530206738 are incorrect due to the incorrect Gz(30) = 5876236938019300 and Gz(31) = 22974847474172100. - Wang Kai, Jun 05 2016

Examples

			For n = 3, there are 4 different graphic sequences possible: 0 0 0; 1 1 0; 2 1 1; 2 2 2. - Daan van Berkel (daan.v.berkel.1980(AT)gmail.com), Jun 25 2010
From _Gus Wiseman_, Dec 31 2020: (Start)
The a(0) = 1 through a(4) = 11 sorted degree sequences:
  ()  (0)  (0,0)  (0,0,0)  (0,0,0,0)
           (1,1)  (0,1,1)  (0,0,1,1)
                  (1,1,2)  (0,1,1,2)
                  (2,2,2)  (0,2,2,2)
                           (1,1,1,1)
                           (1,1,1,3)
                           (1,1,2,2)
                           (1,2,2,3)
                           (2,2,2,2)
                           (2,2,3,3)
                           (3,3,3,3)
For example, the graph {{2,3},{2,4}} has degrees (0,2,1,1), so (0,1,1,2) is counted under a(4).
(End)
		

References

  • R. A. Brualdi and H. J. Ryser, Combinatorial Matrix Theory, Cambridge Univ. Press, 1992.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • P. R. Stein, On the number of graphical partitions, pp. 671-684 of Proc. 9th S-E Conf. Combinatorics, Graph Theory, Computing, Congr. Numer. 21 (1978).

Crossrefs

Counting the positive partitions by sum gives A000569, ranked by A320922.
The version with half-loops is A029889, with covering case A339843.
The covering case (no zeros) is A095268.
Covering simple graphs are ranked by A309356 and A320458.
Non-graphical partitions are counted by A339617 and ranked by A339618.
The version with loops is A339844, with covering case A339845.
A006125 counts simple graphs, with covering case A006129.
A027187 counts partitions of even length, ranked by A028260.
A058696 counts partitions of even numbers, ranked by A300061.
A320921 counts connected graphical partitions.
A322353 counts factorizations into distinct semiprimes.
A339659 counts graphical partitions of 2n into k parts.
A339661 counts factorizations into distinct squarefree semiprimes.

Programs

  • Mathematica
    Table[Length[Union[Sort[Table[Count[Join@@#,i],{i,n}]]&/@Subsets[Subsets[Range[n],{2}]]]],{n,0,5}] (* Gus Wiseman, Dec 31 2020 *)

Formula

G.f. = 1 + x + 2*x^2 + 4*x^3 + 11*x^4 + 31*x^5 + 102*x^6 + 342*x^7 + 1213*x^8 + ...
a(n) ~ c * 4^n / n^(3/4) for some constant c > 0. Computational estimates suggest c ≈ 0.099094. - Tom Johnston, Jan 18 2023

Extensions

More terms from Torsten Sillke, torsten.sillke(AT)lhsystems.com, using Cor. 6.3.3, Th. 6.3.6, Cor. 6.2.5 of Brualdi-Ryser.
a(19) from Herman Jamke (hermanjamke(AT)fastmail.fm), May 19 2007
a(20)-a(23) from Nathann Cohen, Jul 09 2011
a(24)-a(29) from Antal Iványi, Nov 15 2011
a(30) and a(31) corrected by Wang Kai, Jun 05 2016

A071356 Expansion of (1 - 2*x - sqrt(1 - 4*x - 4*x^2))/(4*x^2).

Original entry on oeis.org

1, 2, 6, 20, 72, 272, 1064, 4272, 17504, 72896, 307648, 1312896, 5655808, 24562176, 107419264, 472675072, 2091206144, 9296612352, 41507566592, 186045061120, 836830457856, 3776131489792, 17089399689216, 77548125675520, 352766964908032
Offset: 0

Views

Author

N. J. A. Sloane, Jun 12 2002

Keywords

Comments

Number of underdiagonal lattice paths from (0,0) to the line x=n, using only steps R=(1,0), V=(0,1) and D=(1,2). Also number of Motzkin paths of length n in which both the "up" and the "level" steps come in two colors. E.g., a(2)=6 because we have RR, RVR, RRV, RD, RVRV and RRVV. - Emeric Deutsch, Dec 21 2003
Inverse binomial transform of little Schroeder numbers 1,3,11,... (A001003 with first term deleted). - David Callan, Feb 07 2004
a(n) is the number of planar trees satisfying: 1) Every internal node has at least two children, 2) Among the children of a node, only the leftmost and the rightmost children can be leaves, 3) The tree has n+1 leaves. For instance, a(3)=6. - Marcelo Aguiar (maguiar(AT)math.tamu.edu), Oct 14 2005
Hankel transform is A006125(n+1)=2^C(n+1,2). - Paul Barry, Jan 08 2008
Equals binomial transform of A025235: (1, 1, 3, 7, 21, 61, 191, ...). - Gary W. Adamson, Sep 03 2010
Conjecturally, the number of sequences (e(1), ..., e(n+1)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i) > e(j) <= e(k). [Martinez and Savage, 2.19] - Eric M. Schmidt, Jul 17 2017
Let s denote West's stack-sorting map, and let Av_n(tau_1, ..., tau_r) denote the set of permutations of [n] that avoid the patterns tau_1, ..., tau_r. It is conjectured that a(n) = |s^{-1}(Av_{n+1}(132, 231))| = |s^{-1}(Av_{n+1}(132, 312))| = |s^{-1}(Av_{n+1}(231, 312))|. Only the last of these equalities is known. - Colin Defant, Sep 16 2018

Examples

			a(3) = 20 = sum of top row terms in M^3 = (9 + 7 + 3 + 1).
		

Crossrefs

A036774(n) = a(n-1) * n! / 2^(n-1).
Row sums of A071943.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!((1 - 2*x - Sqrt(1 - 4*x - 4*x^2))/(4*x^2))); // Vincenzo Librandi, Jan 21 2020
  • Mathematica
    CoefficientList[Series[(1-2*x-Sqrt[1-4*x-4*x^2])/(4*x^2), {x, 0, 20}], x] (* Vaclav Kotesovec, Sep 24 2013 *)
    a[n_] := 2^n Hypergeometric2F1[(1-n)/2, -n/2, 2, 2];
    Table[a[n], {n, 0, 24}] (* Peter Luschny, May 30 2021 *)
  • PARI
    a(n)=if(n<0,0,n++; polcoeff(serreverse(x/(1+2*x+2*x^2)+x*O(x^n)),n))
    
  • PARI
    {a(n)= if(n<1, n==0, polcoeff( 2/(1 -2*x +sqrt(1 -4*x -4*x^2 +x*O(x^n))), n))}
    
  • PARI
    {a(n)= local(A); if(n<0, 0, A= x*O(x^n); n!*simplify(polcoeff( exp(2*x +A)* besseli(1, 2*x* quadgen(8) +A), n)))} /* Michael Somos, Mar 31 2007 */
    
  • Sage
    def A071356_list(n):  # n>=1
        T = [0]*(n+1); R = [1]
        for m in (1..n-1):
            a,b,c = 1,0,0
            for k in range(m,-1,-1):
                r = a + 2*(b + c)
                if k < m : T[k+2] = u;
                a,b,c = T[k-1],a,b
                u = r
            T[1] = u; R.append(u)
        return R
    A071356_list(25)  # Peter Luschny, Nov 01 2012
    

Formula

G.f. A(x) satisfies 2x^2*A(x)^2+(2x-1)*A(x)+1=0 and A(x)=1/(1-2x-2x^2/A(x)). - Michael Somos, Sep 06 2003
a(n) = Sum_{k=0..floor(n/2)} C(n, 2k)C(k)2^(n-2k)*2^k. - Paul Barry, May 18 2005
G.f.: (1 - 2*x - sqrt(1 - 4*x - 4*x^2) )/(4*x^2) = 2/(1 - 2*x +sqrt(1 - 4*x - 4*x^2)).
Moment representation is a(n) = (1/(4*Pi))*int(x^n*sqrt(4-4x-x^2), x, -2*sqrt(2)-2, 2*sqrt(2)-2). - Paul Barry, Jan 08 2008
G.f.: 1/(1-2x-2x^2/(1-2x-2x^2/(1-2x-2x^2/(1-2x-2x^2/(1-2x-2x^2/.... (continued fraction). - Paul Barry, Dec 06 2008
From Gary W. Adamson, Jul 22 2011: (Start)
a(n) = sum of top row terms of M^n, M = an infinite square production matrix as follows:
1, 1, 0, 0, 0, 0, ...
2, 1, 1, 0, 0, 0, ...
2, 2, 1, 1, 0, 0, ...
2, 2, 2, 1, 1, 0, ...
2, 2, 2, 2, 1, 1, ...
2, 2, 2, 2, 2, 1, ... (End)
E.g.f.: a(n) = n!* [x^n] exp(2*x)*BesselI(1, 2*sqrt(2)*x)/(sqrt(2)*x). - Peter Luschny, Aug 25 2012
D-finite with recurrence: (n+2)*a(n) +2*(-2*n-1)*a(n-1) +4*(-n+1)*a(n-2)=0. - R. J. Mathar, Dec 02 2012 (Formula verified and used for computations. - Fung Lam, Feb 24 2014)
a(n) ~ 2^(n - 1/4) * (1+sqrt(2))^(n + 3/2) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Sep 24 2013, simplified Jan 26 2019
a(n) = A179190(n+2)/4. - R. J. Mathar, Jan 20 2020
a(n) = 2^n * hypergeom((1 - n)/2, -n/2, 2, 2). - Peter Luschny, May 30 2021
a(n) = (-2*î)^(n+2) * (Legendre_P(n+2, i) - Legendre_P(n, i))/(4*(2*n + 3)). - Peter Bala, May 06 2024
From Emanuele Munarini, Jun 13 2024: (Start)
a(n) = Sum_{k=0..floor(n/2)} binomial(n, k)*binomial(n-k, k)*2^(n-k)/(k+1).
a(n) = Sum_{k=0..floor((n+2)/3)} binomial(n-2k+2, 2k)*Catalan(n-2k+1).
a(n) = Sum_{k=0..floor((n+2)/4)} binomial(n-2k+1, 2k+1)*Catalan(n-2k). (End)

A327069 Triangle read by rows where T(n,k) is the number of labeled simple graphs with n vertices and spanning edge-connectivity k.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 4, 3, 1, 0, 26, 28, 9, 1, 0, 296, 475, 227, 25, 1, 0, 6064, 14736, 10110, 1782, 75, 1, 0
Offset: 0

Views

Author

Gus Wiseman, Aug 23 2019

Keywords

Comments

The spanning edge-connectivity of a graph is the minimum number of edges that must be removed (without removing incident vertices) to obtain a disconnected or empty graph.
We consider a graph with one vertex and no edges to be disconnected.

Examples

			Triangle begins:
    1
    1   0
    1   1   0
    4   3   1   0
   26  28   9   1   0
  296 475 227  25   1   0
		

Crossrefs

Row sums are A006125.
Column k = 0 is A054592, if we assume A054592(1) = 1.
Column k = 1 is A327071.
Column k = 2 is A327146.
The unlabeled version (except with offset 1) is A263296.

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    spanEdgeConn[vts_,eds_]:=Length[eds]-Max@@Length/@Select[Subsets[eds],Union@@#!=vts||Length[csm[#]]!=1&];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],spanEdgeConn[Range[n],#]==k&]],{n,0,5},{k,0,n}]

Extensions

a(21)-a(27) from Robert Price, May 25 2021
Previous Showing 41-50 of 361 results. Next