login
Positive integers that are divisible by neither 8k-1 nor 8k+1, for all k > 0.
1

%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