Hemos demostrado que las únicas soluciones son $x=2+\sqrt{2}$ y $x=2-\sqrt{2}$.
Nota. Si procedemos directamente simplificando la ecuación inicial y tenemos en cuenta que $x+1\neq 0$, llegamos a la ecuación de cuarto grado $x^4-6x^3+14x^2-20x+8=0$. Esta se puede factorizar sobre los enteros como producto de dos ecuaciones de segundo grado $(x^2-4x+2)(x^2-2x+4)=0$, de donde también se deduce la solución.
El caso de tener $17$ filas y $13$ columnas es completamente similar ya que la fórmula $(\star)$ se demuestra de la misma manera (ahora para $1\leq i\leq 17$ y $1\leq j\leq 13$), lo que equivale a cambiar filas por columnas (es decir, cambiar $d$ por $d_1$) en el cálculo anterior (¿sabrías justificarlo?).
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$.