A020555 -id:A020555 - cp's OEIS frontend

A316980 Number of non-isomorphic strict multiset partitions of weight n.

1, 1, 3, 8, 23, 63, 197, 588, 1892, 6140, 20734, 71472, 254090, 923900, 3446572, 13149295, 51316445, 204556612, 832467052, 3455533022, 14621598811, 63023667027, 276559371189, 1234802595648, 5606647482646, 25875459311317, 121324797470067, 577692044073205

Offset: 0

Gus Wiseman, Jul 18 2018

nonn

Also the number of nonnegative integer n X n matrices with sum of elements equal to n, under row and column permutations, with no equal rows (or alternatively, with no equal columns).

Also the number of non-isomorphic multiset partitions of weight n with no equivalent vertices. 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.

			Non-isomorphic representatives of the a(3) = 8 multiset partitions with no equivalent vertices (first column) and with no equal blocks (second column):
      (111) <-> (111)
      (122) <-> (1)(11)
    (1)(11) <-> (122)
    (1)(22) <-> (1)(22)
    (2)(12) <-> (2)(12)
  (1)(1)(1) <-> (123)
  (1)(2)(2) <-> (1)(23)
  (1)(2)(3) <-> (1)(2)(3)

Andrew Howroyd, Table of n, a(n) for n = 0..50

Cf. A000009, A001055, A007716, A007717, A020555, A045778, A130091.
Cf. A316974, A316978, A316979, A316981, A316983.
Cf. A049311, A059201, A319558, A319560, A319567.

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
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}
K(q, t, k)={EulerT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
a(n)={if(n==0, 1, my(s=0); forpart(q=n, my(p=sum(t=1, n, subst(x*Ser(K(q, t, n\t))/t, x, x^t))); s+=permcount(q)*polcoef(exp(p-subst(p,x,x^2)), n)); s/n!)} \\ Andrew Howroyd, Jan 21 2023

Euler transform of A319557. - Gus Wiseman, Sep 23 2018

a(7)-a(10) from Gus Wiseman, Sep 23 2018

Terms a(11) and beyond from Andrew Howroyd, Jan 19 2023

A316983 Number of non-isomorphic self-dual multiset partitions of weight n.

Original entry on oeis.org

1, 1, 2, 4, 9, 17, 36, 72, 155, 319, 677, 1429, 3094, 6648, 14518, 31796, 70491, 156818, 352371, 795952, 1813580, 4155367, 9594425, 22283566, 52122379, 122631874, 290432439, 691831161, 1658270316, 3997272089, 9692519896, 23631827354, 57943821449, 142834652193

Offset: 0

Gus Wiseman, Jul 18 2018

nonn

Also the number of nonnegative integer square symmetric matrices with sum of elements equal to n, under row and column permutations.

The dual of a multiset partition has, for each vertex, one block consisting of the indices (or positions) of the blocks containing that vertex, counted with multiplicity.

			Non-isomorphic representatives of the a(4) = 9 self-dual multiset partitions:
  (1111),
  (1)(222), (2)(122), (11)(22), (12)(12),
  (1)(1)(23), (1)(2)(33), (1)(3)(23),
  (1)(2)(3)(4).
The a(4) = 9 square symmetric matrices:
. [4]
.
. [3 0]  [2 0]  [2 1]  [1 1]
. [0 1]  [0 2]  [1 0]  [1 1]
.
. [2 0 0]  [1 1 0]  [0 1 1]
. [0 1 0]  [1 0 0]  [1 0 0]
. [0 0 1]  [0 0 1]  [1 0 0]
.
. [1 0 0 0]
. [0 1 0 0]
. [0 0 1 0]
. [0 0 0 1]

Andrew Howroyd, Table of n, a(n) for n = 0..50

Row sums of A320796.
Main diagonal of A318805.
Cf. A000009, A001055, A007716, A007717, A020555, A045778.
Cf. A316974, A316978, A316979, A316980, A316981.

vector(25, n, n--; T(n,n)) \\ T(n,k) defined in A318805. - Andrew Howroyd, Jan 16 2024

Terms a(9) and beyond from Andrew Howroyd, Sep 03 2018

A007717 Number of symmetric polynomial functions of degree n of a symmetric matrix (of indefinitely large size) under joint row and column permutations. Also number of multigraphs with n edges (allowing loops) on an infinite set of nodes.

Original entry on oeis.org

1, 2, 7, 23, 79, 274, 1003, 3763, 14723, 59663, 250738, 1090608, 4905430, 22777420, 109040012, 537401702, 2723210617, 14170838544, 75639280146, 413692111521, 2316122210804, 13261980807830, 77598959094772, 463626704130058, 2826406013488180, 17569700716557737

Offset: 0

Colin Mallows

nonn

Euler transform of A007719.
Also the number of non-isomorphic multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}. - Gus Wiseman, Jul 18 2018
Number of distinct n X 2n matrices with integer entries and rows sums 2, up to row and column permutations. - Andrew Howroyd, Sep 06 2018
a(n) is the number of unlabeled loopless multigraphs with n edges rooted at one vertex. - Andrew Howroyd, Nov 22 2020

			a(2) = 7 (here - denotes an edge, = denotes a pair of parallel edges and o is a loop):
  oo
  o o
  o-
  o -
  =
  --
  - -
From _Gus Wiseman_, Jul 18 2018: (Start)
Non-isomorphic representatives of the a(2) = 7 multiset partitions of {1, 1, 2, 2}:
  (1122),
  (1)(122), (11)(22), (12)(12),
  (1)(1)(22), (1)(2)(12),
  (1)(1)(2)(2).
(End)
From _Gus Wiseman_, Jan 08 2024: (Start)
Non-isomorphic representatives of the a(1) = 1 through a(3) = 7 rooted loopless multigraphs (root shown as singleton):
  {{1}}  {{1},{1,2}}  {{1},{1,2},{1,2}}
         {{1},{2,3}}  {{1},{1,2},{1,3}}
                      {{1},{1,2},{2,3}}
                      {{1},{1,2},{3,4}}
                      {{1},{2,3},{2,3}}
                      {{1},{2,3},{2,4}}
                      {{1},{2,3},{4,5}}
(End)

Huaien Li and David C. Torney, Enumerations of Multigraphs, 2002.

Andrew Howroyd, Table of n, a(n) for n = 0..50
Mateo Díaz, Dmitriy Drusvyatskiy, Jack Kendrick, and Rekha R. Thomas, Invariant Kernels: Rank Stabilization and Generalization Across Dimensions, arXiv:2502.01886 [math.OC], 2025. See p. 43.
Huaien Li and David C. Torney, Enumeration of unlabelled multigraphs, Ars Combin. 75 (2005) 171-188. MR2133219.
R. J. Mathar, Statistics on Small Graphs, arXiv:1709.09000 [math.CO] (2017) table 67.

Row n=2 of A331485.

Cf. A000664, A002620, A007716, A007719, A020555, A050531, A050532, A050535, A052171, A053418, A053419, A094574, A316972, A316974, A318951, A339065.

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!];
a[n_] := RowSumMats[n, 2n, 2];
Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 0, 25}] (* Jean-François Alcover, Oct 27 2018, after Andrew Howroyd *)

\\ See A318951 for RowSumMats
a(n)=RowSumMats(n, 2*n, 2); \\ Andrew Howroyd, Sep 06 2018

\\ See A339065 for G.
seq(n)={my(A=O(x*x^n)); Vec(G(2*n, x+A, [1]))} \\ Andrew Howroyd, Nov 22 2020

More terms from Vladeta Jovovic, Jan 26 2000

a(0)=1 prepended and a(16)-a(25) added by Max Alekseyev, Jun 21 2011

A050535 Number of loopless multigraphs on infinite set of nodes with n edges.

Original entry on oeis.org

1, 1, 3, 8, 23, 66, 212, 686, 2389, 8682, 33160, 132277, 550835, 2384411, 10709827, 49782637, 238998910, 1182772364, 6023860266, 31525780044, 169316000494, 932078457785, 5253664040426, 30290320077851, 178480713438362, 1073918172017297

Offset: 0

Vladeta Jovovic, Dec 29 1999

nonn

Also, a(n) is the number of n-rowed binary matrices with all row sums equal to 2, up to row and column permutation (see Jovovic's formula). Also, a(n) is the limit of A192517(m,n) as m grows. - Max Alekseyev, Oct 18 2017
Row sums of the triangle defined by the Multiset Transformation of A076864,
;
1;
2 1;
5 2 1;
12 8 2 1;
33 22 8 2 1;
103 72 26 8 2 1;
333 229 87 26 8 2 1;
1183 782 295 92 26 8 2 1;
4442 2760 1036 315 92 26 8 2 1;
17576 10270 3735 1129 321 92 26 8 2 1;
72810 39770 13976 4117 1154 321 92 26 8 2 1;
314595 160713 54132 15547 4237 1161 321 92 26 8 2 1;
- R. J. Mathar, Jul 18 2017
Also the number of non-isomorphic set multipartitions (multisets of sets) of {1, 1, 2, 2, 3, 3, ..., n, n}. - Gus Wiseman, Jul 18 2018

			From _Gus Wiseman_, Jul 18 2018: (Start)
Non-isomorphic representatives of the a(3) = 8 set multipartitions of {1, 1, 2, 2, 3, 3}:
  (123)(123)
  (1)(23)(123)
  (12)(13)(23)
  (1)(1)(23)(23)
  (1)(2)(3)(123)
  (1)(2)(13)(23)
  (1)(1)(2)(3)(23)
  (1)(1)(2)(2)(3)(3)
(End)

Frank Harary and Edgar M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 88, Eq. (4.1.18).

Andrew Howroyd, Table of n, a(n) for n = 0..50
George Barnes, Sanjaye Ramgoolam, and Michael Stephanou, Permutation invariant Gaussian matrix models for financial correlation matrices, arXiv:2306.04569 [q-fin.ST], 2023.
Frank Harary, The number of linear, directed, rooted, and connected graphs, Trans. Am. Math. Soc. 78 (1955) 445-463, eq. (24).
Vladeta Jovovic, Number of m-rowed binary matrices with all row sums equal to n, up to row and column permutation
Patrick T. Komiske, Eric M. Metodiev, and Jesse Thaler, Energy flow polynomials: A complete linear basis for jet substructure, arXiv:1712.07124 [hep-ph], 2017.
Tsuyoshi Miezaki, Akihiro Munemasa, Yusaku Nishimura, Tadashi Sakuma, and Shuhei Tsujie, Universal graph series, chromatic functions, and their index theory, arXiv:2403.09985 [math.CO], 2024. See p. 23.

Cf. A001399, A003082, A014395, A014396, A014397, A014398.
Cf. A058389, A050913, A058783, A058390, A058784, A058785, A058391, A058392, A001501, A058528.
Cf. A007716, A007717, A020555, A050535, A076864 (inverse Euler transf.), A076867 (Euler transform) A094574, A316974.

seq[n_] := G[2n, x+O[x]^n, {}] // CoefficientList[#, x]&;
seq[15] (* Jean-François Alcover, Dec 02 2020, using Andrew Howroyd's code for G in A339065 *)

a(n) = A192517(2*n,n) = A192517(m,n) for any m>=2*n. - Max Alekseyev, Oct 18 2017

Euler transform of A076864. - Andrew Howroyd, Oct 23 2019

More terms from Sean A. Irvine, Oct 02 2011

A316978 Number of factorizations of n into factors > 1 with no equivalent primes.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 4, 1, 1, 1, 5, 1, 4, 1, 4, 1, 1, 1, 7, 2, 1, 3, 4, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 4, 4, 1, 1, 12, 2, 4, 1, 4, 1, 7, 1, 7, 1, 1, 1, 7, 1, 1, 4, 11, 1, 1, 1, 4, 1, 1, 1, 16, 1, 1, 4, 4, 1, 1, 1, 12, 5, 1, 1, 7, 1, 1

Offset: 1

Gus Wiseman, Jul 18 2018

nonn

In a factorization, two primes are equivalent if each factor has in its prime factorization the same multiplicity of both primes.

			The a(36) = 7 factorizations are (2*2*3*3), (2*2*9), (2*3*6), (3*3*4), (2*18), (3*12), (4*9). Missing from this list are (6*6) and (36).

Cf. A001055, A007716, A007717, A020555, A045778, A162247, A281116, A303386.

Cf. A316974, A316979, A316980, A316981.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
dual[eds_]:=Table[First/@Position[eds,x],{x,Union@@eds}];
Table[Length[Select[facs[n],UnsameQ@@dual[primeMS/@#]&]],{n,100}]

a(prime^n) = A000041(n).

a(squarefree) = 1.

A278990 Number of loopless linear chord diagrams with n chords.

Original entry on oeis.org

1, 0, 1, 5, 36, 329, 3655, 47844, 721315, 12310199, 234615096, 4939227215, 113836841041, 2850860253240, 77087063678521, 2238375706930349, 69466733978519340, 2294640596998068569, 80381887628910919255, 2976424482866702081004, 116160936719430292078411

Offset: 0

N. J. A. Sloane, Dec 07 2016

nonn

See the signed version of these numbers, A000806, for much more information about these numbers.
From Gus Wiseman, Feb 27 2019: (Start)
Also the number of 2-uniform set partitions of {1..2n} containing no two successive vertices in the same block. For example, the a(3) = 5 set partitions are:
  {{1,3},{2,5},{4,6}}
  {{1,4},{2,5},{3,6}}
  {{1,4},{2,6},{3,5}}
  {{1,5},{2,4},{3,6}}
  {{1,6},{2,4},{3,5}}
(End)
From Gus Wiseman, Jul 05 2020: (Start)
Also the number of permutations of the multiset {1,1,2,2,...,n,n} with no two consecutive terms equal and where the first i appears before the first j for i < j. For example, the a(3) = 5 permutations are the following.
  (1,2,3,1,2,3)
  (1,2,3,1,3,2)
  (1,2,3,2,1,3)
  (1,2,3,2,3,1)
  (1,2,1,3,2,3)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..404 (terms 0..200 from Gheorghe Coserea)
Dmitry Efimov, The hafnian of Toeplitz matrices of a special type, perfect matchings and Bessel polynomials, arXiv:1904.08651 [math.CO], 2019.
H. Eriksson and A. Martin, Enumeration of Carlitz multipermutations, arXiv:1702.04177 [math.CO], 2017.
E. Krasko, I. Labutin, and A. Omelchenko, Enumeration of labelled and unlabelled Hamiltonian Cycles in complete k-partite graphs, arXiv:1709.03218 [math.CO], 2017, Table 1.
E. Krasko and A. Omelchenko, Enumeration of Chord Diagrams without Loops and Parallel Chords, arXiv:1601.05073 [math.CO], 2016.
E. Krasko and A. Omelchenko, Enumeration of Chord Diagrams without Loops and Parallel Chords, The Electronic Journal of Combinatorics, 24(3) (2017), #P3.43.
Gus Wiseman, The a(4) = 36 loopless linear chord diagrams.
Donovan Young, Counting Bubbles in Linear Chord Diagrams, arXiv:2311.01569 [math.CO], 2023.
Donovan Young, Bubbles in Linear Chord Diagrams: Bridges and Crystallized Diagrams, arXiv:2408.17232 [math.CO], 2024.

Column k=0 of A079267.
Column k=2 of A293157.
Row n=2 of A322013.
Cf. A000110, A000699 (topologically connected 2-uniform), A000806, A001147 (2-uniform), A003436 (cyclical version), A005493, A170941, A190823 (distance 3+ version), A322402, A324011, A324172.
Anti-run compositions are A003242.
Separable partitions are A325534.
Cf. A007716, A292884, A333489.
Other sequences involving the multiset {1,1,2,2,...,n,n}: A001147, A007717, A020555, A094574, A316972.

[n le 2 select 2-n else (2*n-3)*Self(n-1) + Self(n-2): n in [1..30]]; // G. C. Greubel, Sep 26 2023

RecurrenceTable[{a[n]== (2n-1)a[n-1] +a[n-2], a[0]==1, a[1]==0}, a, {n,0,20}] (* Vaclav Kotesovec, Sep 15 2017 *)
FullSimplify[Table[-I*(BesselI[1/2+n,-1] BesselK[3/2,1] - BesselI[3/2,-1] BesselK[1/2+ n,1]), {n,0,20}]] (* Vaclav Kotesovec, Sep 15 2017 *)
Table[(2 n-1)!! Hypergeometric1F1[-n,-2 n,-2], {n,0,20}] (* Eric W. Weisstein, Nov 14 2018 *)
Table[Sqrt[2/Pi]/E ((-1)^n Pi BesselI[1/2+n,1] +BesselK[1/2+n,1]), {n,0,20}] // FunctionExpand // FullSimplify (* Eric W. Weisstein, Nov 14 2018 *)
twouniflin[{}]:={{}};twouniflin[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@twouniflin[Complement[set,s]]]/@Table[{i,j},{j,Select[set,#>i+1&]}];
Table[Length[twouniflin[Range[n]]],{n,0,14,2}] (* Gus Wiseman, Feb 27 2019 *)

seq(N) = {
  my(a = vector(N)); a[1] = 0; a[2] = 1;
  for (n = 3, N, a[n] = (2*n-1)*a[n-1] + a[n-2]);
  concat(1, a);
};
seq(20) \\ Gheorghe Coserea, Dec 09 2016

def A278990_list(prec):
    P. = PowerSeriesRing(QQ, prec)
    return P( exp(-1+sqrt(1-2*x))/sqrt(1-2*x) ).egf_to_ogf().list()
A278990_list(30) # G. C. Greubel, Sep 26 2023

From Gheorghe Coserea, Dec 09 2016: (Start)
D-finite with recurrence a(n) = (2*n-1)*a(n-1) + a(n-2), with a(0) = 1, a(1) = 0.
E.g.f. y satisfies: 0 = (1-2*x)*y'' - 3*y' - y.
a(n) - a(n-1) = A003436(n) for all n >= 2. (End)
From Vaclav Kotesovec, Sep 15 2017: (Start)
a(n) = sqrt(2)*exp(-1)*(BesselK(1/2 + n, 1)/sqrt(Pi) - i*sqrt(Pi)*BesselI(1/2 + n, -1)), where i is the imaginary unit.
a(n) ~ 2^(n+1/2) * n^n  / exp(n+1). (End)
a(n) = A114938(n)/n! - Gus Wiseman, Jul 05 2020 (from Alexander Burstein's formula at A114938).
From G. C. Greubel, Sep 26 2023: (Start)
a(n) = (-1)^n * (i/e)*Sqrt(2/Pi) * BesselK(n + 1/2, -1).
G.f.: sqrt(Pi/(2*x)) * exp(-(1+x)^2/(2*x)) * Erfi((1+x)/sqrt(2*x)).
E.g.f.: exp(-1 + sqrt(1-2*x))/sqrt(1-2*x). (End)

a(0)=1 prepended by Gheorghe Coserea, Dec 09 2016

A114938 Number of permutations of the multiset {1,1,2,2,...,n,n} with no two consecutive terms equal.

Original entry on oeis.org

1, 0, 2, 30, 864, 39480, 2631600, 241133760, 29083420800, 4467125013120, 851371260364800, 197158144895712000, 54528028997584665600, 17752366094818747392000, 6720318485119046923315200, 2927066537906697348594432000, 1453437879238150456164433920000

Offset: 0

Hugo Pfoertner, Jan 08 2006

a(n) is also the number of (0,1)-matrices A=(a_ij) of size n X 2n such that each row has exactly two 1's and each column has exactly one 1 and with the restriction that no 1 stands on the line from a_11 to a_22. - Shanzhen Gao, Feb 24 2010
a(n) is the number of permutations of the multiset {1,1,2,2,...,n,n} with no fixed points. - Alexander Burstein, May 16 2020
Also the number of 2-uniform ordered set partitions of {1...2n} containing no two successive vertices in the same block. - Gus Wiseman, Jul 04 2020

			a(2) = 2 because there are two permutations of {1,1,2,2} avoiding equal consecutive terms: 1212 and 2121.

R. P. Stanley, Enumerative Combinatorics Volume I, Cambridge University Press, 1997. Chapter 2, Sieve Methods, Example 2.2.3, page 68.

Seiichi Manyama, Table of n, a(n) for n = 0..238 (terms 1..100 from Andrew Woods)
H. Eriksson and A. Martin, Enumeration of Carlitz multipermutations, arXiv:1702.04177 [math.CO], 2017.

Cf. A114939 = preferred seating arrangements of n couples.
Cf. A007060 = arrangements of n couples with no adjacent spouses; A007060(n) = 2^n * A114938(n) (this sequence).
Cf. A104548, A193638.
Cf. A278990 = number of loopless linear chord diagrams with n chords.
Cf. A000806 = Bessel polynomial y_n(-1).
The version for multisets with prescribed multiplicities is A335125.
The version for prime indices is A335452.
Anti-run compositions are counted by A003242.
Anti-run compositions are ranked by A333489.
Inseparable partitions are counted by A325535.
Inseparable partitions are ranked by A335448.
Separable partitions are counted by A325534.
Separable partitions are ranked by A335433.
Cf. A007716, A128695, A292884, A335126, A335434, A335451.
Other sequences involving the multiset {1,1,2,2,...,n,n}: A001147, A007717, A020555, A094574, A316972.
Row n=2 of A322093.

[1] cat [n le 2 select 2*(n-1) else n*(2*n-1)*Self(n-1) + (n-1)*n*Self(n-2): n in [1..20]]; // Vincenzo Librandi, Aug 10 2015

Table[Sum[Binomial[n,i](2n-i)!/2^(n-i) (-1)^i,{i,0,n}],{n,0,20}]  (* Geoffrey Critzer, Jan 02 2013, and adapted to the extension by Stefano Spezia, Nov 15 2018 *)
Table[Length[Select[Permutations[Join[Range[n],Range[n]]],!MatchQ[#,{_,x_,x_,_}]&]],{n,0,5}] (* Gus Wiseman, Jul 04 2020 *)
A114938[n_] := ((2 n)! Hypergeometric1F1[-n, -2 n, -2]) / 2^n;
Array[A114938, 17, 0]  (* Peter Luschny, Sep 04 2025 *)

A114938(n)=sum(k=0, n, binomial(n, k)*(-1)^(n-k)*(n+k)!/2^k);
vector(20, n, A114938(n-1)) \\ Michel Marcus, Aug 10 2015

def A114938(n): return (-1)^n*sum(binomial(n,k)*factorial(n+k)//(-2)^k for k in range(n+1))
[A114938(n) for n in range(31)] # G. C. Greubel, Sep 26 2023

a(n) = Sum_{k=0..n} ((binomial(n, k)*(-1)^(n-k)*(n+k)!)/2^k).
a(n) = (-1)^n * n! * A000806(n), n>0. - Vladeta Jovovic, Nov 19 2009
a(n) = n*(2*n-1)*a(n-1) + (n-1)*n*a(n-2). - Vaclav Kotesovec, Aug 07 2013
a(n) ~ 2^(n+1)*n^(2*n)*sqrt(Pi*n)/exp(2*n+1). - Vaclav Kotesovec, Aug 07 2013
a(n) = n! * A278990(n). - Alexander Burstein, May 16 2020
From G. C. Greubel, Sep 26 2023: (Start)
a(n) = (-1)^n * (i/e)*sqrt(2/Pi) * n! * BesselK(n+1/2, -1).
a(n) = [n! * (1/x) * p_{n+1}(x)]|A104548%20for%20p">{x= -1} (See A104548 for p{n}(x)).
E.g.f.: sqrt(Pi/(2*x)) * exp(-(1+x)^2/(2*x)) * erfi((1+x)/sqrt(2*x)).
Sum_{n >= 0} a(n)*x^n/(n!)^2 = exp(sqrt(1-2*x))/sqrt(1-2*x).
Sum_{n >= 0} a(n)*x^n/(n!*(n+1)!) = ( 1 - exp(-1 + sqrt(1-2*x)) )/x. (End)
a(n) = ((2*n)!/2^n) * hypergeom([-n], [-2*n], -2]) = A007060(n) / 2^n. - Peter Luschny, Sep 04 2025

a(0)=1 prepended by Seiichi Manyama, Nov 15 2018

A014500 Number of graphs with unlabeled (non-isolated) nodes and n labeled edges.

Original entry on oeis.org

1, 1, 2, 9, 70, 794, 12055, 233238, 5556725, 158931613, 5350854707, 208746406117, 9315261027289, 470405726166241, 26636882237942128, 1678097862705130667, 116818375064650241036, 8932347052564257212796, 746244486452472386213939, 67796741482683128375533560

Offset: 0

Simon Plouffe, Gilbert Labelle (gilbert(AT)lacim.uqam.ca)

nonn

G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004.

Alois P. Heinz, Table of n, a(n) for n = 0..100
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Peter Cameron, Thomas Prellberg, Dudley Stark, Asymptotic enumeration of 2-covers and line graphs, Discrete Math. 310 (2010), no. 2, 230-240 (see u_n).
G. Labelle, Counting enriched multigraphs according to the number of their edges (or arcs), Discrete Math., 217 (2000), 237-248.
G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004. [Cached copy, with permission]

Row n=2 of A331126.

Cf. A020554, A020555, A014501, A060053.

read("transforms") ;
A020556 := proc(n) local k; add((-1)^(n+k)*binomial(n, k)*combinat[bell](n+k), k=0..n) end proc:
A014500 := proc(n) local i,gexp,lexp;
gexp := [seq(1/2^i/i!,i=0..n+1)] ;
lexp := add( A020556(i)*((log(1+x))/2)^i/i!,i=0..n+1) ;
lexp := taylor(lexp,x=0,n+1) ;
lexp := gfun[seriestolist](lexp,'ogf') ;
CONV(gexp,lexp) ; op(n+1,%)*n! ; end proc:
seq(A014500(n),n=0..20) ; # R. J. Mathar, Jul 03 2011

max = 20; A020556[n_] := Sum[(-1)^(n+k)*Binomial[n, k]*BellB[n+k], {k, 0, n}]; egf = Exp[x/2]*Sum[A020556[n]*(Log[1+x]/2)^n/n!, {n, 0, max}] + O[x]^max; CoefficientList[egf, x]*Range[0, max-1]! (* Jean-François Alcover, Feb 19 2017, after Vladeta Jovovic *)

\\ here egf1 is A020556 as e.g.f.
egf1(n)={my(bell=serlaplace(exp(exp(x + O(x^(2*n+1)))-1))); sum(i=0, n, sum(k=0, i, (-1)^k*binomial(i,k)*polcoef(bell, 2*i-k))*x^i/i!) + O(x*x^n)}
seq(n)={my(B=egf1(n), L=log(1+x + O(x*x^n))/2); Vec(serlaplace(exp(x/2 + O(x*x^n))*sum(k=0, n, polcoef(B ,k)*L^k)))} \\ Andrew Howroyd, Jan 13 2020

E.g.f.: exp(-1+x/2)*Sum((1+x)^binomial(n, 2)/n!, n=0..infinity) [probably in the Labelle paper]. - Vladeta Jovovic, Apr 27 2004
E.g.f.: exp(x/2)*Sum(A020556(n)*(log(1+x)/2)^n/n!, n=0..infinity). - Vladeta Jovovic, May 02 2004
Binomial transform of A060053.
The e.g.f.'s of A020554 (S(x)) and A014500 (U(x)) are related by S(x) = U(e^x-1).
The e.g.f.'s of A014500 (U(x)) and A060053 (V(x)) are related by U(x) = e^x*V(x).

A094574 Number of (<=2)-covers of an n-set.

Original entry on oeis.org

1, 1, 5, 40, 457, 6995, 136771, 3299218, 95668354, 3268445951, 129468914524, 5868774803537, 301122189141524, 17327463910351045, 1109375488487304027, 78484513540137938209, 6098627708074641312182, 517736625823888411991202, 47791900951140948275632148

Offset: 0

Goran Kilibarda, Vladeta Jovovic, May 12 2004

nonn

Also the number of strict multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}. For example, the a(2) = 5 strict multiset partitions of {1, 1, 2, 2} are (1122), (1)(122), (2)(112), (11)(22), (1)(2)(12). - Gus Wiseman, Jul 18 2018

			From _Gus Wiseman_, Sep 02 2019: (Start)
These are set-systems covering {1..n} with vertex-degrees <= 2. For example, the a(3) = 40 covers are:
  {123}  {1}{23}    {1}{2}{3}     {1}{2}{3}{12}
         {2}{13}    {1}{2}{13}    {1}{2}{3}{13}
         {3}{12}    {1}{2}{23}    {1}{2}{3}{23}
         {1}{123}   {1}{3}{12}    {1}{2}{13}{23}
         {12}{13}   {1}{3}{23}    {1}{2}{3}{123}
         {12}{23}   {2}{3}{12}    {1}{3}{12}{23}
         {13}{23}   {2}{3}{13}    {2}{3}{12}{13}
         {2}{123}   {1}{12}{23}
         {3}{123}   {1}{13}{23}
         {12}{123}  {1}{2}{123}
         {13}{123}  {1}{3}{123}
         {23}{123}  {2}{12}{13}
                    {2}{13}{23}
                    {2}{3}{123}
                    {3}{12}{13}
                    {3}{12}{23}
                    {12}{13}{23}
                    {1}{23}{123}
                    {2}{13}{123}
                    {3}{12}{123}
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..300

Row n=2 of A219585. - Alois P. Heinz, Nov 23 2012
Dominated by A003465.
Graphs with vertex-degrees <= 2 are A136281.
Cf. A002718, A007716, A020554, A020555, A050535, A094574, A136284, A316974, A327104, A327106, A327229.
Main diagonal of A346517.

facs[n_]:=facs[n]=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[Array[Prime,n,1,Times]^2],UnsameQ@@#&]],{n,0,6}] (* Gus Wiseman, Jul 18 2018 *)
m = 20;
a094577[n_] := Sum[Binomial[n, k]*BellB[2 n - k], {k, 0, n}];
egf = Exp[(1 - Exp[x])/2]*Sum[a094577[n]*(x/2)^n/n!, {n, 0, m}] + O[x]^m;
CoefficientList[egf + O[x]^m, x]*Range[0, m-1]! (* Jean-François Alcover, May 13 2019 *)

Row sums of A094573.

E.g.f: exp(-1-1/2*(exp(x)-1))*Sum(exp(x*binomial(n+1, 2))/n!, n=0..infinity) or exp((1-exp(x))/2)*Sum(A094577 (n)*(x/2)^n/n!, n=0..infinity).

A219727 Number A(n,k) of k-partite partitions of {n}^k into k-tuples; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 5, 9, 3, 1, 1, 15, 66, 31, 5, 1, 1, 52, 712, 686, 109, 7, 1, 1, 203, 10457, 27036, 6721, 339, 11, 1, 1, 877, 198091, 1688360, 911838, 58616, 1043, 15, 1, 1, 4140, 4659138, 154703688, 231575143, 26908756, 476781, 2998, 22, 1

Offset: 0

Alois P. Heinz, Nov 26 2012

A(n,k) is the number of factorizations of m^n where m is a product of k distinct primes. A(2,2) = 9: (2*3)^2 = 36 has 9 factorizations: 36, 3*12, 4*9, 3*3*4, 2*18, 6*6, 2*3*6, 2*2*9, 2*2*3*3.

A(n,k) is the number of (n*k) X k matrices with nonnegative integer entries and column sums n up to permutation of rows. - Andrew Howroyd, Dec 10 2018

			A(1,3) = 5: [(1,1,1)], [(1,1,0),(0,0,1)], [(1,0,1),(0,1,0)], [(1,0,0),(0,1,0),(0,0,1)], [(0,1,1),(1,0,0)].
A(2,2) = 9: [(2,2)], [(2,1),(0,1)], [(2,0),(0,2)], [(2,0),(0,1),(0,1)], [(1,2),(1,0)], [(1,1),(1,1)], [(1,1),(1,0),(0,1)], [(1,0),(1,0),(0,2)], [(1,0),(1,0),(0,1),(0,1)].
Square array A(n,k) begins:
  1,   1,    1,      1,        1,         1,         1,       1, ...
  1,   1,    2,      5,       15,        52,       203,     877, ...
  1,   2,    9,     66,      712,     10457,    198091, 4659138, ...
  1,   3,   31,    686,    27036,   1688360, 154703688, ...
  1,   5,  109,   6721,   911838, 231575143, ...
  1,   7,  339,  58616, 26908756, ...
  1,  11, 1043, 476781, ...
  1,  15, 2998, ...

Andrew Howroyd, Table of n, a(n) for n = 0..209

Columns k=0..3 give: A000012, A000041, A002774, A219678.
Rows n=0..4 give: A000012, A000110, A020555, A322487, A358781.
Main diagonal gives A322488.
Cf. A188392, A219585 (partitions of {n}^k into distinct k-tuples), A256384, A318284, A318951.

EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); EulerT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)}
T(n, k)={my(m=n*k, q=Vec(exp(O(x*x^m) + intformal((x^n-1)/(1-x)))/(1-x))); if(n==0, 1, sum(j=0, m, my(s=0); forpart(p=j, s+=D(p,n,k), [1,n]); s*q[#q-j]))} \\ Andrew Howroyd, Dec 11 2018

cp's OEIS Frontend

A316980 Number of non-isomorphic strict multiset partitions of weight n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

Extensions

A316983 Number of non-isomorphic self-dual multiset partitions of weight n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Extensions

A007717 Number of symmetric polynomial functions of degree n of a symmetric matrix (of indefinitely large size) under joint row and column permutations. Also number of multigraphs with n edges (allowing loops) on an infinite set of nodes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Extensions

A050535 Number of loopless multigraphs on infinite set of nodes with n edges.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A316978 Number of factorizations of n into factors > 1 with no equivalent primes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A278990 Number of loopless linear chord diagrams with n chords.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

Extensions

A114938 Number of permutations of the multiset {1,1,2,2,...,n,n} with no two consecutive terms equal.

Views

Author

Keywords

Comments