División de los casos según el rango de $f$. Si el rango de $f$ tiene cuatro elementos, sólo tenemos la función identidad. Si el rango de $f$ tiene tres elementos $a<b<c$, $f$ tiene tres tipos diferentes, de acuerdo a $(1,2,3,4)$ que se asigna a $(a,a,b,c),(a,b,b,c)$ o $(a,b,c,c)$. De ello se sigue que tenemos $12$ débilmente creciente en función de$S$$S$, de tal manera que $|f(S)|=3$. Si el rango de $f$ tiene sólo dos elementos $a<b$, tenemos tres tipos, de acuerdo a $(1,2,3,4)$ que se asigna a $(a,a,a,b),(a,a,b,b)$ o $(a,b,b,b)$, por lo tanto $18$ funciones. Si el rango de $f$ tiene un solo elemento, tenemos cuatro funciones. Resumiendo, la respuesta es dada por
$$ 1+12+18+4 = \color{red}{35} = \binom{7}{3}.$$
Este número es también la suma de los coeficientes de $1,x,x^2,x^3$$(1+x+x^2+x^3)^4$, no por casualidad. ¿Puede usted imaginar ¿por qué?