A140135 Product of largest semiprime <= n and smallest semiprime >= n.

16, 24, 36, 54, 54, 81, 100, 140, 140, 140, 196, 225, 315, 315, 315, 315, 315, 441, 484, 550, 550, 625, 676, 858, 858, 858, 858, 858, 858, 1089, 1156, 1225, 1330, 1330, 1444, 1521, 1794, 1794, 1794, 1794, 1794, 1794, 2116, 2254, 2254, 2401, 2499, 2601, 2805

Offset

4,1

Comments

This is to A030664 as semiprimes A001358 are to primes A000040. Subset of A014613.

Links

Table of n, a(n) for n=4..52.

Formula

a(n) = MAX{j in A001358 and j <= n} * MIN{j in A001358 and j >= n}

Example

a(10) = 100 because the largest semiprime <= 10 is 10, the smallest semiprime >= 10 is 10 and 10*10=100.

Maple

isA001358 := proc(n) RETURN( numtheory[bigomega](n) = 2) ; end: A001358 := proc(n) option remember ; local a; if n = 1 then 4; else for a from A001358(n-1)+1 do if isA001358(a) then RETURN(a) ; fi ; od: fi ; end: prevsemiprime := proc(n) local a; for a from n to 4 by -1 do if isA001358(a) then RETURN(a) ; fi ; od: RETURN(-1) ; end: nextsemiprime := proc(n) local a; for a from n do if isA001358(a) then RETURN(a) ; fi ; od: RETURN(-1) ; end: A140135 := proc(n) prevsemiprime(n)*nextsemiprime(n) ; end: seq(A140135(n), n=4..80) ; # R. J. Mathar, May 11 2008

Mathematica

ls[n_]:=Module[{i=0}, While[PrimeOmega[n+i]!=2, i++]; n+i]; ss[n_]:=Module[ {i=0}, While[PrimeOmega[n-i]!=2, i++]; n-i]; Table[ls[n]*ss[n], {n, 4, 60}] (* Harvey P. Dale, Oct 18 2013 *)

Crossrefs

Cf. A001358, A014613, A030664.

Sequence in context: A046370 A103248 A286195 * A349241 A120142 A110228

Adjacent sequences: A140132 A140133 A140134 * A140136 A140137 A140138

Keyword

easy,nonn

Author

Jonathan Vos Post, May 09 2008

Extensions

More terms from R. J. Mathar, May 11 2008

Status

approved