Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%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