login
Numbers k with A033676(k)=2, where A033676 is the largest divisor <= sqrt(k).
53

%I #58 Aug 30 2024 10:42:57

%S 4,6,8,10,14,22,26,34,38,46,58,62,74,82,86,94,106,118,122,134,142,146,

%T 158,166,178,194,202,206,214,218,226,254,262,274,278,298,302,314,326,

%U 334,346,358,362,382,386,394,398,422,446,454,458,466,478,482,502,514

%N Numbers k with A033676(k)=2, where A033676 is the largest divisor <= sqrt(k).

%C Define a sieve operation with parameter s that eliminates integers of the form s^2 + s*i (i >= 0) from the set A000027 of natural numbers. The sequence lists those natural numbers that are eliminated by the sieve s=2 and cannot be eliminated by any sieve s >= 3. - _R. J. Mathar_, Jun 24 2009

%C After a(3)=8 all terms are 2*prime; for n > 3, a(n) = 2*prime(n-1) = 2*A000040(n-1). - _Zak Seidov_, Jul 18 2009

%C From _Omar E. Pol_, Jul 18 2009: (Start)

%C A classification of the natural numbers A000027.

%C =============================================================

%C Numbers k whose largest divisor <= sqrt(k) equals j

%C =============================================================

%C j Sequence Comment

%C =============================================================

%C 1 ..... A008578 1 together with the prime numbers

%C 2 ..... A161344 This sequence

%C 3 ..... A161345

%C 4 ..... A161424

%C 5 ..... A161835

%C 6 ..... A162526

%C 7 ..... A162527

%C 8 ..... A162528

%C 9 ..... A162529

%C 10 .... A162530

%C 11 .... A162531

%C 12 .... A162532

%C ... And so on. (End)

%C The numbers k whose largest divisor <= sqrt(k) is j are exactly those numbers j*m where m is either a prime >= k or one of the numbers in row j of A163925. - _Franklin T. Adams-Watters_, Aug 06 2009

%C See also A163280, the main entry for this sequence. - _Omar E. Pol_, Oct 24 2009

%C Also A100484 UNION 8. - _Omar E. Pol_, Nov 29 2012 (after Seidov and Hasler)

%H Omar E. Pol, <a href="http://www.polprimos.com">Determinacion geometrica de los numeros primos y perfectos</a>

%H Omar E. Pol, <a href="http://www.polprimos.com/imagenespub/polprdipi.jpg">Illustration: Divisors and pi(x)</a>

%H Omar E. Pol, <a href="http://www.polprimos.com/imagenespub/poldiv08.jpg">Illustration of initial terms</a>

%H Omar E. Pol. <a href="http://www.polprimos.com/imagenespub/poldiv06.jpg">Illustration of initial terms of A008578, A161344, A161345, A161424</a>

%F Equals 2*A000040 union {8}. - _M. F. Hasler_, Nov 27 2012

%F a(n) = 2*A046022(n+1) = 2*A175787(n). - _Omar E. Pol_, Nov 27 2012

%p isA := proc(n,s) if n mod s <> 0 then RETURN(false); fi; if n/s-s >= 0 then RETURN(true); else RETURN(false); fi; end: isA161344 := proc(n) for s from 3 to n do if isA(n,s) then RETURN(false); fi; od: isA(n,2) ; end: for n from 1 to 3000 do if isA161344(n) then printf("%d,",n) ; fi; od; # _R. J. Mathar_, Jun 24 2009

%t a[n_] := If[n <= 3, 2n+2, 2*Prime[n-1]]; Table[a[n], {n, 1, 56}] (* _Jean-François Alcover_, Nov 26 2012, after _Zak Seidov_ *)

%o (PARI) a(n)=if(n>3,prime(n-1),n+1)*2 \\ _M. F. Hasler_, Nov 27 2012

%Y Cf. A000005, A018253, A160811, A160812, A161205, A161346, A033676, A008578, A161345, A161424, A161835, A162526, A162527, A162528, A162529, A162530, A162531, A162532, A163925.

%Y Second column of array in A163280. Also, second row of array in A163990.

%K easy,nonn

%O 1,1

%A _Omar E. Pol_, Jun 20 2009

%E More terms from _R. J. Mathar_, Jun 24 2009

%E Definition added by _R. J. Mathar_, Jun 28 2009