%I #13 Sep 08 2022 08:45:33
%S 4,27,250,1029,6655,13182,39304,61731,133837,341446,446865,911754,
%T 1378420,1669647,2387929,3870802,5955991,6809430,9925179,12526885,
%U 14004612,19228521,23443267,31018636,43808304,51515050,55729077
%N a(n) = (prime(n)^4 - prime(n)^3)/2.
%H Vincenzo Librandi, <a href="/A138423/b138423.txt">Table of n, a(n) for n = 1..200</a>
%F a(n) = A138403(n)/2. - _R. J. Mathar_, Oct 15 2017
%t a = {}; Do[p = Prime[n]; AppendTo[a, (p^4 - p^3)/2], {n, 1, 50}]; a
%o (Magma) [(NthPrime((n))^4 - NthPrime((n))^3)/2: n in [1..30] ]; // _Vincenzo Librandi_, Jun 18 2011
%o (PARI) forprime(p=2,1e3,print1((p^4-p^3)/2", ")) \\ _Charles R Greathouse IV_, Jul 15 2011
%K nonn,easy
%O 1,1
%A _Artur Jasinski_, Mar 19 2008
|