En definitiva, hemos probado que $f(n)$ es la suma de los dígitos de $n$ en base $2$. Dado $k\in\mathbb{N}$, está claro entonces que el menor número con $f(n)=k$ es $n=2^k-1$ (que tiene $k$ dígitos en base $2$, todos ellos iguales a $1$). Como $2^{10}\lt 2001\lt 2^{11}-1$, deducimos que el máximo de $f(n)$ cuando $n\leq 2001$ es $10$. También deducimos que el menor natural $n$ tal que $f(n)=2001$ es $2^{2001}-1$.
Nota. Observemos que el propio cálculo anterior nos dice que $\tan(\alpha+\beta)=-1$, luego $\alpha=\arctan(2)$ y $\beta=\arctan(3)$ suman $135^\circ$. Esto nos lleva directamente a la famosa identidad \[\arctan(1)+\arctan(2)+\arctan(3)=180^\circ.\]
Nota. Las exponenciales pueden ocultar la aplicación de la desigualdad entre las medias aritmética-geométrica, pero una solución similar se tiene aplicando la desigualdade de Jensen a la función convexa $f(t)=3^t$. ¿Sabrías escribir los detalles?
La única posibilidad que nos queda es $x=2$, que nos da $y=f(x)=f(2)=2$ y $z=f(y)=f(2)=2$. En consecuencia, $x=y=z=2$ es la única solución del sistema.