A138409 a(n) = prime(n)^6 - prime(n)^2.

60, 720, 15600, 117600, 1771440, 4826640, 24137280, 47045520, 148035360, 594822480, 887502720, 2565725040, 4750102560, 6321361200, 10779213120, 22164358320, 42180530160, 51520370640, 90458377680, 128100278880

Offset: 1

Author

Artur Jasinski, Mar 19 2008

Keywords

• nonn
• easy

Links

• Vincenzo Librandi, Table of n, a(n) for n = 1..200
• Index to sequences related to prime powers

Crossrefs

Cf. A030516, A001248, A138441, A030514.

Programs

Magma

[NthPrime((n))^6 - NthPrime((n))^2: n in [1..30] ]; // Vincenzo Librandi, Jun 17 2011

Mathematica

a = {}; Do[p = Prime[n]; AppendTo[a, p^6 - p^2], {n, 1, 50}]; a
#^6-#^2&/@Prime[Range[20]] (* Harvey P. Dale, Jul 04 2023 *)

PARI

forprime(p=2,1e3,print1(p^6-p^2", ")) \\ Charles R Greathouse IV, Jun 16 2011

Formula

From Elmo R. Oliveira, Jan 17 2023: (Start)

a(n) = A030516(n) - A001248(n).

a(n) = 2 * A138441(n).

a(n) = A001248(n) * (A030514(n) - 1). (End)