cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 27 results. Next

A368252 Numbers with a record value of number of uniform divisors (A327527).

Original entry on oeis.org

1, 2, 4, 6, 12, 24, 30, 60, 120, 180, 210, 420, 840, 1260, 2310, 4620, 9240, 13860, 27720, 30030, 60060, 120120, 180180, 360360, 510510, 1021020, 2042040, 3063060, 6126120, 9699690, 19399380, 38798760, 58198140, 116396280, 223092870, 446185740, 892371480, 1338557220
Offset: 1

Views

Author

Amiram Eldar, Dec 19 2023

Keywords

Comments

First differs from A346016 at n = 48, and from A282472 at n = 155.

Examples

			The sequence A327527 starts with 1, 2, 2, 3, 2, 4, 2, 4, 3, 4, 2, 5, ... . The record values, 1, 2, 3, 4 and 5, occur at indices 1, 2, 4, 6, 12, the first 5 terms of this sequence.
		

Crossrefs

Subsequence of A025487.

Programs

  • Mathematica
    d[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, 1 + Total[2^Accumulate[Count[e, #] & /@ Range[Max[e], 1, -1]] - 1]]; d[1] = 1;
    v = Cases[Import["https://oeis.org/A025487/b025487.txt", "Table"], {, }][[;; , 2]];
    seq = {}; dm = 0; Do[If[(dk = d[v[[k]]]) > dm, dm = dk; AppendTo[seq, v[[k]]]], {k, 1, Length[v]}]; seq

A336416 Number of perfect-power divisors of n!.

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 7, 7, 11, 18, 36, 36, 47, 47, 84, 122, 166, 166, 221, 221, 346, 416, 717, 717, 1001, 1360, 2513, 2942, 4652, 4652, 5675, 5675, 6507, 6980, 13892, 17212, 20408, 20408, 39869, 45329, 51018, 51018, 68758, 68758, 105573, 138617, 284718, 284718, 338126, 421126
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

A number is a perfect power iff it is 1 or its prime exponents (signature) are not relatively prime.

Examples

			The a(1) = 0 through a(9) = 18 divisors:
       1: 1
       2: 1
       6: 1
      24: 1,4,8
     120: 1,4,8
     720: 1,4,8,9,16,36,144
    5040: 1,4,8,9,16,36,144
   40320: 1,4,8,9,16,32,36,64,128,144,576
  362880: 1,4,8,9,16,27,32,36,64,81,128,144,216,324,576,1296,1728,5184
		

Crossrefs

The maximum among these divisors is A090630, with quotient A251753.
The version for distinct prime exponents is A336414.
The uniform version is A336415.
Replacing factorials with Chernoff numbers (A006939) gives A336417.
Prime powers are A000961.
Perfect powers are A001597, with complement A007916.
Prime power divisors are counted by A022559.

Programs

  • Mathematica
    perpouQ[n_]:=Or[n==1,GCD@@FactorInteger[n][[All,2]]>1];
    Table[Length[Select[Divisors[n!],perpouQ]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, (d==1) || ispower(d)); \\ Michel Marcus, Aug 19 2020
    
  • PARI
    addhelp(val, "exponent of prime p in n!")
    val(n, p) = my(r=0); while(n, r+=n\=p);r
    a(n) = {if(n<=3, return(1)); my(pr = primes(primepi(n\2)), v = vector(#pr, i, val(n, pr[i])), res = 1, cv); for(i = 2, v[1], if(issquarefree(i), cv = v\i; res-=(prod(i = 1, #cv, cv[i]+1)-1)*(-1)^omega(i) ) ); res } \\ David A. Corneth, Aug 19 2020

Formula

a(p) = a(p-1) for prime p. - David A. Corneth, Aug 19 2020

Extensions

a(26)-a(34) from Jinyuan Wang, Aug 19 2020
a(35)-a(49) from David A. Corneth, Aug 19 2020

A062319 Number of divisors of n^n, or of A000312(n).

Original entry on oeis.org

1, 1, 3, 4, 9, 6, 49, 8, 25, 19, 121, 12, 325, 14, 225, 256, 65, 18, 703, 20, 861, 484, 529, 24, 1825, 51, 729, 82, 1653, 30, 29791, 32, 161, 1156, 1225, 1296, 5329, 38, 1521, 1600, 4961, 42, 79507, 44, 4005, 4186, 2209, 48, 9457, 99, 5151, 2704, 5565, 54
Offset: 0

Views

Author

Jason Earls, Jul 05 2001

Keywords

Comments

From Gus Wiseman, May 02 2021: (Start)
Conjecture: The number of divisors of n^n equals the number of pairwise coprime ordered n-tuples of divisors of n. Confirmed up to n = 30. For example, the a(1) = 1 through a(5) = 6 tuples are:
(1) (1,1) (1,1,1) (1,1,1,1) (1,1,1,1,1)
(1,2) (1,1,3) (1,1,1,2) (1,1,1,1,5)
(2,1) (1,3,1) (1,1,1,4) (1,1,1,5,1)
(3,1,1) (1,1,2,1) (1,1,5,1,1)
(1,1,4,1) (1,5,1,1,1)
(1,2,1,1) (5,1,1,1,1)
(1,4,1,1)
(2,1,1,1)
(4,1,1,1)
The unordered case (pairwise coprime n-multisets of divisors of n) is counted by A343654.
(End)

Examples

			From _Gus Wiseman_, May 02 2021: (Start)
The a(1) = 1 through a(5) = 6 divisors:
  1  1  1   1    1
     2  3   2    5
     4  9   4    25
        27  8    125
            16   625
            32   3125
            64
            128
            256
(End)
		

Crossrefs

Number of divisors of A000312(n).
Taking Omega instead of sigma gives A066959.
Positions of squares are A173339.
Diagonal n = k of the array A343656.
A000005 counts divisors.
A059481 counts k-multisets of elements of {1..n}.
A334997 counts length-k strict chains of divisors of n.
A343658 counts k-multisets of divisors.
Pairwise coprimality:
- A018892 counts coprime pairs of divisors.
- A084422 counts pairwise coprime subsets of {1..n}.
- A100565 counts pairwise coprime triples of divisors.
- A225520 counts pairwise coprime sets of divisors.
- A343652 counts maximal pairwise coprime sets of divisors.
- A343653 counts pairwise coprime non-singleton sets of divisors > 1.
- A343654 counts pairwise coprime sets of divisors > 1.

Programs

  • Magma
    [NumberOfDivisors(n^n): n in  [0..60]]; // Vincenzo Librandi, Nov 09 2014
    
  • Mathematica
    A062319[n_IntegerQ]:=DivisorSigma[0,n^n]; (* Enrique Pérez Herrero, Nov 09 2010 *)
    Join[{1},DivisorSigma[0,#^#]&/@Range[60]] (* Harvey P. Dale, Jun 06 2024 *)
  • PARI
    je=[]; for(n=0,200,je=concat(je,numdiv(n^n))); je
    
  • PARI
    { for (n=0, 1000, write("b062319.txt", n, " ", numdiv(n^n)); ) } \\ Harry J. Smith, Aug 04 2009
    
  • PARI
    a(n)=local(fm);fm=factor(n);prod(k=1,matsize(fm)[1],fm[k,2]*n+1) \\ Franklin T. Adams-Watters, May 03 2011
    
  • PARI
    a(n) = if(n==0, 1, sumdiv(n, d, n^omega(d))); \\ Seiichi Manyama, May 12 2021
    
  • Python
    from math import prod
    from sympy import factorint
    def A062319(n): return prod(n*d+1 for d in factorint(n).values()) # Chai Wah Wu, Jun 03 2021

Formula

a(n) = A000005(A000312(n)). - Enrique Pérez Herrero, Nov 09 2010
a(2^n) = A002064(n). - Gus Wiseman, May 02 2021
a(prime(n)) = prime(n) + 1. - Gus Wiseman, May 02 2021
a(n) = Product_{i=1..s} (1 + n * m_i) where (m_1,...,m_s) is the sequence of prime multiplicities (prime signature) of n. - Gus Wiseman, May 02 2021
a(n) = Sum_{d|n} n^omega(d) for n > 0. - Seiichi Manyama May 12 2021

A009998 Triangle in which j-th entry in i-th row is (j+1)^(i-j).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 4, 1, 1, 16, 27, 16, 5, 1, 1, 32, 81, 64, 25, 6, 1, 1, 64, 243, 256, 125, 36, 7, 1, 1, 128, 729, 1024, 625, 216, 49, 8, 1, 1, 256, 2187, 4096, 3125, 1296, 343, 64, 9, 1, 1, 512, 6561, 16384, 15625, 7776, 2401, 512, 81, 10, 1
Offset: 0

Views

Author

Keywords

Comments

Read as a square array this is the Hilbert transform of triangle A123125 (see A145905 for the definition of this term). For example, the fourth row of A123125 is (0,1,4,1) and the expansion (x + 4*x^2 + x^3)/(1-x)^4 = x + 8*x^2 + 27*x^3 + 64*x^4 + ... generates the entries in the fourth row of this array read as a square. - Peter Bala, Oct 28 2008

Examples

			Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  3,  1;
  1,  8,  9,  4,  1;
  1, 16, 27, 16,  5,  1;
  1, 32, 81, 64, 25,  6,  1;
  ...
From _Gus Wiseman_, May 01 2021: (Start)
The rows of the triangle are obtained by reading antidiagonals upward in the following table of A(k,n) = n^k, with offset k = 0, n = 1:
         n=1:     n=2:     n=3:     n=4:     n=5:     n=6:
   k=0:   1        1        1        1        1        1
   k=1:   1        2        3        4        5        6
   k=2:   1        4        9       16       25       36
   k=3:   1        8       27       64      125      216
   k=4:   1       16       81      256      625     1296
   k=5:   1       32      243     1024     3125     7776
   k=6:   1       64      729     4096    15625    46656
   k=7:   1      128     2187    16384    78125   279936
   k=8:   1      256     6561    65536   390625  1679616
   k=9:   1      512    19683   262144  1953125 10077696
  k=10:   1     1024    59049  1048576  9765625 60466176
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 24.

Crossrefs

Row sums give A026898.
Column n = 2 of the array is A000079.
Column n = 3 of the array is A000244.
Row k = 2 of the array is A000290.
Row k = 3 of the array is A000578.
Diagonal n = k of the array is A000312.
Diagonal n = k + 1 of the array is A000169.
Diagonal n = k + 2 of the array is A000272.
The transpose of the array is A009999.
The numbers of divisors of the entries are A343656 (row sums: A343657).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Haskell
    a009998 n k = (k + 1) ^ (n - k)
    a009998_row n = a009998_tabl !! n
    a009998_tabl = map reverse a009999_tabl
    -- Reinhard Zumkeller, Feb 02 2014
    
  • Maple
    E := (n,x) -> `if`(n=0,1,x*(1-x)*diff(E(n-1,x),x)+E(n-1,x)*(1+(n-1)*x));
    G := (n,x) -> E(n,x)/(1-x)^(n+1);
    A009998 := (n,k) -> coeff(series(G(n-k,x),x,18),x,k);
    seq(print(seq(A009998(n,k),k=0..n)),n=0..6);
    # Peter Luschny, Aug 02 2010
  • Mathematica
    Flatten[Table[(j+1)^(i-j),{i,0,20},{j,0,i}]] (* Harvey P. Dale, Dec 25 2012 *)
  • PARI
    T(i,j)=(j+1)^(i-j) \\ Charles R Greathouse IV, Feb 06 2017

Formula

T(n,n) = 1; T(n,k) = (k+1)*T(n-1,k) for k=0..n-1. - Reinhard Zumkeller, Feb 02 2014
T(n,m) = (m+1)*Sum_{k=0..n-m}((n+1)^(k-1)*(n-m)^(n-m-k)*(-1)^(n-m-k)*binomial(n-m-1,k-1)). - Vladimir Kruchinin, Sep 12 2015

Extensions

a(62) corrected to 512 by T. D. Noe, Dec 20 2007

A334997 Array T read by ascending antidiagonals: T(n, k) = Sum_{d divides n} T(d, k-1) with T(n, 0) = 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 3, 3, 4, 1, 1, 2, 6, 4, 5, 1, 1, 4, 3, 10, 5, 6, 1, 1, 2, 9, 4, 15, 6, 7, 1, 1, 4, 3, 16, 5, 21, 7, 8, 1, 1, 3, 10, 4, 25, 6, 28, 8, 9, 1, 1, 4, 6, 20, 5, 36, 7, 36, 9, 10, 1, 1, 2, 9, 10, 35, 6, 49, 8, 45, 10, 11, 1, 1, 6, 3, 16, 15, 56, 7, 64, 9, 55, 11, 12, 1
Offset: 1

Views

Author

Stefano Spezia, May 19 2020

Keywords

Comments

T(n, k) is called the generalized divisor function (see Beekman).
As an array with offset n=1, k=0, T(n,k) is the number of length-k chains of divisors of n. For example, the T(4,3) = 10 chains are: 111, 211, 221, 222, 411, 421, 422, 441, 442, 444. - Gus Wiseman, Aug 04 2022

Examples

			From _Gus Wiseman_, Aug 04 2022: (Start)
Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=1:  1   1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8   9
  n=3:  1   2   3   4   5   6   7   8   9
  n=4:  1   3   6  10  15  21  28  36  45
  n=5:  1   2   3   4   5   6   7   8   9
  n=6:  1   4   9  16  25  36  49  64  81
  n=7:  1   2   3   4   5   6   7   8   9
  n=8:  1   4  10  20  35  56  84 120 165
The T(4,5) = 21 chains:
  (1,1,1,1,1)  (4,2,1,1,1)  (4,4,2,2,2)
  (2,1,1,1,1)  (4,2,2,1,1)  (4,4,4,1,1)
  (2,2,1,1,1)  (4,2,2,2,1)  (4,4,4,2,1)
  (2,2,2,1,1)  (4,2,2,2,2)  (4,4,4,2,2)
  (2,2,2,2,1)  (4,4,1,1,1)  (4,4,4,4,1)
  (2,2,2,2,2)  (4,4,2,1,1)  (4,4,4,4,2)
  (4,1,1,1,1)  (4,4,2,2,1)  (4,4,4,4,4)
The T(6,3) = 16 chains:
  (1,1,1)  (3,1,1)  (6,2,1)  (6,6,1)
  (2,1,1)  (3,3,1)  (6,2,2)  (6,6,2)
  (2,2,1)  (3,3,3)  (6,3,1)  (6,6,3)
  (2,2,2)  (6,1,1)  (6,3,3)  (6,6,6)
The triangular form T(n-k,k) gives the number of length k chains of divisors of n - k. It begins:
  1
  1  1
  1  2  1
  1  2  3  1
  1  3  3  4  1
  1  2  6  4  5  1
  1  4  3 10  5  6  1
  1  2  9  4 15  6  7  1
  1  4  3 16  5 21  7  8  1
  1  3 10  4 25  6 28  8  9  1
  1  4  6 20  5 36  7 36  9 10  1
  1  2  9 10 35  6 49  8 45 10 11  1
(End)
		

References

  • Richard Beekman, An Introduction to Number-Theoretic Combinatorics, Lulu Press 2017.

Crossrefs

Cf. A000217 (4th row), A000290 (6th row), A000292 (8th row), A000332 (16th row), A000389 (32nd row), A000537 (36th row), A000578 (30th row), A002411 (12th row), A002417 (24th row), A007318, A027800 (48th row), A335078, A335079.
Column k = 2 of the array is A007425.
Column k = 3 of the array is A007426.
Column k = 4 of the array is A061200.
The transpose of the array is A077592.
The subdiagonal n = k + 1 of the array is A163767.
The version counting all multisets of divisors (not just chains) is A343658.
The strict case is A343662 (row sums: A337256).
Diagonal n = k of the array is A343939.
Antidiagonal sums of the array (or row sums of the triangle) are A343940.
A067824(n) counts strict chains of divisors starting with n.
A074206(n) counts strict chains of divisors from n to 1.
A146291 counts divisors by Omega.
A251683(n,k) counts strict length k + 1 chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict length k chains of divisors from n to 1.
A337255(n,k) counts strict length k chains of divisors starting with n.

Programs

  • Mathematica
    T[n_,k_]:=If[n==1,1,Product[Binomial[Extract[Extract[FactorInteger[n],i],2]+k,k],{i,1,Length[FactorInteger[n]]}]]; Table[T[n-k,k],{n,1,13},{k,0,n-1}]//Flatten
  • PARI
    T(n, k) = if (k==0, 1, sumdiv(n, d, T(d, k-1)));
    matrix(10, 10, n, k, T(n, k-1)) \\ to see the array for n>=1, k >=0; \\ Michel Marcus, May 20 2020

Formula

T(n, k) = Sum_{d divides n} T(d, k-1) with T(n, 0) = 1 (see Theorem 3 in Beekman's article).
T(i*j, k) = T(i, k)*T(j, k) if i and j are coprime positive integers (see Lemma 1 in Beekman's article).
T(p^m, k) = binomial(m+k, k) for every prime p (see Lemma 2 in Beekman's article).

Extensions

Duplicate term removed by Stefano Spezia, Jun 03 2020

A336414 Number of divisors of n! with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 2, 3, 7, 10, 20, 27, 48, 86, 147, 195, 311, 390, 595, 1031, 1459, 1791, 2637, 3134, 4747, 7312, 10766, 12633, 16785, 26377, 36142, 48931, 71144, 82591, 112308, 128023, 155523, 231049, 304326, 459203, 568095, 642446, 812245, 1137063, 1441067, 1612998, 2193307, 2429362
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

A number has distinct prime multiplicities iff its prime signature is strict.

Examples

			The first and second columns below are the a(6) = 20 counted divisors of 6! together with their prime signatures. The third column shows the A000005(6!) - a(6) = 10 remaining divisors.
      1: ()      20: (2,1)    |    6: (1,1)
      2: (1)     24: (3,1)    |   10: (1,1)
      3: (1)     40: (3,1)    |   15: (1,1)
      4: (2)     45: (2,1)    |   30: (1,1,1)
      5: (1)     48: (4,1)    |   36: (2,2)
      8: (3)     72: (3,2)    |   60: (2,1,1)
      9: (2)     80: (4,1)    |   90: (1,2,1)
     12: (2,1)  144: (4,2)    |  120: (3,1,1)
     16: (4)    360: (3,2,1)  |  180: (2,2,1)
     18: (1,2)  720: (4,2,1)  |  240: (4,1,1)
		

Crossrefs

Perfect-powers are A001597, with complement A007916.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
The maximum divisor with distinct prime multiplicities is A327498.
Divisors of n! with equal prime multiplicities are counted by A336415.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n!],UnsameQ@@Last/@FactorInteger[#]&]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); #vecsort(ex,,8) == #ex); \\ Michel Marcus, Jul 24 2020

Formula

a(n) = A181796(n!).

Extensions

a(21)-a(41) from Alois P. Heinz, Jul 24 2020

A077592 Table by antidiagonals of tau_k(n), the k-th Piltz function (see A007425), or n-th term of the sequence resulting from applying the inverse Möbius transform (k-1) times to the all-ones sequence.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 6, 2, 1, 1, 6, 5, 10, 3, 4, 1, 1, 7, 6, 15, 4, 9, 2, 1, 1, 8, 7, 21, 5, 16, 3, 4, 1, 1, 9, 8, 28, 6, 25, 4, 10, 3, 1, 1, 10, 9, 36, 7, 36, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 49, 6, 35, 10, 9, 2, 1, 1, 12, 11, 55, 9, 64, 7, 56, 15, 16, 3, 6, 1
Offset: 1

Views

Author

Henry Bottomley, Nov 08 2002

Keywords

Comments

As an array with offset n=0, k=1, also the number of length n chains of divisors of k. - Gus Wiseman, Aug 04 2022

Examples

			T(6,3) = 9 because we have: 1*1*6, 1*2*3, 1*3*2, 1*6*1, 2*1*3, 2*3*1, 3*1*2, 3*2*1, 6*1*1. - _Geoffrey Critzer_, Feb 16 2015
From _Gus Wiseman_, May 03 2021: (Start)
Array begins:
       k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=0:  1   1   1   1   1   1   1   1
  n=1:  1   2   2   3   2   4   2   4
  n=2:  1   3   3   6   3   9   3  10
  n=3:  1   4   4  10   4  16   4  20
  n=4:  1   5   5  15   5  25   5  35
  n=5:  1   6   6  21   6  36   6  56
  n=6:  1   7   7  28   7  49   7  84
  n=7:  1   8   8  36   8  64   8 120
  n=8:  1   9   9  45   9  81   9 165
The triangular form T(n,k) = A(n-k,k) gives the number of length n - k chains of divisors of k. It begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  6  2  1
  1  6  5 10  3  4  1
  1  7  6 15  4  9  2  1
  1  8  7 21  5 16  3  4  1
  1  9  8 28  6 25  4 10  3  1
  1 10  9 36  7 36  5 20  6  4  1
  1 11 10 45  8 49  6 35 10  9  2  1
(End)
		

Crossrefs

Columns include (with multiplicity and some offsets) A000012, A000027, A000027, A000217, A000027, A000290, A000027, A000292, A000217, A000290, A000027, A002411, A000027, A000290, A000290, A000332 etc.
Cf. A077593.
Row n = 2 of the array is A007425.
Row n = 3 of the array is A007426.
Row n = 4 of the array is A061200.
The diagonal n = k of the array (central column of the triangle) is A163767.
The transpose of the array is A334997.
Diagonal n = k of the array is A343939.
Antidiagonal sums of the array (or row sums of the triangle) are A343940.
A067824(n) counts strict chains of divisors starting with n.
A074206(n) counts strict chains of divisors from n to 1.
A146291(n,k) counts divisors of n with k prime factors (with multiplicity).
A251683(n,k) counts strict length k + 1 chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict length k chains of divisors from n to 1.
A337255(n,k) counts strict length k chains of divisors starting with n.

Programs

  • Maple
    with(numtheory):
    A:= proc(n,k) option remember; `if`(k=1, 1,
          add(A(d, k-1), d=divisors(n)))
        end:
    seq(seq(A(n, 1+d-n), n=1..d), d=1..14);  # Alois P. Heinz, Feb 25 2015
  • Mathematica
    tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[tau[n - k + 1, k], {n, 14}, {k, n, 1, -1}] // Flatten (* Robert G. Wilson v *)
    tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#] + k - 1, k - 1] & /@ FactorInteger[n]); Table[tau[k, n - k + 1], {n, 1, 13}, {k, 1, n}] // Flatten (* Amiram Eldar, Sep 13 2020 *)
    Table[Length[Select[Tuples[Divisors[k],n-k],And@@Divisible@@@Partition[#,2,1]&]],{n,12},{k,1,n}] (* TRIANGLE, Gus Wiseman, May 03 2021 *)
    Table[Length[Select[Tuples[Divisors[k],n-1],And@@Divisible@@@Partition[#,2,1]&]],{n,6},{k,6}] (* ARRAY, Gus Wiseman, May 03 2021 *)

Formula

If n = Product_i p_i^e_i, then T(n,k) = Product_i C(k+e_i-1, e_i). T(n,k) = Sum_d{d|n} T(n-1,d) = A077593(n,k) - A077593(n-1,k).
Columns are multiplicative.
Dirichlet g.f. for column k: Zeta(s)^k. - Geoffrey Critzer, Feb 16 2015
A(n,k) = A334997(k,n). - Gus Wiseman, Aug 04 2022

Extensions

Typo in formula fixed by Geoffrey Critzer, Feb 16 2015

A343656 Array read by antidiagonals where A(n,k) is the number of divisors of n^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 5, 2, 1, 1, 6, 5, 7, 3, 4, 1, 1, 7, 6, 9, 4, 9, 2, 1, 1, 8, 7, 11, 5, 16, 3, 4, 1, 1, 9, 8, 13, 6, 25, 4, 7, 3, 1, 1, 10, 9, 15, 7, 36, 5, 10, 5, 4, 1, 1, 11, 10, 17, 8, 49, 6, 13, 7, 9, 2, 1, 1, 12, 11, 19, 9, 64, 7, 16, 9, 16, 3, 6, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 28 2021

Keywords

Comments

First differs from A343658 at A(4,2) = 5, A343658(4,2) = 6.
As a triangle, T(n,k) = number of divisors of k^(n-k).

Examples

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7
  n=1:  1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8
  n=3:  1   2   3   4   5   6   7   8
  n=4:  1   3   5   7   9  11  13  15
  n=5:  1   2   3   4   5   6   7   8
  n=6:  1   4   9  16  25  36  49  64
  n=7:  1   2   3   4   5   6   7   8
  n=8:  1   4   7  10  13  16  19  22
  n=9:  1   3   5   7   9  11  13  15
Triangle begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  5  2  1
  1  6  5  7  3  4  1
  1  7  6  9  4  9  2  1
  1  8  7 11  5 16  3  4  1
  1  9  8 13  6 25  4  7  3  1
  1 10  9 15  7 36  5 10  5  4  1
  1 11 10 17  8 49  6 13  7  9  2  1
  1 12 11 19  9 64  7 16  9 16  3  6  1
  1 13 12 21 10 81  8 19 11 25  4 15  2  1
For example, row n = 8 counts the following divisors:
  1  64  243  256  125  36  7  1
     32  81   128  25   18  1
     16  27   64   5    12
     8   9    32   1    9
     4   3    16        6
     2   1    8         4
     1        4         3
              2         2
              1         1
		

Crossrefs

Columns k=1..9 of the array give A000005, A048691, A048785, A344327, A344328, A344329, A343526, A344335, A344336.
Row n = 6 of the array is A000290.
Diagonal n = k of the array is A062319.
Array antidiagonal sums (row sums of the triangle) are A343657.
Dominated by A343658.
A000312 = n^n.
A007318 counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Mathematica
    Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]
  • PARI
    A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

Formula

A(n,k) = A000005(A009998(n,k)), where A009998(n,k) = n^k is the interpretation as an array.
A(n,k) = Sum_{d|n} k^omega(d). - Seiichi Manyama, May 15 2021

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

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

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.

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n!],SameQ@@Last/@FactorInteger[#]&]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); (#ex==0) || (vecmin(ex) == vecmax(ex))); \\ Michel Marcus, Jul 24 2020
    
  • PARI
    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

Formula

a(n) = A327527(n!).

Extensions

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

A163767 a(n) = tau_{n}(n) = number of ordered n-factorizations of n.

Original entry on oeis.org

1, 2, 3, 10, 5, 36, 7, 120, 45, 100, 11, 936, 13, 196, 225, 3876, 17, 3078, 19, 4200, 441, 484, 23, 62400, 325, 676, 3654, 11368, 29, 27000, 31, 376992, 1089, 1156, 1225, 443556, 37, 1444, 1521, 459200, 41, 74088, 43, 43560, 46575, 2116, 47, 11995200, 1225
Offset: 1

Views

Author

Paul D. Hanna, Aug 04 2009

Keywords

Comments

Also the number of length n - 1 chains of divisors of n. - Gus Wiseman, May 07 2021

Examples

			Successive Dirichlet self-convolutions of the all 1's sequence begin:
(1),1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,... (A000012)
1,(2),2,3,2,4,2,4,3,4,2,6,2,4,4,5,... (A000005)
1,3,(3),6,3,9,3,10,6,9,3,18,3,9,9,15,... (A007425)
1,4,4,(10),4,16,4,20,10,16,4,40,4,16,16,35,... (A007426)
1,5,5,15,(5),25,5,35,15,25,5,75,5,25,25,70,... (A061200)
1,6,6,21,6,(36),6,56,21,36,6,126,6,36,36,126,... (A034695)
1,7,7,28,7,49,(7),84,28,49,7,196,7,49,49,210,... (A111217)
1,8,8,36,8,64,8,(120),36,64,8,288,8,64,64,330,... (A111218)
1,9,9,45,9,81,9,165,(45),81,9,405,9,81,81,495,... (A111219)
1,10,10,55,10,100,10,220,55,(100),10,550,10,100,... (A111220)
1,11,11,66,11,121,11,286,66,121,(11),726,11,121,... (A111221)
1,12,12,78,12,144,12,364,78,144,12,(936),12,144,... (A111306)
...
where the main diagonal forms this sequence.
From _Gus Wiseman_, May 07 2021: (Start)
The a(1) = 1 through a(5) = 5 chains of divisors:
  ()  (1)  (1/1)  (1/1/1)  (1/1/1/1)
      (2)  (3/1)  (2/1/1)  (5/1/1/1)
           (3/3)  (2/2/1)  (5/5/1/1)
                  (2/2/2)  (5/5/5/1)
                  (4/1/1)  (5/5/5/5)
                  (4/2/1)
                  (4/2/2)
                  (4/4/1)
                  (4/4/2)
                  (4/4/4)
(End)
		

Crossrefs

Main diagonal of A077592.
Diagonal n = k + 1 of the array A334997.
The version counting all multisets of divisors (not just chains) is A343935.
A000005 counts divisors.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A146291 counts divisors of n with k prime factors (with multiplicity).
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts nonempty strict chains of divisors of n.
A251683/A334996 count strict nonempty length-k divisor chains from n to 1.
A337255 counts strict length-k chains of divisors starting with n.
A339564 counts factorizations with a selected factor.
A343662 counts strict length-k chains of divisors (row sums: A337256).
Cf. A060690.

Programs

  • Mathematica
    Table[Times@@(Binomial[#+n-1,n-1]&/@FactorInteger[n][[All,2]]),{n,1,50}] (* Enrique Pérez Herrero, Dec 25 2013 *)
  • PARI
    {a(n,m=n)=if(n==1,1,if(m==1,1,sumdiv(n,d,a(d,1)*a(n/d,m-1))))}
    
  • Python
    from math import prod, comb
    from sympy import factorint
    def A163767(n): return prod(comb(n+e-1,e) for e in factorint(n).values()) # Chai Wah Wu, Jul 05 2024

Formula

a(p) = p for prime p.
a(n) = n^k when n is the product of k distinct primes (conjecture).
a(n) = n-th term of the n-th Dirichlet self-convolution of the all 1's sequence.
a(2^n) = A060690(n). - Alois P. Heinz, Jun 12 2024
Showing 1-10 of 27 results. Next