Assume that $ n $ is a given positive integer. Determine all functions $ f : \mathbb R \to \mathbb R $ such that $$ f ( x + y ) = f ( x ) + f ( y ) \tag 0 \label 0 $$ forall $ x , y \in \mathbb R $, and $$ f \left( x ^ { n + 1 } \right) = x ^ n f ( x ) \tag 1 \label 1 $$ for all $ x \in \mathbb R $.
I've recently encountered two posts (here and here) where the problem was somehow similar to the above one; one was in fact equivalent to the above question for $ n = 3 $, and the other one was similar to the case $ n = 1 $. While those posts asked for solutions with additional regularity assumptions like continuity and differentiability, I found out that \eqref{0} and \eqref{1} themselves were enough for showing that all solutions are of the form $ f ( x ) = a x $ for some constant $ a \in \mathbb R $, and posted answers to them not using the additional regularity conditions.
Let's consider the case $ n = 2 $. Letting $ a = f ( 1 ) $ and using \eqref{1} together with \eqref{0} and its well-known consequences, we have $$ \left( x ^ 2 + 2 x + 1 \right) \big( f ( x ) + a \big) = ( x + 1 ) ^ 2 f ( x + 1 ) = f \left( ( x + 1 ) ^ 3 \right) \\ = f \left( x ^ 3 \right) + f \left( 3 x ^ 2 \right) + f ( 3 x ) + f ( 1 ) = 3 f \left( x ^ 2 \right) + \left( x ^ 2 + 3 \right) f ( x ) + a \text . \tag 2 \label 2 $$ Similarly, by calculating $ f \left( ( x - 1 ) ^ 3 \right) $ in two different ways, we get $$ \left( x ^ 2 - 2 x + 1 \right) \big( f ( x ) - a \big) = - 3 f \left( x ^ 2 \right) + \left( x ^ 2 + 3 \right) f ( x ) - a \text . \tag 3 \label 3 $$ Adding \eqref{2} and \eqref{3} (and dividing by $ 2 $) we get $$ \left( x ^ 2 + 1 \right) f ( x ) + 2 a x = \left( x ^ 2 + 3 \right) f ( x ) \text , $$ which shows that $ f ( x ) = a x $.
For $ n = 1 $, the problem can be solved more easily, as we only need to consider $ f \left( ( x + 1 ) ^ 2 \right) $ (see my answer to the related post). For $ n = 3 $, the "$ f \left( ( x \pm 1 ) ^ 2 \right) $ technique" above gives us an equation containing $ f \left( x ^ 2 \right) $ and $ f ( x ) $, and then we can use this equation and continue similar to the case of $ n = 1 $ (see my answer to the related post).
It seems that as \eqref{1} takes the multiplicative structure of $ \mathbb R $ into account in a very specific way, it could be true than for any $ n $ the solutions are linear. But the problem with my above approach for $ n \ge 4 $ is that using the expressions for $ f \left( ( x \pm 1 ) ^ { n + 1 } \right) $, there will be more than one term of the form $ f ( x ^ m ) $ ($ m > 1 $) in the equation, and I don't know how to handle that. Is there any other way to get around this problem? Or is there in fact a nonlinear solution for larger $ n $?