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-4 of 4 results.

A001142 a(n) = Product_{k=1..n} k^(2k - 1 - n).

Original entry on oeis.org

1, 1, 2, 9, 96, 2500, 162000, 26471025, 11014635520, 11759522374656, 32406091200000000, 231627686043080250000, 4311500661703860387840000, 209706417310526095716965894400, 26729809777664965932590782608648192
Offset: 0

Views

Author

Keywords

Comments

Absolute value of determinant of triangular matrix containing binomial coefficients.
These are also the products of consecutive horizontal rows of the Pascal triangle. - Jeremy Hehn (ROBO_HEN5000(AT)rose.net), Mar 29 2007
Limit_{n->oo} a(n)*a(n+2)/a(n+1)^2 = e, as follows from lim_{n->oo} (1 + 1/n)^n = e. - Harlan J. Brothers, Nov 26 2009
A000225 gives the positions of odd terms. - Antti Karttunen, Nov 02 2014
Product of all unreduced fractions h/k with 1 <= k <= h <= n. - Jonathan Sondow, Aug 06 2015
a(n) is a product of the binomial coefficients from the n-th row of the Pascal triangle, for n= 0, 1, 2, ... For n > 0, a(n) means the number of all maximum chains in the poset formed by the n-dimensional Boolean cube {0,1}^n and the relation "precedes by weight". This relation is defined over {0,1}^n as follows: for arbitrary vectors u, v of {0,1}^n we say that "u precedes by weight v" if wt(u) < wt(v) or if u = v, where wt(u) denotes the (Hamming) weight of u. For more details, see the sequence A051459. - Valentin Bakoev, May 17 2019

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. 828.
  • 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

Cf. also A004788, A056606 (squarefree kernel), A256113.

Programs

  • GAP
    List([0..15], n-> Product([0..n], k-> Binomial(n,k) )); # G. C. Greubel, May 23 2019
    
  • Haskell
    a001142 = product . a007318_row -- Reinhard Zumkeller, Mar 16 2015
    
  • Magma
    [(&*[Binomial(n,k): k in [0..n]]): n in [0..15]]; // G. C. Greubel, May 23 2019
    
  • Maple
    a:=n->mul(binomial(n,k), k=0..n): seq(a(n), n=0..14); # Zerinvary Lajos, Jan 22 2008
  • Mathematica
    Table[Product[k^(2*k - 1 - n), {k, n}], {n, 0, 20}] (* Harlan J. Brothers, Nov 26 2009 *)
    Table[Hyperfactorial[n]/BarnesG[n+2], {n, 0, 20}] (* Peter Luschny, Nov 29 2015 *)
    Table[Product[(n - k + 1)^(n - 2 k + 1), {k, 1, n}], {n, 0, 20}] (* Harlan J. Brothers, Aug 26 2023 *)
  • Maxima
    a(n):= prod(binomial(n,k),k,0,n); n : 15; for i from 0 thru n do print (a(i)); /* Valentin Bakoev, May 17 2019 */
    
  • PARI
    for(n=0,16,print(prod(m=1,n,binomial(n,m))))
    
  • PARI
    A001142(n) = prod(k=1, n, k^((k+k)-1-n)); \\ Antti Karttunen, Nov 02 2014
    
  • Python
    from math import factorial, prod
    from fractions import Fraction
    def A001142(n): return prod(Fraction((k+1)**k,factorial(k)) for k in range(1,n)) # Chai Wah Wu, Jul 15 2022
  • Sage
    a = lambda n: prod(k^k/factorial(k) for k in (1..n))
    [a(n) for n in range(20)] # Peter Luschny, Nov 29 2015
    
  • Scheme
    (define (A001142 n) (mul (lambda (k) (expt k (+ k k -1 (- n)))) 1 n))
    (define (mul intfun lowlim uplim) (let multloop ((i lowlim) (res 1)) (cond ((> i uplim) res) (else (multloop (+ 1 i) (* res (intfun i)))))))
    ;; Antti Karttunen, Oct 28 2014
    

Formula

a(n) = C(n, 0)*C(n, 1)* ... *C(n, n).
From Harlan J. Brothers, Nov 26 2009: (Start)
a(n) = Product_{j=1..n-2} Product_{k=1..j} (1+1/k)^k, n >= 3.
a(1) = a(2) = 1, a(n) = a(n-1) * Product_{k=1..n-2} (1+1/k)^k. (End)
a(n) = hyperfactorial(n)/superfactorial(n) = A002109(n)/A000178(n). - Peter Luschny, Jun 24 2012
a(n) ~ A^2 * exp(n^2/2 + n - 1/12) / (n^(n/2 + 1/3) * (2*Pi)^((n+1)/2)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Jul 10 2015
a(n) = Product_{i=1..n} Product_{j=1..i} (i/j). - Pedro Caceres, Apr 06 2019
a(n) = Product_{k=1..n} (n - k + 1)^(n - 2*k + 1). - Harlan J. Brothers, Aug 26 2023

Extensions

More terms from James Sellers, May 01 2000
Better description from Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 30 2001

A249428 Numbers n such that A249151(2n+1) = n.

Original entry on oeis.org

1, 2, 4, 6, 10, 11, 12, 16, 18, 22, 28, 30, 36, 40, 42, 46, 52, 58, 60, 66, 70, 72, 78, 82, 88, 96, 100, 102, 106, 108, 112, 126, 130, 136, 138, 148, 150, 156, 162, 166, 172, 178, 180, 190, 192, 196, 198, 210, 222, 226, 228, 232, 238, 240, 250, 256, 262, 268, 270, 276, 280, 282, 292, 306, 310, 312, 316, 330, 336, 346, 348
Offset: 1

Views

Author

Antti Karttunen, Oct 28 2014

Keywords

Comments

Seems to be A006093 with at least one additional term, 11 at a(6).

Crossrefs

Cf. A006093 (seems to be a subsequence).
Cf. also A056077, A249151, A249424.

Programs

Formula

a(n) = A249151(A249424(n)).

A132858 Composite "antimutinous" numbers. An antimutinous number is an integer m > 1 where m/p^k < p, where p is the largest prime divisor of m and p^k is the largest power of p dividing m.

Original entry on oeis.org

4, 6, 8, 9, 10, 14, 15, 16, 18, 20, 21, 22, 25, 26, 27, 28, 32, 33, 34, 35, 38, 39, 42, 44, 46, 49, 50, 51, 52, 54, 55, 57, 58, 62, 64, 65, 66, 68, 69, 74, 75, 76, 77, 78, 81, 82, 85, 86, 87, 88, 91, 92, 93, 94, 95, 98, 99, 100, 102, 104, 106, 110, 111, 114, 115, 116, 117, 118
Offset: 1

Views

Author

Leroy Quet, Nov 21 2007

Keywords

Comments

{a(k)-1} is the complement of sequence A056077. In other words, {a(k)} contains precisely those positive integers m where A001142(m-1) (= product{k=1 to m-1} k^(2k-m)) is not divisible by all primes <= m-1.

Crossrefs

Programs

  • Mathematica
    antiQ[n_] := Module[{f = FactorInteger[n], p, k}, p = f[[-1, 1]]; k = f[[-1, 2]]; n/p^k < p]; Select[Range[118], CompositeQ[#] && antiQ[#] &] (* Amiram Eldar, Feb 24 2020 *)

Extensions

Extended by Ray Chandler, Nov 17 2008

A329947 Integers k such that the radical of the cumulative product of k^k/k! equals its predecessor.

Original entry on oeis.org

1, 12, 30, 36, 40, 60, 70, 72, 96, 108, 112, 126, 150, 175, 180, 192, 198, 210, 224, 240, 270, 280, 306, 312, 324, 330, 336, 350, 351, 352, 378, 384, 396, 399, 400, 408, 418, 420, 432, 440, 441, 442, 448, 455, 456, 462, 475, 490, 494, 520, 539, 540, 546, 560
Offset: 1

Views

Author

Peter Luschny, Dec 21 2019

Keywords

Comments

No prime numbers appear in this sequence.

Examples

			Consider the rows 11 and 12 of Pascal's triangle.
P11 = [1, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11, 1].
P12 = [1, 12, 66, 220, 495, 792, 924, 792, 495, 220, 66, 12, 1].
lcm(P11) = 2310 and radical(2310) = 2310.
lcm(P12) = 27720 and radical(27720) = 2310.
Since radical(lcm(P11)) = radical(lcm(P12)) 12 is in this sequence.
Also: 1 is in this sequence because radical(lcm(P0)) = radical(lcm([1])) = radical(1) = 1 = radical(lcm([1, 1])) = radical(lcm(P1)).
		

Crossrefs

Programs

  • Maple
    h := n -> mul(k^k/factorial(k), k=0..n):
    rad := n -> mul(k, k = numtheory[factorset](n)):
    g := proc(n) option remember; rad(h(n)) end:
    isA329947 := n -> g(n) = g(n-1): select(isA329947, [$1..560]);
  • Mathematica
    h[n_] := Product[k^k/k!, {k, 1, n}];
    rad[n_] := Times @@ FactorInteger[n][[All, 1]];
    g[n_] := g[n] = rad[h[n]];
    isA329947[n_] := g[n] == g[n-1];
    Select[Range[560], isA329947] (* Jean-François Alcover, Feb 28 2024, after Maple code *)
Showing 1-4 of 4 results.