|
%I #23 Nov 22 2022 02:34:52
%S 12,72,600,2352,14520,28392,83232,129960,279312,706440,922560,1872792,
%T 2824080,3416952,4877472,7887672,12113880,13842120,20146632,25406640,
%U 28392912,38943840,47451432,62734320,88519872,104050200,112540272
%N a(n) = (n-th prime)^4-(n-th prime)^2.
%H Vincenzo Librandi, <a href="/A138402/b138402.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 Product_{n>=1} (1 - 1/a(n)) = A065471.
%F From _Amiram Eldar_, Nov 22 2022: (Start)
%F a(n) = A001248(n) * A084920(n).
%F a(n) = A036689(n) * A036690(n). (End)
%t a = {}; Do[p = Prime[n]; AppendTo[a, p^4 - p^2], {n, 1, 50}]; a
%t #^4-#^2&/@Prime[Range[30]] (* _Harvey P. Dale_, Sep 19 2018 *)
%o (PARI) forprime(p=2,1e3,print1(p^4-p^2", ")) \\ _Charles R Greathouse IV_, Jun 16 2011
%o (PARI) apply(p->p^4-p^2, primes(100)) \\ _Charles R Greathouse IV_, Apr 17 2015
%o (Magma) [NthPrime((n))^4 - NthPrime((n))^2: n in [1..30] ]; // _Vincenzo Librandi_, Jun 17 2011
%Y Cf. A001248, A036689, A036690, A065471, A084920.
%K nonn,easy
%O 1,1
%A _Artur Jasinski_, Mar 19 2008
|
