%I #15 May 12 2023 09:33:19
%S 4,10,16,22,25,34,40,46,55,58,64,82,85,88,94,100,106,115,118,121,136,
%T 142,145,154,160,166,178,184,187,202,205,214,220,226,232,235,238,250,
%U 253,256,262,265,274,289,295,298,319,322,328,334,340,346,352,355,358,376
%N Numbers of the form 3k+1 with greatest prime divisor of the form 3m-1.
%C All numbers of the form 2p, where p==2(mod 3) is prime, are in the sequence.
%H Amiram Eldar, <a href="/A176262/b176262.txt">Table of n, a(n) for n = 1..10000</a>
%t Select[3 Range@ 120 + 1, Mod[#, 3] == 2 &[FactorInteger[#][[-1, 1]]] &] (* _Michael De Vlieger_, Feb 07 2016 *)
%o (PARI) isok(n) = ((n % 3) == 1) && (n != 1) && ((vecmax(factor(n)[,1]) % 3) == 2); \\ _Michel Marcus_, Feb 07 2016
%Y Cf. A176258, A176257, A176256, A176255, A002148, A002145, A016813, A004767.
%K nonn
%O 1,1
%A _Vladimir Shevelev_, Apr 13 2010
%E Corrected and extended by _Michel Marcus_, Feb 07 2016
|