%I #14 Sep 08 2022 08:46:04
%S 73,1801,3529,10369,20809,103969,115201,426889,649801,2080801,2205001,
%T 2654209,3266569,3328201,4428289,5171329,10017289,10672201,11347849,
%U 14709889,21780001,22177801,28395649,29675809,30701449,32320801,35583049,40176649,41368609
%N Primes p such that p+1 is the sum of the squares of a pair of twin primes.
%C Primes in A184417.
%C Obviously, no prime has the form q^2+(q+2)^2+1, where q and q+2 are twin primes.
%H Bruno Berselli, <a href="/A224505/b224505.txt">Table of n, a(n) for n = 1..1000</a>
%e 3529 (prime) is in the sequence because 3529+1 = 41^2+43^2, where 41 and 43 are twin primes.
%p A224505:=proc(q) local a,n;
%p for n from 1 to q do
%p if ithprime(n+1)-ithprime(n)=2 then a:=ithprime(n+1)^2+ithprime(n)^2-1;
%p if isprime(a) then print(a); fi; fi;
%p od; end: A224505(10^6); # _Paolo P. Lava_, Apr 17 2013
%t Select[(#[[1]]^2 + #[[2]]^2 - 1) & /@ Select[Partition[Prime[Range[700]], 2, 1], #[[2]] - #[[1]] == 2 &], PrimeQ]
%o (Magma) [p: r in PrimesUpTo(5000) | IsPrime(r+2) and IsPrime(p) where p is 2*r^2+4*r+3];
%Y Cf. A063533 (sums of the squares of a pair of twin primes), A118072 (primes which are sum of a pair of twin primes minus 1), A184417.
%K nonn
%O 1,1
%A _Bruno Berselli_, Apr 08 2013
|