Encontrar todos los números naturales $n$ tal que % $ $$n\mid 1^n+2^n+\cdots+n^n.$
Sabemos a través de la fórmula de Faulhaber, %#% $ de #% donde $$\sum_{k=1}^{n}k^n=\frac1{n+1}\sum_{k=0}^n\binom{n+1}{k}B_k n^{n+1-k},$ es un número de Bernoulli. Revisé algunos docena valores de $B_k$ y parece que los números impares sólo son soluciones. ¿Cualquier ideas sobre cómo proceder de aquí?