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-10 of 14 results. Next

A191646 Triangle read by rows: T(n,k) = number of connected multigraphs with n >= 0 edges and 1 <= k <= n+1 vertices, with no loops allowed.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 5, 3, 0, 1, 4, 11, 11, 6, 0, 1, 6, 22, 34, 29, 11, 0, 1, 7, 37, 85, 110, 70, 23, 0, 1, 9, 61, 193, 348, 339, 185, 47, 0, 1, 11, 95, 396, 969, 1318, 1067, 479, 106, 0, 1, 13, 141, 771, 2445, 4457, 4940, 3294, 1279, 235
Offset: 0

Views

Author

Alberto Tacchella, Jul 04 2011

Keywords

Examples

			Triangle T(n,k) (with rows n >= 0 and columns k >= 1) begins as follows:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 2,  2;
  0, 1, 3,  5,  3;
  0, 1, 4, 11, 11,  6;
  0, 1, 6, 22, 34, 29, 11;
  ...
		

Crossrefs

Row sums give A076864. Diagonal is A000055.
Cf. A034253, A054923, A192517, A253186 (column k=3), A290778 (column k=4).

Programs

  • PARI
    EulerT(v)={my(p=exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1); Vec(p/x,-#v)}
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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,x)={sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i],v[j])); g*x^(v[i]*v[j]/g))) + sum(i=1, #v, my(t=v[i]); ((t-1)\2)*x^t + if(t%2,0,x^(t/2)))}
    G(n,m)={my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(edges(p,x) + O(x*x^m), -m))); s/n!}
    R(n)={Mat(apply(p->Col(p+O(y^n),-n), InvEulerMT(vector(n, k, 1 + y*Ser(G(k,n-1), y)))))}
    { my(A=R(10)); for(n=1, #A, for(k=1, n, print1(A[n,k], ", "));print) } \\ Andrew Howroyd, May 14 2018

Formula

T(n,k=3) = A253186(n) = A034253(n,k=2) for n >= 1. - Petros Hadjicostas, Oct 02 2019

A054923 Triangle read by rows: number of connected graphs with k >= 0 edges and n nodes (1<=n<=k+1).

Original entry on oeis.org

1, 0, 1, 0, 0, 1, 0, 0, 1, 2, 0, 0, 0, 2, 3, 0, 0, 0, 1, 5, 6, 0, 0, 0, 1, 5, 13, 11, 0, 0, 0, 0, 4, 19, 33, 23, 0, 0, 0, 0, 2, 22, 67, 89, 47, 0, 0, 0, 0, 1, 20, 107, 236, 240, 106, 0, 0, 0, 0, 1, 14, 132, 486, 797, 657, 235, 0, 0, 0, 0, 0, 9, 138, 814, 2075, 2678, 1806, 551, 0, 0, 0, 0, 0, 5, 126, 1169, 4495, 8548, 8833, 5026, 1301
Offset: 0

Views

Author

Keywords

Comments

The diagonal n = k+1 is A000055(n). - Jonathan Vos Post, Aug 10 2008

Examples

			Triangle begins:
  1;
  0, 1;
  0, 0, 1;
  0, 0, 1, 2;
  0, 0, 0, 2, 3;
  0, 0, 0, 1, 5   6;
  0, 0, 0, 1, 5, 13,  11;
  0, 0, 0, 0, 4, 19,  33,  23;
  0, 0, 0, 0, 2, 22,  67,  89,  47;
  0, 0, 0, 0, 1, 20, 107, 236, 240, 106;
  ... (so with 5 edges there's 1 graph with 4 nodes, 5 with 5 nodes and 6 with 6 nodes). [Typo corrected by Anders Haglund, Jul 08 2008]
		

References

  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 93, Table 4.2.2; p. 241, Table A2.

Crossrefs

Main diagonal is A000055.
Subsequent diagonals give the number of connected unlabeled graphs with n nodes and n+k edges for k=0..2: A001429, A001435, A001436.
Cf. A002905 (row sums), A001349 (column sums), A008406, A046751 (transpose), A054924 (transpose), A046742 (w/o left column), A343088 (labeled).

Programs

  • PARI
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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,t) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i],v[j])); t(v[i]*v[j]/g)^g )) * prod(i=1, #v, my(c=v[i]); t(c)^((c-1)\2)*if(c%2, 1, t(c/2)))}
    G(n, x)={my(s=0); forpart(p=n, s+=permcount(p)*edges(p,i->1+x^i)); s/n!}
    T(n)={Mat([Col(p+O(y^n), -n) | p<-InvEulerMT(vector(n, k, G(k, y + O(y^n))))])}
    {my(A=T(10)); for(n=1, #A, print(A[n,1..n]))} \\ Andrew Howroyd, Oct 23 2019

Extensions

a(83)-a(89) corrected by Andrew Howroyd, Oct 24 2019

A322114 Regular triangle read by rows where T(n,k) is the number of unlabeled connected graphs with loops with n edges and k vertices, 1 <= k <= n+1.

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 0, 1, 3, 2, 0, 0, 3, 6, 3, 0, 0, 2, 11, 14, 6, 0, 0, 1, 13, 35, 33, 11, 0, 0, 0, 10, 61, 112, 81, 23, 0, 0, 0, 5, 75, 262, 347, 204, 47, 0, 0, 0, 2, 68, 463, 1059, 1085, 526, 106, 0, 0, 0, 1, 49, 625, 2458, 4091, 3348, 1376, 235
Offset: 0

Views

Author

Gus Wiseman, Nov 26 2018

Keywords

Examples

			Triangle begins:
   1
   1   1
   0   1   1
   0   1   3   2
   0   0   3   6   3
   0   0   2  11  14   6
   0   0   1  13  35  33  11
Non-isomorphic representatives of the graphs counted in row 4:
  {{2}{3}{12}{13}}   {{4}{12}{23}{34}}   {{13}{24}{35}{45}}
  {{2}{3}{13}{23}}   {{4}{13}{23}{34}}   {{14}{25}{35}{45}}
  {{3}{12}{13}{23}}  {{4}{13}{24}{34}}   {{15}{25}{35}{45}}
                     {{4}{14}{24}{34}}
                     {{12}{13}{24}{34}}
                     {{14}{23}{24}{34}}
		

Crossrefs

Row sums are A191970. Last column is A000055.

Programs

  • PARI
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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,t) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i],v[j])); t(v[i]*v[j]/g)^g )) * prod(i=1, #v, my(c=v[i]); t(c)^((c+1)\2)*if(c%2, 1, t(c/2)))}
    G(n, x)={my(s=0); forpart(p=n, s+=permcount(p)*edges(p,i->1+x^i)); s/n!}
    T(n)={Mat([Col(p+O(y^n), -n) | p<-InvEulerMT(vector(n, k, G(k, y + O(y^n))))])}
    {my(A=T(10)); for(n=1, #A, print(A[n,1..n]))} \\ Andrew Howroyd, Oct 22 2019

Extensions

Terms a(28) and beyond from Andrew Howroyd, Oct 22 2019

A076864 Number of connected loopless multigraphs with n edges.

Original entry on oeis.org

1, 1, 2, 5, 12, 33, 103, 333, 1183, 4442, 17576, 72810, 314595, 1410139, 6541959, 31322474, 154468852, 783240943, 4077445511, 21765312779, 118999764062, 665739100725, 3807640240209, 22246105114743, 132672322938379, 807126762251748
Offset: 0

Views

Author

N. J. A. Sloane, Nov 23 2002

Keywords

Comments

Inverse Euler transform of A050535.

Crossrefs

Programs

  • Mathematica
    A050535 = Cases[Import["https://oeis.org/A050535/b050535.txt", "Table"], {, }][[All, 2]];
    (* EulerInvTransform is defined in A022562 *)
    Join[{1}, EulerInvTransform[A050535 // Rest]] (* Jean-François Alcover, Feb 11 2020, updated Mar 17 2020 *)

Extensions

More terms from Sean A. Irvine, Oct 02 2011
Name and comment swapped by Gus Wiseman, Nov 28 2018
a(0)=1 prepended by Andrew Howroyd, Oct 23 2019

A007719 Number of independent polynomial invariants of symmetric matrix of order n.

Original entry on oeis.org

1, 2, 4, 11, 30, 95, 328, 1211, 4779, 19902, 86682, 393072, 1847264, 8965027, 44814034, 230232789, 1213534723, 6552995689, 36207886517, 204499421849, 1179555353219, 6942908667578, 41673453738272, 254918441681030, 1588256152307002, 10073760672179505
Offset: 0

Views

Author

Keywords

Comments

Also, number of connected multigraphs with n edges (allowing loops) and any number of nodes.
Also the number of non-isomorphic connected multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}. - Gus Wiseman, Jul 18 2018

Examples

			From _Gus Wiseman_, Jul 18 2018: (Start)
Non-isomorphic representatives of the a(3) = 11 connected multiset partitions of {1, 1, 2, 2, 3, 3}:
  (112233),
  (1)(12233), (12)(1233), (112)(233), (123)(123),
  (1)(2)(1233), (1)(12)(233), (1)(23)(123), (12)(13)(23),
  (1)(2)(3)(123), (1)(2)(13)(23).
(End)
		

Crossrefs

Programs

  • Mathematica
    mob[m_, n_] := If[Mod[m, n] == 0, MoebiusMu[m/n], 0];
    EULERi[b_] := Module[{a, c, i, d}, c = {}; For[i = 1, i <= Length[b], i++,
      c = Append[c, i*b[[i]] - Sum[c[[d]]*b[[i - d]], {d, 1, i - 1}]]]; a = {};
      For[i = 1, i <= Length[b], i++, a = Append[a, (1/i)*Sum[mob[i, d]*c[[d]], {d, 1, i}]]]; Return[a]];
    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];
    Kq[q_, t_, k_] := SeriesCoefficient[1/Product[g = GCD[t, q[[j]]]; (1 - x^(q[[j]]/g))^g, {j, 1, Length[q]}], {x, 0, k}];
    RowSumMats[n_, m_, k_] := Module[{s = 0}, Do[s += permcount[q]* SeriesCoefficient[ Exp[Sum[Kq[q, t, k]/t x^t, {t, 1, n}]], {x, 0, n}], {q, IntegerPartitions[m]}]; s/m!];
    A007717 = Table[Print[n]; RowSumMats[n, 2 n, 2], {n, 0, 20}];
    Join[{1}, EULERi[Rest[A007717]]] (* Jean-François Alcover, Oct 29 2018, using Andrew Howroyd's code for A007717 *)

Formula

Inverse Euler transform of A007717.

Extensions

a(0)=1 added by Alberto Tacchella, Jun 20 2011
a(7)-a(25) from Franklin T. Adams-Watters, Jun 21 2011

A053419 Number of graphs with loops (symmetric relations) with n edges.

Original entry on oeis.org

1, 2, 5, 14, 38, 107, 318, 972, 3111, 10410, 36371, 132656, 504636, 1998361, 8224448, 35112342, 155211522, 709123787, 3342875421, 16234342515, 81102926848, 416244824068, 2192018373522, 11831511359378, 65387590986455, 369661585869273, 2135966349269550, 12604385044890628
Offset: 0

Views

Author

Vladeta Jovovic, Jan 10 2000

Keywords

Comments

In a multiset partition, two vertices are equivalent if in every block the multiplicity of the first is equal to the multiplicity of the second. a(n) is the number of non-isomorphic multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n} with no equivalent vertices. For example, non-isomorphic representatives of the a(2) = 5 multiset partitions are (1)(122), (11)(22), (1)(1)(22), (1)(2)(12), (1)(1)(2)(2). - Gus Wiseman, Jul 18 2018
a(n) is the number of unlabeled simple graphs with n edges rooted at one vertex. - Andrew Howroyd, Nov 22 2020

Crossrefs

Programs

Formula

Euler transform of A191970. - Andrew Howroyd, Oct 22 2019

Extensions

a(16)-a(24) from Max Alekseyev, Jan 22 2010
Terms a(25) and beyond from Andrew Howroyd, Oct 22 2019

A322115 Triangle read by rows where T(n,k) is the number of unlabeled connected multigraphs with loops with n edges and k vertices.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 4, 2, 1, 6, 11, 9, 3, 1, 9, 25, 34, 20, 6, 1, 12, 52, 104, 99, 49, 11, 1, 16, 94, 274, 387, 298, 118, 23, 1, 20, 162, 645, 1295, 1428, 881, 300, 47, 1, 25, 263, 1399, 3809, 5803, 5088, 2643, 765, 106, 1, 30, 407, 2823, 10187, 20645, 24606, 17872, 7878, 1998, 235
Offset: 0

Views

Author

Gus Wiseman, Nov 26 2018

Keywords

Examples

			Triangle begins:
  1
  1   1
  1   2   1
  1   4   4   2
  1   6  11   9   3
  1   9  25  34  20   6
  1  12  52 104  99  49  11
		

Crossrefs

Row sums are A007719. Diagonal k = n-1 is A000055.

Programs

  • PARI
    EulerT(v)={my(p=exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1); Vec(p/x,-#v)}
    InvEulerMT(u)={my(n=#u, p=log(1+x*Ser(u)), vars=variables(p)); Vec(serchop( sum(i=1, n, moebius(i)*substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i), 1))}
    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,x)={sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i],v[j])); g*x^(v[i]*v[j]/g))) + sum(i=1, #v, my(t=v[i]); ((t+1)\2)*x^t + if(t%2, 0, x^(t/2)))}
    G(n,m)={my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(edges(p,x) + O(x*x^m), -m))); s/n!}
    R(n)={Mat(apply(p->Col(p+O(y^n), -n), InvEulerMT(vector(n, k, 1 + y*Ser(G(k,n-1), y)))))}
    { my(T=R(10)); for(n=1, #T, print(T[n, 1..n])) } \\ Andrew Howroyd, Nov 30 2018

Extensions

Terms a(28) and beyond from Andrew Howroyd, Nov 30 2018

A321911 Number of distinct chromatic symmetric functions of simple connected graphs with n vertices.

Original entry on oeis.org

1, 1, 2, 6, 20, 103, 759
Offset: 1

Views

Author

Gus Wiseman, Nov 21 2018

Keywords

Comments

A stable partition of a graph is a set partition of the vertices where no edge has both ends in the same block. The chromatic symmetric function is given by X_G = Sum_p m(t(p)) where the sum is over all stable partitions p of G, t(p) is the integer partition whose parts are the block-sizes of p, and m is augmented monomial symmetric functions (see A321895).

Examples

			The a(4) = 6 connected chromatic symmetric functions (m is the augmented monomial symmetric function basis):
                    m(1111)
           m(211) + m(1111)
          2m(211) + m(1111)
  m(22) + 2m(211) + m(1111)
  m(22) + 3m(211) + m(1111)
  m(31) + 3m(211) + m(1111)
		

Crossrefs

Programs

  • Mathematica
    spsu[,{}]:={{}};spsu[foo,set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@spsu[Select[foo,Complement[#,Complement[set,s]]=={}&],Complement[set,s]]]/@Cases[foo,{i,_}];
    chromSF[g_]:=Sum[m[Sort[Length/@stn,Greater]],{stn,spsu[Select[Subsets[Union@@g],Select[DeleteCases[g,{_}],Function[ed,Complement[ed,#]=={}]]=={}&],Union@@g]}];
    simpleSpans[n_]:=simpleSpans[n]=If[n==0,{{}},Union@@Table[If[#=={},Union[ine,{{n}}],Union[Complement[ine,List/@#],{#,n}&/@#]]&/@Subsets[Range[n-1]],{ine,simpleSpans[n-1]}]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Union[chromSF/@Select[simpleSpans[n],Length[csm[#]]==1&]]],{n,6}]

A322151 Number of labeled connected graphs with loops with n edges (the vertices are {1,2,...,k} for some k).

Original entry on oeis.org

1, 2, 5, 27, 216, 2311, 30988, 499919, 9431026, 203743252, 4960335470, 134382267082, 4009794148101, 130668970606412, 4617468180528235, 175867725701333896, 7182126650899080024, 313063334893103361130, 14507460736615554141354, 712192629608088061633746
Offset: 0

Views

Author

Gus Wiseman, Nov 28 2018

Keywords

Crossrefs

Row sums of A322147. The unlabeled version is A191970.

Programs

  • Mathematica
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[multsubs[Range[n+1],2],{n}],And[Union@@#==Range[Max@@Union@@#],Length[csm[#]]==1]&]],{n,5}]
  • PARI
    Connected(v)={my(u=vector(#v)); for(n=1, #u, u[n]=v[n] - sum(k=1, n-1, binomial(n-1, k)*v[k]*u[n-k])); u}
    seq(n)={Vec(vecsum(Connected(vector(2*n, j, (1 + x + O(x*x^n))^binomial(j+1,2)))))} \\ Andrew Howroyd, Nov 28 2018

Extensions

Terms a(7) and beyond from Andrew Howroyd, Nov 28 2018

A322152 Number of labeled connected multigraphs with loops with n edges (the vertices are {1,2,...,k} for some k).

Original entry on oeis.org

1, 2, 7, 39, 314, 3359, 45000, 725269, 13670256, 295099184, 7179749707, 194399095705, 5797793490859, 188855813757729, 6671188010874785, 254007814638737649, 10370334196814589256, 451923738493729293016, 20937747226064522726151, 1027666505638118490940059
Offset: 0

Views

Author

Gus Wiseman, Nov 28 2018

Keywords

Crossrefs

Row sums of A322148. The unlabeled version is A007719.

Programs

  • Mathematica
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[multsubs[multsubs[Range[n+1],2],n],And[Union@@#==Range[Max@@Union@@#],Length[csm[#]]==1]&]],{n,5}]
  • PARI
    Connected(v)={my(u=vector(#v)); for(n=1, #u, u[n]=v[n] - sum(k=1, n-1, binomial(n-1, k)*v[k]*u[n-k])); u}
    seq(n)={Vec(vecsum(Connected(vector(2*n, j, 1/(1 - x + O(x*x^n))^binomial(j+1,2)))))} \\ Andrew Howroyd, Nov 28 2018

Extensions

Terms a(7) and beyond from Andrew Howroyd, Nov 28 2018
Showing 1-10 of 14 results. Next