var('n')
f(n) = floor(n/floor(sqrt(n))) show(f)
N=300 coords = [(j,f(j)) for j in range(1,N)] P = plot(point(coords,size=2)) P+=plot(sqrt(x),(x,1,N),color='red') P