login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A181619 Numbers k such that k^2+1 = 2p,(k+1)^2+1 = 5q, (k+2)^2+1 = 10r where p, q, and r are primes. 1
11, 51, 61, 101, 221, 261, 571, 2271, 2821, 2871, 5071, 5651, 5761, 6561, 6951, 9751, 10461, 10851, 11621, 11711, 14961, 15911, 16551, 17171, 17601, 18511, 19071, 19551, 23151, 25261, 27351, 27751 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

a(n) == 1 (mod 10).

LINKS

Harvey P. Dale, Table of n, a(n) for n = 1..1000

EXAMPLE

a(2) = 51 because 51^2+1 = 2*1301, 52^2+1 = 5*541, 53^2+1 = 10*281.

MAPLE

with(numtheory):for n from 1 to 30000 do : x:=n^2+1:y:=(n+1)^2+1:z:=(n+2)^2+1:x1:=factorset(x):y1:=factorset(y):z1:=factorset(z):n1:=bigomega(x):n2:=bigomega(y):n3:=bigomega(z):if

  x1[1]=2 and n1=2 and y1[1]=5 and n2 = 2 and z1[1]=2 and z1[2]=5 and n3=3 then

  printf(`%d, `, n):else fi:od:

MATHEMATICA

ksQ[k_]:=And@@PrimeQ[{(k^2+1)/2, ((k+1)^2+1)/5, ((k+2)^2+1)/10}]; Select[ Range[30000], ksQ] (* Harvey P. Dale, Sep 01 2013 *)

PROG

(PARI) forstep(k=1, 1e5, 10, if(isprime(k^2\2+1)&isprime((k+1)^2\5+1)&isprime((k+2)^2\10+1), print1(k", ")))

CROSSREFS

Cf. A002144, A005574, A002731.

Sequence in context: A026618 A266034 A026684 * A067983 A175360 A226451

Adjacent sequences:  A181616 A181617 A181618 * A181620 A181621 A181622

KEYWORD

nonn

AUTHOR

Michel Lagneau, Jan 31 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
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 16 22:04 EDT 2021. Contains 343955 sequences. (Running on oeis4.)