OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Zak Seidov)
FORMULA
EXAMPLE
For k=15: prime(15) = 47 and prime(16) = 53, 47 + 53 = 100 = 10^2.
MATHEMATICA
lst={}; Do[p1=Prime[n]; p2=Prime[n+1]; q=(p1+p2)^0.5; If[q==IntegerPart[q], AppendTo[lst, n]], {n, 1, 9!}]; lst (* Vladimir Joseph Stephan Orlovsky, Sep 02 2008 *)
PROG
(PARI) j=[]; for(n=1, 30000, x=prime(n)+prime(n+1); if(issquare(x), j=concat(j, n))); j
(PARI) { n=0; default(primelimit, 8500000); for (m=1, 10^9, if (issquare(prime(m) + prime(m + 1)), write("b064397.txt", n++, " ", m); if (n==175, break)) ) } \\ Harry J. Smith, Sep 13 2009
(Magma) [n: n in [0..50000]| IsSquare(NthPrime(n) +NthPrime(n+1))]; // Vincenzo Librandi, Apr 06 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Sep 29 2001
STATUS
approved