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 74 results. Next

A084859 Binomial transform of Cullen numbers A002064.

Original entry on oeis.org

1, 4, 16, 62, 232, 842, 2980, 10334, 35248, 118610, 394684, 1301126, 4255624, 13825658, 44657428, 143521838, 459230560, 1463719586, 4649308012, 14722502870, 46491507256, 146447041994, 460259735236, 1443537130622
Offset: 0

Views

Author

Paul Barry, Jun 12 2003

Keywords

Crossrefs

Cf. A004277.

Programs

  • Mathematica
    LinearRecurrence[{8,-21,18},{1,4,16},30] (* Harvey P. Dale, Mar 31 2015 *)

Formula

a(n)=2^n+2n*3^(n-1); a(n)=sum{k=0..n, C(n, k)(1+k*2^k) }.
O.g.f.: (1-4x+5x^2)/[(1-2x)(-1+3x)^2]. - R. J. Mathar, Apr 02 2008
a(0)=1, a(1)=4, a(2)=16, a(n)=8*a(n-1)-21*a(n-2)+18*a(n-3). - Harvey P. Dale, Mar 31 2015

A056820 Euler totient function of Cullen numbers (A002064), n*2^n +1.

Original entry on oeis.org

1, 2, 6, 20, 48, 132, 240, 528, 1364, 4180, 7560, 20784, 42768, 70992, 151632, 446208, 986880, 1688400, 4006800, 6080256, 13843440, 42125380, 79567488, 152712000, 292573440, 526051328, 797636448, 2942476992, 7257012896
Offset: 0

Views

Author

Robert G. Wilson v, Aug 31 2000

Keywords

Crossrefs

Programs

  • Mathematica
    Table[ EulerPhi[ n*2^n + 1 ], {n, 1, 30} ]

Formula

a(n) = A000010(A002064(n)). - Amiram Eldar, Nov 29 2019

A001787 a(n) = n*2^(n-1).

Original entry on oeis.org

0, 1, 4, 12, 32, 80, 192, 448, 1024, 2304, 5120, 11264, 24576, 53248, 114688, 245760, 524288, 1114112, 2359296, 4980736, 10485760, 22020096, 46137344, 96468992, 201326592, 419430400, 872415232, 1811939328, 3758096384, 7784628224, 16106127360, 33285996544
Offset: 0

Views

Author

Keywords

Comments

Number of edges in an n-dimensional hypercube.
Number of 132-avoiding permutations of [n+2] containing exactly one 123 pattern. - Emeric Deutsch, Jul 13 2001
Number of ways to place n-1 nonattacking kings on a 2 X 2(n-1) chessboard for n >= 2. - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), May 22 2001
Arithmetic derivative of 2^n: a(n) = A003415(A000079(n)). - Reinhard Zumkeller, Feb 26 2002
(-1) times the determinant of matrix A_{i,j} = -|i-j|, 0 <= i,j <= n.
a(n) is the number of ones in binary numbers 1 to 111...1 (n bits). a(n) = A000337(n) - A000337(n-1) for n = 2,3,... . - Emeric Deutsch, May 24 2003
The number of 2 X n 0-1 matrices containing n+1 1's and having no zero row or column. The number of spanning trees of the complete bipartite graph K(2,n). This is the case m = 2 of K(m,n). See A072590. - W. Edwin Clark, May 27 2003
Binomial transform of 0,1,2,3,4,5,... (A001477). Without the initial 0, binomial transform of odd numbers.
With an additional leading zero, [0,0,1,4,...] this is the binomial transform of the integers repeated A004526. Its formula is then (2^n*(n-1) + 0^n)/4. - Paul Barry, May 20 2003
Number of zeros in all different (n+1)-bit integers. - Ralf Stephan, Aug 02 2003
From Lekraj Beedassy, Jun 03 2004: (Start)
Final element of a summation table (as opposed to a difference table) whose first row consists of integers 0 through n (or first n+1 nonnegative integers A001477); illustrating the case n=5:
0 1 2 3 4 5
1 3 5 7 9
4 8 12 16
12 20 28
32 48
80
and the final element is a(5)=80. (End)
This sequence and A001871 arise in counting ordered trees of height at most k where only the rightmost branch at the root actually achieves this height and the count is by the number of edges, with k = 3 for this sequence and k = 4 for A001871.
Let R be a binary relation on the power set P(A) of a set A having n = |A| elements such that for all elements x,y of P(A), xRy if x is a proper subset of y and there are no z in P(A) such that x is a proper subset of z and z is a proper subset of y. Then a(n) = |R|. - Ross La Haye, Sep 21 2004
Number of 2 X n binary matrices avoiding simultaneously the right-angled numbered polyomino patterns (ranpp) (00;1) and (10;1). An occurrence of a ranpp (xy;z) in a matrix A=(a(i,j)) is a triple (a(i1,j1), a(i1,j2), a(i2,j1)) where i1 < i2, j1 < j2 and these elements are in same relative order as those in the triple (x,y,z). - Sergey Kitaev, Nov 11 2004
Number of subsequences 00 in all binary words of length n+1. Example: a(2)=4 because in 000,001,010,011,100,101,110,111 the sequence 00 occurs 4 times. - Emeric Deutsch, Apr 04 2005
If you expand the n-factor expression (a+1)*(b+1)*(c+1)*...*(z+1), there are a(n) variables in the result. For example, the 3-factor expression (a+1)*(b+1)*(c+1) expands to abc+ab+ac+bc+a+b+c+1 with a(3) = 12 variables. - David W. Wilson, May 08 2005
An inverse Chebyshev transform of n^2, where g(x)->(1/sqrt(1-4*x^2))*g(x*c(x^2)), c(x) the g.f. of A000108. - Paul Barry, May 13 2005
Sequences A018215 and A058962 interleaved. - Graeme McRae, Jul 12 2006
The number of never-decreasing positive integer sequences of length n with a maximum value of 2*n. - Ben Paul Thurston, Nov 13 2006
Total size of all the subsets of an n-element set. For example, a 2-element set has 1 subset of size 0, 2 subsets of size 1 and 1 of size 2. - Ross La Haye, Dec 30 2006
Convolution of the natural numbers [A000027] and A045623 beginning [0,1,2,5,...]. - Ross La Haye, Feb 03 2007
If M is the matrix (given by rows) [2,1;0,2] then the sequence gives the (1,2) entry in M^n. - Antonio M. Oller-Marcén, May 21 2007
If X_1,X_2,...,X_n is a partition of a 2n-set X into 2-blocks then, for n > 0, a(n) is equal to the number of (n+1)-subsets of X intersecting each X_i (i=1,2,...,n). - Milan Janjic, Jul 21 2007
Number of n-permutations of 3 objects u,v,w, with repetition allowed, containing exactly one u. Example: a(2)=4 because we have uv, vu, uw and wu. - Zerinvary Lajos, Dec 27 2007
A member of the family of sequences defined by a(n) = n*[c(1)*...*c(r)]^(n-1); c(i) integer. This sequence has c(1)=2, A027471 has c(1)=3. - Ctibor O. Zizka, Feb 23 2008
a(n) is the number of ways to split {1,2,...,n-1} into two (possibly empty) complementary intervals {1,2,...,i} and {i+1,i+2,...,n-1} and then select a subset from each interval. - Geoffrey Critzer, Jan 31 2009
Equals the Jacobsthal sequence A001045 convolved with A003945: (1, 3, 6, 12, ...). - Gary W. Adamson, May 23 2009
Starting with offset 1 = A059570: (1, 2, 6, 14, 34, ...) convolved with (1, 2, 2, 2, ...). - Gary W. Adamson, May 23 2009
Equals the first left hand column of A167591. - Johannes W. Meijer, Nov 12 2009
The number of tatami tilings of an n X n square with n monomers is n*2^(n-1). - Frank Ruskey, Sep 25 2010
Under T. D. Noe's variant of the hypersigma function, this sequence gives hypersigma(2^n): a(n) = A191161(A000079(n)). - Alonso del Arte, Nov 04 2011
Number of Dyck (n+2)-paths with exactly one valley at height 1 and no higher valley. - David Scambler, Nov 07 2011
Equals triangle A059260 * A016777 as a vector, where A016777 = (3n + 1): [1, 4, 7, 10, 13, ...]. - Gary W. Adamson, Mar 06 2012
Main transitions in systems of n particles with spin 1/2 (see A212697 with b=2). - Stanislav Sykora, May 25 2012
Let T(n,k) be the triangle with (first column) T(n,1) = 2*n-1 for n >= 1, otherwise T(n,k) = T(n,k-1) + T(n-1,k-1), then a(n) = T(n,n). - J. M. Bergot, Jan 17 2013
Sum of all parts of all compositions (ordered partitions) of n. The equivalent sequence for partitions is A066186. - Omar E. Pol, Aug 28 2013
Starting with a(1)=1: powers of 2 (A000079) self-convolved. - Bob Selcoe, Aug 05 2015
Coefficients of the series expansion of the normalized Schwarzian derivative -S{p(x)}/6 of the polynomial p(x) = -(x-x1)*(x-x2) with x1 + x2 = 1 (cf. A263646). - Tom Copeland, Nov 02 2015
a(n) is the number of North-East lattice paths from (0,0) to (n+1,n+1) that have exactly one east step below y = x-1 and no east steps above y = x+1. Details can be found in Pan and Remmel's link. - Ran Pan, Feb 03 2016
Also the number of maximal and maximum cliques in the n-hypercube graph for n > 0. - Eric W. Weisstein, Dec 01 2017
Let [n]={1,2,...,n}; then a(n-1) is the total number of elements missing in proper subsets of [n] that contain n to form [n]. For example, for n = 3, a(2) = 4 since the proper subsets of [3] that contain 3 are {3}, {1,3}, {2,3} and the total number of elements missing in these subsets to form [3] is 4: 2 in the first subset, 1 in the second, and 1 in the third. - Enrique Navarrete, Aug 08 2020
Number of 3-permutations of n elements avoiding the patterns 132, 231. See Bonichon and Sun. - Michel Marcus, Aug 19 2022

Examples

			a(2)=4 since 2314, 2341,3124 and 4123 are the only 132-avoiding permutations of 1234 containing exactly one increasing subsequence of length 3.
x + 4*x^2 + 12*x^3 + 32*x^4 + 80*x^5 + 192*x^6 + 448*x^7 + ...
a(5) = 1*0 + 5*1 + 10*2 + 10*3 + 5*4 + 1*5 = 80, with 1,5,10,10,5,1 the 5th row of Pascal's triangle. - _J. M. Bergot_, Apr 29 2014
		

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. 796.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 131.
  • Clifford A. Pickover, The Math Book, From Pythagoras to the 57th Dimension, 250 Milestones in the History of Mathematics, Sterling Publ., NY, 2009, page 282.
  • 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).

Crossrefs

Three other versions, essentially identical, are A085750, A097067, A118442.
Partial sums of A001792.
A058922(n+1) = 4*A001787(n).
Equals A090802(n, 1).
Column k=1 of A038207.
Row sums of A003506, A322427, A322428.

Programs

  • Haskell
    a001787 n = n * 2 ^ (n - 1)
    a001787_list = zipWith (*) [0..] $ 0 : a000079_list
    -- Reinhard Zumkeller, Jul 11 2014
    
  • Magma
    [n*2^(n-1): n in [0..40]]; // Vincenzo Librandi, Feb 04 2016
    
  • Maple
    spec := [S, {B=Set(Z, 0 <= card), S=Prod(Z, B, B)}, labeled]: seq(combstruct[count](spec, size=n), n=0..29); # Zerinvary Lajos, Oct 09 2006
    A001787:=1/(2*z-1)^2; # Simon Plouffe in his 1992 dissertation, dropping the initial zero
  • Mathematica
    Table[Sum[Binomial[n, i] i, {i, 0, n}], {n, 0, 30}] (* Geoffrey Critzer, Mar 18 2009 *)
    f[n_] := n 2^(n - 1); f[Range[0, 40]] (* Vladimir Joseph Stephan Orlovsky, Feb 09 2011 *)
    Array[# 2^(# - 1) &, 40, 0] (* Harvey P. Dale, Jul 26 2011 *)
    Join[{0}, Table[n 2^(n - 1), {n, 20}]] (* Eric W. Weisstein, Dec 01 2017 *)
    Join[{0}, LinearRecurrence[{4, -4}, {1, 4}, 20]] (* Eric W. Weisstein, Dec 01 2017 *)
    CoefficientList[Series[x/(-1 + 2 x)^2, {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
  • PARI
    {a(n) = if( n<0, 0, n * 2^(n-1))}
    
  • PARI
    concat(0, Vec(x/(1-2*x)^2 + O(x^50))) \\ Altug Alkan, Nov 03 2015
    
  • Python
    def A001787(n): return n*(1<Chai Wah Wu, Nov 14 2022

Formula

a(n) = Sum_{k=1..n} k*binomial(n, k). - Benoit Cloitre, Dec 06 2002
E.g.f.: x*exp(2x). - Paul Barry, Apr 10 2003
G.f.: x/(1-2*x)^2.
G.f.: x / (1 - 4*x / (1 + x / (1 - x))). - Michael Somos, Apr 07 2012
A108666(n) = Sum_{k=0..n} binomial(n, k)^2 * a(n). - Michael Somos, Apr 07 2012
PSumSIGN transform of A053220. PSumSIGN transform is A045883. Binomial transform is A027471(n+1). - Michael Somos, Jul 10 2003
Starting at a(1)=1, INVERT transform is A002450, INVERT transform of A049072, MOBIUS transform of A083413, PSUM transform is A000337, BINOMIAL transform is A081038, BINOMIAL transform of A005408. - Michael Somos, Apr 07 2012
a(n) = 2*a(n-1)+2^(n-1).
a(2*n) = n*4^n, a(2*n+1) = (2*n+1)4^n.
G.f.: x/det(I-x*M) where M=[1,i;i,1], i=sqrt(-1). - Paul Barry, Apr 27 2005
Starting 1, 1, 4, 12, ... this is 0^n + n2^(n-1), the binomial transform of the 'pair-reversed' natural numbers A004442. - Paul Barry, Jul 24 2003
Convolution of [1, 2, 4, 8, ...] with itself. - Jon Perry, Aug 07 2003
The signed version of this sequence, n(-2)^(n-1), is the inverse binomial transform of n(-1)^(n-1) (alternating sign natural numbers). - Paul Barry, Aug 20 2003
a(n-1) = (Sum_{k=0..n} 2^(n-k-1)*C(n-k, k)*C(1,(k+1)/2)*(1-(-1)^k)/2) - 0^n/4. - Paul Barry, Oct 15 2004
a(n) = Sum_{k=0..floor(n/2)} binomial(n, k)(n-2k)^2. - Paul Barry, May 13 2005
a(n+2) = A049611(n+2) - A001788(n).
a(n) = n! * Sum_{k=0..n} 1/((k - 1)!(n - k)!). - Paul Barry, Mar 26 2003
a(n+1) = Sum_{k=0..n} 4^k * A109466(n,k). - Philippe Deléham, Nov 13 2006
Row sums of A130300 starting (1, 4, 12, 32, ...). - Gary W. Adamson, May 20 2007
Equals row sums of triangle A134083. Equals A002064(n) + (2^n - 1). - Gary W. Adamson, Oct 07 2007
a(n) = 4*a(n-1) - 4*a(n-2), a(0)=0, a(1)=1. - Philippe Deléham, Nov 16 2008
Sum_{n>0} 1/a(n) = 2*log(2). - Jaume Oliver Lafont, Feb 10 2009
a(n) = A000788(A000225(n)) = A173921(A000225(n)). - Reinhard Zumkeller, Mar 04 2010
a(n) = n * A011782(n). - Omar E. Pol, Aug 28 2013
a(n-1) = Sum_{t_1+2*t_2+...+n*t_n=n} (t_1+t_2+...+t_n-1)*multinomial(t_1+t_2 +...+t_n,t_1,t_2,...,t_n). - Mircea Merca, Dec 06 2013
a(n+1) = Sum_{r=0..n} (2*r+1)*C(n,r). - J. M. Bergot, Apr 07 2014
a(n) = A007283(n)*n/6. - Enxhell Luzhnica, Apr 16 2016
a(n) = (A000225(n) + A000337(n))/2. - Anton Zakharov, Sep 17 2016
Sum_{n>0} (-1)^(n+1)/a(n) = 2*log(3/2) = 2*A016578. - Ilya Gutkovskiy, Sep 17 2016
a(n) = Sum_{k=0..n-1} Sum_{i=0..n-1} (i+1) * C(k,i). - Wesley Ivan Hurt, Sep 21 2017
a(n) = Sum_{i=1..n} Sum_{j=1..n} phi(i)*binomial(n, i*j). - Ridouane Oudra, Feb 17 2024

A000325 a(n) = 2^n - n.

Original entry on oeis.org

1, 1, 2, 5, 12, 27, 58, 121, 248, 503, 1014, 2037, 4084, 8179, 16370, 32753, 65520, 131055, 262126, 524269, 1048556, 2097131, 4194282, 8388585, 16777192, 33554407, 67108838, 134217701, 268435428, 536870883, 1073741794, 2147483617
Offset: 0

Views

Author

Rosario Salamone (Rosario.Salamone(AT)risc.uni-linz.ac.at)

Keywords

Comments

Number of permutations of degree n with at most one fall; called "Grassmannian permutations" by Lascoux and Schützenberger. - Axel Kohnert (Axel.Kohnert(AT)uni-bayreuth.de)
Number of different permutations of a deck of n cards that can be produced by a single shuffle. [DeSario]
Number of Dyck paths of semilength n having at most one long ascent (i.e., ascent of length at least two). Example: a(4)=12 because among the 14 Dyck paths of semilength 4, the only paths that have more than one long ascent are UUDDUUDD and UUDUUDDD (each with two long ascents). Here U = (1, 1) and D = (1, -1). Also number of ordered trees with n edges having at most one branch node (i.e., vertex of outdegree at least two). - Emeric Deutsch, Feb 22 2004
Number of {12,1*2*,21*}-avoiding signed permutations in the hyperoctahedral group.
Number of 1342-avoiding circular permutations on [n+1].
2^n - n is the number of ways to partition {1, 2, ..., n} into arithmetic progressions, where in each partition all the progressions have the same common difference and have lengths at least 1. - Marty Getz (ffmpg1(AT)uaf.edu) and Dixon Jones (fndjj(AT)uaf.edu), May 21 2005
if b(0) = x and b(n) = b(n-1) + b(n-1)^2*x^(n-2) for n > 0, then b(n) is a polynomial of degree a(n). - Michael Somos, Nov 04 2006
The chromatic invariant of the Mobius ladder graph M_n for n >= 2. - Jonathan Vos Post, Aug 29 2008
Dimension sequence of the dual alternative operad (i.e., associative and satisfying the identity xyz + yxz + zxy + xzy + yzx + zyx = 0) over the field of characteristic 3. - Pasha Zusmanovich, Jun 09 2009
An elephant sequence, see A175654. For the corner squares six A[5] vectors, with decimal values between 26 and 176, lead to this sequence (without the first leading 1). For the central square these vectors lead to the companion sequence A168604. - Johannes W. Meijer, Aug 15 2010
a(n+1) is also the number of order-preserving and order-decreasing partial isometries (of an n-chain). - Abdullahi Umar, Jan 13 2011
A040001(n) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, May 12 2012
A130103(n+1) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, May 12 2012
The number of labeled graphs with n vertices whose vertex set can be partitioned into a clique and a set of isolated points. - Alex J. Best, Nov 20 2012
For n > 0, a(n) is a B_2 sequence. - Thomas Ordowski, Sep 23 2014
See coefficients of the linear terms of the polynomials of the table on p. 10 of the Getzler link. - Tom Copeland, Mar 24 2016
Consider n points lying on a circle, then for n>=2 a(n-1) is the maximum number of ways to connect two points with non-intersecting chords. - Anton Zakharov, Dec 31 2016
Also the number of cliques in the (n-1)-triangular honeycomb rook graph. - Eric W. Weisstein, Jul 14 2017
From Eric M. Schmidt, Jul 17 2017: (Start)
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i) != e(j) < e(k). [Martinez and Savage, 2.7]
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i), e(j), e(k) pairwise distinct. [Martinez and Savage, 2.7]
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(j) >= e(k) and e(i) != e(k) pairwise distinct. [Martinez and Savage, 2.7]
(End)
Number of F-equivalence classes of Łukasiewicz paths. Łukasiewicz paths are F-equivalent iff the positions of pattern F are identical in these paths. - Sergey Kirgizov, Apr 08 2018
From Gus Wiseman, Feb 10 2019: (Start)
Also the number of connected partitions of an n-cycle. For example, the a(1) = 1 through a(4) = 12 connected partitions are:
{{1}} {{12}} {{123}} {{1234}}
{{1}{2}} {{1}{23}} {{1}{234}}
{{12}{3}} {{12}{34}}
{{13}{2}} {{123}{4}}
{{1}{2}{3}} {{124}{3}}
{{134}{2}}
{{14}{23}}
{{1}{2}{34}}
{{1}{23}{4}}
{{12}{3}{4}}
{{14}{2}{3}}
{{1}{2}{3}{4}}
(End)
Number of subsets of n-set without the single-element subsets. - Yuchun Ji, Jul 16 2019
For every prime p, there are infinitely many terms of this sequence that are divisible by p (see IMO Compendium link and Doob reference). Corresponding indices n are: for p = 2, even numbers A299174; for p = 3, A047257; for p = 5, A349767. - Bernard Schott, Dec 10 2021
Primes are in A081296 and corresponding indices in A048744. - Bernard Schott, Dec 12 2021

Examples

			G.f. = 1 + x + 2*x^2 + 5*x^3 + 12*x^4 + 27*x^5 + 58*x^6 + 121*x^7 + ...
		

References

  • Michael Doob, The Canadian Mathematical Olympiad & L'Olympiade Mathématique du Canada 1969-1993, Canadian Mathematical Society & Société Mathématique du Canada, Problem 4, 1983, page 158, 1993.

Crossrefs

Column 1 of triangle A008518.
Row sum of triangles A184049 and A184050.

Programs

  • Haskell
    a000325 n = 2 ^ n - n
    a000325_list = zipWith (-) a000079_list [0..]
    -- Reinhard Zumkeller, Jul 17 2012
    
  • Magma
    [2^n - n: n in [0..35]]; // Vincenzo Librandi, May 13 2011
    
  • Maple
    A000325 := proc(n) option remember; if n <=1 then n+1 else 2*A000325(n-1)+n-1; fi; end;
    g:=1/(1-2*z): gser:=series(g, z=0, 43): seq(coeff(gser, z, n)-n, n=0..31); # Zerinvary Lajos, Jan 09 2009
  • Mathematica
    Table[2^n - n, {n, 0, 39}] (* Alonso del Arte, Sep 15 2014 *)
    LinearRecurrence[{4, -5, 2}, {1, 2, 5}, {0, 20}] (* Eric W. Weisstein, Jul 14 2017 *)
  • PARI
    {a(n) = 2^n - n}; /* Michael Somos, Nov 04 2006 */
    
  • Python
    def A000325(n): return (1<Chai Wah Wu, Jan 11 2023

Formula

a(n+1) = 2*a(n) + n - 1, a(0) = 1. - Reinhard Zumkeller, Apr 12 2003
Binomial transform of 1, 0, 1, 1, 1, .... The sequence starting 1, 2, 5, ... has a(n) = 1 + n + 2*Sum_{k=2..n} binomial(n, k) = 2^(n+1) - n - 1. This is the binomial transform of 1, 1, 2, 2, 2, 2, .... a(n) = 1 + Sum_{k=2..n} C(n, k). - Paul Barry, Jun 06 2003
G.f.: (1-3x+3x^2)/((1-2x)*(1-x)^2). - Emeric Deutsch, Feb 22 2004
A107907(a(n+2)) = A000051(n+2) for n > 0. - Reinhard Zumkeller, May 28 2005
a(n+1) = sum of n-th row of the triangle in A109128. - Reinhard Zumkeller, Jun 20 2005
Row sums of triangle A133116. - Gary W. Adamson, Sep 14 2007
G.f.: 1 / (1 - x / (1 - x / ( 1 - x / (1 + x / (1 - 2*x))))). - Michael Somos, May 12 2012
First difference is A000225. PSUM transform is A084634. - Michael Somos, May 12 2012
a(n) = [x^n](B(x)^n-B(x)^(n-1)), n>0, a(0)=1, where B(x) = (1+2*x+sqrt(1+4*x^2))/2. - Vladimir Kruchinin, Mar 07 2014
E.g.f.: (exp(x) - x)*exp(x). - Ilya Gutkovskiy, Aug 07 2016
a(n) = A125128(n) - A000225(n) + 1. - Miquel Cerda, Aug 12 2016
a(n) = 2*A125128(n) - A095151(n) + 1. - Miquel Cerda, Aug 12 2016
a(n) = A079583(n-1) - A000225(n-1). - Miquel Cerda, Aug 15 2016
a(n)^2 - 4*a(n-1)^2 = (n-2)*(a(n)+2*a(n-1)). - Yuchun Ji, Jul 13 2018
a(n) = 2^(-n) * A186947(n) = 2^n * A002064(-n) for all n in Z. - Michael Somos, Jul 18 2018
a(2^n) = (2^a(n) - 1)*2^n. - Lorenzo Sauras Altuzarra, Feb 01 2022

A036289 a(n) = n*2^n.

Original entry on oeis.org

0, 2, 8, 24, 64, 160, 384, 896, 2048, 4608, 10240, 22528, 49152, 106496, 229376, 491520, 1048576, 2228224, 4718592, 9961472, 20971520, 44040192, 92274688, 192937984, 402653184, 838860800, 1744830464, 3623878656, 7516192768, 15569256448, 32212254720
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Right side of the binomial sum Sum_{i = 0..n} (n-2*i)^2 * binomial(n, i) = n*2^n. - Yong Kong (ykong(AT)curagen.com), Dec 28 2000
Let W be a binary relation on the power set P(A) of a set A having n = |A| elements such that for all elements x, y of P(A), xRy if x is a proper subset of y and there are no z in P(A) such that x is a proper subset of z and z is a proper subset of y, or y is a proper subset of x and there are no z in P(A) such that y is a proper subset of z and z is a proper subset of x. Then a(n) = |W|. - Ross La Haye, Sep 26 2007
Partial sums give A036799. - Vladimir Joseph Stephan Orlovsky, Jul 09 2011
a(n) = n with the bits shifted to the left by n places (new bits on the right hand side are zeros). - Indranil Ghosh, Jan 05 2017
Satisfies Benford's law [Theodore P. Hill, Personal communication, Feb 06, 2017]. - N. J. A. Sloane, Feb 08 2017
Also the circumference of the n-cube connected cycle graph. - Eric W. Weisstein, Sep 03 2017
a(n) is also the number of derangements in S_{n+3} with a descent set of {i, i+1} such that i ranges from 1 to n-2. - Isabella Huang, Mar 17 2018
a(n-1) is also the number of multiplications required to compute the permanent of general n X n matrices using Glynn's formula (see Theorem 2.1 in Glynn). - Stefano Spezia, Oct 27 2021

References

  • Arno Berger and Theodore P. Hill. An Introduction to Benford's Law. Princeton University Press, 2015.
  • A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992, Eq. (4.2.2.29)

Crossrefs

Equals 2*A001787. Equals A003261(n) + 1.

Programs

  • Haskell
    a036289 n = n * 2 ^ n
    a036289_list = zipWith (*) [0..] a000079_list
    -- Reinhard Zumkeller, Mar 05 2012
    
  • Maple
    g:=1/(1-2*z): gser:=series(g, z=0, 43): seq(coeff(gser, z, n)*n, n=0..34); # Zerinvary Lajos, Jan 11 2009
  • Mathematica
    Table[n*2^n, {n, 0, 50}] (* Vladimir Joseph Stephan Orlovsky, Mar 18 2010 *)
    LinearRecurrence[{4,-4},{0,2},40] (* Harvey P. Dale, Mar 02 2018 *)
  • PARI
    a(n)=n<Charles R Greathouse IV, Jun 15 2011
    
  • Python
    a=lambda n: n<Indranil Ghosh, Jan 05 2017

Formula

Main diagonal of array (A085454) defined by T(i, 1) = i, T(1, j) = 2j, T(i, j) = T(i-1, j) + T(i-1, j-1). - Benoit Cloitre, Aug 05 2003
Binomial transform of A005843, the even numbers. - Joshua Zucker, Jan 13 2006
G.f.: 2*x/(1-2*x)^2. - R. J. Mathar, Nov 21 2007
a(n) = A000079(n)*n. - Omar E. Pol, Dec 21 2008
E.g.f.: 2*x exp(2*x). - Geoffrey Critzer, Oct 03 2011
a(n) = A002064(n) - 1. - Reinhard Zumkeller, Mar 16 2013
From Vaclav Kotesovec, Feb 14 2015: (Start)
Sum_{n>=1} 1/a(n) = log(2).
Sum_{n>=1} (-1)^(n+1)/a(n) = log(3/2).
(End)

A057711 a(0)=0, a(1)=1, a(n) = n*2^(n-2) for n >= 2.

Original entry on oeis.org

0, 1, 2, 6, 16, 40, 96, 224, 512, 1152, 2560, 5632, 12288, 26624, 57344, 122880, 262144, 557056, 1179648, 2490368, 5242880, 11010048, 23068672, 48234496, 100663296, 209715200, 436207616, 905969664, 1879048192, 3892314112, 8053063680
Offset: 0

Views

Author

Bernhard Wolf (wolf(AT)cs.tu-berlin.de), Oct 24 2000

Keywords

Comments

Number of states in the planning domain FERRY, when n-3 cars are at one of two shores while the (n-2)nd car may be on the ferry or at one of the shores.
If the ferry could board any number of cars (instead of only one), the number of states would form the Pisot sequence P(2,6) (A008776). In addition, if k shores existed, the sequence would form the Pisot sequence P(k,k(k+1)). This corresponds to the BRIEFCASE planning domain.
a(i) is the number of occurrences of the number 1 in all palindromic compositions of n = 2*(i+1). - Silvia Heubach (sheubac(AT)calstatela.edu), Jan 10 2003. E.g., there are 5 palindromic compositions of 6, namely 111111 11211 2112 1221 141, containing a total of 16 1's.
Number of occurrences of 00's in all circular binary words of length n. Example: a(3)=6 because in the circular binary words 000, 001, 010, 011, 100, 101, 110 and 111 we have a total of 3+1+1+0+1+0+0+0=6 occurrences of 00. a(n) = Sum_{k=0..n} k*A119458(n,k). - Emeric Deutsch, May 20 2006
a(n) is the number of permutations on [n] for which the entries of each left factor form a circular subinterval of [n]. A subset I of [n] forms a circular subinterval of [n] if it is an ordinary interval [a,b] or has the form [1,a]-union-[b,n] for 1 <= a < b <= n. For example, (5,4,2) is a left factor of the permutation (5,4,2,1,3) which does not form a circular subinterval of [5] and a(4)=16 counts all 24 permutations of [4] except the eight whose first two entries are 1,3 (in either order) or 2,4. - David Callan, Mar 30 2007
a(n) is the total number of runs in all Boolean (n-1)-strings. For example, the 8 Boolean 3-strings, 000, 001, 010, 011, 100, 101, 110, 111 have 1, 2, 3, 2, 2, 3, 2, 1 runs respectively. - David Callan, Jul 22 2008
From Gary W. Adamson, Jul 31 2010: (Start)
Starting with "1" = (1, 2, 4, 8, ...) convolved with (1, 0, 2, 4, 8, ...).
Example: a(6) = 96 = (32, 16, 8, 4, 2, 1) dot (1, 0, 2, 4, 8, 16) = (32 + 0 + 16 + 16 + 16, + 16) = 32 + 4*16 (End)
An elephant sequence, see A175654. For the corner squares 24 A[5] vectors, with decimal values between 27 and 432, lead to this sequence (without the leading 0). For the central square these vectors lead to the companion sequence A087447 (without the first leading 1). - Johannes W. Meijer, Aug 15 2010
Starting with 1 = (1, 1, 2, 4, 8, 16, ...) convolved with (1, 1, 3, 7, 15, 31, ...). - Gary W. Adamson, Oct 26 2010
a(n) is the number of ways to draw simple polygonal chains for n vertices lying on a circle. - Anton Zakharov, Dec 31 2016
Also the number of edges, maximal cliques, and maximum cliques in the n-folded cube graph for n > 3. - Eric W. Weisstein, Dec 01 2017 and Mar 21 2018
Number of pairs of compositions of n corresponding to a seaweed algebra of index n-2 for n > 2. - Nick Mayers, Jun 25 2018
Starting with 1, 2, 6, 16, ..., number of permutations of length n>0 avoiding the partially ordered pattern (POP) {1>2, 1>3} of length 4. That is, number of length n permutations having no subsequences of length 4 in which the first element is larger than the second and third elements. - Sergey Kitaev, Dec 08 2020

Examples

			a(1)=6 because the palindromic compositions of n=4 are 4, 1+2+1, 1+1+1+1 and 2+2 and they contain 6 ones. - Silvia Heubach (sheubac(AT)calstatela.edu), Jan 10 2003
		

Crossrefs

Pisot sequence P(2, 6) (A008776), Pisot sequence P(k, k(k+1))
Cf. A119458.

Programs

  • Magma
    [Ceiling(n*2^(n-2)) : n in [0..40]]; // Vincenzo Librandi, Sep 22 2011
    
  • Mathematica
    Join[{0, 1}, Table[n 2^(n - 2), {n, 2, 30}]] (* Eric W. Weisstein, Dec 01 2017 *)
    Join[{0, 1}, LinearRecurrence[{4, -4}, {2, 6}, 20]] (* Eric W. Weisstein, Dec 01 2017 *)
    CoefficientList[Series[x (1 - 2 x + 2 x^2)/(1 - 2 x)^2, {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
  • PARI
    a(n)=ceil(n*2^(n-2)) \\ Charles R Greathouse IV, Oct 31 2011
    
  • PARI
    x='x+O('x^50); concat(0, Vec(x*(1-2*x+2*x^2)/(1-2*x)^2)) \\ Altug Alkan, Nov 01 2015

Formula

a(n) = ceiling(n*2^(n-2)).
Binomial transform of (0, 1, 0, 3, 0, 5, 0, 7, ...).
From Paul Barry, Apr 06 2003: (Start)
a(0)=0, a(n) = n*(0^(n-1) + 2^(n-1))/2, n > 0.
a(n) = Sum_{k=0..n} binomial(n, 2k+1)*(2k+1).
E.g.f.: x*exp(x)*cosh(x). (End)
The sequence 1, 1, 6, 16, ... is the binomial transform of A016813 with interpolated zeros. - Paul Barry, Jul 25 2003
For n > 1, a(n) = Sum_{k=0..n} (k-n/2)^2 C(n, k). (n+1)*a(n) = A001788(n). - Mario Catalani (mario.catalani(AT)unito.it), Nov 26 2003
From Paul Barry, May 07 2004: (Start)
a(n) = n*2^(n-2) - Sum_{k=0..n} binomial(n, k)*k*(-1)^k.
G.f.: x*(1-2*x+2*x^2)/(1-2*x)^2. (End)
a(n+1) = ceiling(binomial(n+1,1)*2^(n-1)). - Zerinvary Lajos, Nov 01 2006
a(n+1) = Sum_{k=0..n} A196389(n,k)*2^k. - Philippe Deléham, Oct 31 2011
a(0)=0, a(1)=1, a(2)=2, a(3)=6, a(n+1) = 4*a(n)-4*a(n-1) for n >= 3. - Philippe Deléham, Feb 20 2013
a(n) = A002064(n-1) - A002064(n-2), for n >= 2. - Ivan N. Ianakiev, Dec 29 2013
From Amiram Eldar, Aug 05 2020: (Start)
Sum_{n>=1} 1/a(n) = 4*log(2) - 1.
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(3/2) - 1. (End)

A003261 Woodall (or Riesel) numbers: n*2^n - 1.

Original entry on oeis.org

1, 7, 23, 63, 159, 383, 895, 2047, 4607, 10239, 22527, 49151, 106495, 229375, 491519, 1048575, 2228223, 4718591, 9961471, 20971519, 44040191, 92274687, 192937983, 402653183, 838860799, 1744830463, 3623878655, 7516192767, 15569256447, 32212254719, 66571993087
Offset: 1

Views

Author

Keywords

Comments

For n>1, a(n) is base at which zero is reached for the function "write f(j) in base j, read as base j+1 and then subtract 1 to give f(j+1)" starting from f(n) = n^2 - 1. - Henry Bottomley, Aug 06 2000
Sequence corresponds also to the maximum chain length of the classic puzzle whereby, under agreed commercial terms, an asset of unringed golden chain, when judiciously fragmented into as few as n pieces and n-1 opened links (through n-1 cuts), might be used to settle debt sequentially, with a golden link covering for unit cost. Here beside the n-1 opened links, the n fragmented pieces have lengths n, 2*n, 4*n, ..., 2^(n-1)*n. For instance, the chain of original length a(5)=159, if segregated by 4 cuts into 5+1+10+1+20+1+40+1+80, may be used to pay sequentially, i.e., a link-cost at a time, for an equivalent cost up to 159 links, to the same creditor. - Lekraj Beedassy, Feb 06 2003

Examples

			G.f. = x + 7*x^2 + 23*x^3 + 63*x^4 + 159*x^5 + 383*x^6 + 895*x^7 + ... - _Michael Somos_, Nov 04 2018
		

References

  • A. Brousseau, Number Theory Tables. Fibonacci Association, San Jose, CA, 1973, p. 159.
  • K. R. Bhutani and A. B. Levin, "The Problem of Sawing a Chain", Journal of Recreational Mathematics 2002-3 31(1) 32-35.
  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.
  • M. Gardner, Martin Gardner's Sixth Book of Mathematical Diversions from Scientific American, "Gold Links", Problem 4, pp. 50-51; 57-58, University of Chicago Press, 1983.
  • O. O'Shea, Mathematical Brainteasers with Surprising Solutions, Problem 76, pp. 183-185, Prometheus Books, Guilford, Connecticut, 2020.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 241.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = A036289(n) - 1 = A002064(n) - 2.
Cf. A133653.

Programs

Formula

G.f.: x*(-1-2*x+4*x^2) / ( (x-1)*(-1+2*x)^2 ). - Simon Plouffe in his 1992 dissertation
Binomial transform of A133653 and double binomial transform of [1, 5, -1, 1, -1, 1, ...]. - Gary W. Adamson, Sep 19 2007
a(n) = -(2)^n * A006127(-n) for all n in Z. - Michael Somos, Nov 04 2018
E.g.f.: 1 + exp(x)*(2*exp(x)*x - 1). - Stefano Spezia, Nov 24 2024

A004277 1 together with positive even numbers.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, 124, 126, 128, 130, 132
Offset: 0

Views

Author

Keywords

Comments

Also number of non-attacking bishops on n X n board. - Koksal Karakus (karakusk(AT)hotmail.com), May 27 2002
Engel expansion of e^(1/2) (see A006784 for definition) [when offset by 1]. - Henry Bottomley, Dec 18 2000
Numbers n such that a 2n-group (i.e., a group of order 2n) has subgroup C_2. - Lekraj Beedassy, Oct 14 2004
Image of 1/(1-2x) under the mapping g(x)->g(x/(1+x^2)). - Paul Barry, Jan 16 2005
Position of n in A113322: A113322(a(n-1)) = n for n>0. - Reinhard Zumkeller, Oct 26 2005
Incrementally largest terms in the continued fraction for e. - Nick Hobson, Jan 11 2007
Conjecturally, the differences of two consecutive primes (without repetition). - Juri-Stepan Gerasimov, Nov 09 2009
Equals (1, 2, 2, 2, ...) convolved with (1, 0, 2, 0, 2, 0, 2, ...). - Gary W. Adamson, Mar 03 2010
a(n) is the number of 0-dimensional elements (vertices) in an n-cross polytope. - Patrick J. McNab, Jul 06 2015
Numbers k such that in the symmetric representation of sigma(k) there is no pair bars as its ends (Cf. A237593). - Omar E. Pol, Sep 28 2018
Also, the coordination sequence of the L-lattice (see A332419). - Sean A. Irvine, Jul 29 2020

Crossrefs

INVERT transformation yields A098182 without A098182(0). - R. J. Mathar, Sep 11 2008

Programs

Formula

G.f.: (1+x^2)/(1-x)^2. - Paul Barry, Feb 28 2003
Inverse binomial transform of Cullen numbers A002064. a(n)=2n+0^n. - Paul Barry, Jun 12 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(n-k-1)*(-1)^k*2^(n-2k). - Paul Barry, Jan 16 2005
Equals binomial transform of [1, 1, 1, -1, 1, -1, 1, ...]. - Gary W. Adamson, Jul 15 2008
E.g.f.: 1+x*sinh(x) (aerated sequence). - Paul Barry, Oct 11 2009
a(n) = 0^n + 2*n = A000007(n) + A005843(n). - Reinhard Zumkeller, Jan 11 2012

Extensions

Corrected by Charles R Greathouse IV, Mar 18 2010

A005849 Indices of prime Cullen numbers: numbers k such that k*2^k + 1 is prime.

Original entry on oeis.org

1, 141, 4713, 5795, 6611, 18496, 32292, 32469, 59656, 90825, 262419, 361275, 481899, 1354828, 6328548, 6679881
Offset: 1

Views

Author

Keywords

Comments

From Amiram Eldar, Jun 05 2021: (Start)
The terms were found by:
a(1) - Cullen (1905). He found that there are no other terms up to 100 with the possible exception of 53. Cunningham (1906) showed that the 53rd Cullen number is composite and that the only possible term up to 200 is 141.
a(2) - Robinson (1958).
a(3)-a(6) - Keller (1995).
a(7)-a(8) - Masakatu Morii (1997).
a(9)-a(10) - Jeffrey Young (1997).
a(11)-a(12) - Darren Smith (1998).
a(13) - Masakatu Morii (1998).
a(14) - Mark Rodenkirch (2005).
a(15) - Dennis R. Gesker (2009).
a(16) - Magnus Bergman (2009). (End)

References

  • A. J. Cunningham, Solution of question 15897, Math. Quest. Educ. Times, Vol. 10 (1906), pp. 44-47.
  • Jean-Marie De Koninck, Ces nombres qui nous fascinent, Entry 141, p. 48, Ellipses, Paris 2008.
  • Harvey Dubner, Generalized Cullen numbers, J. Rec. Math., Vol. 21, No. 3 (1989), pp. 190-191.
  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section B20.
  • Paulo Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 283.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A002064, A002234, A050920, A173474 (complement).

Programs

  • Mathematica
    Select[Range[1000], PrimeQ[# 2^# + 1] &] (* Alonso del Arte, Jul 30 2017 *)
  • PARI
    is(n)=isprime(n<Charles R Greathouse IV, Feb 06 2017

Extensions

a(14) = 1354828 from old Proth Search pages by Mohammed Bouayoun (mohammed.bouayoun(AT)sanef.com), Apr 20 2006
The term 1467763 was added in error and has now been deleted; Jens Kruse Andersen, Nov 28 2007, remarks that 1467763 * 2^1467763 - 1 is a Woodall prime, but 3 divides the Cullen number 1467763 * 2^1467763 + 1.
6328548 from John Blazek, May 14 2009. He later reports that the search of the range from 6300000 to 6328548 was completed on May 28 2009.
Added a(16) = 6679881 from Caldwell's page, fixed broken link. - M. F. Hasler, Jan 18 2015
Name edited by Andrey Zabolotskiy and Felix Fröhlich, May 28 2021

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
Showing 1-10 of 74 results. Next