Let $r$ be an non-negative integer. Let $\lambda_k = \alpha_k + \beta_k i$ be any
set of roots of the equation:
$$ 1 + 2^{-\lambda} + 3^{-\lambda} + 4^{-\lambda} + \cdots + N^{-\lambda} = 0\tag{*1}$$
subject to the constraint $\alpha_i > r$ and $\beta_i > 0$.
For any real numbers $A_k^{\pm}$ and $\delta_k^{\pm}$, consider the function $g$ defined by:
$$g(x) = \begin{cases}
\sum_{k} A_k^{+} \Re\left(e^{i\delta_k^{+}} x^{\lambda_k} \right)\\
\sum_{k} A_k^{-} \Re\left(e^{i\delta_k^{-}} |x|^{\lambda_k} \right)
\end{cases}
=
\begin{cases}
\sum_{k} A_k^{+} x^{\alpha_k}\cos(\beta_k \log x + \delta_k^{+}), &\quad\text{ for } x \ge 0\\
\sum_{k} A_k^{-} |x|^{\alpha_k}\cos(\beta_k \log |x| + \delta_k^{-}),&\quad\text{ for } x \le 0
\end{cases}
$$
It is easy to check $f(x) = g(x) + \frac{1}{N}$ satisfy the functional equation
$$1 = f(x) + f(\frac{x}{2}) + \cdots + f(\frac{x}{N})\tag{*2}$$
provided the sum over $k$ converges. Furthermore, if $A_k^{\pm}$ decrease fast enough,
the constraint $\alpha_i > r$ will force $g$ and its first $r^{th}$ derivatives vanishes and continuous at $0$.
In short,
if the set of roots of $(*1)$ is non-empty, then there are non-trivial
$C^{r}$ solution for $(*2)$.
For an example, consider the case $N = 4$ and $r = 0$. It seems there are infinite many
solution for $(*1)$. The one with smallest $\beta$ is give by:
$$\lambda_1 = \alpha_1 + \beta_1 i \sim 0.62597108186373 + 3.127120203586539 i$$
Which leads to a non-trivial continuous solution for $(*2)$:
$$f(x) = \frac{1}{4} + |x|^{\alpha_1} \cos(\beta_1 \log|x|) \sim \frac14 + |x|^{0.62597108186373} \cos( 3.127120203586539 \log |x|)$$
Please note that for $N = 4$, it seems all the complex roots of $(*1)$ has
$\Re(\lambda) < 1$. It is not clear whether $(*2)$ has any non-trivial $C^{1}$ solution at all.