login
a(n) = prime(n)^6 - prime(n).
4

%I #19 Feb 13 2023 03:36:30

%S 62,726,15620,117642,1771550,4826796,24137552,47045862,148035866,

%T 594823292,887503650,2565726372,4750104200,6321363006,10779215282,

%U 22164361076,42180533582,51520374300,90458382102,128100283850

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

%H Vincenzo Librandi, <a href="/A138408/b138408.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 a(n) = A030516(n) - A000040(n). - _Elmo R. Oliveira_, Jan 28 2023

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

%t #^6-#&/@Prime[Range[20]] (* _Harvey P. Dale_, Jun 09 2013 *)

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

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

%Y Cf. A000040, A030516.

%K nonn,easy

%O 1,1

%A _Artur Jasinski_, Mar 19 2008