%I #18 Sep 08 2022 08:45:20
%S 9,16,36,64,144,196,324,400,576,900,1024,1444,1764,1936,2304,2916,
%T 3600,3844,4624,5184,5476,6400,7056,8100,9604,10404,10816,11664,12100,
%U 12996,16384,17424,19044,19600,22500,23104,24964,26896,28224,30276,32400,33124,36864
%N a(n) = (prime(n)+1)^2.
%H Vincenzo Librandi, <a href="/A110833/b110833.txt">Table of n, a(n) for n = 1..1000</a>
%F From _Amiram Eldar_, Jan 23 2021: (Start)
%F a(n) = A008864(n)^2.
%F Product_{n>=1} (1 + 1/a(n)) = A065486.
%F Product_{n>=1} (1 - 1/a(n)) = A065472. (End)
%t Table[(Prime[n] + 1)^2, {n, 200}] (* _Vincenzo Librandi_, Mar 27 2014 *)
%o (Magma) [(p+1)^2: p in PrimesUpTo(200)]; // _Vincenzo Librandi_, Mar 27 2014
%o (Python)
%o from sympy import primerange
%o print([(p+1)**2 for p in primerange(1, 192)]) # _Michael S. Branicky_, Sep 16 2021
%Y Cf. A008864, A065472, A065486.
%K nonn,easy
%O 1,1
%A _Giovanni Teofilatto_, Sep 18 2005
%E Corrected and extended by _Ray Chandler_, Oct 08 2005
|