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

A045779 Number of factorizations of n into distinct factors for some n (image of A045778).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 17, 18, 19, 21, 22, 25, 27, 31, 32, 33, 34, 38, 40, 42, 43, 44, 46, 52, 54, 55, 56, 57, 59, 61, 64, 67, 70, 74, 76, 80, 83, 88, 89, 91, 93, 100, 104, 110, 111, 112, 116, 117, 120, 122, 123, 132, 137, 140, 141, 142, 143, 148
Offset: 1

Views

Author

Keywords

Comments

We may use A045778(k*m) >= A045778(k) for any k, m >= 1 to disprove presence of some positive integer in this sequence. - David A. Corneth, Oct 24 2024

Examples

			From _David A. Corneth_, Oct 24 2024: (Start)
5 is a term as 24 has five factorizations into distinct divisors of 24 namely 24 = 2 * 12 = 3 * 8 = 4 * 6 = 2 * 3 * 4 which is five such factorizations.
11 is not a term. From terms in A025487 only the numbers 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 72, 96, 128, 256, 512, 1024 have no more than 11 such factorizations. Any multiple of these numbers in A025487 that is not already listed has more than 11 such factorizations which proves 11 is not in this sequence. (End)
		

Crossrefs

Factorizations are A001055, with image A045782, with complement A330976.
Strict factorizations are A045778 with image A045779 and complement A330975.
The least number with A045779(n) strict factorizations is A045780(n).
The least number with n strict factorizations is A330974(n).

Extensions

Name edited by Gus Wiseman, Jan 11 2020

A331023 Numerator: factorizations divided by strict factorizations A001055(n)/A045778(n).

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, 9, 1, 1, 1, 7, 1, 1, 1, 4, 4, 1, 1, 12, 2, 4, 1, 4, 1, 7, 1, 7, 1, 1, 1, 11, 1, 1, 4, 11, 1, 1, 1, 4, 1, 1, 1, 16, 1, 1, 4, 4, 1, 1, 1, 12, 5, 1, 1, 11, 1, 1, 1, 7, 1, 11, 1, 4, 1, 1, 1, 19, 1, 4, 4, 9, 1, 1, 1, 7, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 08 2020

Keywords

Comments

A factorization of n is a finite, nondecreasing sequence of positive integers > 1 with product n. It is strict if the factors are all different. Factorizations and strict factorizations are counted by A001055 and A045778 respectively.

Crossrefs

Positions of 1's are A005117.
Positions of 2's appear to be A001248.
The denominators are A331024.
The rounded quotients are A331048.
The same for integer partitions is A330994.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[facs[n]]/Length[Select[facs[n],UnsameQ@@#&]],{n,100}]//Numerator
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
    A045778(n, m=n) = ((n<=m) + sumdiv(n, d, if((d>1)&&(d<=m)&&(dA045778(n/d, d-1))));
    A331023(n) = numerator(A001055(n)/A045778(n)); \\ Antti Karttunen, May 27 2021

Formula

a(2^n) = A330994(n).

Extensions

More terms from Antti Karttunen, May 27 2021

A331024 Denominator: factorizations divided by strict factorizations A001055(n)/A045778(n).

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, 5, 1, 1, 1, 5, 1, 1, 1, 3, 3, 1, 1, 7, 1, 3, 1, 3, 1, 5, 1, 5, 1, 1, 1, 9, 1, 1, 3, 4, 1, 1, 1, 3, 1, 1, 1, 9, 1, 1, 3, 3, 1, 1, 1, 7, 2, 1, 1, 9, 1, 1, 1, 5, 1, 9, 1, 3, 1, 1, 1, 10, 1, 3, 3, 5, 1, 1, 1, 5, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 08 2020

Keywords

Comments

A factorization of n is a finite, nondecreasing sequence of positive integers > 1 with product n. It is strict if the factors are all different. Factorizations and strict factorizations are counted by A001055 and A045778 respectively.

Crossrefs

Positions of 1's include all elements of A001248 as well as A005117. The first position of a 1 that is not in A167207 is 128.
The numerators are A331023.
The rounded quotients are A331048.
The same for integer partitions is A330995.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[facs[n]]/Length[Select[facs[n],UnsameQ@@#&]],{n,100}]//Denominator
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
    A045778(n, m=n) = ((n<=m) + sumdiv(n, d, if((d>1)&&(d<=m)&&(dA045778(n/d, d-1))));
    A331024(n) = denominator(A001055(n)/A045778(n)); \\ Antti Karttunen, May 27 2021

Formula

a(2^n) = A330995(n).

Extensions

More terms from Antti Karttunen, May 27 2021

A331048 Nearest integer to A001055(n)/A045778(n), where A001055 is factorizations and A045778 is strict factorizations.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 10 2020

Keywords

Comments

A factorization of n is a finite, nondecreasing sequence of positive integers > 1 with product n. It is strict if the factors are all different.

Crossrefs

The exact quotient is A331023/A331024.
The same for integer partitions is A330996 ~ A330994/A330995.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Round[Length[facs[n]]/Length[Select[facs[n],UnsameQ@@#&]]],{n,100}]

A050323 Number of factorizations into distinct factors indexed by prime signatures: A045778(A025487).

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 2, 5, 5, 3, 5, 7, 9, 4, 9, 10, 16, 5, 14, 18, 14, 15, 17, 25, 6, 21, 34, 19, 31, 27, 38, 8, 31, 57, 25, 59, 42, 40, 56, 10, 67, 44, 70, 46, 91, 33, 100, 64, 80, 80, 12, 117, 61, 52, 141, 74, 140, 43, 161, 93, 143, 111, 15, 193, 83, 120, 253, 116, 167, 208
Offset: 1

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Examples

			a(8) = 5 as A025487(8) = 24 and there are 5 ways to write 24 as a product of distinct divisors of 24 namely 24 = 2 * 12 = 3 * 8 = 4 * 6 = 2 * 3 * 4. - _David A. Corneth_, Oct 24 2024
		

Crossrefs

Programs

A331198 Numbers n with exactly three times as many factorizations (A001055) as strict factorizations (A045778).

Original entry on oeis.org

128, 2187, 10368, 34992, 78125, 80000, 307328, 823543, 1250000, 1366875, 1874048, 3655808, 5250987, 6328125, 10690688, 13176688, 16681088, 19487171, 32019867, 35819648, 62462907, 62748517, 66706983, 90531968, 118210688, 182660427, 187578125, 239892608, 285012027
Offset: 1

Views

Author

Gus Wiseman, Jan 12 2020

Keywords

Comments

Contains p^7 for all primes p.

Examples

			The 15 factorizations and 5 strict factorizations of 2187:
  (2187)           (2187)
  (27*81)          (27*81)
  (3*729)          (3*729)
  (9*243)          (9*243)
  (3*9*81)         (3*9*81)
  (9*9*27)
  (3*27*27)
  (3*3*243)
  (3*9*9*9)
  (3*3*3*81)
  (3*3*9*27)
  (3*3*3*9*9)
  (3*3*3*3*27)
  (3*3*3*3*3*9)
  (3*3*3*3*3*3*3)
		

Crossrefs

Factorizations are A001055.
Strict factorizations are A045778.
Taking "twice" instead of "three times" gives A001248.

Programs

  • Mathematica
    facsm[n_]:=facsm[n]=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsm[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Select[Range[100000],3==Length[facsm[#]]/Length[Select[facsm[#],UnsameQ@@#&]]&]

Extensions

a(7)-(10) from Alois P. Heinz, Jan 17 2020
a(11)-a(29) from Giovanni Resta, Jan 20 2020

A001055 The multiplicative partition function: number of ways of factoring n with all factors greater than 1 (a(1) = 1 by convention).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 5, 1, 4, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 5, 1, 7, 2, 2, 2, 9, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 11, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 5, 2, 1, 11, 2, 2, 2, 7, 1, 11, 2, 4, 2, 2, 2, 19, 1, 4, 4, 9, 1, 5, 1
Offset: 1

Views

Author

Keywords

Comments

From David W. Wilson, Feb 28 2009: (Start)
By a factorization of n we mean a multiset of integers > 1 whose product is n.
For example, 6 is the product of 2 such multisets, {2, 3} and {6}, so a(6) = 2.
Similarly 8 is the product of 3 such multisets, {2, 2, 2}, {2, 4} and {8}, so a(8) = 3.
1 is the product of 1 such multiset, namely the empty multiset {}, whose product is by definition the multiplicative identity 1. Hence a(1) = 1. (End)
a(n) = # { k | A064553(k) = n }. - Reinhard Zumkeller, Sep 21 2001; Benoit Cloitre and N. J. A. Sloane, May 15 2002
Number of members of A025487 with n divisors. - Matthew Vandermast, Jul 12 2004
See sequence A162247 for a list of the factorizations of n and a program for generating the factorizations for any n. - T. D. Noe, Jun 28 2009
So a(n) gives the number of different prime signatures that can be found among the integers that have n divisors. - Michel Marcus, Nov 11 2015
For n > 0, also the number of integer partitions of n with product n, ranked by A301987. For example, the a(12) = 4 partitions are: (12), (6,2,1,1,1,1), (4,3,1,1,1,1,1), (3,2,2,1,1,1,1,1). See also A380218. In general, A379666(m,n) = a(n) for any m >= n. - Gus Wiseman, Feb 07 2025

Examples

			1: 1, a(1) = 1
2: 2, a(2) = 1
3: 3, a(3) = 1
4: 4 = 2*2, a(4) = 2
6: 6 = 2*3, a(6) = 2
8: 8 = 2*4 = 2*2*2, a(8) = 3
etc.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 292-295.
  • Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 1.4.
  • 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).
  • G. Tenenbaum, Introduction to analytic and probabilistic number theory, Cambridge University Press, 1995, p. 198, exercise 9 (in the third edition 2015, p. 296, exercise 211).

Crossrefs

A045782 gives the range of a(n).
For records see A033833, A033834.
Row sums of A316439 (for n>1).
Cf. A096276 (partial sums).
The additive version is A000041 (integer partitions), strict A000009.
Row sums of A318950.
A002865 counts partitions into parts > 1.
A069016 counts distinct sums of factorizations.
A319000 counts partitions by product and sum, row sums A319916.
A379666 (transpose A380959) counts partitions by sum and product, without 1's A379668, strict A379671.

Programs

  • Haskell
    a001055 = (map last a066032_tabl !!) . (subtract 1)
    -- Reinhard Zumkeller, Oct 01 2012
    
  • Java
    public class MultiPart {
        public static void main(String[] argV) {
            for (int i=1;i<=100;++i) System.out.println(1+getDivisors(2,i));
        }
        public static int getDivisors(int min,int n) {
            int total = 0;
            for (int i=min;i=i) { ++total; if (n/i>i) total+=getDivisors(i,n/i); }
            return total;
        }
    } \\ Scott R. Shannon, Aug 21 2019
  • Maple
    with(numtheory):
    T := proc(n::integer, m::integer)
            local A, summe, d:
            if isprime(n) then
                    if n <= m then
                            return 1;
                    end if:
                    return 0 ;
            end if:
            A := divisors(n) minus {n, 1}:
            for d in A do
                    if d > m then
                            A := A minus {d}:
                    end if:
            end do:
            summe := add(T(n/d,d),d=A) ;
            if n <=m then
                    summe := summe + 1:
            end if:
            summe ;
    end proc:
    A001055 := n -> T(n, n):
    [seq(A001055(n), n=1..100)]; # Reinhard Zumkeller and Ulrich Schimke (ulrschimke(AT)aol.com)
  • Mathematica
    c[1, r_] := c[1, r]=1; c[n_, r_] := c[n, r] = Module[{ds, i}, ds = Select[Divisors[n], 1 < # <= r &]; Sum[c[n/ds[[i]], ds[[i]]], {i, 1, Length[ds]}]]; a[n_] := c[n, n]; a/@Range[100] (* c[n, r] is the number of factorizations of n with factors <= r. - Dean Hickerson, Oct 28 2002 *)
    T[, 1] = T[1, ] = 1;
    T[n_, m_] := T[n, m] = DivisorSum[n, Boole[1 < # <= m] * T[n/#, #]&];
    a[n_] := T[n, n];
    a /@ Range[100] (* Jean-François Alcover, Jan 03 2020 *)
  • PARI
    /* factorizations of n with factors <= m (n,m positive integers) */
    fcnt(n,m) = {local(s);s=0;if(n == 1,s=1,fordiv(n,d,if(d > 1 & d <= m,s=s+fcnt(n/d,d))));s}
    A001055(n) = fcnt(n,n) \\ Michael B. Porter, Oct 29 2009
    
  • PARI
    \\ code using Dirichlet g.f., based on Somos's code for A007896
    {a(n) = my(A, v, w, m);
    if(
    n<1, 0,
    \\ define unit vector v = [1, 0, 0, ...] of length n
    v = vector(n, k, k==1);
    for(k=2, n,
    m = #digits(n, k) - 1;
    \\ expand 1/(1-x)^k out far enough
    A = (1 - x)^ -1 + x * O(x^m);
    \\ w = zero vector of length n
    w = vector(n);
    \\ convert A to a vector
    for(i=0, m, w[k^i] = polcoeff(A, i));
    \\ build the answer
    v = dirmul(v, w)
    );
    v[n]
    )
    };
    \\ produce the sequence
    vector(100,n,a(n)) \\ N. J. A. Sloane, May 26 2014
    
  • PARI
    v=vector(100, k, k==1); for(n=2, #v, v+=dirmul(v, vector(#v, k, (k>1) && n^valuation(k,n)==k)) ); v \\ Max Alekseyev, Jul 16 2014
    
  • Python
    from sympy import divisors, isprime
    def T(n, m):
        if isprime(n): return 1 if n<=m else 0
        A=filter(lambda d: d<=m, divisors(n)[1:-1])
        s=sum(T(n//d, d) for d in A)
        return s + 1 if n<=m else s
    def a(n): return T(n, n)
    print([a(n) for n in range(1, 106)]) # Indranil Ghosh, Aug 19 2017
    

Formula

The asymptotic behavior of this sequence was studied by Canfield, Erdős & Pomerance and Luca, Mukhopadhyay, & Srinivas. - Jonathan Vos Post, Jul 07 2008
Dirichlet g.f.: Product_{k>=2} 1/(1 - 1/k^s).
If n = p^k for a prime p, a(n) = partitions(k) = A000041(k).
Since the sequence a(n) is the right diagonal of A066032, the given recursive formula for A066032 applies (see Maple program). - Reinhard Zumkeller and Ulrich Schimke (ulrschimke(AT)aol.com)
a(A002110(n)) = A000110(n).
a(p^k*q^k) = A002774(k) if p and q are distinct primes. - R. J. Mathar, Jun 06 2024
a(n) = A028422(n) + 1. - Gus Wiseman, Feb 07 2025

Extensions

Incorrect assertion about asymptotic behavior deleted by N. J. A. Sloane, Jun 08 2009

A050320 Number of ways n is a product of squarefree numbers > 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 3, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 6, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 3, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 6, 2, 2, 2, 2, 1, 6, 2, 2, 2, 2, 2, 2, 1, 2, 2, 3, 1, 5, 1
Offset: 1

Views

Author

Christian G. Bower, Sep 15 1999

Keywords

Comments

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).
Broughan shows (Theorem 8) that the average value of a(n) is k exp(2*sqrt(log n)/sqrt(zeta(2)))/log(n)^(3/4) where k is about 0.18504. - Charles R Greathouse IV, May 21 2013
From Gus Wiseman, Aug 20 2020: (Start)
Also the number of set multipartitions (multisets of sets) of the multiset of prime indices of n. For example, the a(n) set multipartitions for n = 2, 6, 36, 60, 360 are:
{1} {12} {12}{12} {1}{123} {1}{12}{123}
{1}{2} {1}{2}{12} {12}{13} {12}{12}{13}
{1}{1}{2}{2} {1}{1}{23} {1}{1}{12}{23}
{1}{2}{13} {1}{1}{2}{123}
{1}{3}{12} {1}{2}{12}{13}
{1}{1}{2}{3} {1}{3}{12}{12}
{1}{1}{1}{2}{23}
{1}{1}{2}{2}{13}
{1}{1}{2}{3}{12}
{1}{1}{1}{2}{2}{3}
(End)

Examples

			For n = 36 we have three choices as 36 = 2*2*3*3 = 6*6 = 2*3*6 (but no factorizations with factors 4, 9, 12, 18 or 36 are allowed), thus a(36) = 3. - _Antti Karttunen_, Oct 21 2017
		

Crossrefs

Cf. A001055, A005117, A050325. a(p^k)=1. a(A002110)=A000110.
a(n!)=A103774(n).
Cf. A206778.
Differs from A259936 for the first time at n=36.
A050326 is the strict case.
A045778 counts strict factorizations.
A089259 counts set multipartitions of integer partitions.
A116540 counts normal set multipartitions.

Programs

  • Haskell
    a050320 n = h n $ tail $ a206778_row n where
       h 1 _          = 1
       h _ []         = 0
       h m fs'@(f:fs) =
         if f > m then 0 else if r > 0 then h m fs else h m' fs' + h m fs
         where (m', r) = divMod m f
    -- Reinhard Zumkeller, Dec 16 2013
  • Mathematica
    sub[w_, e_] := Block[{v = w}, v[[e]]--; v]; ric[w_, k_] := If[Max[w] == 0, 1, Block[{e, s, p = Flatten@Position[Sign@w, 1]}, s = Select[Prepend[#, First@p] & /@ Subsets[Rest@p], Total[1/2^#] <= k &]; Sum[ric[sub[w, e], Total[1/2^e]], {e, s}]]]; sig[w_] := sig[w] = ric[w, 1];  a[n_] := sig@ Sort[Last /@ FactorInteger[n]]; Array[a, 103] (* Giovanni Resta, May 21 2013 *)
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sqfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]]
    Table[Length[sqfacs[n]],{n,100}] (* Gus Wiseman, Aug 20 2020 *)

Formula

Dirichlet g.f.: Product_{n is squarefree and > 1} (1/(1-1/n^s)).
a(n) = A050325(A101296(n)). - R. J. Mathar, May 26 2017
a(n!) = A103774(n); a(A006939(n)) = A337072(n). - Gus Wiseman, Aug 20 2020

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
Showing 1-10 of 332 results. Next