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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A180252 Numbers where all prime divisors are of the form k^2+1. 5
 1, 2, 4, 5, 8, 10, 16, 17, 20, 25, 32, 34, 37, 40, 50, 64, 68, 74, 80, 85, 100, 101, 125, 128, 136, 148, 160, 170, 185, 197, 200, 202, 250, 256, 257, 272, 289, 296, 320, 340, 370, 394, 400, 401, 404, 425, 500, 505, 512, 514 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 LINKS Ivan Neretin, Table of n, a(n) for n = 1..10000 FORMULA Sum_{n>=1} 1/a(n) = Product_{p in A002496} p/(p-1) = Product_{k in A005574} (1 + 1/k^2) = 2.809865... - Amiram Eldar, Sep 27 2020 EXAMPLE a(17) = 74 because 74 = 2*37 = (1^2+1)*(6^2+1). MAPLE with(numtheory):T:=array(1..50):U:=array(1..1000):k:=1:for m from 1 to 300   do:x:=m^2+1:if type(x, prime)=true then T[k]:=x:k:=k+1:else fi:od:for x from   2 to 2000 do: B:=factorset(x):yy:=nops(B):A:=convert(T, set):if A intersect   B = B then printf(`%d, `, x):else fi:od: MATHEMATICA Select[Range@520, And @@ IntegerQ /@ Sqrt[FactorInteger[#][[All, 1]] - 1] &] (* Ivan Neretin, Aug 31 2016 *) CROSSREFS Cf. A002496, A002522, A005574. Sequence in context: A032937 A126684 A089653 * A191203 A216686 A114652 Adjacent sequences:  A180249 A180250 A180251 * A180253 A180254 A180255 KEYWORD nonn AUTHOR Michel Lagneau, Jan 20 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.

Last modified June 15 07:12 EDT 2021. Contains 345043 sequences. (Running on oeis4.)