OFFSET
1,1
COMMENTS
Subset of A002315 (Numbers k such that k^2 + 1 = 2*q^2).
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..22
EXAMPLE
a(2) = 41 because 41^2+1 = 2*29^2.
MAPLE
with(numtheory):for n from 1 to 1000000 do : p:=ithprime(n):x:=2*p^2: y:=sqrt(x-1):if
y=floor(y) then print(y):else fi:od:
PROG
(PARI) list(lim)=my(v=List(), w=3+quadgen(32), k, n); while((k=imag((1+w)*w^n++))<=lim, if(ispseudoprime(sqrtint((k^2+1)/2)), listput(v, k))); Vec(v) \\ Charles R Greathouse IV, Sep 14 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Feb 01 2011
EXTENSIONS
More terms from Charles R Greathouse IV, Feb 01 2011
STATUS
approved