login
a(n) = prime(n)^6 - prime(n)^2.
8

%I #25 Jul 04 2023 15:34:51

%S 60,720,15600,117600,1771440,4826640,24137280,47045520,148035360,

%T 594822480,887502720,2565725040,4750102560,6321361200,10779213120,

%U 22164358320,42180530160,51520370640,90458377680,128100278880

%N a(n) = prime(n)^6 - prime(n)^2.

%H Vincenzo Librandi, <a href="/A138409/b138409.txt">Table of n, a(n) for n = 1..200</a>

%H <a href="/index/Pri#prime_powers">Index to sequences related to prime powers</a>

%F From _Elmo R. Oliveira_, Jan 17 2023: (Start)

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

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

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

%t a = {}; Do[p = Prime[n]; AppendTo[a, p^6 - p^2], {n, 1, 50}]; a

%t #^6-#^2&/@Prime[Range[20]] (* _Harvey P. Dale_, Jul 04 2023 *)

%o (PARI) forprime(p=2,1e3,print1(p^6-p^2", ")) \\ _Charles R Greathouse IV_, Jun 16 2011

%o (Magma) [NthPrime((n))^6 - NthPrime((n))^2: n in [1..30] ]; // _Vincenzo Librandi_, Jun 17 2011

%Y Cf. A030516, A001248, A138441, A030514.

%K nonn,easy

%O 1,1

%A _Artur Jasinski_, Mar 19 2008