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!)
A286658 Primes of the form p*b^b + 1, where p is a prime and b>1. 1
13, 29, 53, 149, 173, 269, 293, 317, 389, 509, 557, 653, 769, 773, 797, 1109, 1229, 1493, 1637, 1733, 1949, 1997, 2309, 2477, 2693, 2837, 2909, 2957, 3329, 3413, 3533, 3677, 3989, 4133, 4157, 4253, 4349, 4373, 4493, 4517, 5189, 5309, 5693, 5717, 5813, 6173 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

Robert Israel, Table of n, a(n) for n = 1..10000

EXAMPLE

a(1) = 3*(2^2)+1 = 13.

a(2) = 7*(2^2)+1 = 29.

a(3) = 13*(2^2)+1 = 53.

MAPLE

N:= 10000: # for all terms <= N

Res:= NULL:

P:= select(isprime, [2, seq(i, i=3..N/4, 2)]):

for b from 2  do

  q:= b^b; if q > N/2 then break fi;

  for i from 1 to nops(P) do

     x:= P[i]*q+1;

     if x > N then break fi;

     if isprime(x) then Res:= Res, x fi;

od od:

sort(convert({Res}, list)); # Robert Israel, Nov 12 2019

MATHEMATICA

nmax=10^4; pimax=PrimePi[nmax]; bmax=1; While[(bmax+1)^(bmax+1)<=nmax, bmax++]; Select[Union@Flatten@Table[Prime[pi] b^b+1, {b, 2, bmax}, {pi, pimax}], PrimeQ[#]&&#<=nmax&]

CROSSREFS

Cf. A175768, A285015.

Sequence in context: A010337 A244637 A162579 * A090866 A098062 A094481

Adjacent sequences:  A286655 A286656 A286657 * A286659 A286660 A286661

KEYWORD

nonn,easy

AUTHOR

Vincenzo Librandi, May 12 2017

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 September 24 03:26 EDT 2021. Contains 347623 sequences. (Running on oeis4.)