|
|
A209255
|
|
Integers n such that (prime(n)*prime(n+1)+1)/2 is prime.
|
|
1
|
|
|
91, 327, 404, 653, 802, 804, 885, 999, 1003, 1045, 1208, 1362, 1365, 1452, 1461, 1541, 1551, 1717, 1738, 1759, 1820, 1877, 1883, 1908, 1939, 1943, 2064, 2135, 2146, 2155, 2214, 2313, 2316, 2325, 2498, 2507, 2508, 2555, 2613, 2623, 2686, 2736, 2750, 2905
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Corresponding primes: 111847, 2400169, 3872527, 11921827, 18954307, 19127041, 23646547, 31307767, 31545607, 3478606, 47990587, 63472627, 63765907, 73580509, 74823127, 83812387, 84773149, 107501767, 11006826, 113507227, 121882867, 130136827, 131300941, 135548041, 141103129, 141742267.
|
|
LINKS
|
|
|
EXAMPLE
|
a(1) = 91 because prime(91) = 467, prime(92) = 479 and (467*479+1)/2 = 111847 (prime).
|
|
MATHEMATICA
|
Select[Range[2000], PrimeQ[(Prime[#]*Prime[# + 1] + 1)/2] &] (* G. C. Greubel, Jan 05 2018 *)
Position[Partition[Prime[Range[3000]], 2, 1], _?(PrimeQ[(#[[1]]#[[2]]+1)/2]&), 1, Heads->False]//Flatten (* Harvey P. Dale, Jun 30 2022 *)
|
|
PROG
|
(PARI) for(n=2, 2000, if(isprime((prime(n)*prime(n+1)+1)/2), print1(n", ")))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|