login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A114574
a(n) = p*p! where p = prime(n).
2
4, 18, 600, 35280, 439084800, 80951270400, 6046686277632000, 2311256907767808000, 594596384994354462720000, 256411097818451356681764864000000, 254907998279515607349492449280000000, 509258864375374766713691244518493388800000000, 1371553591139716091434972544191070818271232000000000
OFFSET
1,1
REFERENCES
Originally published on Nov 09 2004 by indiejade(AT)gmail.com
FORMULA
a(n) = A000040(n)*A039716(n). - Michel Marcus, Feb 08 2015
a(n) = A001563(A000040(n)). - Michel Marcus, Aug 04 2015
EXAMPLE
2 * 2! = 4, 3 * 3! = 18, 5 * 5! = 600.
MATHEMATICA
Array[Prime[#] Prime[#]! &, 13] (* Michael De Vlieger, Aug 04 2015 *)
#*#!&/@Prime[Range[15]] (* Harvey P. Dale, Sep 21 2024 *)
PROG
(Magma) [p*Factorial(p): p in PrimesUpTo(50)]; // Vincenzo Librandi, Jun 09 2013
(Python)
from sympy import factorial, prime
def a(n): p = prime(n); return p * factorial(p)
print([a(n) for n in range(1, 14)]) # Michael S. Branicky, Apr 18 2021
CROSSREFS
Cf. A000040 (prime(n)), A039716 (prime(n)!).
Sequence in context: A086400 A082022 A354303 * A295223 A173222 A342467
KEYWORD
nonn
AUTHOR
L. Shawnee Cook (shawneecook(AT)gmail.com), Feb 17 2006
STATUS
approved