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 466 results. Next

A321390 Third Moebius transform of A007716. Number of non-isomorphic aperiodic multiset partitions of weight n whose parts have relatively prime periods and whose dual is also an aperiodic multiset partition.

Original entry on oeis.org

1, 1, 1, 7, 24, 88, 265, 907, 2929, 9918, 33931, 119366, 428314, 1574221, 5913415, 22699536, 88994103, 356058537, 1453049451, 6044132791, 25612496016, 110503624870, 485160989937, 2166488899639, 9835208617114, 45370059225048
Offset: 0

Views

Author

Gus Wiseman, Nov 08 2018

Keywords

Comments

The Moebius transform c of a sequence b is c(n) = Sum_{d|n} mu(d) * b(n/d).
Also the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns where the multiset of rows and the multiset of columns are both aperiodic and the nonzero entries are relatively prime, up to row and column permutations.
A multiset is aperiodic if its multiplicities are relatively prime. The period of a multiset is the GCD of its multiplicities.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 24 multiset partitions:
  {{1}}  {{1},{2}}  {{1,2,2}}      {{1,2,2,2}}
                    {{1},{1,1}}    {{1,2,3,3}}
                    {{1},{2,2}}    {{1},{1,1,1}}
                    {{1},{2,3}}    {{1},{1,2,2}}
                    {{2},{1,2}}    {{1},{2,2,2}}
                    {{1},{2},{2}}  {{1,2},{2,2}}
                    {{1},{2},{3}}  {{1},{2,3,3}}
                                   {{1,2},{3,3}}
                                   {{1},{2,3,4}}
                                   {{1,3},{2,3}}
                                   {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{1,1}}
                                   {{1},{1},{2,2}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}
		

Crossrefs

A001970 Functional determinants; partitions of partitions; Euler transform applied twice to all 1's sequence.

Original entry on oeis.org

1, 1, 3, 6, 14, 27, 58, 111, 223, 424, 817, 1527, 2870, 5279, 9710, 17622, 31877, 57100, 101887, 180406, 318106, 557453, 972796, 1688797, 2920123, 5026410, 8619551, 14722230, 25057499, 42494975, 71832114, 121024876, 203286806, 340435588, 568496753, 946695386
Offset: 0

Views

Author

Keywords

Comments

a(n) = number of partitions of n, when for each k there are p(k) different copies of part k. E.g., let the parts be 1, 2a, 2b, 3a, 3b, 3c, 4a, 4b, 4c, 4d, 4e, ... Then the a(4) = 14 partitions of 4 are: 4 = 4a = 4b = ... = 4e = 3a+1 = 3b+1 = 3c+1 = 2a+2a = 2a+2b = 2b+2b = 2a+1 = 2b+1 = 1+1+1+1.
Equivalently (Cayley), a(n) = number of 2-dimensional partitions of n. E.g., for n = 4 we have:
4 31 3 22 2 211 21 2 2 1111 111 11 11 1
1 2 1 11 1 1 11 1 1
1 1 1
1
Also total number of different species of singularity for conjugate functions with n letters (Sylvester).
According to [Belmans], this sequence gives "[t]he number of Segre symbols for the intersection of two quadrics in a fixed dimension". - Eric M. Schmidt, Sep 02 2017
From Gus Wiseman, Jul 30 2022: (Start)
Also the number of non-isomorphic multiset partitions of weight n with all constant blocks. The strict case is A089259. For example, non-isomorphic representatives of the a(1) = 1 through a(3) = 6 multiset partitions are:
{{1}} {{1,1}} {{1,1,1}}
{{1},{1}} {{1},{1,1}}
{{1},{2}} {{1},{2,2}}
{{1},{1},{1}}
{{1},{2},{2}}
{{1},{2},{3}}
A000688 counts factorizations into prime powers.
A007716 counts non-isomorphic multiset partitions by weight.
A279784 counts twice-partitions of type PPR, factorizations A295935.
Constant partitions are ranked by prime-powers: A000961, A023894, A054685, A246655, A355743.
(End)

Examples

			G.f. = 1 + x + 3*x^2 + 6*x^3 + 15*x^4 + 28*x^5 + 66*x^6 + 122*x^7 + ...
a(3) = 6 because we have (111) = (111) = (11)(1) = (1)(1)(1), (12) = (12) = (1)(2), (3) = (3).
The a(4)=14 multiset partitions whose total sum of parts is 4 are:
((4)),
((13)), ((1)(3)),
((22)), ((2)(2)),
((112)), ((1)(12)), ((2)(11)), ((1)(1)(2)),
((1111)), ((1)(111)), ((11)(11)), ((1)(1)(11)), ((1)(1)(1)(1)). - _Gus Wiseman_, Dec 19 2016
		

References

  • A. Cayley, Recherches sur les matrices dont les termes sont des fonctions linéaires d'une seule indéterminée, J. Reine angew. Math., 50 (1855), 313-317; Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 2, p. 219.
  • V. A. Liskovets, Counting rooted initially connected directed graphs. Vesci Akad. Nauk. BSSR, ser. fiz.-mat., No 5, 23-32 (1969), MR44 #3927.
  • 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).
  • J. J. Sylvester, An Enumeration of the Contacts of Lines and Surfaces of the Second Order, Phil. Mag. 1 (1851), 119-140. Reprinted in Collected Papers, Vol. 1. See p. 239, where one finds a(n)-2, but with errors.
  • J. J. Sylvester, Note on the 'Enumeration of the Contacts of Lines and Surfaces of the Second Order', Phil. Mag., Vol. VII (1854), pp. 331-334. Reprinted in Collected Papers, Vol. 2, pp. 30-33.

Crossrefs

Related to A001383 via generating function.
The multiplicative version (factorizations) is A050336.
The ordered version (sequences of partitions) is A055887.
Row-sums of A061260.
Main diagonal of A055885.
We have A271619(n) <= a(n) <= A063834(n).
Column k=3 of A290353.
The strict case is A316980.
Cf. A089300.

Programs

  • Haskell
    Following Vladeta Jovovic:
    a001970 n = a001970_list !! (n-1)
    a001970_list = 1 : f 1 [1] where
       f x ys = y : f (x + 1) (y : ys) where
                y = sum (zipWith (*) ys a061259_list) `div` x
    -- Reinhard Zumkeller, Oct 31 2015
    
  • Maple
    with(combstruct); SetSetSetU := [T, {T=Set(S), S=Set(U,card >= 1), U=Set(Z,card >=1)},unlabeled];
    # second Maple program:
    with(numtheory): with(combinat):
    a:= proc(n) option remember; `if`(n=0, 1, add(add(d*
          numbpart(d), d=divisors(j))*a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..35);  # Alois P. Heinz, Dec 19 2016
  • Mathematica
    m = 32; f[x_] = Product[1/(1-x^k)^PartitionsP[k], {k, 1, m}]; CoefficientList[ Series[f[x], {x, 0, m-1}], x] (* Jean-François Alcover, Jul 19 2011, after g.f. *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / prod(k=1, n, 1 - numbpart(k) * x^k + x * O(x^n)), n))}; /* Michael Somos, Dec 20 2016 */
    
  • Python
    from sympy.core.cache import cacheit
    from sympy import npartitions, divisors
    @cacheit
    def a(n): return 1 if n == 0 else sum([sum([d*npartitions(d) for d in divisors(j)])*a(n - j) for j in range(1, n + 1)]) / n
    [a(n) for n in range(51)]  # Indranil Ghosh, Aug 19 2017, after Maple code
    # (Sage) # uses[EulerTransform from A166861]
    b = BinaryRecurrenceSequence(0, 1, 1)
    a = EulerTransform(EulerTransform(b))
    print([a(n) for n in range(36)]) # Peter Luschny, Nov 17 2022

Formula

G.f.: Product_{k >= 1} 1/(1-x^k)^p(k), where p(k) = number of partitions of k = A000041. [Cayley]
a(n) = (1/n)*Sum_{k = 1..n} a(n-k)*b(k), n > 1, a(0) = 1, b(k) = Sum_{d|k} d*numbpart(d), where numbpart(d) = number of partitions of d, cf. A061259. - Vladeta Jovovic, Apr 21 2001
Logarithmic derivative yields A061259 (equivalent to above formula from Vladeta Jovovic). - Paul D. Hanna, Sep 05 2012
a(n) = Sum_{k=1..A000041(n)} A001055(A215366(n,k)) = number of factorizations of Heinz numbers of integer partitions of n. - Gus Wiseman, Dec 19 2016
a(n) = |{m>=1 : n = Sum_{k=1..A001222(m)} A056239(A112798(m,k)+1)}| = number of normalized twice-prime-factored multiset partitions (see A275024) whose total sum of parts is n. - Gus Wiseman, Dec 19 2016

Extensions

Additional comments from Valery A. Liskovets
Sylvester references from Barry Cipra, Oct 07 2003

A302242 Total weight of the n-th multiset multisystem. Totally additive with a(prime(n)) = Omega(n).

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 2, 0, 2, 1, 1, 1, 2, 2, 2, 0, 1, 2, 3, 1, 3, 1, 2, 1, 2, 2, 3, 2, 2, 2, 1, 0, 2, 1, 3, 2, 3, 3, 3, 1, 1, 3, 2, 1, 3, 2, 2, 1, 4, 2, 2, 2, 4, 3, 2, 2, 4, 2, 1, 2, 3, 1, 4, 0, 3, 2, 1, 1, 3, 3, 3, 2, 2, 3, 3, 3, 3, 3, 2, 1, 4, 1, 1, 3, 2, 2, 3, 1, 4
Offset: 1

Views

Author

Gus Wiseman, Apr 03 2018

Keywords

Comments

A multiset multisystem is a finite multiset of finite multisets of positive integers. The n-th multiset multisystem is constructed by factoring n into prime numbers and then factoring each prime index into prime numbers and taking their prime indices. This produces a unique multiset multisystem for each n, and every possible multiset multisystem is so constructed as n ranges over all positive integers.

Examples

			Sequence of finite multisets of finite multisets of positive integers begins: (), (()), ((1)), (()()), ((2)), (()(1)), ((11)), (()()()), ((1)(1)), (()(2)), ((3)), (()()(1)), ((12)), (()(11)), ((1)(2)), (()()()()), ((4)), (()(1)(1)), ((111)), (()()(2)).
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    a:= n-> add(add(j[2], j=ifactors(pi(i[1]))[2])*i[2], i=ifactors(n)[2]):
    seq(a(n), n=1..100);  # Alois P. Heinz, Sep 07 2018
  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Total[PrimeOmega/@primeMS[n]],{n,100}]
  • PARI
    a(n,f=factor(n))=sum(i=1,#f~, bigomega(primepi(f[i,1]))*f[i,2]) \\ Charles R Greathouse IV, Nov 10 2021

A283877 Number of non-isomorphic set-systems of weight n.

Original entry on oeis.org

1, 1, 2, 4, 9, 18, 44, 98, 244, 605, 1595, 4273, 12048, 34790, 104480, 322954, 1031556, 3389413, 11464454, 39820812, 141962355, 518663683, 1940341269, 7424565391, 29033121685, 115921101414, 472219204088, 1961177127371, 8298334192288, 35751364047676, 156736154469354
Offset: 0

Views

Author

Gus Wiseman, Mar 17 2017

Keywords

Comments

A set-system is a finite set of finite nonempty sets. The weight of a set-system is the sum of cardinalities of its elements.

Examples

			Non-isomorphic representatives of the a(4)=9 set-systems are:
((1234)),
((1)(234)), ((3)(123)), ((12)(34)), ((13)(23)),
((1)(2)(12)), ((1)(2)(34)), ((1)(3)(23)),
((1)(2)(3)(4)).
		

Crossrefs

Programs

  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(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)={WeighT(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(g=sum(t=1, n, subst(x*Ser(K(q, t, n\t)/t),x,x^t) )); s+=permcount(q)*polcoef(exp(g - subst(g,x,x^2)), n)); s/n!)} \\ Andrew Howroyd, Jan 16 2024

Formula

Euler transform of A300913.

Extensions

a(0) = 1 prepended and terms a(11) and beyond from Andrew Howroyd, Sep 01 2019

A049311 Number of (0,1) matrices with n ones and no zero rows or columns, up to row and column permutations.

Original entry on oeis.org

1, 3, 6, 16, 34, 90, 211, 558, 1430, 3908, 10725, 30825, 90156, 273234, 848355, 2714399, 8909057, 30042866, 103859678, 368075596, 1335537312, 4958599228, 18820993913, 72980867400, 288885080660, 1166541823566, 4802259167367, 20141650236664
Offset: 1

Views

Author

Keywords

Comments

Also the number of bipartite graphs with n edges, no isolated vertices and a distinguished bipartite block, up to isomorphism.
The EULERi transform (A056156) is also interesting.
a(n) is also the number of non-isomorphic set multipartitions (multisets of sets) of weight n. - Gus Wiseman, Mar 17 2017

Examples

			E.g. a(2) = 3: two ones in same row, two ones in same column, or neither.
a(3) = 6 is coefficient of x^3 in (1/36)*((1 + x)^9 + 6*(1 + x)^3*(1 + x^2)^3 + 8*(1 + x^3)^3 + 9*(1 + x)*(1 + x^2)^4 + 12*(1 + x^3)*(1 + x^6))=1 + x + 3*x^2 + 6*x^3 + 7*x^4 + 7*x^5 + 6*x^6 + 3*x^7 + x^8 + x^9.
There are a(3) = 6 binary matrices with 3 ones, with no zero rows or columns, up to row and column permutation:
  [1 0 0] [1 1 0] [1 0] [1 1] [1 1 1] [1]
  [0 1 0] [0 0 1] [1 0] [1 0] ....... [1].
  [0 0 1] ....... [0 1] ............. [1]
Non-isomorphic representatives of the a(3)=6 set multipartitions are: ((123)), ((1)(23)), ((2)(12)), ((1)(1)(1)), ((1)(2)(2)), ((1)(2)(3)). - _Gus Wiseman_, Mar 17 2017
		

Crossrefs

Programs

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

Formula

Calculate number of connected bipartite graphs + number of connected bipartite graphs with no duality automorphism, then apply EULER transform.
a(n) is the coefficient of x^n in the cycle index Z(S_n X S_n; 1+x, 1+x^2, ...), where S_n X S_n is Cartesian product of symmetric groups S_n of degree n.

Extensions

More terms and formula from Vladeta Jovovic, Jul 29 2000
a(19)-a(28) from Max Alekseyev, Jul 22 2009
a(29)-a(102) from Aliaksandr Siarhei, Dec 13 2013
Name edited by Gus Wiseman, Dec 18 2018

A000612 Number of P-equivalence classes of switching functions of n or fewer variables, divided by 2.

Original entry on oeis.org

1, 2, 6, 40, 1992, 18666624, 12813206169137152, 33758171486592987164087845043830784, 1435913805026242504952006868879460423834904914948818373264705576411070464
Offset: 0

Views

Author

Keywords

Comments

Also number of nonisomorphic sets of nonempty subsets of an n-set.
Equivalently, number of nonisomorphic fillings of a Venn diagram of n sets. - Joerg Arndt, Mar 24 2020
Number of hypergraphs on n unlabeled nodes. - Charles R Greathouse IV, Apr 06 2021

Examples

			Non-isomorphic representatives of the a(2) = 6 set-systems are 0, {1}, {12}, {1}{2}, {1}{12}, {1}{2}{12}. - _Gus Wiseman_, Aug 07 2018
		

References

  • M. A. Harrison, Introduction to Switching and Automata Theory. McGraw Hill, NY, 1965, p. 153.
  • S. Muroga, Threshold Logic and Its Applications. Wiley, NY, 1971, p. 38 Table 2.3.2. - Row 5.
  • 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

Programs

  • Maple
    a:= n-> add(1/(p-> mul((c-> j^c*c!)(coeff(p, x, j)), j=1..degree(p)))(
            add(x^i, i=l))*2^((w-> add(mul(2^igcd(t, l[i]), i=1..nops(l)),
            t=1..w)/w)(ilcm(l[]))), l=combinat[partition](n))/2:
    seq(a(n), n=0..9);  # Alois P. Heinz, Aug 12 2019
  • Mathematica
    sysnorm[{}] := {};sysnorm[m_]:=If[Union@@m!=Range[Max@@Flatten[m]],sysnorm[m/.Rule@@@Table[{(Union@@m)[[i]],i},{i,Length[Union@@m]}]],First[Sort[sysnorm[m,1]]]];sysnorm[m_,aft_]:=If[Length[Union@@m]<=aft,{m},With[{mx=Table[Count[m,i,{2}],{i,Select[Union@@m,#>=aft&]}]},Union@@(sysnorm[#,aft+1]&/@Union[Table[Map[Sort,m/.{par+aft-1->aft,aft->par+aft-1},{0,1}],{par,First/@Position[mx,Max[mx]]}]])]];
    Table[Length[Union[sysnorm/@Subsets[Rest[Subsets[Range[n]]]]]],{n,4}] (* Gus Wiseman, Aug 07 2018 *)
    a[n_] := Sum[1/Function[p, Product[Function[c, j^c*c!][Coefficient[p, x, j]], {j, 1, Exponent[p, x]}]][Total[x^l]]*2^(Function[w, Sum[Product[2^GCD[t, l[[i]]], {i, 1, Length[l]}], {t, 1, w}]/w][If[l=={}, 1, LCM @@ l]]), {l, IntegerPartitions[n]}]/2;
    a /@ Range[0, 9] (* Jean-François Alcover, Feb 04 2020, after Alois P. Heinz *)
  • Python
    def partition(n, I=1):
      yield () if n==0 else (n,)
      for i in range(I, n//2 + 1):
        for p in partition(n-i, i):
          yield (i,) + p
    def a(n):
      import math, operator, functools
      fracs = [(1<<(sum(functools.reduce(operator.mul, (1<Gregory Morse, Dec 23 2024

Formula

a(n) = A003180(n)/2.

Extensions

More terms from Vladeta Jovovic, Feb 23 2000

A007718 Number of independent polynomial invariants of matrix of order n.

Original entry on oeis.org

1, 1, 3, 6, 17, 40, 125, 354, 1159, 3774, 13113, 46426, 171027, 644038, 2493848, 9867688, 39922991, 164747459, 693093407, 2968918400, 12940917244, 57353242370, 258306634422, 1181572250326, 5486982683756, 25856584485254
Offset: 0

Views

Author

Keywords

Comments

Also the number of non-isomorphic connected multiset partitions of weight n. The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices. - Gus Wiseman, Sep 23 2018

Examples

			From _Gus Wiseman_, Sep 24 2018: (Start)
Non-isomorphic representatives of the a(1) = 1 through a(4) = 17 connected multiset partitions:
  {{1}}  {{1,1}}    {{1,1,1}}      {{1,1,1,1}}
         {{1,2}}    {{1,2,2}}      {{1,1,2,2}}
        {{1},{1}}   {{1,2,3}}      {{1,2,2,2}}
                   {{1},{1,1}}     {{1,2,3,3}}
                   {{2},{1,2}}     {{1,2,3,4}}
                  {{1},{1},{1}}   {{1},{1,1,1}}
                                  {{1},{1,2,2}}
                                  {{2},{1,2,2}}
                                  {{3},{1,2,3}}
                                  {{1,1},{1,1}}
                                  {{1,2},{1,2}}
                                  {{1,2},{2,2}}
                                  {{1,3},{2,3}}
                                 {{1},{1},{1,1}}
                                 {{1},{2},{1,2}}
                                 {{2},{2},{1,2}}
                                {{1},{1},{1},{1}}
(End)
		

Crossrefs

Formula

Inverse Euler transform of A007716.

Extensions

a(7)-a(25) from Franklin T. Adams-Watters, Jun 21 2011
a(0)=1 prepended by Andrew Howroyd, Jan 15 2023

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

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 18 2018

Keywords

Comments

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.

Examples

			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)
		

Crossrefs

Programs

  • PARI
    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

Formula

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

Extensions

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

Views

Author

Gus Wiseman, Jul 18 2018

Keywords

Comments

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.

Examples

			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]
		

Crossrefs

Row sums of A320796.
Main diagonal of A318805.

Programs

Extensions

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

A255906 Number of collections of nonempty multisets with a total of n objects having color set {1,...,k} for some k<=n.

Original entry on oeis.org

1, 1, 4, 16, 76, 400, 2356, 15200, 106644, 806320, 6526580, 56231024, 513207740, 4941362512, 50013751812, 530481210672, 5880285873060, 67954587978448, 816935340368068, 10196643652651664, 131904973822724540, 1765645473517011568, 24420203895517396180
Offset: 0

Views

Author

Alois P. Heinz, Mar 10 2015

Keywords

Comments

Number of multiset partitions of normal multisets of size n, where a multiset is normal if it spans an initial interval of positive integers. - Gus Wiseman, Jul 30 2018

Examples

			a(0) = 1: {}.
a(1) = 1: {{1}}.
a(2) = 4: {{1},{1}}, {{1,1}}, {{1},{2}}, {{1,2}}.
a(3) = 16: {{1},{1},{1}}, {{1},{1,1}}, {{1,1,1}}, {{1},{1},{2}}, {{1},{2},{2}}, {{1},{1,2}}, {{1},{2,2}}, {{2},{1,1}}, {{2},{1,2}}, {{1,1,2}}, {{1,2,2}}, {{1},{2},{3}}, {{1},{2,3}}, {{2},{1,3}}, {{3},{1,2}}, {{1,2,3}}.
		

Crossrefs

Row sums of A255903. Also row sums of A317532.

Programs

  • Mathematica
    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]]]];
    allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    Table[Length[Join@@mps/@allnorm[n]],{n,6}] (* Gus Wiseman, Jul 30 2018 *)
  • PARI
    R(n, k)={Vec(-1 + 1/prod(j=1, n, (1 - x^j + O(x*x^n))^binomial(k+j-1, j) ))}
    seq(n) = {concat([1], sum(k=1, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) ))} \\ Andrew Howroyd, Sep 23 2023

Formula

a(n) = Sum_{k=0..n} A255903(n,k).
Showing 1-10 of 466 results. Next