A325272 -id:A325272 - cp's OEIS frontend

A336415 Number of divisors of n! with equal prime multiplicities.

1, 1, 2, 4, 6, 10, 13, 21, 24, 28, 33, 49, 53, 85, 94, 100, 104, 168, 173, 301, 307, 317, 334, 590, 595, 603, 636, 642, 652, 1164, 1171, 2195, 2200, 2218, 2283, 2295, 2301, 4349, 4478, 4512, 4519, 8615, 8626, 16818, 16836, 16844, 17101, 33485, 33491, 33507, 33516, 33582

Offset: 0

Gus Wiseman, Jul 22 2020

nonn

A number k has "equal prime multiplicities" (or is "uniform") iff its prime signature is constant, meaning that k is a power of a squarefree number.

			The a(n) uniform divisors of n for n = 1, 2, 6, 8, 30, 36 are the columns:
  1  2  6  8  30  36
     1  3  6  15  30
        2  4  10  16
        1  3   8  15
           2   6  10
           1   5   9
               4   8
               3   6
               2   5
               1   4
                   3
                   2
                   1
In 20!, the multiplicity of the third prime (5) is 4 but the multiplicity of the fourth prime (7) is 2. Hence there are 2^3 - 1 = 3 divisors with all exponents 3 (we subtract |{1}| = 1 from that count as 1 has no exponent 3). - _David A. Corneth_, Jul 27 2020

David A. Corneth, Table of n, a(n) for n = 0..9999
Jon Maiga, Computer-generated formulas for A336415, Sequence Machine.

The version for distinct prime multiplicities is A336414.
The version for nonprime perfect powers is A336416.
Uniform partitions are counted by A047966.
Uniform numbers are A072774, with nonprime terms A182853.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
Maximum divisor with distinct prime multiplicities is A327498.
Uniform divisors are counted by A327527.
Maximum uniform divisor is A336618.
1st differences are given by A048675.
Cf. A000005, A000961, A001222, A001597, A007916, A112798, A124010, A327499.
Factorial numbers: A000142, A007489, A022559, A027423, A048656, A071626, A108731, A325272, A325273, A325617.

Table[Length[Select[Divisors[n!],SameQ@@Last/@FactorInteger[#]&]],{n,0,15}]

a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); (#ex==0) || (vecmin(ex) == vecmax(ex))); \\ Michel Marcus, Jul 24 2020

a(n) = {if(n<2, return(1)); my(f = primes(primepi(n)), res = 1, t = #f); f = vector(#f, i, val(n, f[i])); for(i = 1, f[1], while(f[t] < i, t--; ); res+=(1<David A. Corneth, Jul 27 2020

a(n) = A327527(n!).

Terms a(31) and onwards from David A. Corneth, Jul 27 2020

A325276 Irregular triangle read by rows where row n is the omega-sequence of n!.

Original entry on oeis.org

1, 2, 2, 1, 4, 2, 2, 1, 5, 3, 2, 2, 1, 7, 3, 3, 1, 8, 4, 3, 2, 2, 1, 11, 4, 3, 2, 2, 1, 13, 4, 3, 2, 2, 1, 15, 4, 4, 1, 16, 5, 4, 2, 2, 1, 19, 5, 4, 2, 2, 1, 20, 6, 4, 2, 2, 1, 22, 6, 4, 2, 1, 24, 6, 5, 2, 2, 1, 28, 6, 5, 2, 2, 1, 29, 7, 5, 2, 2, 1

Offset: 0

Gus Wiseman, Apr 18 2019

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).

			Triangle begins:
  {}
  {}
   1
   2  2  1
   4  2  2  1
   5  3  2  2  1
   7  3  3  1
   8  4  3  2  2  1
  11  4  3  2  2  1
  13  4  3  2  2  1
  15  4  4  1
  16  5  4  2  2  1
  19  5  4  2  2  1
  20  6  4  2  2  1
  22  6  4  2  1
  24  6  5  2  2  1
  28  6  5  2  2  1
  29  7  5  2  2  1
  32  7  5  2  2  1
  33  8  5  2  2  1
  36  8  5  2  2  1
  38  8  5  2  2  1
  40  8  6  2  2  1
  41  9  6  2  2  1
  45  9  6  2  2  1
  47  9  6  2  2  1
  49  9  6  3  2  2  1
  52  9  6  3  2  2  1
  55  9  6  3  2  2  1
  56 10  6  3  2  2  1
  59 10  6  3  2  2  1

Row lengths are A325272. Row sums are A325274. Row n is row A325275(n) of A112798. Second-to-last column is A325273. Column k = 1 is A022559. Column k = 2 is A000720. Column k = 3 is A071626.
Cf. A000142, A006939, A081401, A303555, A323023, A325238, A325275, A325277.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).

omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
Table[omseq[n!],{n,0,30}]

A325278 Smallest number with adjusted frequency depth n.

Original entry on oeis.org

1, 2, 4, 6, 12, 60, 2520, 1286485200, 35933692027611398678865941374040400000

Offset: 0

Gus Wiseman, Apr 17 2019

nonn

The adjusted frequency depth of a positive integer n is 0 if n = 1, and otherwise it is 1 plus the number of times one must apply A181819 to reach a prime number, where A181819(k = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of k. For example, 180 has adjusted frequency depth 5 because we have: 180 -> 18 -> 6 -> 4 -> 3.

Differs from A182857 in having 2 instead of 3.

A subsequence of A325238.
Cf. A011784, A056239, A112798, A118914, A181819, A181821, A182857, A225485, A323023, A325258, A325272, A325277, A325278, A325280.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number).

nn=10000;
fd[n_]:=Switch[n,1,0,?PrimeQ,1,,1+fd[Times@@Prime/@Last/@FactorInteger[n]]];
fds=fd/@Range[nn];
Sort[Table[Position[fds,x][[1,1]],{x,Union[fds]}]]

A335407 Number of anti-run permutations of the prime indices of n!.

Original entry on oeis.org

1, 1, 1, 2, 0, 2, 3, 54, 0, 30, 105, 6090, 1512, 133056, 816480, 127209600, 0, 10090080, 562161600, 69864795000, 49989139200, 29593652088000, 382147120555200, 41810689605484800, 4359985823793600, 3025062801079038720, 49052072750637116160, 25835971971637227375360

Offset: 0

Gus Wiseman, Jul 01 2020

nonn

An anti-run is a sequence with no adjacent equal parts.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Conjecture: Only vanishes at n = 4 and n = 8.
a(16) = 0. Proof: 16! = 2^15 * m where bigomega(m) = A001222(m) = 13. We can't separate 15 1's with 13 other numbers. - David A. Corneth, Jul 04 2020

			The a(0) = 1 through a(6) = 3 anti-run permutations:
  ()  ()  (1)  (1,2)  .  (1,2,1,3,1)  (1,2,1,2,1,3,1)
               (2,1)     (1,3,1,2,1)  (1,2,1,3,1,2,1)
                                      (1,3,1,2,1,2,1)

Andrew Howroyd, Table of n, a(n) for n = 0..200

The version for Mersenne numbers is A335432.
Anti-run compositions are A003242.
Anti-run patterns are counted by A005649.
Permutations of prime indices are A008480.
Anti-runs are ranked by A333489.
Separable partitions are ranked by A335433.
Inseparable partitions are ranked by A335448.
Anti-run permutations of prime indices are A335452.
Strict permutations of prime indices are A335489.
Cf. A056239, A106351, A112798, A114938, A278990, A292884, A325535, A335125.
Factorial numbers: A000142, A001222, A002982, A007489, A022559, A027423, A054991, A108731, A181819, A181821, A325272, A325273, A325617.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n!]],!MatchQ[#,{_,x_,x_,_}]&]],{n,0,10}]

\\ See A335452 for count.
a(n)={count(factor(n!)[,2])} \\ Andrew Howroyd, Feb 03 2021

a(n) = A335452(A000142(n)). - Andrew Howroyd, Feb 03 2021

Terms a(14) and beyond from Andrew Howroyd, Feb 03 2021

A325509 Number of factorizations of n! into factorial numbers > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 2, 2, 3, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1

Offset: 0

Gus Wiseman, May 08 2019

nonn

			n = 10:
  (6*120*5040)
  (720*5040)
  (3628800)
n = 16:
  (2*2*2*2*1307674368000)
  (2*120*87178291200)
  (20922789888000)
n = 24:
  (2*2*6*25852016738884976640000)
  (24*25852016738884976640000)
  (620448401733239439360000)

Cf. A000142, A011371, A022559, A034876, A034878, A076934, A109100, A109101, A109102, A115627, A135291, A322583, A325272, A325273, A325276, A325508, A325510, A325511.

facs[n_,u_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d,u],Min@@#>=d&]],{d,Intersection[u,Rest[Divisors[n]]]}]];
Table[Length[facs[n!,Rest[Array[#!&,n]]]],{n,15}]

a(n) = 1 + A034876(n).

More terms from Alois P. Heinz, May 08 2019

A325274 Sum of the omega-sequence of n!.

Original entry on oeis.org

0, 0, 1, 5, 9, 13, 14, 20, 23, 25, 24, 30, 33, 35, 35, 40, 44, 46, 49, 51, 54, 56, 59, 61, 65, 67, 72, 75, 78, 80, 83, 85, 90, 90, 95, 97, 101, 103, 105, 106, 110, 112, 115, 117, 122, 125, 127, 129, 134, 136, 139, 140, 143, 145, 149, 153, 157, 159, 160, 162

Offset: 0

Gus Wiseman, Apr 18 2019

nonn

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1), with sum 13.

a(n) = A056239(A325275(n)).
Cf. A000142, A006939, A303555, A323023, A325238, A325272, A325273, A325276, A325277.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).

omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
Table[Total[omseq[n!]],{n,0,100}]

A325275 Heinz number of the omega-sequence of n!.

Original entry on oeis.org

1, 1, 2, 18, 126, 990, 850, 11970, 19530, 25830, 4606, 73458, 92862, 116298, 43134, 229086, 275418, 366894, 440946, 515394, 568062, 613206, 769158, 963378, 1060254, 1135602, 6108570, 6431490, 6915870, 8923590, 9398610, 10191870, 11352510, 3139866, 16458210

Offset: 0

Gus Wiseman, Apr 18 2019

nonn

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

A001222(a(n)) = A325272.
A055396(a(n)/2) = A325273.
A056239(a(n)) = A325274.
Row n of A325276 is row a(n) of A112798.
Cf. A000142, A006939, A056239, A112798, A303555, A323023, A325238, A325277.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).

omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
Table[Times@@Prime/@omseq[n!],{n,30}]

A336425 Number of ways to choose a divisor with distinct prime exponents of a divisor with distinct prime exponents of n!.

Original entry on oeis.org

1, 1, 3, 5, 24, 38, 132, 195, 570, 1588, 4193, 6086, 14561, 19232, 37142, 106479, 207291, 266871, 549726, 674330, 1465399, 3086598, 5939574, 7182133, 12324512, 28968994, 46819193, 82873443, 165205159, 196666406, 350397910, 406894074, 593725529, 1229814478, 1853300600, 4024414209, 6049714096, 6968090487, 9700557121, 16810076542, 26339337285

Offset: 0

Gus Wiseman, Aug 06 2020

nonn

			The a(4) = 24 divisors of divisors:
  1/1  2/1  3/1  4/1  8/1  12/1   24/1
       2/2  3/3  4/2  8/2  12/2   24/2
                 4/4  8/4  12/3   24/3
                      8/8  12/4   24/4
                           12/12  24/8
                                  24/12
                                  24/24

Max Alekseyev, Table of n, a(n) for n = 0..85

A336422 is the non-factorial generalization.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A327526 gives the maximum divisor of n with equal prime exponents.
A327498 gives the maximum divisor of n with distinct prime exponents.
A336414 counts divisors of n! with distinct prime exponents.
A336415 counts divisors of n! with equal prime exponents.
A336423 counts chains in A130091, with maximal version A336569.
Cf. A000005, A000110, A098859, A118914, A124010, A336424, A336500, A336568, A336570, A336571, A336865, A336866, A336869.
Factorial numbers: A000142, A022559, A027423, A048656, A048742, A071626, A325272, A325273, A325617, A327499, A336416, A336418, A336617.

strsigQ[n_]:=UnsameQ@@Last/@FactorInteger[n];
Table[Total[Cases[Divisors[n!],d_?strsigQ:>Count[Divisors[d],e_?strsigQ]]],{n,0,20}]

Terms a(21) onward from Max Alekseyev, Nov 07 2024

A336498 Irregular triangle read by rows where T(n,k) is the number of divisors of n! with k prime factors, counted with multiplicity.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 2, 2, 1, 1, 3, 4, 4, 3, 1, 1, 3, 5, 6, 6, 5, 3, 1, 1, 4, 8, 11, 12, 11, 8, 4, 1, 1, 4, 8, 11, 12, 12, 12, 12, 11, 8, 4, 1, 1, 4, 8, 12, 16, 19, 20, 20, 19, 16, 12, 8, 4, 1, 1, 4, 9, 15, 21, 26, 29, 30, 30, 29, 26, 21, 15, 9, 4, 1

Offset: 0

Gus Wiseman, Aug 03 2020

Row n is row n! of A146291. Row lengths are A022559(n) + 1.

			Triangle begins:
  1
  1
  1  1
  1  2  1
  1  2  2  2  1
  1  3  4  4  3  1
  1  3  5  6  6  5  3  1
  1  4  8 11 12 11  8  4  1
  1  4  8 11 12 12 12 12 11  8  4  1
  1  4  8 12 16 19 20 20 19 16 12  8  4  1
Row n = 6 counts the following divisors:
  1  2   4   8  16   48  144  720
     3   6  12  24   72  240
     5   9  18  36   80  360
        10  20  40  120
        15  30  60  180
            45  90
Row n = 7 counts the following divisors:
  1  2   4    8   16   48   144   720  5040
     3   6   12   24   72   240  1008
     5   9   18   36   80   336  1680
     7  10   20   40  112   360  2520
        14   28   56  120   504
        15   30   60  168   560
        21   42   84  180   840
        35   45   90  252  1260
             63  126  280
             70  140  420
            105  210  630
                 315

A000720 is column k = 1.
A008302 is the version for superprimorials.
A022559 gives row lengths minus one.
A027423 gives row sums.
A146291 is the generalization to non-factorials.
A336499 is the restriction to divisors in A130091.
A000142 lists factorial numbers.
A336415 counts uniform divisors of n!.
Cf. A000005, A001222, A118914, A124010, A181796, A327526, A336420.
Factorial numbers: A002982, A007489, A048656, A054991, A071626, A325272, A325617, A336414, A336415, A336416, A336418.

Table[Length[Select[Divisors[n!],PrimeOmega[#]==k&]],{n,0,10},{k,0,PrimeOmega[n!]}]

A336616 Maximum divisor of n! with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 2, 3, 24, 40, 720, 1008, 8064, 72576, 3628800, 5702400, 68428800, 80870400, 317011968, 118879488000, 1902071808000, 2487324672000, 44771844096000, 50039119872000, 1000782397440000, 21016430346240000, 5085976143790080000, 6156707963535360000

Offset: 0

Gus Wiseman, Jul 29 2020

nonn

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

			The sequence of terms together with their prime signatures begins:
             1: ()
             1: ()
             2: (1)
             3: (1)
            24: (3,1)
            40: (3,1)
           720: (4,2,1)
          1008: (4,2,1)
          8064: (7,2,1)
         72576: (7,4,1)
       3628800: (8,4,2,1)
       5702400: (8,4,2,1)
      68428800: (10,5,2,1)
      80870400: (10,5,2,1)
     317011968: (11,5,2,1)
  118879488000: (11,6,3,2,1)

David A. Corneth, Table of n, a(n) for n = 0..589
Gus Wiseman, Sequences counting and encoding certain classes of multisets

A327498 is the version not restricted to factorials, with quotient A327499.
A336414 counts these divisors.
A336617 is the quotient n!/a(n).
A336618 is the version for equal prime multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A327526 gives the maximum divisor of n with equal prime multiplicities.
A336415 counts divisors of n! with equal prime multiplicities.
Cf. A000005, A000110, A098859, A118914, A124010, A336424, A336500, A336568.
Factorial numbers: A000142, A007489, A022559, A027423, A048656, A048742, A071626, A325272, A325273, A325617, A336416.

Table[Max@@Select[Divisors[n!],UnsameQ@@Last/@If[#==1,{},FactorInteger[#]]&],{n,0,15}]

a(n) = { if(n < 2, return(1)); my(pr = primes(primepi(n)), res = pr[#pr]); for(i = 1, #pr, pr[i] = [pr[i], val(n, pr[i])] ); forstep(i = #pr, 2, -1, if(pr[i][2] < pr[i-1][2], res*=pr[i-1][1]^pr[i-1][2] ) ); res }
val(n, p) = my(r=0); while(n, r+=n\=p); r \\ David A. Corneth, Aug 25 2020

a(n) = A327498(n!).

cp's OEIS Frontend

A336415 Number of divisors of n! with equal prime multiplicities.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A325276 Irregular triangle read by rows where row n is the omega-sequence of n!.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A325278 Smallest number with adjusted frequency depth n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A335407 Number of anti-run permutations of the prime indices of n!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A325509 Number of factorizations of n! into factorial numbers > 1.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Formula

Extensions

A325274 Sum of the omega-sequence of n!.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A325275 Heinz number of the omega-sequence of n!.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A336425 Number of ways to choose a divisor with distinct prime exponents of a divisor with distinct prime exponents of n!.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica