|
|
A014574
|
|
Average of twin prime pairs.
|
|
384
|
|
|
4, 6, 12, 18, 30, 42, 60, 72, 102, 108, 138, 150, 180, 192, 198, 228, 240, 270, 282, 312, 348, 420, 432, 462, 522, 570, 600, 618, 642, 660, 810, 822, 828, 858, 882, 1020, 1032, 1050, 1062, 1092, 1152, 1230, 1278, 1290, 1302, 1320, 1428, 1452, 1482, 1488, 1608
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
With an initial 1 added, this is the complement of the closure of {2} under a*b+1 and a*b-1. - Franklin T. Adams-Watters, Jan 11 2006
Also the square root of the product of twin prime pairs + 1. Two consecutive odd numbers can be written as 2k+1,2k+3. Then (2k+1)(2k+3)+1 = 4(k^2+2k+1) = 4(k+1)^2, a perfect square. Since twin prime pairs are two consecutive odd numbers, the statement is true for all twin prime pairs. - Cino Hilliard, May 03 2006
Or, single (or isolated) composites. Nonprimes k such that neither k-1 nor k+1 is nonprime. - Juri-Stepan Gerasimov, Aug 11 2009
Aside from the first term in the sequence, all remaining terms have digital root 3, 6, or 9. - J. W. Helkenberg, Jul 24 2013
Every term but the first is a multiple of 6. - Harvey P. Dale, Mar 31 2023
|
|
REFERENCES
|
Archimedeans Problems Drive, Eureka, 30 (1967).
|
|
LINKS
|
|
|
FORMULA
|
|
|
MAPLE
|
P := select(isprime, [$1..1609]): map(p->p+1, select(p->member(p+2, P), P)); # Peter Luschny, Mar 03 2011
A014574 := proc(n) option remember; local p ; if n = 1 then 4 ; else p := nextprime( procname(n-1) ) ; while not isprime(p+2) do p := nextprime(p) ; od ; return p+1 ; end if ; end proc: # R. J. Mathar, Jun 11 2011
|
|
MATHEMATICA
|
Select[Table[Prime[n] + 1, {n, 260}], PrimeQ[ # + 1] &] (* Ray Chandler, Oct 12 2005 *)
Mean/@Select[Partition[Prime[Range[300]], 2, 1], Last[#]-First[#]==2&] (* Harvey P. Dale, Jan 16 2014 *)
|
|
PROG
|
if n = 1 then
return(4),
for k : 2 step 2 do (
if primep(p+k-1) and primep(p+k+1) then
return(p+k)
)
(Haskell)
a014574 n = a014574_list !! (n-1)
a014574_list = [x | x <- [2, 4..], a010051 (x-1) == 1, a010051 (x+1) == 1]
(GAP) a:=1+Filtered([1..2000], p->IsPrime(p) and IsPrime(p+2)); # Muniru A Asiru, May 20 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|