login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A183064 Numbers k such that k^2+1 = 2*p^2, p prime. 3
7, 41, 8119, 47321, 63018038201, 2470433131948081, 96845919575610633161, 19175002942688032928599, 5834531641231893991002972081099601, 6733044458057842709277507685523012161, 228725309250740208744750893347264645481 (list; graph; refs; listen; history; text; internal format)
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

Cf. A002315, A005574, A002731, A174492.

Sequence in context: A057006 A144747 A176090 * A290045 A188066 A225327

Adjacent sequences:  A183061 A183062 A183063 * A183065 A183066 A183067

KEYWORD

nonn

AUTHOR

Michel Lagneau, Feb 01 2011

EXTENSIONS

More terms from Charles R Greathouse IV, Feb 01 2011

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified February 20 09:02 EST 2018. Contains 299384 sequences. (Running on oeis4.)