A130133 - OEIS

%I #16 Nov 02 2020 07:56:23

%S 1,2,4,8,248,731,1333,3503,17608,35003,50963,62611,82603,90148,94643,

%T 106978,201295,231311,253828,335723,364808,374573,425323,490915,

%U 592595,628015,725203,984343,1031803,1112023,1136195,1376903,1411343,1430003,1642798,1926088

%N Terms in A015922 not divisible by 3.

%C Intersection of A015922 and A001651. - _Michel Marcus_, Oct 13 2013

%H Alois P. Heinz, <a href="/A130133/b130133.txt">Table of n, a(n) for n = 1..100</a>

%p a:= proc(n) option remember; local k;

%p       for k from 1+`if`(n=1, 0, a(n-1)) while

%p       irem(k, 3)=0 or 2&^k mod k <> 8 mod k do od; k

%p     end:

%p seq(a(n), n=1..30);  # _Alois P. Heinz_, Jun 04 2014

%t {1, 2, 4, 8} ~Join~ Select[Range[2 10^6], PowerMod[2, #, #] == 8 && !Divisible[#, 3]&] (* _Jean-François Alcover_, Nov 02 2020 *)

%o (PARI) isok(n) = (n % 3) && (Mod(2^n, n) == Mod(8, n)); \\ _Michel Marcus_, Oct 13 2013

%Y Cf. A001651, A015922, A130134.

%Y Intersection with A033553 gives A277344.

%K nonn

%O 1,2

%A _Zak Seidov_, May 12 2007

%E a(17)-a(28) from _Michel Marcus_, Oct 13 2013

%E a(29)-a(36) from _Alois P. Heinz_, Jun 04 2014