This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A004625 #16 Oct 26 2023 00:45:24 %S A004625 1,17,41,73,89,97,113,137,193,233,241,257,281,289,313,337,353,401,409, %T A004625 433,449,457,521,569,577,593,601,617,641,673,697,761,769,809,857,881, %U A004625 929,937,953,977,1009,1033,1049 %N A004625 Numbers divisible only by primes congruent to 1 mod 8. %H A004625 Vincenzo Librandi, <a href="/A004625/b004625.txt">Table of n, a(n) for n = 1..1000</a> %t A004625 ok[n_] := And @@ (#1[[2]] == 1 && Mod[#1[[1]], 8] == 1 & ) /@ FactorInteger[n]; Select[ Range[2000], ok] (* _Jean-François Alcover_, Mar 09 2012 *) %o A004625 (Magma) [n: n in [1..2000] | forall{d: d in PrimeDivisors(n) | d mod 8 eq 1}]; // _Vincenzo Librandi_, Aug 20 2012 %o A004625 (PARI) is(n)=is(n%8!=1, return(0)); my(f=factor(n)[,1]); for(i=1,#f~, if(f[i]%8!=1, return(0))); 1 \\ _Charles R Greathouse IV_, Jul 02 2016 %K A004625 nonn %O A004625 1,2 %A A004625 _N. J. A. Sloane_ %E A004625 Offset changed from 0 to 1 by _Vincenzo Librandi_, Aug 20 2012