A001142 -id:A001142 - cp's OEIS frontend

A323444 Sum of exponents in prime-power factorization of Product_{k=0..n} binomial(n,k) (A001142).

0, 0, 1, 2, 6, 6, 11, 10, 23, 28, 33, 28, 45, 38, 44, 50, 86, 74, 96, 82, 106, 110, 114, 96, 147, 150, 153, 182, 211, 184, 215, 186, 281, 280, 279, 278, 347, 308, 306, 304, 380, 336, 374, 328, 368, 408, 403, 352, 489, 482, 524, 516, 559, 498, 596, 586, 686, 674

Offset: 0

Daniel Suteu, Jan 15 2019

nonn

Also sum of exponents in prime-power factorization of hyperfactorial(n) / superfactorial(n).

			a(4) = 6 because C(4,0)*C(4,1)*C(4,2)*C(4,3)*C(4,4) = 2^5 * 3^1 and 5 + 1 = 6, where C(n,k) is the binomial coefficient.

Amiram Eldar, Table of n, a(n) for n = 0..10000
Jeffrey C. Lagarias and Harsh Mehta, Products of binomial coefficients and unreduced Farey fractions, International Journal of Number Theory, Vol. 12, No. 1 (2016), pp. 57-91; arXiv preprint, arXiv:1409.4145 [math.NT], 2014-2015.
Eric Weisstein's World of Mathematics, Hyperfactorial.
Eric Weisstein's World of Mathematics, Superfactorial.
Index entries for sequences computed from exponents in factorization of n.
Index entries for sequences related to factorial numbers.

Cf. A000178, A001142, A001222, A002109, A004788, A007318, A022559, A303279, A303281.

Array[Sum[PrimeOmega@ Binomial[#, k], {k, 0, #}] &, 57] (* Michael De Vlieger, Jan 19 2019 *)

a(n) = sum(k=0, n, bigomega(binomial(n, k)));

a(n) = my(t=0); sum(k=1, n, my(b=bigomega(k)); t+=b; k*b-t);

first(n) = my(res = List([0]), r=0, t=0, b=0); for(k=1, n, b=bigomega(k); t += b; r += k*b-t; listput(res, r)); res \\ adapted from Daniel Suteu \\ David A. Corneth, Jan 16 2019

a(n) = A303281(n) - A303279(n), for n > 0.

a(n) = A001222(A001142(n)).

A092592 a(n) = A001142(n)/A002944(n), i.e., the product of C(n,j) binomial coefficients (for j=0..n) is divided by the least common multiple of them.

Original entry on oeis.org

1, 1, 3, 8, 250, 2700, 252105, 39337984, 46664771328, 12859560000000, 100271725559775000, 155537541908508672000, 8147102459616398435002560, 1112629444624748831692923019008, 198709463856994388404947123486328125

Offset: 1

Labos Elemer, Mar 10 2004

nonn

Cf. A001142, A002944.

a(n) = my(v = vector(n+1, j, binomial(n, j-1))); prod(k=1, #v, v[k])/ lcm(v); \\ Michel Marcus, Jul 29 2017

A000312 a(n) = n^n; number of labeled mappings from n points to themselves (endofunctions).

Original entry on oeis.org

1, 1, 4, 27, 256, 3125, 46656, 823543, 16777216, 387420489, 10000000000, 285311670611, 8916100448256, 302875106592253, 11112006825558016, 437893890380859375, 18446744073709551616, 827240261886336764177, 39346408075296537575424, 1978419655660313589123979

Offset: 0

N. J. A. Sloane

Also number of labeled pointed rooted trees (or vertebrates) on n nodes.
For n >= 1 a(n) is also the number of n X n (0,1) matrices in which each row contains exactly one entry equal to 1. - Avi Peretz (njk(AT)netvision.net.il), Apr 21 2001
Also the number of labeled rooted trees on (n+1) nodes such that the root is lower than its children. Also the number of alternating labeled rooted ordered trees on (n+1) nodes such that the root is lower than its children. - Cedric Chauve (chauve(AT)lacim.uqam.ca), Mar 27 2002
With p(n) = the number of integer partitions of n, p(i) = the number of parts of the i-th partition of n, d(i) = the number of different parts of the i-th partition of n, p(j, i) = the j-th part of the i-th partition of n, m(i, j) = multiplicity of the j-th part of the i-th partition of n, one has: a(n) = Sum_{i=1..p(n)} (n!/(Product_{j=1..p(i)} p(i, j)!)) * ((n!/(n - p(i)))!/(Product_{j=1..d(i)} m(i, j)!)). - Thomas Wieder, May 18 2005
All rational solutions to the equation x^y = y^x, with x < y, are given by x = A000169(n+1)/A000312(n), y = A000312(n+1)/A007778(n), where n = 1, 2, 3, ... . - Nick Hobson, Nov 30 2006
a(n) is the total number of leaves in all (n+1)^(n-1) trees on {0,1,2,...,n} rooted at 0. For example, with edges directed away from the root, the trees on {0,1,2} are {0->1,0->2},{0->1->2},{0->2->1} and contain a total of a(2)=4 leaves. - David Callan, Feb 01 2007
Limit_{n->infinity} A000169(n+1)/a(n) = exp(1). Convergence is slow, e.g., it takes n > 74 to get one decimal place correct and n > 163 to get two of them. - Alonso del Arte, Jun 20 2011
Also smallest k such that binomial(k, n) is divisible by n^(n-1), n > 0. - Michel Lagneau, Jul 29 2013
For n >= 2 a(n) is represented in base n as "one followed by n zeros". - R. J. Cano, Aug 22 2014
Number of length-n words over the alphabet of n letters. - Joerg Arndt, May 15 2015
Number of prime parking functions of length n+1. - Rui Duarte, Jul 27 2015
The probability density functions p(x, m=q, n=q, mu=1) = A000312(q)*E(x, q, q) and p(x, m=q, n=1, mu=q) = (A000312(q)/A000142(q-1))*x^(q-1)*E(x, q, 1), with q >= 1, lead to this sequence, see A163931, A274181 and A008276. - Johannes W. Meijer, Jun 17 2016
Satisfies Benford's law [Miller, 2015]. - N. J. A. Sloane, Feb 12 2017
A signed version of this sequence apart from the first term (1, -4, -27, 256, 3125, -46656, ...), has the following property: for every prime p == 1 (mod 2n), (-1)^(n(n-1)/2)*n^n = A057077(n)*a(n) is always a 2n-th power residue modulo p. - Jianing Song, Sep 05 2018
From Juhani Heino, May 07 2019: (Start)
n^n is both Sum_{i=0..n} binomial(n,i)*(n-1)^(n-i)
  and Sum_{i=0..n} binomial(n,i)*(n-1)^(n-i)*i.
The former is the familiar binomial distribution of a throw of n n-sided dice, according to how many times a required side appears, 0 to n. The latter is the same but each term is multiplied by its amount. This means that if the bank pays the player 1 token for each die that has the chosen side, it is always a fair game if the player pays 1 token to enter - neither bank nor player wins on average.
Examples:
2-sided dice (2 coins): 4 = 1 + 2 + 1 = 1*0 + 2*1 + 1*2 (0 omitted from now on);
3-sided dice (3 long triangular prisms): 27 = 8 + 12 + 6 + 1 = 12*1 + 6*2 + 1*3;
4-sided dice (4 long square prisms or 4 tetrahedrons): 256 = 81 + 108 + 54 + 12 + 1 = 108*1 + 54*2 + 12*3 + 1*4;
5-sided dice (5 long pentagonal prisms): 3125 = 1024 + 1280 + 640 + 160 + 20 + 1 = 1280*1 + 640*2 + 160*3 + 20*4 + 1*5;
6-sided dice (6 cubes): 46656 = 15625 + 18750 + 9375 + 2500 + 375 + 30 + 1 = 18750*1 + 9375*2 + 2500*3 + 375*4 + 30*5 + 1*6.
(End)
For each n >= 1 there is a graph on a(n) vertices whose largest independent set has size n and whose independent set sequence is constant (specifically, for each k=1,2,...,n, the graph has n^n independent sets of size k). There is no graph of smaller order with this property (Ball et al. 2019). - David Galvin, Jun 13 2019
For n >= 2 and 1 <= k <= n, a(n)*(n + 1)/4 + a(n)*(k - 1)*(n + 1 - k)/2*n is equal to the sum over all words w = w(1)...w(n) of length n over the alphabet {1, 2, ..., n} of the following quantity: Sum_{i=1..w(k)} w(i). Inspired by Problem 12432 in the AMM (see links). - Sela Fried, Dec 10 2023
Also, dimension of the unique cohomology group of the smallest interval containing the poset of partitions decorated by Perm, i.e. the poset of pointed partitions. - Bérénice Delcroix-Oger, Jun 25 2025

			G.f. = 1 + x + 4*x^2 + 27*x^3 + 256*x^4 + 3125*x^5 + 46656*x^6 + 823543*x^7 + ...

F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge, 1998, pp. 62, 63, 87.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 173, #39.
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.37)
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).

Kenny Lau, Table of n, a(n) for n = 0..385 [First 100 terms computed by T. D. Noe]
Taylor Ball, David Galvin, Katie Hyry, and Kyle Weingartner, Independent set and matching permutations, arXiv:1901.06579 [math.CO], 2019.
Arthur T. Benjamin and Fritz Juhnke, Another way of counting n^n, SIAM J. Discrete Math., Vol. 5, No. 3 (1992), pp. 377-379. - _N. J. A. Sloane_, Jun 09 2011
H. Bottomley, Illustration of initial terms.
H. J. Brothers and J. A. Knox, New closed-form approximations to the logarithmic constant e, The Mathematical Intelligencer, Vol. 20 (4), 1998, pp. 25-29. (Sequence appears as formula in Eq. (8))
C. Chauve, S. Dulucq and O. Guibert, Enumeration of some labeled trees, Proceedings of FPSAC/SFCA 2000 (Moscow), Springer, pp. 146-157.
Bérénice Delcroix-Oger and Clément Dupont, Lie-operads and operadic modules from poset cohomology, arXiv:2505.06094 [math.CO], 2025. See p. 34.
Frank Ellermann, Illustration of binomial transforms.
José María Grau and Antonio M. Oller-Marcén, On the last digit and the last non-zero digit of n^n in base b, Bulletin of the Korean Mathematical Society, Vol. 51, No. 5 (2014), pp. 1325-1337; arXiv preprint, arXiv:1203.4066 [math.NT], 2012.
Nick Hobson, Solution to puzzle 48: Exponential equation.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 36.
Steven J. Miller (ed.), Exercises to "The Theory and Applications of Benford's Law", Princeton University Press, 2015.
Mustafa Obaid et al., The number of complete exceptional sequences for a Dynkin algebra, arXiv preprint arXiv:1307.7573 [math.RT], 2013.
Franck Ramaharo, A generating polynomial for the pretzel knot, arXiv:1805.10680 [math.CO], 2018.
E. Vigren (Proposer), Problem 12432, Amer. Math. Monthly 130 (2023), p. 953.
Elena L. Wang and Guoce Xin, On Ward Numbers and Increasing Schröder Trees, arXiv:2507.15654 [math.CO], 2025. See pp. 12-13.
Eric Weisstein's World of Mathematics, Hadamard's Maximum Determinant Problem.
Eric Weisstein's World of Mathematics, Hankel Matrix.
Dimitri Zvonkine, An algebra of power series..., arXiv:math/0403092 [math.AG], 2004.
Index entries for "core" sequences
Index entries for sequences related to rooted trees
Index entries for sequences related to Benford's law

Cf. A000107, A000169, A000272, A001372, A007778, A007830, A008785-A008791, A019538, A048993, A008279, A085741, A062206, A212333.
First column of triangle A055858. Row sums of A066324.
Cf. A001923 (partial sums), A002109 (partial products), A007781 (first differences), A066588 (sum of digits).
Cf. A056665, A081721, A130293, A168658, A275549-A275558 (various classes of endofunctions).
Cf. A174824, A204688.
Cf. A055137, A083648.

a000312 n = n ^ n
a000312_list = zipWith (^) [0..] [0..]  -- Reinhard Zumkeller, Jul 07 2012

A000312 := n->n^n: seq(A000312(n), n=0..17);

Array[ #^# &, 16] (* Vladimir Joseph Stephan Orlovsky, May 01 2008 *)
Table[Sum[StirlingS2[n, i] i! Binomial[n, i], {i, 0, n}], {n, 0, 20}] (* Geoffrey Critzer, Mar 17 2009 *)
a[ n_] := If[ n < 1, Boole[n == 0], n^n]; (* Michael Somos, May 24 2014 *)
a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ 1 / (1 + LambertW[-x]), {x, 0, n}]]; (* Michael Somos, May 24 2014 *)
a[ n_] := If[n < 0, 0, n! SeriesCoefficient[ Nest[ 1 / (1 - x / (1 - Integrate[#, x])) &, 1 + O[x], n], {x, 0, n}]]; (* Michael Somos, May 24 2014 *)
a[ n_] := If[ n < 0, 0, With[{m = n + 1}, m! SeriesCoefficient[ InverseSeries[ Series[ (x - 1) Log[1 - x], {x, 0, m}]], m]]]; (* Michael Somos, May 24 2014 *)

A000312[n]:=if n=0 then 1 else n^n$
makelist(A000312[n],n,0,30); /* Martin Ettl, Oct 29 2012 */

PARI
```
{a(n) = n^n};
```

is(n)=my(b,k=ispower(n,,&b));if(k,for(e=1,valuation(k,b), if(k/b^e == e, return(1)))); n==1 \\ Charles R Greathouse IV, Jan 14 2013

{a(n) = my(A = 1 + O(x)); if( n<0, 0, for(k=1, n, A = 1 / (1 - x / (1 - intformal( A)))); n! * polcoeff( A, n))}; /* Michael Somos, May 24 2014 */

def A000312(n): return n**n # Chai Wah Wu, Nov 07 2022

a(n-1) = -Sum_{i=1..n} (-1)^i*i*n^(n-1-i)*binomial(n, i). - Yong Kong (ykong(AT)curagen.com), Dec 28 2000
E.g.f.: 1/(1 + W(-x)), W(x) = principal branch of Lambert's function.
a(n) = Sum_{k>=0} binomial(n, k)*Stirling2(n, k)*k! = Sum_{k>=0} A008279(n,k)*A048993(n,k) = Sum_{k>=0} A019538(n,k)*A007318(n,k). - Philippe Deléham, Dec 14 2003
E.g.f.: 1/(1 - T), where T = T(x) is Euler's tree function (see A000169).
a(n) = A000169(n+1)*A128433(n+1,1)/A128434(n+1,1). - Reinhard Zumkeller, Mar 03 2007
Comment on power series with denominators a(n): Let f(x) = 1 + Sum_{n>=1} x^n/n^n. Then as x -> infinity, f(x) ~ exp(x/e)*sqrt(2*Pi*x/e). - Philippe Flajolet, Sep 11 2008
E.g.f.: 1 - exp(W(-x)) with an offset of 1 where W(x) = principal branch of Lambert's function. - Vladimir Kruchinin, Sep 15 2010
a(n) = (n-1)*a(n-1) + Sum_{i=1..n} binomial(n, i)*a(i-1)*a(n-i). - Vladimir Shevelev, Sep 30 2010
With an offset of 1, the e.g.f. is the compositional inverse ((x - 1)*log(1 - x))^(-1) = x + x^2/2! + 4*x^3/3! + 27*x^4/4! + .... - Peter Bala, Dec 09 2011
a(n) = denominator((1 + 1/n)^n) for n > 0. - Jean-François Alcover, Jan 14 2013
a(n) = A089072(n,n) for n > 0. - Reinhard Zumkeller, Mar 18 2013
a(n) = (n-1)^(n-1)*(2*n) + Sum_{i=1..n-2} binomial(n, i)*(i^i*(n-i-1)^(n-i-1)), n > 1, a(0) = 1, a(1) = 1. - Vladimir Kruchinin, Nov 28 2014
log(a(n)) = lim_{k->infinity} k*(n^(1+1/k) - n). - Richard R. Forberg, Feb 04 2015
From Ilya Gutkovskiy, Jun 18 2016: (Start)
Sum_{n>=1} 1/a(n) = 1.291285997... = A073009.
Sum_{n>=1} 1/a(n)^2 = 1.063887103... = A086648.
Sum_{n>=1} n!/a(n) = 1.879853862... = A094082. (End)
A000169(n+1)/a(n) -> e, as n -> oo. - Daniel Suteu, Jul 23 2016
a(n) = n!*Product_{k=1..n} binomial(n, k)/Product_{k=1..n-1} binomial(n-1, k) = n!*A001142(n)/A001142(n-1). - Tony Foster III, Sep 05 2018
a(n-1) = abs(p_n(2-n)), for n > 2, the single local extremum of the n-th row polynomial of A055137 with Bagula's sign convention. - Tom Copeland, Nov 15 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = A083648. - Amiram Eldar, Jun 25 2021
Limit_{n->oo} (a(n+1)/a(n) - a(n)/a(n-1)) = e (see Brothers/Knox link). - Harlan J. Brothers, Oct 24 2021
Conjecture: a(n) = Sum_{i=0..n} A048994(n, i) * A048993(n+i, n) for n >= 0; proved by Mike Earnest, see link at A354797. - Werner Schulte, Jun 19 2022

A036987 Fredholm-Rueppel sequence.

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 0

N. J. A. Sloane

Binary representation of the Kempner-Mahler number Sum_{k>=0} 1/2^(2^k) = A007404.
Also a(n) = A(n) mod 2 where A is any of A001700, A005573, A007854, A026641, A049027, A064063, A064088, A064090, A064092, A064325, A064327, A064329, A064331, A064613, A076026, A105523, A123273, A126694, A126930, A126931, A126982, A126983, A126987, A127016, A127053, A127358, A127360, A127361, A127363. - Philippe Deléham, May 26 2007
a(n) = (product of digits of n; n in binary notation) mod 2. This sequence is a transformation of the Thue-Morse sequence (A010060), since there exists a function f such that f(sum of digits of n) = (product of digits of n). - Ctibor O. Zizka, Feb 12 2008
a(n-1), n >= 1, the characteristic sequence for powers of 2, A000079, is the unique solution of the following formal product and formal power series identity: Product_{j>=1} (1 + a(j-1)*x^j) = 1 + Sum_{k>=1} x^k = 1/(1-x). The product is therefore Product_{l>=1} (1 + x^(2^l)). Proof. Compare coefficients of x^n and use the binary representation of n. Uniqueness follows from the recurrence relation given for the general case under A147542. - Wolfdieter Lang, Mar 05 2009
a(n) is also the number of orbits of length n for the map x -> 1-cx^2 on [-1,1] at the Feigenbaum critical value c=1.401155... . - Thomas Ward, Apr 08 2009
A054525 (Mobius transform) * A001511 = A036987 = A047999^(-1) * A001511 = the inverse of Sierpiński's gasket * the ruler sequence. - Gary W. Adamson, Oct 26 2009 [Of course this is only vaguely correct depending on how the fuzzy indexing in these formulas is made concrete. - R. J. Mathar, Jun 20 2014]
Characteristic function of A000225. - Reinhard Zumkeller, Mar 06 2012
Also parity of the Catalan numbers A000108. - Omar E. Pol, Jan 17 2012
For n >= 2, also the largest exponent k >= 0 such that n^k in binary notation does not contain both 0 and 1. Unlike for the decimal version of this sequence, A062518, where the terms are only conjectural, for this sequence the values of a(n) can be proved to be the characteristic function of A000225, as follows: n^k will contain both 0 and 1 unless n^k = 2^r-1 for some r. But this is a special case of Catalan's equation x^p = y^q-1, which was proved by Preda Mihăilescu to have no nontrivial solution except 2^3 = 3^2 - 1. - Christopher J. Smyth, Aug 22 2014
Image, under the coding a,b -> 1; c -> 0, of the fixed point, starting with a, of the morphism a -> ab, b -> cb, c -> cc. - Jeffrey Shallit, May 14 2016
Number of nonisomorphic Boolean algebras of order n+1. - Jianing Song, Jan 23 2020

			G.f. = 1 + x + x^3 + x^7 + x^15 + x^31 + x^63 + x^127 + x^255 + x^511 + ...
a(7) = 1 since 7 = 2^3 - 1, while a(10) = 0 since 10 is not of the form 2^k - 1 for any integer k.

Antti Karttunen, Table of n, a(n) for n = 0..65537
D. Bailey et al., On the binary expansions of algebraic numbers, Journal de Théorie des Nombres de Bordeaux 16 (2004), 487-518.
Paul Barry, Some observations on the Rueppel sequence and associated Hankel determinants, arXiv:2005.04066 [math.CO], 2020.
Paul Barry, Conjectures and results on some generalized Rueppel sequences, arXiv:2107.00442 [math.CO], 2021.
Daniele A. Gewurz and Francesca Merola, Sequences realized as Parker vectors of oligomorphic permutation groups, J. Integer Seqs., Vol. 6, 2003.
Mats Granvik, Mathematica program for computing Fredholm Rueppel sequences
D. Kohel, S. Ling and C. Xing, Explicit Sequence Expansions, in Sequences and their Applications, C. Ding, T. Helleseth, and H. Niederreiter, eds., Proceedings of SETA'98 (Singapore, 1998), 308-317, 1999.
Preda Mihăilescu, Primary Cyclotomic Units and a Proof of Catalan's Conjecture, J. Reine angew. Math. 572 (2004): 167-195. doi:10.1515/crll.2004.048. MR 2076124.
H. Niederreiter and M. Vielhaber, Tree Complexity and a Doubly Exponential Gap between Structured and Random Sequences, J. Complexity, 12 (1996), 187-198.
Apisit Pakapongpun and Thomas Ward, Functorial orbit counting, Journal of Integer Sequences, 12 (2009) Article 09.2.4. [From _Thomas Ward_, Apr 08 2009]
Eric Rowland and Reem Yassawi, Profinite automata, arXiv:1403.7659 [math.DS], 2014. See p. 8.
E. Sheppard, net.math post (1985)
Stephen Wolfram, [Page 1092] A New Kind of Science | Online.
Index entries for characteristic functions

Cf. A007404, A043545, A062518, A078885, A078585, A078886, A078887, A078888, A078889, A078890.
The first row of A073346. Occurs for first time in A073202 as row 6 (and again as row 8).
Congruent to any of the sequences A000108, A007460, A007461, A007463, A007464, A061922, A068068 reduced modulo 2. Characteristic function of A000225.
If interpreted with offset=1 instead of 0 (i.e., a(1)=1, a(2)=1, a(3)=0, a(4)=1, ...) then this is the characteristic function of 2^n (A000079) and as such occurs as the first row of A073265. Also, in that case the INVERT transform will produce A023359.
This is Guy Steele's sequence GS(1, 3), also GS(3, 1) (see A135416).
Cf. A054525, A047999. - Gary W. Adamson, Oct 26 2009
Cf. A043529, A127802.

a036987 n = ibp (n+1) where
   ibp 1 = 1
   ibp n = if r > 0 then 0 else ibp n' where (n',r) = divMod n 2
a036987_list = 1 : f [0,1] where f (x:y:xs) = y : f (x:xs ++ [x,x+y])
-- Same list generator function as for a091090_list, cf. A091090.
-- Reinhard Zumkeller, May 19 2015, Apr 13 2013, Mar 13 2013

A036987:= n-> `if`(2^ilog2(n+1) = n+1, 1, 0):
seq(A036987(n), n=0..128);

RealDigits[ N[ Sum[1/10^(2^n), {n, 0, Infinity}], 110]][[1]]
(* Recurrence: *)
t[n_, 1] = 1; t[1, k_] = 1;
t[n_, k_] := t[n, k] =
  If[n < k, If[n > 1 && k > 1, -Sum[t[k - i, n], {i, 1, n - 1}], 0],
   If[n > 1 && k > 1, Sum[t[n - i, k], {i, 1, k - 1}], 0]];
Table[t[n, k], {k, n, n}, {n, 104}]
(* Mats Granvik, Jun 03 2011 *)
mb2d[n_]:=1 - Module[{n2 = IntegerDigits[n, 2]}, Max[n2] - Min[n2]]; Array[mb2d, 120, 0] (* Vincenzo Librandi, Jul 19 2019 *)
Table[PadRight[{1},2^k,0],{k,0,7}]//Flatten (* Harvey P. Dale, Apr 23 2022 *)

{a(n) =( n++) == 2^valuation(n, 2)}; /* Michael Somos, Aug 25 2003 */

a(n) = !bitand(n, n+1); \\ Ruud H.G. van Tol, Apr 05 2023

from sympy import catalan
def a(n): return catalan(n)%2 # Indranil Ghosh, May 25 2017

def A036987(n): return int(not(n&(n+1))) # Chai Wah Wu, Jul 06 2022

1 followed by a string of 2^k - 1 0's. Also a(n)=1 iff n = 2^m - 1.
a(n) = a(floor(n/2)) * (n mod 2) for n>0 with a(0)=1. - Reinhard Zumkeller, Aug 02 2002 [Corrected by Mikhail Kurkov, Jul 16 2019]
Sum_{n>=0} 1/10^(2^n) = 0.110100010000000100000000000000010...
1 if n=0, floor(log_2(n+1)) - floor(log_2(n)) otherwise. G.f.: (1/x) * Sum_{k>=0} x^(2^k) = Sum_{k>=0} x^(2^k-1). - Ralf Stephan, Apr 28 2003
a(n) = 1 - A043545(n). - Michael Somos, Aug 25 2003
a(n) = -Sum_{d|n+1} mu(2*d). - Benoit Cloitre, Oct 24 2003
Dirichlet g.f. for right-shifted sequence: 2^(-s)/(1-2^(-s)).
a(n) = A000108(n) mod 2 = A001405(n) mod 2. - Paul Barry, Nov 22 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n,k)*Sum_{j=0..k} binomial(k, 2^j-1). - Paul Barry, Jun 01 2006
A000523(n+1) = Sum_{k=1..n} a(k). - Mitch Harris, Jul 22 2011
a(n) = A209229(n+1). - Reinhard Zumkeller, Mar 07 2012
a(n) = Sum_{k=1..n} A191898(n,k)*cos(Pi*(n-1)*(k-1))/n; (conjecture). - Mats Granvik, Mar 04 2013
a(n) = A000035(A000108(n)). - Omar E. Pol, Aug 06 2013
a(n) = 1 iff n=2^k-1 for some k, 0 otherwise. - M. F. Hasler, Jun 20 2014
a(n) = ceiling(log_2(n+2)) - ceiling(log_2(n+1)). - Gionata Neri, Sep 06 2015
From John M. Campbell, Jul 21 2016: (Start)
a(n) = (A000168(n-1) mod 2).
a(n) = (A000531(n+1) mod 2).
a(n) = (A000699(n+1) mod 2).
a(n) = (A000891(n) mod 2).
a(n) = (A000913(n-1) mod 2), for n>1.
a(n) = (A000917(n-1) mod 2), for n>0.
a(n) = (A001142(n) mod 2).
a(n) = (A001246(n) mod 2).
a(n) = (A001246(n) mod 4).
a(n) = (A002057(n-2) mod 2), for n>1.
a(n) = (A002430(n+1) mod 2). (End)
a(n) = 2 - A043529(n). - Antti Karttunen, Nov 19 2017
a(n) = floor(1+log(n+1)/log(2)) - floor(log(2n+1)/log(2)). - Adriano Caroli, Sep 22 2019
This is also the decimal expansion of -Sum_{k>=1} mu(2*k)/(10^k - 1), where mu is the Möbius function (A008683). - Amiram Eldar, Jul 12 2020

Edited by M. F. Hasler, Jun 20 2014

A002944 a(n) = LCM(1,2,...,n) / n.

Original entry on oeis.org

1, 1, 2, 3, 12, 10, 60, 105, 280, 252, 2520, 2310, 27720, 25740, 24024, 45045, 720720, 680680, 12252240, 11639628, 11085360, 10581480, 232792560, 223092870, 1070845776, 1029659400, 2974571600, 2868336900, 80313433200, 77636318760, 2329089562800

Offset: 1

N. J. A. Sloane

Equals LCM of all numbers of (n-1)-st row of Pascal's triangle [Montgomery-Breusch]. - J. Lowell, Apr 16 2014. Corrected by N. J. A. Sloane, Sep 04 2019

Williams proves that a(n+1) = A034386(n) for n=2, 11 and 23 only. This is trivially the case for n=0 and 1, too. - Michel Marcus, Apr 16 2020

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 = 1..500
Bakir Farhi, An identity involving the least common multiple of binomial coefficients and its application, arXiv:0906.2295 [math.NT], 2009.
Bakir Farhi, An identity involving the least common multiple of binomial coefficients and its application, Amer. Math. Monthly, 116 (2009), 836-839.
Peter L. Montgomery (proposer) and Robert Breusch (solver), LCM of Binomial Coefficients, Problem E2686, American Mathematical Monthly, Vol. 84 (1977), p. 820.
Peter L. Montgomery (proposer) and Robert Breusch (solver), LCM of Binomial Coefficients, Solution to Problem E2686, American Mathematical Monthly, Vol. 86 (1979), p. 131.
Ian S. Williams, On a problem of Kurt Mahler concerning binomial coefficents (sic), Bulletin of the Australian Mathematical Society, Volume 14, Issue 2, April 1976, pp. 299-302.
Index entries for sequences related to lcm's

Cf. A025527 and A025537.
Cf. A100561, A003418, A001142, A001405. - Franklin T. Adams-Watters, Jul 05 2009
Cf. A056606 (squarefree kernel).

a002944 n = a003418 n `div` n  -- Reinhard Zumkeller, Mar 16 2015

A003418 := n-> lcm(seq(i,i=1..n)); f:=n->A003418(n)/n;
BB:=n->sum(1/sqrt(k), k=1..n): a:=n->floor(denom(BB(n))/n): seq(a(n), n=1..29); # Zerinvary Lajos, Mar 29 2007

Table[Apply[LCM,Range[n]]/n,{n,1,30}]  (* Geoffrey Critzer, Feb 10 2013 *)

a(n) = lcm(vector(n, i, i))/n; \\ Michel Marcus, Apr 16 2014

a(n) = A003418(n) / n.
a(n) = LCM of C(n-1, 0), C(n-1, 1), ..., C(n-1, n-1). [Montgomery-Breusch] [Corrected by N. J. A. Sloane, Jun 11 2008]
Equally, a(n+1) = LCM_{k=0..n} binomial(n,k). - Franklin T. Adams-Watters, Jul 05 2009

More terms from Jud McCranie, Jan 17 2000

Edited by N. J. A. Sloane, Jun 11 2008 and Sep 04 2019

A055775 a(n) = floor(n^n / n!).

Original entry on oeis.org

1, 1, 2, 4, 10, 26, 64, 163, 416, 1067, 2755, 7147, 18613, 48638, 127463, 334864, 881657, 2325750, 6145596, 16263866, 43099804, 114356611, 303761260, 807692034, 2149632061, 5726042115, 15264691107, 40722913454, 108713644516

Offset: 0

Henry Bottomley, Jul 12 2000

Stirling's approximation for n! suggests that this should be about e^n/sqrt(pi*2n). Bill Gosper has noted that e^n/sqrt(pi*(2n+1/3)) is significantly better.
n^n/n! = A001142(n)/A001142(n-1), where A001142(n) is product{k=0 to n} C(n,k) (where C() is a binomial coefficient). - Leroy Quet, May 01 2004
There are n^n distinct functions from [n] to [n] or sequences on n symbols of length n, the number of those sequences having n distinct symbols is n!. So the probability P(n) of bijection is n!/n^n. The expected value of the number of functions that we pick until we found a bijection is the reciprocal of P(n), or n^n/n!. - Washington Bomfim, Mar 05 2012

			a(5)=26 since 5^5=3125, 5!=120, 3125/120=26.0416666...

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Washington Bomfim, A method to find bijections from a set of n integers to {0,1, ... ,n-1}
Eric Weisstein's World of Mathematics, Stirling's Approximation for n!

Cf. A073225, A094082, A053042, A036679, A061711, A152170, A209081, A208846, A208847.

[Floor((n^n)/Factorial(n)): n in [0..30]]; // Vincenzo Librandi, Aug 22 2011

Join[{1}, Table[Floor[n^n/n!], {n, 30}]] (* Vladimir Joseph Stephan Orlovsky, Jan 15 2009 *)

a(n)=n^n\n! \\ Charles R Greathouse IV, Apr 17 2012

a(n) = floor(A000312(n)/A000142(n)).

More terms from James Sellers, Jul 13 2000

A007685 a(n) = Product_{k=1..n} binomial(2*k,k).

Original entry on oeis.org

1, 2, 12, 240, 16800, 4233600, 3911846400, 13425456844800, 172785629592576000, 8400837310791045120000, 1552105098192510332190720000, 1094904603628138948657963991040000, 2960792853328653706847125274154762240000, 30794022150329995743434211126374020153344000000

Offset: 0

N. J. A. Sloane

H. W. Gould, A class of binomial sums and a series transform, Utilitas Math., 45 (1994), 71-83.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

G. C. Greubel, Table of n, a(n) for n = 0..50
H. W. Gould, A class of binomial sums and a series transform, Utilitas Math., 45 (1994), 71-83. (Annotated scanned copy)
Bernd C. Kellner, Asymptotic products of binomial and multinomial coefficients revisited, Integers 24 (2024), Article #A59, 10 pp.; arXiv:2312.11369 [math.CO], 2023.

Cf. A000984, A001142, A112332, A203471, A268196, A296589, A296590, A338550.

[seq(mul(binomial(2*k,k),k=1..n),n=0..16)];

Table[Product[Binomial[2*k, k], {k, 1, n}], {n, 0, 50}] (* G. C. Greubel, Feb 02 2017 *)

a(n) = prod(k=1,n, binomial(2*k, k)); \\ Michel Marcus, Sep 18 2015

a(0) = 1, a(n) = (2^(2*n)*a(n - 1)*Gamma(n + 1/2))/(sqrt(Pi)*Gamma(n + 1)). - Ilya Gutkovskiy, Sep 18 2015
a(n) = (2^(n^2 + n - 1/24)*A^(3/2)*Pi^(-n/2 - 1/4)*BarnesG(n + 3/2))/(e^(1/8)*BarnesG(n + 2)), where A is the Glaisher-Kinkelin constant (A074962), BarnesG is the Barnes G-function. - Ilya Gutkovskiy, Sep 18 2015
a(n) ~ A^(3/2) * 2^(n^2 + n - 7/24) * exp(n/2 - 1/8) / (Pi^((n+1)/2) * n^(n/2 + 3/8)), where A = A074962 is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Apr 16 2016
For n>0, a(n) = 2^((n+1)/2) * sqrt(BarnesG(2*n)) * Gamma(2*n) / (n * BarnesG(n)^2 * Gamma(n)^(7/2)). - Vaclav Kotesovec, Apr 20 2024
Product_{1 <= j <= i <= n}  (i + j)/(i - j + 1). - Peter Bala, Oct 25 2024

A268196 a(n) = Product_{k=0..n} binomial(3*k,k).

Original entry on oeis.org

1, 3, 45, 3780, 1871100, 5618913300, 104309506501200, 12129109415959536000, 8920608231265175901456000, 41809329673499408044341517200000, 1256161937180234817183361549396758000000, 243113461110708695347467432844366521953760000000

Offset: 0

Vaclav Kotesovec, Apr 16 2016

Harvey P. Dale, Table of n, a(n) for n = 0..49
Eric Weisstein's World of Mathematics, Barnes G-Function.
Wikipedia, Barnes G-function.

Cf. A001142, A005809, A007685, A086205, A112332, A268504, A262261.

Table[Product[Binomial[3k,k],{k,0,n}],{n,0,12}]
FoldList[Times,Table[Binomial[3n,n],{n,0,15}]] (* Harvey P. Dale, Apr 23 2018 *)

a(n) = A^(7/6) * Gamma(1/3)^(1/3) * 3^(3*n^2/2 + 2*n + 11/36)* BarnesG(n + 4/3) * BarnesG(n + 5/3) / (exp(7/72) * 2^(n^2 + 2*n + 5/8) * Pi^(n/2 + 5/12) * BarnesG(n + 3/2) * BarnesG(n + 2)), where A = A074962 is the Glaisher-Kinkelin constant.
a(n) ~ A^(7/6) * Gamma(1/3)^(1/3) * 3^(11/36 + 2*n + 3*n^2/2) * exp(n/2 - 7/72) / (2^(n^2 + 2*n + 7/8) * Pi^(n/2 + 2/3) * n^(n/2 + 25/72)), where A = A074962 is the Glaisher-Kinkelin constant.
a(n) = A268504(n) / (A000178(n) * A098694(n)).

A249434 Integers m such that m! divides the product of elements on row m of Pascal's triangle.

Original entry on oeis.org

0, 1, 2, 4, 6, 10, 12, 16, 18, 22, 28, 30, 35, 36, 39, 40, 42, 46, 52, 58, 60, 62, 66, 70, 72, 78, 79, 82, 83, 88, 89, 96, 100, 102, 104, 106, 107, 108, 112, 126, 130, 131, 136, 138, 143, 148, 149, 150, 153, 156, 159, 162, 164, 166, 167, 172, 174, 175, 178, 179, 180, 181, 190, 192, 194, 196, 197, 198, 199, 207, 209, 210, 219, 222, 226, 228, 232, 238, 240, 250, 256

Offset: 1

Antti Karttunen, Nov 02 2014

nonn

Integers m such that A249151(m) >= m.
Equally: Integers m such that A249431(m) is nonnegative.
It seems that A006093 gives all those k for which A249151(k) = k. If that is true, then this is a disjoint union of A006093 and A249429.

			0! = 1 divides the product of binomial coefficients on row 0 of A007318, namely {1}, thus a(1) = 0.
1! = 1 divides the product of row 1 (1*1), thus a(2) = 1.
2! = 2 divides the product of row 2 (1*2*1), thus a(3) = 2.
3! = 6 does not divide the product of row 3 (1*3*3*1), but 4! = 24 divides the product of row 4 (1*4*6*4*1), as 96 = 4*24, thus a(4) = 4.

Antti Karttunen, Table of n, a(n) for n = 1..1421

Complement: A249433.
Subsequences: A006093 (conjectured), A249429, A249430, A249432.
Cf. A000142, A001142, A007318, A249151, A249431.

A249421 A(n,k) = exponent of the largest power of n-th prime which divides the product of the elements on row (k-1) of Pascal's triangle; a square array read by antidiagonals.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 5, 2, 0, 0, 0, 2, 1, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 4, 4, 0, 0, 0, 0, 0, 17, 2, 3, 0, 0, 0, 0, 0, 0, 10, 0, 2, 0, 0, 0, 0, 0, 0, 0, 12, 14, 1, 6, 0, 0, 0, 0, 0, 0, 0, 4, 10, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 18, 6, 8, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8, 13, 6, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 8, 4, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 1

Antti Karttunen, Oct 28 2014

Square array A(n,k), where n = row, k = column, read by antidiagonals: A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ....

A(n,k) is A000040(n)-adic valuation of A001142(k-1).

			The top left corner of the array:
0, 0, 1, 0, 5, 2, 4, 0, 17, 10, 12,  4, 18,  8, 11,  0, 49, 34, 36, 20, 42,
0, 0, 0, 2, 1, 0, 4, 2,  0, 14, 10,  6, 13,  8,  3, 12,  6,  0, 28, 20, 12,
0, 0, 0, 0, 0, 4, 3, 2,  1,  0,  8,  6,  4,  2,  0, 12,  9,  6,  3,  0, 16,
0, 0, 0, 0, 0, 0, 0, 6,  5,  4,  3,  2,  1,  0, 12, 10,  8,  6,  4,  2,  0,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0, 10,  9,  8,  7,  6,  5,  4,  3,  2,  1,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0,  0,  0, 12, 11, 10,  9,  8,  7,  6,  5,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 16, 15, 14, 13,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 18, 17,
...

Jeffrey C. Lagarias, Harsh Mehta, Products of binomial coefficients and unreduced Farey fractions, arXiv:1409.4145 [math.NT], 2014.

Transpose: A249422.
Row 1: A187059, Row 2: A249343, Row 3: A249345, Row 4 A249347. (Cf. also A249346).
Cf. A000040, A001142, A007318, A249344, A249151.

(define (A249421 n) (A249421bi (A002260 n) (A004736 n)))
(define (A249421bi row col) (A249344bi row (A001142 (- col 1)))) ;; Code for A249344bi given in A249344.

A(n, k) = A249344(n, A001142(k-1)).

cp's OEIS Frontend

A323444 Sum of exponents in prime-power factorization of Product_{k=0..n} binomial(n,k) (A001142).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

A092592 a(n) = A001142(n)/A002944(n), i.e., the product of C(n,j) binomial coefficients (for j=0..n) is divided by the least common multiple of them.

Views

Author

Keywords

Crossrefs

Programs

PARI

A000312 a(n) = n^n; number of labeled mappings from n points to themselves (endofunctions).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

PARI

PARI

Python

Formula

A036987 Fredholm-Rueppel sequence.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Python

Formula

Extensions

A002944 a(n) = LCM(1,2,...,n) / n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A055775 a(n) = floor(n^n / n!).

Views

Author

Keywords

Comments

Examples