|
%I #23 Mar 20 2024 05:14:43
%S 2,6,18,20,26,36,44,50,54,68,90,98,114,116,138,140,186,198,228,230,
%T 260,278,300,308,326,354,386,426,438,470,498,516,524,534,546,548,596,
%U 614,644,660,690,714,716,740,746,804,818,870,900,938,966,986,996,998,1044,1068
%N Numbers k in A045954 such that 2*k+1 is prime.
%H Amiram Eldar, <a href="/A045962/b045962.txt">Table of n, a(n) for n = 1..10000</a>
%t seq[max_] := Module[{lst = Range[2, max, 2], i = 2, len}, While[i <= (len = Length@lst) && (k = lst[[i]]) <= len, lst = Drop[lst, {k, len, k}]; i++]; Select[lst, PrimeQ[2*# + 1] &]]; seq[1000] (* _Amiram Eldar_, Mar 20 2024, after _Robert G. Wilson v_ at A045954 *)
%Y Cf. A045954, A045963.
%Y Subsequence of A005097.
%K nonn
%O 1,1
%A _Felice Russo_
%E More terms from _David W. Wilson_
%E Title corrected by _Sean A. Irvine_, Mar 29 2021
|
