A095705 -id:A095705 - cp's OEIS frontend

A035310 Let f(n) = number of ways to factor n = A001055(n); a(n) = sum of f(k) over all terms k in A025487 that have n factors.

1, 4, 12, 47, 170, 750, 3255, 16010, 81199, 448156, 2579626, 15913058, 102488024, 698976419, 4976098729, 37195337408, 289517846210, 2352125666883, 19841666995265, 173888579505200, 1577888354510786, 14820132616197925, 143746389756336173, 1438846957477988926

Offset: 1

Alford Arnold

Ways of partitioning an n-multiset with multiplicities some partition of n.

Number of multiset partitions of strongly normal multisets of size n, where a finite multiset is strongly normal if it covers an initial interval of positive integers with weakly decreasing multiplicities. The (weakly) normal version is A255906. - Gus Wiseman, Dec 31 2019

			a(3) = 12 because there are 3 terms in A025487 with 3 factors, namely 8, 12, 30; and f(8)=3, f(12)=4, f(30)=5 and 3+4+5 = 12.
From _Gus Wiseman_, Dec 31 2019: (Start)
The a(1) = 1 through a(3) = 12 multiset partitions of strongly normal multisets:
  {{1}}  {{1,1}}    {{1,1,1}}
         {{1,2}}    {{1,1,2}}
         {{1},{1}}  {{1,2,3}}
         {{1},{2}}  {{1},{1,1}}
                    {{1},{1,2}}
                    {{1},{2,3}}
                    {{2},{1,1}}
                    {{2},{1,3}}
                    {{3},{1,2}}
                    {{1},{1},{1}}
                    {{1},{1},{2}}
                    {{1},{2},{3}}
(End)

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

Cf. A025487, A000041, A000110, A035098, A080688.
Sequence A035341 counts the ordered cases. Tables A093936 and A095705 distribute the values; e.g. 81199 = 30 + 536 + 3036 + 6181 + 10726 + 11913 + 14548 + 13082 + 21147.
Cf. A035341, A093936, A095705.
Row sums of A317449.
The uniform case is A317584.
The case with empty intersection is A317755.
The strict case is A317775.
The constant case is A047968.
The set-system case is A318402.
The case of strict parts is A330783.
Multiset partitions of integer partitions are A001970.
Unlabeled multiset partitions are A007716.
Cf. A001055, A255906, A269134, A316652, A317654, A318284, A330467, A330471, A330475, A330675.

with(numtheory):
g:= proc(n, k) option remember;
      `if`(n>k, 0, 1) +`if`(isprime(n), 0,
      add(`if`(d>k, 0, g(n/d, d)), d=divisors(n) minus {1, n}))
    end:
b:= proc(n, i, l)
      `if`(n=0, g(mul(ithprime(t)^l[t], t=1..nops(l))$2),
      `if`(i<1, 0, add(b(n-i*j, i-1, [l[], i$j]), j=0..n/i)))
    end:
a:= n-> b(n$2, []):
seq(a(n), n=1..10);  # Alois P. Heinz, May 26 2013

g[n_, k_] := g[n, k] = If[n > k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d > k, 0, g[n/d, d]], {d, Divisors[n] ~Complement~ {1, n}}]]; b[n_, i_, l_] := If[n == 0, g[p = Product[Prime[t]^l[[t]], {t, 1, Length[l]}], p], If[i < 1, 0, Sum[b[n - i*j, i-1, Join[l, Array[i&, j]]], {j, 0, n/i}]]]; a[n_] := b[n, n, {}]; Table[Print[an = a[n]]; an, {n, 1, 13}] (* Jean-François Alcover, Dec 12 2013, after Alois P. Heinz *)

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
D(p, n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=EulerT(v)); Vec(1/prod(k=1, n, 1 - u[k]*x^k + O(x*x^n))-1, -n)/prod(i=1, #v, i^v[i]*v[i]!)}
seq(n)={my(s=0); forpart(p=n, s+=D(p,n)); s} \\ Andrew Howroyd, Dec 30 2020

from sympy.core.cache import cacheit
from sympy import divisors, isprime, prime
from operator import mul
@cacheit
def g(n, k):
    return (0 if n > k else 1) + (0 if isprime(n) else sum(g(n//d, d) for d in divisors(n)[1:-1] if d <= k))
@cacheit
def b(n, i, l):
    if n==0:
        p = reduce(mul, (prime(t + 1)**l[t] for t in range(len(l))))
        return g(p, p)
    else:
        return 0 if i<1 else sum([b(n - i*j, i - 1, l + [i]*j) for j in range(n//i + 1)])
def a(n):
    return b(n, n, [])
for n in range(1, 11): print(a(n)) # Indranil Ghosh, Aug 19 2017, after Maple code

More terms from Erich Friedman.
81199 from Alford Arnold, Mar 04 2008
a(10) from Alford Arnold, Mar 31 2008
a(10) corrected by Alford Arnold, Aug 07 2008
a(11)-a(13) from Alois P. Heinz, May 26 2013
a(14) from Alois P. Heinz, Sep 27 2014
a(15) from Alois P. Heinz, Jan 10 2015
Terms a(16) and beyond from Andrew Howroyd, Dec 30 2020

A035341 Sum of ordered factorizations over all prime signatures with n factors.

Original entry on oeis.org

1, 1, 5, 25, 173, 1297, 12225, 124997, 1492765, 19452389, 284145077, 4500039733, 78159312233, 1460072616929, 29459406350773, 634783708448137, 14613962109584749, 356957383060502945, 9241222160142506097, 252390723655315856437, 7260629936987794508973

Offset: 0

Alford Arnold

Let f(n) = number of ordered factorizations of n (A074206(n)); a(n) = sum of f(k) over all terms k in A025487 that have n factors.
When the unordered spectrum A035310 is so ordered the sequences A000041 A000070 ...A035098 A000110 yield A000079 A001792 ... A005649 A000670 respectively.
Row sums of A095705. - David Wasserman, Feb 22 2008
From Ludovic Schwob, Sep 23 2023: (Start)
a(n) is the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns, with weakly decreasing row sums. The a(3) = 25 matrices:
  [1 1 1] [1 2] [2 1] [3]
  .
  [1 1] [1 1] [1 1 0] [1 0 1] [0 1 1] [2] [0 2] [2 0]
  [1 0] [0 1] [0 0 1] [0 1 0] [1 0 0] [1] [1 0] [0 1]
  .
  [1] [1 0] [0 1] [1 0] [0 1] [1 0 0] [1 0 0] [0 1] [1 0]
  [1] [1 0] [0 1] [0 1] [1 0] [0 1 0] [0 0 1] [1 0] [0 1]
  [1] [0 1] [1 0] [1 0] [0 1] [0 0 1] [0 1 0] [1 0] [0 1]
  .
  [0 1 0] [0 1 0] [0 0 1] [0 0 1]
  [1 0 0] [0 0 1] [1 0 0] [0 1 0]
  [0 0 1] [1 0 0] [0 1 0] [1 0 0] (End)

			a(3) = 25 because there are 3 terms in A025487 with 3 factors, namely 8, 12, 30; and f(8)=4, f(12)=8, f(30)=13 and 4+8+13 = 25.

Alois P. Heinz, Table of n, a(n) for n = 0..250 (first 36 terms from David Wasserman)
Eric Weisstein's World of Mathematics, Perfect Partition

Cf. A005651, A025487, A035310, A255906, A365961.

Row sums of A261719.

b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
      add(b(n-i*j, i-1, k)*binomial(i+k-1, k-1)^j, j=0..n/i)))
    end:
a:= n->add(add(b(n$2, k-i)*(-1)^i*binomial(k, i), i=0..k), k=0..n):
seq(a(n), n=0..25);  # Alois P. Heinz, Aug 29 2015

b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 1, 0, Sum[b[n - i*j, i - 1, k]*If[j == 0, 1, Binomial[i + k - 1, k - 1]^j], {j, 0, n/i}]]];
a[n_] := Sum[Sum[b[n, n, k-i]*(-1)^i*Binomial[k, i], {i, 0, k}], {k, 0, n}];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Oct 26 2015, after Alois P. Heinz, updated Dec 15 2020 *)

R(n,k)=Vec(-1 + 1/prod(j=1, n, 1 - binomial(k+j-1,j)*x^j + O(x*x^n)))
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

a(n) ~ c * n! / log(2)^n, where c = 1/(2*log(2)) * Product_{k>=2} 1/(1-1/k!) = A247551 / (2*log(2)) = 1.8246323... . - Vaclav Kotesovec, Jan 21 2017

More terms from Erich Friedman.

More terms from David Wasserman, Feb 22 2008

A098348 Triangular array read by rows: a(n, k) = number of ordered factorizations of a "hook-type" number with n total prime factors and k distinct prime factors. "Hook-type" means that only one prime can have multiplicity > 1.

Original entry on oeis.org

1, 2, 3, 4, 8, 13, 8, 20, 44, 75, 16, 48, 132, 308, 541, 32, 112, 368, 1076, 2612, 4683, 64, 256, 976, 3408, 10404, 25988, 47293, 128, 576, 2496, 10096, 36848, 116180, 296564, 545835, 256, 1280, 6208, 28480, 120400, 454608, 1469892, 3816548

Offset: 1

Alford Arnold, Sep 04 2004

The first three columns are A000079, A001792 and A098385.
The first two diagonals are A000670 and A005649.
A070175 gives the smallest representative of each hook-type prime signature, so this sequence is a rearrangement of A074206(A070175).

			a(4, 2) = 20 because 24=2*2*2*3 has 20 ordered factorizations and so does any other number with the same prime signature.

Cf. A050324, A070175, A070826, A074206, A095705. A098349 gives the row sums. A098384.

a(n, k) = 1 + (Sum_{i=1..k-1} binomial(k-1, i)*a(i, i)) + (Sum_{j=1..k} Sum_{i=j..j+n-k-1} binomial(k-1, j-1)*a(i, j)) + (Sum_{j=1..k-1} binomial(k-1,j-1)*a(j+n-k, j)). - David Wasserman, Feb 21 2008
a(n, k) = A074206(2^(n+1-k)*A070826(k)). - David Wasserman, Feb 21 2008
The following conjectural formula for the triangle entries agrees with the values listed above: T(n,k) = Sum_{j = 0..n-k} 2^(n-k-j)*binomial(n-k,j)*a(k,j), where a(k,j) = 2^j*Sum_{i = j+1..k+1} binomial(i,j+1)*(i-1)!*Stirling2(k+1,i). See A098384 for related conjectures. - Peter Bala, Apr 20 2012

Edited and extended by David Wasserman, Feb 21 2008

A098349 Sum of ordered factorizations over hook-type prime signatures. (Row sums of A098348).

Original entry on oeis.org

1, 5, 25, 147, 1045, 8883, 88389, 1008723, 12984933, 186065011, 2936648325, 50611122451, 945623586725, 19037956164787, 410846314815941, 9460698944482643, 231538646070940901, 6001404380246870771, 164229726267373422853, 4731532138668375166355, 143154470312498479646245, 4538018949359464639487283, 150411282096757092041331781

Offset: 1

Alford Arnold, Sep 03 2004

The general case for all prime signatures is described by A095705 with row sums A035341.

			The array begins
1
2 3
4 8 13
8 20 44 75
16 48 132 308 541
...
therefore sequence begins 1 5 25 147 1045 ...

Cf. A035341, A095705, A098348.

More terms from David Wasserman, Feb 20 2008

A131420 A tabular sequence of arrays counting ordered factorizations over least prime signatures. The unordered version is described by sequence A129306.

Original entry on oeis.org

1, 2, 3, 4, 8, 13, 8, 20, 44, 75, 26, 16, 48, 132, 308, 541, 76, 176, 32, 112, 368, 1076, 2612, 4683, 208, 604, 1460, 252, 818, 64, 256, 876, 3408, 10404, 25988, 47293, 544, 1888, 5740, 14300, 768, 2316, 3172, 7880, 128, 576, 2496, 10096, 36848, 116180

Offset: 1

Alford Arnold, Jul 10 2007

nonn

The display has 1 2 3 5 7 11 15 ... terms per column. (cf. A000041)
The arrays begin
1.....2.....4......8......16.....32.....64......128
......3.....8.....20......48....112....256......576
...........13.....44.....132....368....976.....2496
..................75.....308...1076...3408....10096
.........................541...2612..10404....36848
...............................4683..25988...116180
.....................................47293...296564
.............................................545835
..................26......76....208....544
.........................176....604...1888
...............................1460...5740
.....................................14300
................................252....768
......................................2316
................................818...3172
......................................7880
with column sums
1....5....25....173....1297....12225....124997 => A035341
Column i corresponds to partitions of i. The rows correspond successively to the partitions {i}, {i-1,1},{i-2,1,1},{i-3,1,1,1}, ..., {i-7,1,1,1,1,1,1,1}, {i-2,2}, {i-3,2,1}, {i-4,2,1,1}, {i-5,2,1,1,1}, {i-3,3}, {i-3,3,1}, {i-4,2,2}, {i-5,2,2,1}. - Roger Lipsett, Feb 26 2016

			36 = 2*2*3*3 and is in A025487. There are 26 ways to factor 36 so a(11) = 26.

Cf. A000041, A000670, A002033, A025487, A035098, A035341, A050324, A074206, A095705, A098348, A104725, A108464, A129306.

gozinta counts ordered factorizations of an integer, and if lst is a partition we have
gozinta[1] = 1;
gozinta[n_] := gozinta[n] = 1 + Sum[gozinta[n/i], {i, Rest@Most@Divisors@n}]
a[lst_] := gozinta[Times @@ (Array[Prime, Length@lst]^lst)] (* Roger Lipsett, Feb 26 2016 *)

Corrected entries in table in comments section - Roger Lipsett, Feb 26 2016

cp's OEIS Frontend

A035310 Let f(n) = number of ways to factor n = A001055(n); a(n) = sum of f(k) over all terms k in A025487 that have n factors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Extensions

A035341 Sum of ordered factorizations over all prime signatures with n factors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A098348 Triangular array read by rows: a(n, k) = number of ordered factorizations of a "hook-type" number with n total prime factors and k distinct prime factors. "Hook-type" means that only one prime can have multiplicity > 1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Formula

Extensions

A098349 Sum of ordered factorizations over hook-type prime signatures. (Row sums of A098348).

Views

Author

Keywords

Comments

Examples

Crossrefs

Extensions

A131420 A tabular sequence of arrays counting ordered factorizations over least prime signatures. The unordered version is described by sequence A129306.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Extensions