login
Primes of the form ((p-1)/2)^2+((p+1)/2), where p is prime.
13

%I #17 Dec 22 2016 23:15:57

%S 3,7,13,31,43,73,211,241,421,463,1123,1723,2551,2971,4831,5701,6163,

%T 8011,8191,9901,11131,12433,14281,17293,19183,20023,23563,24181,28393,

%U 30103,31153,35911,37831,43891,46441,53593,60271,77563,83233,86143

%N Primes of the form ((p-1)/2)^2+((p+1)/2), where p is prime.

%C Subsequence of A002383. - _Charles R Greathouse IV_, Aug 11 2009

%H G. C. Greubel, <a href="/A163418/b163418.txt">Table of n, a(n) for n = 1..1500</a>

%e ((3-1)/2)^2+((3+1)/2)=1+2=3;

%e ((5-1)/2)^2+((5+1)/2)=4+3=7;

%e ((7-1)/2)^2+((7+1)/2)=9+4=13,..

%t f[n_]:=((p-1)/2)^2+((p+1)/2); lst={};Do[p=Prime[n];If[PrimeQ[f[p]],AppendTo[lst,f[p]]],{n,6!}];lst

%t Select[((#-1)/2)^2+((#+1)/2)&/@Prime[Range[200]],PrimeQ] (* _Harvey P. Dale_, Aug 06 2012 *)

%o (PARI) forprime(p=3,10^3,if(isprime(q=((p-1)/2)^2+(p+1)/2),print1(q,", "))) \\ _Derek Orr_, Apr 30 2015

%Y Cf. A162652.

%K nonn,easy

%O 1,1

%A _Vladimir Joseph Stephan Orlovsky_, Jul 27 2009