|
MATHEMATICA
|
Restop = Function[n, pp = Prime[Range[PrimePi[n]]]; Sum[n/pp[[k]] - Floor[n/pp[[k]]], {k, 1, Length[pp]}]]; oops = Floor[Map[Restop, Range[1, 30000]]]; Complement[Range[Max[oops]], Union[oops]]
f[n_] := Floor[Plus @@ FractionalPart[n / Table[ Prime[i], {i, PrimePi[n]}]]]; Do[ a = f[n]; If[ Position[t, a] == {}, AppendTo[t, a]], {n, 81900}]; Complement[ Range[ Sort[t][[ -1]]], t] (from Robert G. Wilson v, Feb 12 2005)
|
