|
%I #16 Jul 05 2023 12:24:43
%S 2,7,30,100,107,135,172,217,322,352,452,562,590,667,707,917,940,975,
%T 1092,1127,1222,1470,1570,1950,2282,2545,2772,2865,2930,3007,3087,
%U 3682,3770,3840,3945,4447,4452,4477,5142,5555,5600,5625,5635,6262,6442,7520,8232
%N Numbers k such that k and k+3 are in A002822.
%e 2 is a term because 2 and 5 are in A002822.
%p isA002822 := proc(n) isprime(6*n-1) and isprime(6*n+1) ; end proc:
%p isA173233 := proc(n) isA002822( n ) and isA002822(n+3 ) ;end proc:
%p for n from 1 to 3000 do if isA173233(n) then printf("%d,",n) ; end if; end do: # _R. J. Mathar_, May 02 2010
%t SequencePosition[Table[If[AllTrue[6n+{1,-1},PrimeQ],1,0],{n,9000}],{1,_,_,1}][[;;,1]] (* _Harvey P. Dale_, Jul 05 2023 *)
%Y Cf. A002822, A173232.
%K nonn
%O 1,1
%A _Juri-Stepan Gerasimov_, Feb 13 2010
%E Definition and sequence corrected (667, 707, 1097 etc inserted) by _R. J. Mathar_, May 02 2010
%E More terms from _Jinyuan Wang_, May 15 2020
|
