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!)
A181436 Numbers k such that the prime divisors of k^2 + 1 are of the form q^2 + 1. 4
1, 2, 3, 4, 6, 7, 10, 13, 14, 16, 20, 24, 26, 36, 38, 40, 43, 54, 56, 66, 68, 74, 84, 90, 94, 110, 116, 117, 120, 124, 126, 130, 134, 146, 150, 156, 160, 170, 176, 180, 183, 184, 204, 206, 210, 224, 230, 236, 240, 250, 256, 260, 264, 270, 280, 284, 293, 300, 306, 314, 326, 327 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Ivan Neretin, Table of n, a(n) for n = 1..1000

EXAMPLE

183 is in the sequence because 183^2 + 1 = 2*5*17*197 and 2 = 1^2 + 1, 5 = 2^2+1, 17 = 4^2+1 and 197 = 14^2 + 1.

MAPLE

with(numtheory):nn:=1000:for n from 1 to nn do: x:=n^2+1:y:=factorset(x):ny:=nops(y):id:=0:for

  q from 1 to ny do: z:=y[q]-1:zz:=sqrt(z):if zz=floor(zz) then id:=id+1:else  fi:od:if id=ny then printf(`%d, `, n):else fi:od:

MATHEMATICA

Select[Range@330, And @@ IntegerQ /@ Sqrt[FactorInteger[#^2 + 1][[All, 1]] - 1] &] (* Ivan Neretin, Aug 31 2016 *)

PROG

(PARI) isok(n) = {fn = factor(n^2+1)[, 1]; for (k=1, #fn, if (!issquare(fn[k]-1), return (0)); ); 1; } \\ Michel Marcus, Sep 01 2016

CROSSREFS

Cf. A180252, A002144, A002522, A005574.

Sequence in context: A239468 A336815 A119793 * A199118 A035941 A039854

Adjacent sequences:  A181433 A181434 A181435 * A181437 A181438 A181439

KEYWORD

nonn

AUTHOR

Michel Lagneau, Jan 29 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 January 23 15:15 EST 2021. Contains 340385 sequences. (Running on oeis4.)