A297707 a(n) = Product_{k=1..n-1} n!k, where n!k is k-tuple factorial of n.

1, 2, 18, 768, 90000, 44789760, 30494620800, 121762322841600, 393644011735296000, 5618427494400000000000, 107587910030480590233600000, 5951222311476064581656248320000, 176804782652901880753915871232000000, 69819090744423637487544223697731584000000

Offset

1,2

Comments

What is the least n > 2 for which a(n) - prevprime(a(n)) is a composite number? If such a number n exists, it is greater than 250.

The least n for which nextprime(a(n)) - a(n) is a composite number is 158.

Links

Michel Marcus, Table of n, a(n) for n = 1..100

Formula

a(n) = Product_{t=1..n-1} (Product_{k=0..floor((n-1)/t)} (n-t*k)).

a(n) = (n^(n-1))*Product_{k=1..n-1} k^tau(n-k).

Example

a(2) = (2!1) = (2*1) = 2;
a(3) = (3!1)*(3!2) = (3*2*1)*(3*1) = 18;
a(4) = (4!1)*(4!2)*(4!3) = (4*3*2*1)*(4*2)*(4*1) = 768;
a(5) = (5!1)*(5!2)*(5!3)*(5!4) = (5*4*3*2*1)*(5*3*1)*(5*2)*(5*1) = 90000.

Maple

b:= proc(n, k) option remember; `if`(n<1, 1, n*b(n-k, k)) end:
a:= n-> mul(b(n, k), k=1..n-1):
seq(a(n), n=1..20);  # Alois P. Heinz, Dec 02 2018

Mathematica

Array[(#^(# - 1)) Product[k^DivisorSigma[0, # - k], {k, # - 1}] &, 13] (* Michael De Vlieger, Jan 04 2018 *)

Prog

(PARI) a(n) = (n^(n-1))*prod(k=1, n-1, k^numdiv(n-k)); \\ Michel Marcus, Dec 02 2018

Crossrefs

Cf. A000142, A006882, A007661, A007662, A085157, A085158, A114799, A114800.

Cf. A114806, A288327, A006990, A033933.

Sequence in context: A013035 A350008 A132520 * A131631 A306789 A015190

Adjacent sequences: A297704 A297705 A297706 * A297708 A297709 A297710

Keyword

nonn

Author

Lechoslaw Ratajczak, Jan 03 2018

Status

approved