cp's OEIS Frontend

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.

A008367 Composite but smallest prime factor >= 17.

This page as a plain text file.
%I A008367 #45 Jul 08 2025 01:03:09
%S A008367 289,323,361,391,437,493,527,529,551,589,629,667,697,703,713,731,779,
%T A008367 799,817,841,851,893,899,901,943,961,989,1003,1007,1037,1073,1081,
%U A008367 1121,1139,1147,1159,1189,1207,1219,1241,1247,1271,1273,1333,1343,1349,1357,1363,1369,1387,1403,1411,1457
%N A008367 Composite but smallest prime factor >= 17.
%C A008367 Composite numbers k such that k^720 mod 30030 = 1. - _Gary Detlefs_, May 02 2012
%C A008367 The asymptotic density of this sequence is 192/1001. - _Amiram Eldar_, Mar 22 2021
%H A008367 T. D. Noe, <a href="/A008367/b008367.txt">Table of n, a(n) for n = 1..10000</a>
%F A008367 For 1 <= n < 107, a(n) = A287391(n+2); then a(107) = 2329, a(108) = 2363 are not in A287391, but again a(n) = A287391(n) for 108 < n < 120. - _M. F. Hasler_, Oct 04 2018
%p A008367 for i from 1 to 2000 do if gcd(i,30030) = 1 and not isprime(i) then print(i); fi; od;
%t A008367 Select[ Range[ 1500 ], (GCD[ #1, 30030 ]==1&&!PrimeQ[ #1 ])& ]
%t A008367 Select[Range[2000], ! PrimeQ[#] && FactorInteger[#][[1, 1]] >= 17 &] (* _T. D. Noe_, Mar 16 2013 *)
%o A008367 (PARI) is(n)={gcd(n,30030)==1 && !ispseudoprime(n)} \\ _M. F. Hasler_, Oct 04 2018
%o A008367 (GAP) Filtered([17..1500],n->PowerMod(n,720,30030)=1 and not IsPrime(n)); # _Muniru A Asiru_, Nov 24 2018
%Y A008367 Intersection of A002808 and A008366.
%Y A008367 Cf. A038511, A084969, A084970.
%Y A008367 Cf. A287391.
%K A008367 nonn
%O A008367 1,1
%A A008367 _N. J. A. Sloane_