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 11-20 of 27 results. Next

A316974 Number of non-isomorphic strict multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}.

Original entry on oeis.org

1, 1, 4, 14, 49, 173, 652, 2494
Offset: 0

Views

Author

Gus Wiseman, Jul 17 2018

Keywords

Comments

Also the number of unlabeled multigraphs with n edges, allowing loops, spanning an initial interval of positive integers with no equivalent vertices (two vertices are equivalent if in every edge the multiplicity of the first is equal to the multiplicity of the second). For example, non-isomorphic representatives of the a(2) = 4 multigraphs are {(1,2),(1,3)}, {(1,1),(1,2)}, {(1,1),(2,2)}, {(1,1),(1,1)}.

Examples

			Non-isomorphic representatives of the a(3) = 14 strict multiset partitions:
  (112233),
  (1)(12233), (11)(2233), (12)(1233), (112)(233),
  (1)(2)(1233), (1)(12)(233), (1)(23)(123), (2)(11)(233), (11)(22)(33), (12)(13)(23),
  (1)(2)(3)(123), (1)(2)(12)(33), (1)(2)(13)(23).
		

Crossrefs

Extensions

a(7) from Andrew Howroyd, Feb 07 2020

A346500 Number A(n,k) of partitions of the (n+k)-multiset {1,2,...,n,1,2,...,k}; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 5, 4, 4, 5, 15, 11, 9, 11, 15, 52, 36, 26, 26, 36, 52, 203, 135, 92, 66, 92, 135, 203, 877, 566, 371, 249, 249, 371, 566, 877, 4140, 2610, 1663, 1075, 712, 1075, 1663, 2610, 4140, 21147, 13082, 8155, 5133, 3274, 3274, 5133, 8155, 13082, 21147
Offset: 0

Views

Author

Alois P. Heinz, Jul 20 2021

Keywords

Comments

Also number A(n,k) of factorizations of Product_{i=1..n} prime(i) * Product_{i=1..k} prime(i); A(2,2) = 9: 2*2*3*3, 3*3*4, 6*6, 2*3*6, 4*9, 2*2*9, 3*12, 2*18, 36.

Examples

			A(2,2) = 9: 1122, 11|22, 12|12, 1|122, 112|2, 11|2|2, 1|1|22, 1|12|2, 1|1|2|2.
Square array A(n,k) begins:
    1,    1,    2,     5,    15,     52,     203,     877, ...
    1,    2,    4,    11,    36,    135,     566,    2610, ...
    2,    4,    9,    26,    92,    371,    1663,    8155, ...
    5,   11,   26,    66,   249,   1075,    5133,   26683, ...
   15,   36,   92,   249,   712,   3274,   16601,   91226, ...
   52,  135,  371,  1075,  3274,  10457,   56135,  325269, ...
  203,  566, 1663,  5133, 16601,  56135,  198091, 1207433, ...
  877, 2610, 8155, 26683, 91226, 325269, 1207433, 4659138, ...
  ...
		

Crossrefs

Columns (or rows) k=0-10 give: A000110, A035098, A322764, A322768, A346881, A346882, A346883, A346884, A346885, A346886, A346887.
Main diagonal gives A020555.
First upper (or lower) diagonal gives A322766.
Second upper (or lower) diagonal gives A322767.
Antidiagonal sums give A346490.
A(2n,n) gives A322769.

Programs

  • Maple
    g:= proc(n, k) option remember; uses numtheory; `if`(n>k, 0, 1)+
         `if`(isprime(n), 0, add(`if`(d>k or max(factorset(n/d))>d, 0,
            g(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    p:= proc(n) option remember; `if`(n=0, 1, p(n-1)*ithprime(n)) end:
    A:= (n, k)-> g(p(n)*p(k)$2):
    seq(seq(A(n, d-n), n=0..d), d=0..10);
    # second Maple program:
    b:= proc(n) option remember; `if`(n=0, 1,
          add(b(n-j)*binomial(n-1, j-1), j=1..n))
        end:
    A:= proc(n, k) option remember; `if`(n
    				
  • Mathematica
    b[n_] := b[n] = If[n == 0, 1, Sum[b[n-j] Binomial[n-1, j-1], {j, 1, n}]];
    A[n_, k_] := A[n, k] = If[n < k, A[k, n],
         If[k == 0, b[n], (A[n + 1, k - 1] + Sum[A[n - k + j, j]*
         Binomial[k - 1, j], {j, 0, k - 1}] + A[n, k - 1])/2]];
    Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 10}] // Flatten (* Jean-François Alcover, Aug 18 2021, after Alois P. Heinz's second program *)

Formula

A(n,k) = A001055(A002110(n)*A002110(k)).
A(n,k) = A(k,n).
A(n,k) = A322765(abs(n-k),min(n,k)).

A020554 Number of multigraphs on n labeled edges (without loops).

Original entry on oeis.org

1, 1, 3, 16, 139, 1750, 29388, 624889, 16255738, 504717929, 18353177160, 769917601384, 36803030137203, 1984024379014193, 119571835094300406, 7995677265437541258, 589356399302126773920, 47609742627231823142029, 4193665147256300117666879
Offset: 0

Views

Author

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

Keywords

Comments

Or, number of bicoverings of an n-set.
Or, number of 2-covers of [1,...,n].
Also the number of set multipartitions (multisets of sets) of {1, 1, 2, 2, 3, 3, ..., n, n}. - Gus Wiseman, Jul 18 2018

Examples

			From _Gus Wiseman_, Jul 18 2018: (Start)
The a(3) = 16 set multipartitions of {1, 1, 2, 2, 3, 3}:
  (123)(123)
  (1)(23)(123) (2)(13)(123) (3)(12)(123) (12)(13)(23)
  (1)(1)(23)(23) (1)(2)(3)(123) (1)(2)(13)(23) (1)(3)(12)(23) (2)(2)(13)(13) (2)(3)(12)(13) (3)(3)(12)(12)
  (1)(1)(2)(3)(23) (1)(2)(2)(3)(13) (1)(2)(3)(3)(12)
  (1)(1)(2)(2)(3)(3)
(End)
		

References

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

Crossrefs

Programs

  • Mathematica
    Ceiling[ CoefficientList[ Series[ Exp[ -1 + (Exp[ z ] - 1)/2 ]Sum[ Exp[ s(s - 1)z/2 ]/s!, {s, 0, 21} ], {z, 0, 9} ], z ] Table[ n!, {n, 0, 9} ] ] (* Mitch Harris, May 01 2004 *)
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Table[Length[Select[mps[Ceiling[Range[1/2,n,1/2]]],And@@UnsameQ@@@#&]],{n,5}] (* Gus Wiseman, Jul 18 2018 *)

Formula

E.g.f.: exp(-3/2+exp(x)/2) * Sum_{n>=0} exp(binomial(n, 2)*x)/n! [Comtet]. - Vladeta Jovovic, Apr 27 2004
E.g.f. (an equivalent version in Maple format): G:=exp(-1+(exp(z)-1)/2)*sum(exp(s*(s-1)*z/2)/s!, s=0..infinity);
E.g.f.: exp((exp(x)-1)/2)*Sum_{n>=0} A020556(n)*(x/2)^n/n!. - Vladeta Jovovic, May 02 2004
Stirling_2 transform of A014500.
The e.g.f.'s of A020554 (S(x)) and A014500 (U(x)) are related by S(x) = U(e^x-1).

A096443 Number of partitions of a multiset whose signature is the n-th partition (in Mathematica order).

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 5, 7, 9, 11, 15, 7, 12, 16, 21, 26, 36, 52, 11, 19, 29, 38, 31, 52, 74, 66, 92, 135, 203, 15, 30, 47, 64, 57, 98, 141, 109, 137, 198, 296, 249, 371, 566, 877, 22, 45, 77, 105, 97, 171, 250, 109, 212, 269, 392, 592, 300, 444, 560, 850, 1315, 712, 1075
Offset: 0

Views

Author

Jon Wild, Aug 11 2004

Keywords

Comments

The signature of a multiset is the partition consisting of the multiplicities of its elements; e.g., {a,a,a,b,c} is represented by [3,1,1]. The Mathematica order for partitions orders by ascending number of total elements, then by descending numerical order of its representation. The list begins:
n.....#elements.....n-th partition
0.....0 elements:....[]
1.....1 element:.....[1]
2.....2 elements:....[2]
3....................[1,1]
4.....3 elements:....[3]
5....................[2,1]
6....................[1,1,1]
7.....4 elements:....[4]
8....................[3,1]
9....................[2,2]
10...................[2,1,1]
11...................[1,1,1,1]
12....5 elements:....[5]
13...................[4,1]
A000041 and A000110 are subsequences for conjugate partitions. A000070 and A035098 are also subsequences for conjugate partitions. - Alford Arnold, Dec 31 2005
A002774 and A020555 is another pair of subsequences for conjugate partitions. - Franklin T. Adams-Watters, May 16 2006

Examples

			The 10th partition is [2,1,1]. The partitions of a multiset whose elements have multiplicities 2,1,1 - for example, {a,a,b,c} - are:
{{a,a,b,c}}
{{a,a,b},{c}}
{{a,a,c},{b}}
{{a,b,c},{a}}
{{a,a},{b,c}}
{{a,b},{a,c}}
{{a,a},{b},{c}}
{{a,b},{a},{c}}
{{a,c},{a},{b}}
{{b,c},{a},{a}}
{{a},{a},{b},{c}}
We see there are 11 partitions of this multiset, so a(10)=11.
Also, a(n) is the number of distinct factorizations of A063008(n). For example, A063008(10) = 60 and 60 has 11 factorizations: 60, 30*2, 20*3, 15*4, 15*2*2, 12*5, 10*6, 10*3*2, 6*5*2, 5*4*3, 5*3*2*2 which confirms that a(10) = 11.
		

Crossrefs

Programs

  • Mathematica
    MultiPartiteP[n : {___Integer?NonNegative}] :=
    Block[{p, $RecursionLimit = 1024, firstPositive},
      firstPositive =
       Compile[{{vv, _Integer, 1}},
        Module[{k = 1}, Do[If[el == 0, k++, Break[]], {el, vv}]; k]];
      p[{0 ...}] := 1;
      p[v_] :=
       p[v] = Module[{len = Length[v], it, k, zeros, sum, pos, gcd},
         it = Array[k, len];
         pos = firstPositive[v];
         zeros = ConstantArray[0, len];
         sum = 0;
         Do[If[it == zeros, Continue[]];
          gcd = GCD @@ it;
          sum += it[[pos]] DivisorSigma[-1, gcd] p[v - it];,
          Evaluate[Sequence @@ Thread[{it, 0, v}]]];
         sum/v[[pos]]];
      p[n]];
    ParallelMap[MultiPartiteP,
    Flatten[Table[IntegerPartitions[k], {k, 0, 8}], 1]]
    (* Oleksandr Pavlyk, Jan 23 2011 *)

Extensions

Edited by Franklin T. Adams-Watters, May 16 2006

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

A322765 Array read by upwards antidiagonals: T(m,n) = number of set partitions of the multiset consisting of one copy each of x_1, x_2, ..., x_m, and two copies each of y_1, y_2, ..., y_n, for m >= 0, n >= 0.

Original entry on oeis.org

1, 1, 2, 2, 4, 9, 5, 11, 26, 66, 15, 36, 92, 249, 712, 52, 135, 371, 1075, 3274, 10457, 203, 566, 1663, 5133, 16601, 56135, 198091, 877, 2610, 8155, 26683, 91226, 325269, 1207433, 4659138, 4140, 13082, 43263, 149410, 537813, 2014321, 7837862, 31638625, 132315780
Offset: 0

Views

Author

N. J. A. Sloane, Dec 30 2018

Keywords

Examples

			The array begins:
    1,    2,     9,     66,      712,     10457,      198091, ...
    1,    4,    26,    249,     3274,     56135,     1207433, ...
    2,   11,    92,   1075,    16601,    325269,     7837862, ...
    5,   36,   371,   5133,    91226,   2014321,    53840640, ...
   15,  135,  1663,  26683,   537813,  13241402,   389498179, ...
   52,  566,  8155, 149410,  3376696,  91914202,  2955909119, ...
  203, 2610, 43263, 894124, 22451030, 670867539, 23456071495, ...
  ...
		

References

  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778.

Crossrefs

Rows include A020555, A322766, A322767.
Columns include A000110, A035098, A322764, A322768.
Main diagonal is A322769.
See A322770 for partitions into distinct parts.

Programs

  • Maple
    B := n -> combinat[bell](n):
    P := proc(m,n) local k; global B; option remember;
    if n = 0 then B(m)  else
    (1/2)*( P(m+2,n-1) + P(m+1,n-1) + add( binomial(n-1,k)*P(m,k), k=0..n-1) ); fi; end; # P(m,n) (which is Knuth's notation) is T(m,n)
  • Mathematica
    P[m_, n_] := P[m, n] = If[n == 0, BellB[m], (1/2)(P[m+2, n-1] + P[m+1, n-1] + Sum[Binomial[n-1, k] P[m, k], {k, 0, n-1}])];
    Table[P[m-n, n], {m, 0, 8}, {n, 0, m}] // Flatten (* Jean-François Alcover, Jan 02 2019, from Maple *)
  • PARI
    {T(n, k) = if(k==0, sum(j=0, n, stirling(n, j, 2)), (T(n+2, k-1)+T(n+1, k-1)+sum(j=0, k-1, binomial(k-1, j)*T(n, j)))/2)} \\ Seiichi Manyama, Nov 21 2020

Formula

Knuth p. 779 gives a recurrence using the Bell numbers A000110 (see Maple program).
From Alois P. Heinz, Jul 21 2021: (Start)
A(n,k) = A001055(A002110(n+k)*A002110(k)).
A(n,k) = A346500(n+k,k). (End)

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

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 3, 1, 3, 1, 1, 1, 5, 1, 1, 2, 3, 1, 1, 1, 3, 1, 1, 1, 4, 1, 1, 1, 5, 1, 1, 1, 3, 3, 1, 1, 7, 1, 3, 1, 3, 1, 5, 1, 5, 1, 1, 1, 6, 1, 1, 3, 4, 1, 1, 1, 3, 1, 1, 1, 9, 1, 1, 3, 3, 1, 1, 1, 7, 2, 1, 1, 6, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Jul 18 2018

Keywords

Comments

In a factorization, two primes are equivalent if each factor has in its prime factorization the same multiplicity of both primes. For example, in 60 = (2*30) the primes {3, 5} are equivalent but {2, 3} and {2, 5} are not.

Examples

			The a(24) = 5 factorizations are (2*3*4), (2*12), (3*8), (4*6), (24).
The a(36) = 4 factorizations are (2*3*6), (2*18), (3*12), (4*9).
		

Crossrefs

Programs

  • Mathematica
    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],And[UnsameQ@@#,UnsameQ@@dual[primeMS/@#]]&]],{n,100}]

Formula

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

A316972 Number of connected multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}.

Original entry on oeis.org

1, 2, 5, 28, 277, 3985, 76117, 1833187, 53756682, 1871041538, 75809298105, 3521419837339, 185235838688677, 10923147890901151, 715989783027216302, 51793686238309903860, 4109310551278549543317, 355667047514571431358297, 33422937748872646130124797
Offset: 0

Views

Author

Gus Wiseman, Jul 17 2018

Keywords

Comments

Note that all connected multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n} are strict except for (123...n)(123...n).

Examples

			The a(2) = 5 connected multiset partitions of {1, 1, 2, 2} are (1122), (1)(122), (2)(112), (12)(12), (1)(2)(12). The multiset partitions (11)(22), (1)(1)(22), (2)(2)(11), (1)(1)(2)(2) are not connected.
		

Crossrefs

Programs

  • Mathematica
    nn=10;
    ser=Exp[-3/2+Exp[x]/2]*Sum[Exp[Binomial[n+1,2]*x]/n!,{n,0,3*nn}];
    Round/@(CoefficientList[Series[1+Log[ser],{x,0,nn}],x]*Array[Factorial,nn+1,0]) (* based on Jean-François Alcover after Vladeta Jovovic *)
    (*second program *)
    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]]]]]]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Length/@Table[Select[mps[Ceiling[Range[1/2,n,1/2]]],Length[csm[#]]==1&],{n,4}]

Formula

Logarithmic transform of A020555.

A014501 Number of graphs with loops, having unlabeled (non-isolated) nodes and n labeled edges.

Original entry on oeis.org

1, 2, 7, 43, 403, 5245, 89132, 1898630, 49209846, 1517275859, 54669946851, 2269075206395, 107199678164289, 5707320919486026, 339510756324234931, 22400182888853554291, 1628654713107465602783, 129754625253841669625051
Offset: 0

Views

Author

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

Keywords

References

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

Crossrefs

Row n=2 of A331161.

Formula

E.g.f.: exp(-1+x/2)*Sum((1+x)^binomial(n+1, 2)/n!, n=0..infinity) [probably in the Labelle paper]. - Vladeta Jovovic, Apr 27 2004
Previous Showing 11-20 of 27 results. Next