A343656 -id:A343656 - cp's OEIS frontend

A002064 Cullen numbers: a(n) = n*2^n + 1.

1, 3, 9, 25, 65, 161, 385, 897, 2049, 4609, 10241, 22529, 49153, 106497, 229377, 491521, 1048577, 2228225, 4718593, 9961473, 20971521, 44040193, 92274689, 192937985, 402653185, 838860801, 1744830465, 3623878657, 7516192769, 15569256449, 32212254721, 66571993089

Offset: 0

N. J. A. Sloane

Binomial transform is A084859. Inverse binomial transform is A004277. - Paul Barry, Jun 12 2003
Let A be the Hessenberg matrix of order n defined by: A[1,j]=1, A[i,i]:=2,(i>1), A[i,i-1] =-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)= (-1)^(n-1)*coeff(charpoly(A,x),x). - Milan Janjic, Jan 26 2010
Indices of primes are listed in A005849. - M. F. Hasler, Jan 18 2015
Add the list of fractions beginning with 1/2 + 3/4 + 7/8 + ... + (2^n - 1)/2^n and take the sums pairwise from left to right. For 1/2 + 3/4 = 5/4, 5 + 4 = 9 = a(2); for 5/4 + 7/8 = 17/8, 17 + 8 = 25 = a(3); for 17/8 + 15/16 = 49/16, 49 + 16 = 65 = a(4); for 49/16 + 31/32 = 129/32, 129 + 32 = 161 = a(5). For each pairwise sum a/b, a + b = n*2^(n+1). - J. M. Bergot, May 06 2015
Number of divisors of (2^n)^(2^n). - Gus Wiseman, May 03 2021
Named after the Irish Jesuit priest James Cullen (1867-1933), who checked the primality of the terms up to n=100. - Amiram Eldar, Jun 05 2021

			G.f. = 1 + 3*x + 9*x^2 + 25*x^3 + 65*x^4 + 161*x^5 + 385*x^6 + 897*x^7 + ... - _Michael Somos_, Jul 18 2018

G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.
R. K. Guy, Unsolved Problems in Number Theory, B20.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 240-242.
W. Sierpiński, Elementary Theory of Numbers. Państ. Wydaw. Nauk., Warsaw, 1964, p. 346.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n=0..300
Ray Ballinger, Cullen Primes: Definition and Status.
Attila Bérczes, István Pink, and Paul Thomas Young, Cullen numbers and Woodall numbers in generalized Fibonacci sequences, J. Num. Theor. (2024) Vol. 262, 86-102.
Yuri Bilu, Diego Marques, and Alain Togbé, Generalized Cullen numbers in linear recurrence sequences, Journal of Number Theory, Vol. 202 (2019), pp. 412-425; arXiv preprint, arXiv:1806.09441 [math.NT], 2018.
Daniel Birmajer, Juan B. Gil, David S. Kenepp, and Michael D. Weiner, Restricted generating trees for weak orderings, arXiv:2108.04302 [math.CO], 2021.
C. K. Caldwell, The Top Twenty: Cullen Primes.
James Cullen, Question 15897, Educational Times, Vol. 58 (December 1905), p. 534.
Orhan Eren and Yüksel Soykan, Gaussian Generalized Woodall Numbers, Arch. Current Res. Int'l (2023) Vol. 23, Iss. 8, Art. No. ACRI.108618, 48-68. See p. 50.
Orhan Eren and Yüksel Soykan, On Dual Hyperbolic Generalized Woodall Numbers, Archives Current Res. Int'l (2024) Vol. 24, Iss. 11, Art. No. ACRI.126420, 398-423. See p. 401.
Jon Grantham and Hester Graves, The abc Conjecture Implies That Only Finitely Many Cullen Numbers Are Repunits, arXiv:2009.04052 [math.NT], 2020.
José María Grau and Florian Luca, Cullen numbers with the Lehmer property, Proceedings of the American Mathematical Society, Vol. 140, No. 1 (2012), pp. 129-134; arXiv preprint, arXiv:1103.3578 [math.NT], Mar 18 2011.
Paul Leyland, Factors of Cullen and Woodall numbers.
Paul Leyland, Generalized Cullen and Woodall numbers.
Diego Marques, On Generalized Cullen and Woodall Numbers That are Also Fibonacci Numbers, Journal of Integer Sequences, Vol. 17 (2014), Article 14.9.4.
Hisanori Mishima, Factorizations of many number sequences, Cullen numbers (n = 1 to 100), (n = 101 to 200), (n = 201 to 300), (n = 301 to 323).
Simon Plouffe, Approximations de Séries Génératrices et Quelques Conjectures, Dissertation, Université du Québec à Montréal, 1992.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Wacław Sierpiński, Elementary Theory of Numbers, Warszawa 1964.
Amelia Carolina Sparavigna, On the generalized sums of Mersenne, Fermat, Cullen and Woodall Numbers, Politecnico di Torino (Italy, 2019).
Amelia Carolina Sparavigna, Composition Operations of Generalized Entropies Applied to the Study of Numbers, International Journal of Sciences, Vol. 8, No. 4 (2019), pp. 87-92.
Amelia Carolina Sparavigna, The groupoids of Mersenne, Fermat, Cullen, Woodall and other Numbers and their representations by means of integer sequences, Politecnico di Torino, Italy (2019), [math.NT].
Amelia Carolina Sparavigna, Some Groupoids and their Representations by Means of Integer Sequences, International Journal of Sciences, Vol. 8, No. 10 (2019).
Eric Weisstein's World of Mathematics, Cullen Number.
Wikipedia, Cullen number.
Index entries for linear recurrences with constant coefficients, signature (5,-8,4).

Cf. A005849, A003261, A050914, A130197, A134081, A001787, A143038, A156708, A181527.
Cf. A000325, A186947.
Diagonal k = n + 1 of A046688.
A000005 counts divisors of n.
A000312 = n^n.
A002109 gives hyperfactorials (sigma: A260146, omega: A303281).
A057156 = (2^n)^(2^n).
A062319 counts divisors of n^n.
A173339 lists positions of squares in A062319.
A188385 gives the highest prime exponent in n^n.
A249784 counts divisors of n^n^n.
Cf. A000169, A000272, A036289, A066959, A176029, A340841, A343656.

a002064 n = n * 2 ^ n + 1
a002064_list = 1 : 3 : (map (+ 1) $ zipWith (-) (tail xs) xs)
   where xs = map (* 4) a002064_list
-- Reinhard Zumkeller, Mar 16 2013

[n*2^n + 1: n in [0..40]]; // Vincenzo Librandi, May 07 2015

A002064:=-(1-2*z+2*z**2)/((z-1)*(-1+2*z)**2); # conjectured by Simon Plouffe in his 1992 dissertation

Table[n*2^n+1,{n,0,2*4!}] (* Vladimir Joseph Stephan Orlovsky, Apr 25 2010 *)
LinearRecurrence[{5,-8,4},{1,3,9},51] (* Harvey P. Dale, Oct 13 2011 *)
CoefficientList[Series[(1 - 2 x + 2 x^2)/((1 - x) (2 x - 1)^2), {x, 0, 50}], x] (* Vincenzo Librandi, May 07 2015 *)

A002064(n)=n*2^n+1  \\ M. F. Hasler, Oct 31 2012

a(n) = 4a(n-1) - 4a(n-2) + 1. - Paul Barry, Jun 12 2003
a(n) = sum of row (n+1) of triangle A130197. Example: a(3) = 25 = (12 + 8 + 4 + 1), row 4 of A130197. - Gary W. Adamson, May 16 2007
Row sums of triangle A134081. - Gary W. Adamson, Oct 07 2007
Equals row sums of triangle A143038. - Gary W. Adamson, Jul 18 2008
Equals row sums of triangle A156708. - Gary W. Adamson, Feb 13 2009
G.f.: -(1-2*x+2*x^2)/((-1+x)*(2*x-1)^2). a(n) = A001787(n+1)+1-A000079(n). - R. J. Mathar, Nov 16 2007
a(n) = 1 + 2^(n + log_2(n)) ~ 1 + A000079(n+A004257(n)). a(n) ~ A000051(n+A004257(n)). - Jonathan Vos Post, Jul 20 2008
a(0)=1, a(1)=3, a(2)=9, a(n) = 5*a(n-1)-8*a(n-2)+4*a(n-3). - Harvey P. Dale, Oct 13 2011
a(n) = A036289(n) + 1 = A003261(n) + 2. - Reinhard Zumkeller, Mar 16 2013
E.g.f.: 2*x*exp(2*x) + exp(x). - Robert Israel, Dec 12 2014
a(n) = 2^n * A000325(n) = 4^n * A186947(-n) for all n in Z. - Michael Somos, Jul 18 2018
a(n) = Sum_{i=0..n-1} a(i) + A000325(n+1). - Ivan N. Ianakiev, Aug 07 2019
a(n) = sigma((2^n)^(2^n)) = A000005(A057156(n)) = A062319(2^n). - Gus Wiseman, May 03 2021
Sum_{n>=0} 1/a(n) = A340841. - Amiram Eldar, Jun 05 2021

Edited by M. F. Hasler, Oct 31 2012

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

Jason Earls, Jul 05 2001

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)

			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)

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Harry J. Smith)

Cf. A046798, A077592, A035116. - Enrique Pérez Herrero, Nov 09 2010
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.
Cf. A000169, A000272, A002064, A002109, A009998, A048691, A143773, A146291, A176029, A327527, A343657.

[NumberOfDivisors(n^n): n in  [0..60]]; // Vincenzo Librandi, Nov 09 2014

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 *)

je=[]; for(n=0,200,je=concat(je,numdiv(n^n))); je

{ for (n=0, 1000, write("b062319.txt", n, " ", numdiv(n^n)); ) } \\ Harry J. Smith, Aug 04 2009

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

a(n) = if(n==0, 1, sumdiv(n, d, n^omega(d))); \\ Seiichi Manyama, May 12 2021

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

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

N. J. A. Sloane

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

			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)

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.

T. D. Noe, Rows n=0..50 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Peter Luschny, Figurate number - a very short introduction

Row sums give A026898.
Cf. A088956, A123125, A179927.
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}.
Cf. A002109, A062319, A066959, A143773, A176029, A204688, A326358, A327527.

a009998 n k = (k + 1) ^ (n - k)
a009998_row n = a009998_tabl !! n
a009998_tabl = map reverse a009999_tabl
-- Reinhard Zumkeller, Feb 02 2014

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

Flatten[Table[(j+1)^(i-j),{i,0,20},{j,0,i}]] (* Harvey P. Dale, Dec 25 2012 *)

T(i,j)=(j+1)^(i-j) \\ Charles R Greathouse IV, Feb 06 2017

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

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

Stefano Spezia, May 19 2020

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

			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)

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

Stefano Spezia, First 150 antidiagonals of the array, flattened
Richard Beekman, A General Solution of the Ferris Wheel Problem.

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.
Cf. A000005, A018892, A051026, A062319, A066959, A176029, A327527, A343652, A343656.

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

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

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).

Duplicate term removed by Stefano Spezia, Jun 03 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

Henry Bottomley, Nov 08 2002

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

			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)

Alois P. Heinz, Antidiagonals n = 1..141, flattened
Adolf Piltz, Ueber das Gesetz, nach welchem die mittlere Darstellbarkeit der natürlichen Zahlen als Produkte einer gegebenen Anzahl Faktoren mit der Grösse der Zahlen wächst, Doctoral Dissertation, Friedrich-Wilhelms-Universität zu Berlin, 1881; the k-th Piltz function tau_k(n) is denoted by phi(n,k) and its recurrence and Dirichlet series appear on p. 6.
Wikipedia, Adolf Piltz.

Rows include: A000012, A000005, A034695, A111217, A111218, A111219, A111220, A111221, A111306.
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.
Cf. A018892, A051026, A062319, A143773, A176029, A327527, A337256, A343656, A343658, A343662.

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

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 *)

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

Typo in formula fixed by Geoffrey Critzer, Feb 16 2015

A343658 Array read by antidiagonals where A(n,k) is the number of ways to choose a multiset of k divisors of n.

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, 10, 2, 1, 1, 8, 7, 21, 5, 20, 3, 4, 1, 1, 9, 8, 28, 6, 35, 4, 10, 3, 1, 1, 10, 9, 36, 7, 56, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 84, 6, 35, 10, 10, 2, 1

Offset: 1

Gus Wiseman, Apr 29 2021

First differs from A343656 at A(4,2) = 6, A343656(4,2) = 5.

As a triangle, T(n,k) = number of ways to choose a multiset of n - k divisors of k.

			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  10  20  35  56  84 120 165
  n=7:  1   2   3   4   5   6   7   8   9
  n=8:  1   4  10  20  35  56  84 120 165
  n=9:  1   3   6  10  15  21  28  36  45
Triangle 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  10   2   1
   1   8   7  21   5  20   3   4   1
   1   9   8  28   6  35   4  10   3   1
   1  10   9  36   7  56   5  20   6   4   1
   1  11  10  45   8  84   6  35  10  10   2   1
For example, row n = 6 counts the following multisets:
  {1,1,1,1,1}  {1,1,1,1}  {1,1,1}  {1,1}  {1}  {}
               {1,1,1,2}  {1,1,3}  {1,2}  {5}
               {1,1,2,2}  {1,3,3}  {1,4}
               {1,2,2,2}  {3,3,3}  {2,2}
               {2,2,2,2}           {2,4}
                                   {4,4}
Note that for n = 6, k = 4 in the triangle, the two multisets {1,4} and {2,2} represent the same divisor 4, so they are only counted once under A343656(4,2) = 5.

Andrew Howroyd, Table of n, a(n) for n = 1..1275 (first 50 antidiagonals)

Row k = 1 of the array is A000005.
Column n = 4 of the array is A000217.
Column n = 6 of the array is A000292.
Row k = 2 of the array is A184389.
The distinct products of these multisets are counted by A343656.
Antidiagonal sums of the array (or row sums of the triangle) are A343661.
A000312 = n^n.
A009998(n,k) = n^k (as an array, offset 1).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.
Cf. A000169, A062319, A067824, A143773, A146291, A176029, A285572, A326077, A327527, A334996, A343652, A343657.

multchoo[n_,k_]:=Binomial[n+k-1,k];
Table[multchoo[DivisorSigma[0,k],n-k],{n,10},{k,n}]

A(n,k) = binomial(numdiv(n) + k - 1, k)
{ for(n=1, 9, for(k=0, 8, print1(A(n,k), ", ")); print ) } \\ Andrew Howroyd, Jan 11 2024

A(n,k) = ((A000005(n), k)) = A007318(A000005(n) + k - 1, k).

T(n,k) = ((A000005(k), n - k)) = A007318(A000005(k) + n - k - 1, n - k).

A343657 Sum of number of divisors of x^y for each x >= 1, y >= 0, x + y = n.

Original entry on oeis.org

1, 2, 4, 7, 12, 18, 27, 39, 56, 77, 103, 134, 174, 223, 283, 356, 445, 547, 666, 802, 959, 1139, 1344, 1574, 1835, 2128, 2454, 2815, 3213, 3648, 4126, 4653, 5239, 5888, 6608, 7407, 8298, 9288, 10385, 11597, 12936, 14408, 16025, 17799, 19746, 21882, 24221

Offset: 1

Gus Wiseman, Apr 29 2021

nonn

			The a(7) = 27 divisors:
  1  32  81  64  25  6  1
     16  27  32  5   3
     8   9   16  1   2
     4   3   8       1
     2   1   4
     1       2
             1

Chai Wah Wu, Table of n, a(n) for n = 1..10000

Antidiagonal row sums (row sums of the triangle) of A343656.
Dominated by A343661.
A000005(n) counts divisors of n.
A000312(n) = n^n.
A007318(n,k) counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481(n,k) counts k-multisets of elements of {1..n}.
A343658(n,k) counts k-multisets of divisors of n.
Cf. A000169, A000272, A002064, A002109, A048691, A062319, A066959, A143773, A146291, A176029, A251683, A282935, A326358, A327527, A334996.

Total/@Table[DivisorSigma[0,k^(n-k)],{n,30},{k,n}]

a(n) = Sum_{k=1..n} A000005(k^(n-k)).

A343661 Sum of numbers of y-multisets of divisors of x for each x >= 1, y >= 0, x + y = n.

Original entry on oeis.org

1, 2, 4, 7, 12, 19, 30, 46, 70, 105, 155, 223, 316, 443, 619, 865, 1210, 1690, 2354, 3263, 4497, 6157, 8368, 11280, 15078, 19989, 26296, 34356, 44626, 57693, 74321, 95503, 122535, 157101, 201377, 258155, 330994, 424398, 544035, 696995, 892104, 1140298, 1455080

Offset: 1

Gus Wiseman, Apr 30 2021

nonn

			The a(5) = 12 multisets of divisors:
  {1,1,1,1}  {1,1,1}  {1,1}  {1}  {}
             {1,1,2}  {1,3}  {2}
             {1,2,2}  {3,3}  {4}
             {2,2,2}

Antidiagonal sums of the array A343658 (or row sums of the triangle).
Dominates A343657.
A000005 counts divisors.
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.
A343656 counts divisors of powers.
Cf. A000169, A000312, A009998, A062319, A067824, A143773, A146291, A176029, A184389, A285572, A326077, A327527, A334996, A343652, A343657.

multchoo[n_,k_]:=Binomial[n+k-1,k];
Table[Sum[multchoo[DivisorSigma[0,k],n-k],{k,n}],{n,10}]

a(n) = Sum_{k=1..n} binomial(sigma(k) + n - k - 1, n - k).

A343939 Number of n-chains of divisors of n.

Original entry on oeis.org

1, 3, 4, 15, 6, 49, 8, 165, 55, 121, 12, 1183, 14, 225, 256, 4845, 18, 3610, 20, 4851, 484, 529, 24, 73125, 351, 729, 4060, 12615, 30, 29791, 32, 435897, 1156, 1225, 1296, 494209, 38, 1521, 1600, 505981, 42, 79507, 44, 46575, 49726, 2209, 48

Offset: 1

Gus Wiseman, May 05 2021

nonn

			The a(1) = 1 through a(5) = 6 chains:
  (1)  (1/1)  (1/1/1)  (1/1/1/1)  (1/1/1/1/1)
       (2/1)  (3/1/1)  (2/1/1/1)  (5/1/1/1/1)
       (2/2)  (3/3/1)  (2/2/1/1)  (5/5/1/1/1)
              (3/3/3)  (2/2/2/1)  (5/5/5/1/1)
                       (2/2/2/2)  (5/5/5/5/1)
                       (4/1/1/1)  (5/5/5/5/5)
                       (4/2/1/1)
                       (4/2/2/1)
                       (4/2/2/2)
                       (4/4/1/1)
                       (4/4/2/1)
                       (4/4/2/2)
                       (4/4/4/1)
                       (4/4/4/2)
                       (4/4/4/4)

Diagonal n = k - 1 of the array A077592.
Chains of length n - 1 are counted by A163767.
Diagonal n = k of the array A334997.
The version counting all multisets of divisors (not just chains) is A343935.
A000005(n) counts divisors of n.
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-1) counts strict k-chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict k-chains of divisors from n to 1.
A337255(n,k) counts strict k-chains of divisors starting with n.
A343658(n,k) counts k-multisets of divisors of n.
A343662(n,k) counts strict k-chains of divisors of n (row sums: A337256).
Cf. A062319, A066959, A143773, A176029, A327527, A343656.

Table[Length[Select[Tuples[Divisors[n],n],OrderedQ[#]&&And@@Divisible@@@Reverse/@Partition[#,2,1]&]],{n,10}]

A046688 Antidiagonals of square array in which k-th row (k>0) is an arithmetic progression of difference 2^(k-1).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 5, 1, 1, 5, 7, 9, 9, 1, 1, 6, 9, 13, 17, 17, 1, 1, 7, 11, 17, 25, 33, 33, 1, 1, 8, 13, 21, 33, 49, 65, 65, 1, 1, 9, 15, 25, 41, 65, 97, 129, 129, 1, 1, 10, 17, 29, 49, 81, 129, 193, 257, 257, 1, 1, 11, 19, 33, 57, 97, 161, 257, 385, 513, 513, 1

Offset: 0

N. J. A. Sloane

			From _Gus Wiseman_, May 08 2021: (Start):
Array A(m,n) = 1 + n*2^(m-1) begins:
       n=0: n=1: n=2: n=3: n=4: n=5: n=6: n=7: n=8: n=9:
  m=0:   1    1    1    1    1    1    1    1    1    1
  m=1:   1    2    3    5    9   17   33   65  129  257
  m=2:   1    3    5    9   17   33   65  129  257  513
  m=3:   1    4    7   13   25   49   97  193  385  769
  m=4:   1    5    9   17   33   65  129  257  513 1025
  m=5:   1    6   11   21   41   81  161  321  641 1281
  m=6:   1    7   13   25   49   97  193  385  769 1537
  m=7:   1    8   15   29   57  113  225  449  897 1793
  m=8:   1    9   17   33   65  129  257  513 1025 2049
  m=9:   1   10   19   37   73  145  289  577 1153 2305
Triangle T(n,k) = 1 + (n-k)*2^(k-1) begins:
   1
   1   1
   1   2   1
   1   3   3   1
   1   4   5   5   1
   1   5   7   9   9   1
   1   6   9  13  17  17   1
   1   7  11  17  25  33  33   1
   1   8  13  21  33  49  65  65   1
   1   9  15  25  41  65  97 129 129   1
   1  10  17  29  49  81 129 193 257 257   1
   1  11  19  33  57  97 161 257 385 513 513   1
(End)

G. H. Hardy, A Theorem Concerning the Infinite Cardinal Numbers, Quart. J. Math., 35 (1904), p. 90 = Collected Papers, Vol. VII, p. 430.

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

Row sums are A000079.
Diagonal n = m + 1 of the array is A002064.
Diagonal n = m of the array is A005183.
Column m = 1 of the array is A094373.
Diagonal n = m - 1 of the array is A131056.
A002109 gives hyperfactorials (sigma: A260146, omega: A303281).
A009998(k,n) = n^k.
A009999(n,k) = n^k.
A057156 = (2^n)^(2^n).
A062319 counts divisors of n^n.
Cf. A000169, A000272, A000312, A036289, A343656, A343658.

Table[If[k==0,1,n*2^(k-1)+1],{n,0,9},{k,0,9}] (* ARRAY, Gus Wiseman, May 08 2021 *)
Table[If[k==0,1,1+(n-k)*2^(k-1)],{n,0,10},{k,0,n}] (* TRIANGLE, Gus Wiseman, May 08 2021 *)

A(m,n)={if(m>0, 1+n*2^(m-1), 1)}
{ for(m=0, 10, for(n=0, 10, print1(A(m,n), ", ")); print) } \\ Andrew Howroyd, Mar 07 2020

A(m,n) = 1 + n*2^(m-1) for m > 1. - Andrew Howroyd, Mar 07 2020

As a triangle, T(n,k) = A(k,n-k) = 1 + (n-k)*2^(k-1). - Gus Wiseman, May 08 2021

More terms from Larry Reeves (larryr(AT)acm.org), Apr 06 2000

cp's OEIS Frontend

A002064 Cullen numbers: a(n) = n*2^n + 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

PARI

PARI

Python

Formula

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A343658 Array read by antidiagonals where A(n,k) is the number of ways to choose a multiset of k divisors of n.

Views