login
Let each integer m (1 <= m <= n) be factorized as m = prime_m(1)*prime_m(2)*...*prime_m(bigomega(m)), with the primes sorted in nonincreasing order. Then a(n) is the number of values of m such that each prime_m(i) <= prime_n(i).
19

%I #22 Nov 07 2024 22:01:27

%S 1,2,3,3,4,5,5,4,6,7,6,7,7,9,9,5,8,9,9,10,12,11,10,9,10,13,10,13,11,

%T 14,12,6,15,15,14,12,13,17,18,13,14,19,15,16,16,19,16,11,15,16,21,19,

%U 17,14,18,17,24,21,18,19,19,23,22,7,22,24,20,22,27,23,21

%N Let each integer m (1 <= m <= n) be factorized as m = prime_m(1)*prime_m(2)*...*prime_m(bigomega(m)), with the primes sorted in nonincreasing order. Then a(n) is the number of values of m such that each prime_m(i) <= prime_n(i).

%C Equivalently, a(n) equals the number of values of m such that each value of A238689 T(m,k) <= A238689 T(n,k). (Since the prime factorization of 1 is the empty factorization, we consider each prime_1(i) not to be greater than prime_n(i) for all positive integers n.)

%C Suppose we say that n "covers" m iff both m and n are factorized as described in the sequence definition and each prime_m(i) <= prime_n(i). At least three sequences (A037019, A108951 and A181821) have the property that a(m) divides a(n) iff n "covers" m. These sequences are also divisibility sequences (i.e., sequences with the property that a(m) divides a(n) if m divides n), since any positive integer "covers" each of its divisors.

%C For any positive integers m and k, the following integer sequences (with n >= 0) are arithmetic progressions:

%C 1. The sequence b(n) = a(m*(2^n)).

%C 2. The sequence b(n) = a(m*(prime(n+k))) if prime(k) >= A006530(m).

%C Also, a(n) = the number of distinct prime signatures that occur among the divisors of any integer m such that A181819(m) = n and/or A238745(m) = n.

%C Number of skew partitions whose numerator has Heinz number n, where a skew partition is a pair y/v of integer partitions such that the diagram of v fits inside the diagram of y. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - _Gus Wiseman_, Feb 24 2018

%H Amiram Eldar, <a href="/A238690/b238690.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A085082(A108951(n)) = A085082(A181821(n)).

%F a(n) = a(A122111(n)).

%F a(prime(n)) = a(2^n) = n+1.

%F a((prime(n))^m) = a((prime(m))^n) = binomial(n+m, n).

%F a(A002110(n)) = A000108(n+1).

%F A000005(n) <= a(n) <= n.

%e The prime factorizations of integers 1 through 9, with prime factors sorted from largest to smallest:

%e 1 - the empty factorization (no prime factors)

%e 2 = 2

%e 3 = 3

%e 4 = 2*2

%e 5 = 5

%e 6 = 3*2

%e 7 = 7

%e 8 = 2*2*2

%e 9 = 3*3

%e To find a(9), we consider 9 = 3*3. There are 6 positive integers (1, 2, 3, 4, 6 and 9) which satisfy the following criteria:

%e 1) The largest prime factor, if one exists, is not greater than 3;

%e 2) The second-largest prime factor, if one exists, is not greater than 3;

%e 3) The total number of prime factors (counting repeated factors) does not exceed 2.

%e Therefore, a(9) = 6.

%e From _Gus Wiseman_, Feb 24 2018: (Start)

%e Heinz numbers of the a(15) = 9 partitions contained within the partition (32) are 1, 2, 3, 4, 5, 6, 9, 10, 15. The a(15) = 9 skew partitions are (32)/(), (32)/(1), (32)/(11), (32)/(2), (32)/(21), (32)/(22), (32)/(3), (32)/(31), (32)/(32).

%e Corresponding diagrams are:

%e o o o . o o . o o . . o . . o . . o . . . . . . . . .

%e o o o o . o o o . o . . o o . o . . (End)

%t undptns[y_]:=Select[Tuples[Range[0,#]&/@y],OrderedQ[#,GreaterEqual]&];

%t primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];

%t Table[Length[undptns[Reverse[primeMS[n]]]],{n,100}] (* _Gus Wiseman_, Feb 24 2018 *)

%Y Rearrangement of A115728, A115729 and A238746. A116473(n) is the number of times n appears in the sequence.

%Y Cf. A000041, A000085, A000720, A056239, A063834, A112798, A122111, A153452, A215366, A238689, A259478, A259480, A296150, A296188, A296561, A297388, A299925, A299926, A299966, A299967.

%K nonn,changed

%O 1,2

%A _Matthew Vandermast_, Apr 28 2014