OFFSET
0,2
COMMENTS
Equivalent to A074139 but using canonical order.
LINKS
Alois P. Heinz, Rows n = 0..30, flattened
S.-H. Cha, E. G. DuCasse, and L. V. Quintas, Graph Invariants Based on the Divides Relation and Ordered by Prime Signatures, arxiv:1405.5283
FORMULA
Trow(n) = List_{p in Partitions(n)} (Product_{t in p}(t + 1)). # Peter Luschny, Dec 11 2023
EXAMPLE
Triangle begins:
1;
2;
3, 4;
4, 6, 8;
5, 8, 9, 12, 16;
6, 10, 12, 16, 18, 24, 32;
7, 12, 15, 20, 16, 24, 32, 27, 36, 48, 64;
...
MAPLE
b:= (n, i)-> `if`(n=0 or i=1, [[1$n]], [map(x->
[i, x[]], b(n-i, min(n-i, i)))[], b(n, i-1)[]]):
T:= n-> map(x-> numtheory[tau](mul(ithprime(i)
^x[i], i=1..nops(x))), b(n$2))[]:
seq(T(n), n=0..9); # Alois P. Heinz, Mar 24 2020
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0 || i == 1, {Table[1, {n}]}, Join[Prepend[#, i]& /@ b[n - i, Min[n - i, i]], b[n, i - 1]]];
T[n_] := DivisorSigma[0, #]&[Product[Prime[i]^#[[i]], {i, 1, Length[#]}]& /@ b[n, n]];
T /@ Range[0, 9] // Flatten (* Jean-François Alcover, Jan 09 2025, after Alois P. Heinz *)
PROG
(PARI) \\ here b(n) is A000005.
b(n)={numdiv(n)}
N(sig)={prod(k=1, #sig, prime(k)^sig[k])}
Row(n)={apply(s->b(N(s)), vecsort([Vecrev(p) | p<-partitions(n)], , 4))}
{ for(n=0, 8, print(Row(n))) } \\ Andrew Howroyd, Mar 24 2020
(SageMath)
def A238963row(n):
return list(product(t + 1 for t in p) for p in Partitions(n))
print([A238963row(n) for n in range(10)]) # Peter Luschny, Dec 11 2023
CROSSREFS
KEYWORD
AUTHOR
Sung-Hyuk Cha, Mar 07 2014
EXTENSIONS
Offset corrected by Andrew Howroyd, Mar 24 2020
STATUS
approved