|
%I
%S 1,17,41,73,89,97,113,137,193,233,241,257,281,289,313,337,353,401,409,
%T 433,449,457,521,569,577,593,601,617,641,673,697,761,769,809,857,881,
%U 929,937,953,977,1009,1033,1049
%N Divisible only by primes congruent to 1 mod 8.
%H Vincenzo Librandi, <a href="/A004625/b004625.txt">Table of n, a(n) for n = 1..1000</a>
%t ok[n_] := And @@ (#1[[2]] == 1 && Mod[#1[[1]], 8] == 1 & ) /@ FactorInteger[n]; Select[ Range[2000], ok] (* From Jean-François Alcover, Mar 09 2012 *)
%o (MAGMA) [n: n in [1..2000] | forall{d: d in PrimeDivisors(n) | d mod 8 eq 1}]; // _Vincenzo Librandi_, Aug 20 2012
%K nonn
%O 1,2
%A _N. J. A. Sloane_.
%E Offset changed from 0 to 1 by _Vincenzo Librandi_, Aug 20 2012
|
