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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A217276 Numbers n such that no prime divisors of n^2 + 1 are of the form a^2 + 1. 2
34, 44, 46, 50, 60, 70, 76, 86, 96, 100, 104, 114, 136, 144, 164, 186, 190, 194, 196, 214, 220, 226, 244, 246, 254, 266, 274, 286, 294, 296, 304, 316, 320, 324, 330, 334, 346, 354, 356, 360, 366, 374, 390, 410, 416, 424, 426, 434, 454, 456, 460, 476, 484, 486 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

Table of n, a(n) for n=1..54.

EXAMPLE

34 is in the sequence because 34^2+1 = 1157 = 13*89 and the prime divisors 13, 89 are not of the form a^2+1.

MAPLE

with(numtheory):for n from 1 to 100 do: x:=factorset(n^2+1):n1:=nops(x):ii:=0:for m from 1 to n1 do:y:=sqrt(x[m]-1):if y=floor(y) then ii:=1:else fi:od:if ii=0 then printf(`%d, `, n):else fi:od:

MATHEMATICA

fQ[n_] := Module[{lst = Transpose[FactorInteger[n^2 + 1]][[1]]}, Length[lst] > 1 && And @@ (Not /@ IntegerQ /@ Sqrt[lst - 1])]; Select[Range[500], fQ] (* T. D. Noe, Oct 01 2012 *)

CROSSREFS

Cf. A002522, A005574, A070303, A181436, A180252, A181436.

Sequence in context: A039417 A043240 A044020 * A076772 A248527 A206262

Adjacent sequences:  A217273 A217274 A217275 * A217277 A217278 A217279

KEYWORD

nonn

AUTHOR

Michel Lagneau, Sep 29 2012

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 18 23:26 EST 2018. Contains 299330 sequences. (Running on oeis4.)