|
%I #9 Apr 19 2019 13:48:59
%S 0,1,1,1,2,1,2,1,1,2,2,2,1,1,1,2,1,3,1,1,3,2,3,1,1,1,3,1,2,3,1,1,2,3,
%T 1,1,4,2,4,1,1,1,4,1,2,4,1,1,2,4,1,1,1,2,4,1,2,2,4,1,1,3,4,1,2,5,1,1,
%U 2,5,1,1,1,2,5,1,2,2,5,1,1,3,5,1,1,2,2
%N Irregular triangle where row n is a list of indices in A002110 with multiplicity whose product is A004394(n).
%C Analogous to A306737.
%C The first 52 terms of a(n) and A306737 are identical, since the first 19 terms of A002182 and A004394 are the same, and the first two terms of row 20 are the same. a(20) = 4,2,1,1,1, while A306737(20) = 4,2,2.
%C Each superabundant number A004394(n) can be expressed as a product of primorials in A002110.
%C Row 1 = {0} by convention.
%C Maximum value in row n = A001221(A004394(n)).
%C Row n in reverse order is the conjugate of the list of the multiplicities of the prime divisors of A004394(n).
%H Michael De Vlieger, <a href="/A307322/b307322.txt">Table of n, a(n) for n = 1..10664</a> (rows 1 <= n <= 1200, flattened).
%H Michael De Vlieger, <a href="/A307322/a307322.txt">Relation between A307322 and A306737</a>.
%e Terms in the first rows n of this sequence, followed by the corresponding primorials whose product = A004394(n):
%e n T(n,k) A002110(T(n,k)) A004394(n)
%e -----------------------------------------------
%e 1: 0; 1 = 1
%e 2: 1; 2 = 2
%e 3: 1, 1; 2 * 2 = 4
%e 4: 2; 6 = 6
%e 5: 1, 2; 2 * 6 = 12
%e 6: 1, 1, 2; 2 * 2 * 6 = 24
%e 7: 2, 2; 6 * 6 = 36
%e 8: 1, 1, 1, 2; 2 * 2 * 2 * 6 = 48
%e 9: 1, 3; 2 * 30 = 60
%e 10: 1, 1, 3; 2 * 2 * 30 = 120
%e 11: 2, 3; 6 * 30 = 180
%e 12: 1, 1, 1, 3; 2 * 2 * 2 * 30 = 240
%e 13: 1, 2, 3; 2 * 6 * 30 = 360
%e 14: 1, 1, 2, 3; 2 * 2 * 6 * 30 = 720
%e 15: 1, 1, 4; 2 * 2 * 210 = 840
%e ...
%e Row 6 = {1,1,2} since A002110(1)*A002110(1)*A002110(2) = 2*2*6 = 24 and A004394(6) = 24. The conjugate of {1,1,2} = {3,1} and 24 = 2^3 * 3^1.
%e Row 10 = {1,1,3} since A002110(1)*A002110(1)*A002110(3) = 2*2*30 = 120 and A004394(10) = 120. The conjugate of {1,1,3} = {3,1,1} and 120 = 2^3 * 3^1 * 5^1.
%t Block[{s = Array[DivisorSigma[1, #]/# &, 10^6]}, Map[Table[LengthWhile[#, # >= i &], {i, Max@ #}] &@ If[# == 1, {0}, Function[f, ReplacePart[Table[0, {PrimePi[f[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, f]]@ FactorInteger@ #] &@ FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]] /. {} -> {0}] // Flatten
%Y Cf. A001221, A002110, A002182, A004394, A306737.
%K nonn,tabf
%O 1,5
%A _Michael De Vlieger_, Apr 02 2019
|
