A072849 Prime(a(n)) + ... + prime(a(n)+3) is a square = A051395(n)^2.
3, 21, 33, 84, 104, 199, 689, 1848, 1923, 3435, 3795, 3985, 4126, 4742, 5968, 6413, 6495, 7649, 8927, 9906, 16885, 17677, 20474, 20996, 22924, 23923, 36902, 38733, 40347, 40654, 41956, 42601, 43047, 44482, 44920, 51608, 52305, 56706, 66032
Offset: 1
Examples
a(1) = 3 because prime(3) + prime(4) + prime(5) + prime(6) = 5+7+11+13 = 36 = 6*6.
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..200 from Harvey P. Dale)
Crossrefs
Programs
-
Mathematica
Flatten[Position[Partition[Prime[Range[70000]],4,1],?(IntegerQ[ Sqrt[ Total[ #]]]&),{1},Heads->False]] (* _Harvey P. Dale, Dec 09 2015 *)
-
PARI
lista(m) = {for (n = 1, m, if (issquare(prime(n) + prime(n+1) + prime(n+2) + prime(n+3)), print1(n, ", ")););} \\ Michel Marcus, Apr 19 2013
Formula
Extensions
Definition corrected by Zak Seidov, Dec 13 2014
Comments