%I #12 Nov 04 2024 18:13:33
%S 1,2,3,4,5,6,8,10,11,12,13,16,19,20,22,24,26,29,32,37,38,40,43,44,48,
%T 52,53,58,59,61,64,67,74,76,80,83,86,88,96,101,104,106,107,109,116,
%U 118,122,128,131,134,139,148,149,152,157,160,163,166,172,173,176,179,181,192
%N Positive integers that are divisible by neither 8k-1 nor 8k+1, for all k > 0.
%H G. C. Greubel, <a href="/A180636/b180636.txt">Table of n, a(n) for n = 1..5000</a>
%t fQ[n_] := Union[ MemberQ[{1, 7}, # ] & /@ Union@ Mod[ Rest@ Divisors@ n, 8]] == {False}; fQ[1] = True; Select[ Range@ 200, fQ] (* _Robert G. Wilson v_, Oct 06 2010 *)
%o (Python) # Works in Python 2 or 3
%o import itertools
%o for n in itertools.count(1):
%o for k in range(1, 2+n//8):
%o if n%(8*k-1)==0 or n%(8*k+1)==0:
%o break
%o else:
%o print(n)
%Y Cf. A047424. - _Robert G. Wilson v_, Oct 06 2010
%K easy,nonn
%O 1,2
%A _Glenn G. Chappell_, Sep 13 2010
%E More terms from _Robert G. Wilson v_, Oct 06 2010